请教关于波通量的计算

MacTaylor

需要算波通量，其中在一个循环部分卡住了，请教各位，部分脚本如下：

'set x 1 720'
'set y 1 361'
'set t 1'

* making basic state (climatology)
'define uclm = ave(UGRDprs,time=jan1979,time=jan1998,1yr)'
'define vclm = ave(VGRDprs,time=jan1979,time=jan1998,1yr)'
'define zclm = ave(HGTprs,time=jan1979,time=jan1998,1yr)'

* anomalies
'set t 1 60'
'define zaa=HGTprs-zclm'

* QG stream function
'define psiaa=g/f*zaa'

* magnitude of basic state wind speed
'define magU = mag(uclm,vclm)'

现在这样等于数据中每个月都和1月的气候态值往下计算，是不合理的，我想要1月至12月分别各有一个气候态，然后1月的值对应1月的气候态计算，以此类推，12月的值对应12月的气候态计算，然后每12个月一个循环，一共5年，请教各位高手这个循环应该怎么做啊？谢谢了！

做个霸气的木头

你看看这个帖子有相关算法讲解http://bbs.06climate.com/forum.php?mod=viewthread&tid=7483

xiaohundanpipi

楼主最后解决问题了吗？也有相同困扰，希望帮助

MacTaylor

xiaohundanpipi 发表于 2013-3-12 16:46

                               
登录/注册后可看大图

楼主最后解决问题了吗？也有相同困扰，希望帮助

最后用的最笨的办法，把每个月的气候态算出来，然后跟每个月的再去计算

xiaohundanpipi

楼主的这个grads脚本最后画得是plumb波通量的散度么？我想画垂直方向的抛面图，请问应该怎么修改这个程序？

MacTaylor

xiaohundanpipi 发表于 2013-3-13 09:24

                               
登录/注册后可看大图

楼主的这个grads脚本最后画得是plumb波通量的散度么？我想画垂直方向的抛面图，请问应该怎么修改这个程序？ ...

这个是单层的，如果要垂直方向的剖面，我的理解就是把每一层算出来，然后把经纬度其中一个固定吧。

xiaohundanpipi

我不太理解它这个程序，d maskout( zaa*abs(f0/f),  abs(lat)-10)这个语句画得是什么？还有d skip(px,6,4);maskout( py , abs(lat)-10)这个呢？按它的脚本公式推出来和我找到的plumb通量的公式有区别... 矢量画出来的是散度么？

MacTaylor

xiaohundanpipi 发表于 2013-3-13 09:55

                               
登录/注册后可看大图

我不太理解它这个程序，d maskout( zaa*abs(f0/f),  abs(lat)-10)这个语句画得是什么？还有d skip(px,6,4); ...

其实我也不太了解，只是之前随便看了下画出来的结果，也没用上，你也是网上找的那个吧？按照他的注释：
* stream-function-like geopotential height
d maskout( zaa*abs(f0/f),  abs(lat)-10)
即流函数
* horizontal wave-activity flux
set arrscl 0.5 20
d skip(px,6,4);maskout( py , abs(lat)-10)
即波通量

xiaohundanpipi

MacTaylor 发表于 2013-3-13 10:09

                               
登录/注册后可看大图

其实我也不太了解，只是之前随便看了下画出来的结果，也没用上，你也是网上找的那个吧？按照他的注释：
...

嗯嗯 好的 非常感谢..

MacTaylor

xiaohundanpipi 发表于 2013-3-13 10:14

                               
登录/注册后可看大图

嗯嗯 好的 非常感谢..

不客气