最后的输出变量是f 

这里面调用了一个关键的子程序，spectral_nudging_filter_fft_2d_ncar 

这是一个谱逼近的二维快速傅里叶变换过滤的程序。

二维的滤波这个还不是很清楚，找找技术文件看下

超详细易懂FFT（快速傅里叶变换）及代码实现_傅立叶变换编程-CSDN博客

快速傅里叶变换（FFT）超详解 - 星夜的文章 - 知乎
https://zhuanlan.zhihu.com/p/347091298

 看起来这个想要把这个FFT实现还是相对困难的,而且还需要和并行计算一起考虑，更加困难。