计算并绘制波作用通量

AXD

用ncl基于NCEP逐月数据计算并绘制TN(2001)和Plumb(1985)水平波作用通量，代码为函数，可以在ncl中调用，调用方法写成了注释。Code主要根据网站http://www.atmos.rcast.u-tokyo.ac.jp/nishii/programs/index.html提供的代码修改。由于不同的软件处理数据可能有所不同，因此绘图结果可能略有不同。网址https://www.jianshu.com/p/97042d9019c0给出了利用Pyhon绘制的脚本，亦可参考。图片为2014年11月TN01波作用通量绘制结果，填色为散度，等值线为位势高异常，箭头为波作用通量。附件为代码和相关计算方法的原始参考文献。图2是最近对比python程序（https://mp.weixin.qq.com/s/Q1nl-n1tXUq3Tw1PvQPbrQ）画的，供大家参考。

注：微信公众号“海气风华”有Code！

张如也

【摸鱼咯】在简书上所给出的python脚本中TN01的Fy计算不对哦，并不是因为不同软件绘图的差异，我已经给他指出了。我认为楼主这个才是正确的。仔细对比两人的结果图【注意看东亚地区】可以看到楼主没有很强的向南波作用通量，这符合日本人给出的结果。

咩咩的羊

感谢大佬！！！！

moerHuang

感谢大佬！！

jksuperman

牛X！！！

44444

您好！首先感谢大神，我有两个小问题想请教
请问程序里的
it_s_c = 198101
it_e_c = 201012
为什么是这两个时间呢？
还有
itime = 0,ntime-13,12
为什么是减13？
期待您的答复

AXD

关于循环程序里：

do itime = 0,ntime-13,12

    cuvar(itime:itime+11,:,:) = buvar(:,:,:)
    cvvar(itime:itime+11,:,:) = bvvar(:,:,:)
end do

有人反馈说如果ntime-13，则最后1年的算不出来，应该是ntime-1。确实如此！
这个循环是以12为间隔，假如ntime = 80（ncl计到79），如果设置为ntime-13，那么最后一循环的结束值是67。如设置为ntime-1，则最后一循环的结束值是79。少一年基本不影响结果，但为了严谨，大家还是按ntime-1计算。


如果大家还有改进或修正的想法，欢迎一起交流。大家也可以与其他人写的代码对照一下。
在此致谢我的师弟whw！

AXD

44444 发表于 2020-7-19 21:49
您好！首先感谢大神，我有两个小问题想请教
请问程序里的
it_s_c = 198101

198101-201012是我们习惯使用的气候态（30年）。目前国际上都这么用，一个不成文的规定吧。减13是隔12一循环，最后肯定到不了ntime，所以得减掉一个循环周期。由于ncl从0开始计数，所以减的是13而不是12。比如总数是80，但算出来的ntime是79。即便你循环到最后一个数，也是ntime-1。至于减多少，这个你也可以自己尝试，然后就理解了

44444

AXD 发表于 2020-7-20 08:21
198101-201012是我们习惯使用的气候态（30年）。目前国际上都这么用，一个不成文的规定吧。减13是隔12一 ...

好的，谢谢大神

biguo

想问一下，do itime=0,ntime-13,12那里，改成do itime=0,ntime-12,12是可行的吗

AXD

biguo 发表于 2020-7-25 07:44
想问一下，do itime=0,ntime-13,12那里，改成do itime=0,ntime-12,12是可行的吗

好好理解一下这段代码是干嘛的，你就知道你的想法可不可行了。