grads如何画walker环流？

姚小娟

我尝试使用NCEP/NCAR再分析资料画walker环流，使用到的资料有uwnd.mon.mean.nc和omega.mon.mean.nc，walker环流的计算思路是：1、求各格点上uwnd和omega的多年平均值（1981-2010年），2、求各月份uwnd和omega的距平值，3、求步骤2结果的5S-5N的平均值，4画出50E-60W的经度-高度风场距平剖面图（grads中输入命令为:d uwnd;-omega.2*1000)。grads处理程序如下：
*************************************************************************************
'reinit'
'sdfopen E:/NCEP/a/E/paper/TC/WS/uwnd.mon.mean.nc'

*-- 写入文件 --*
'set gxout fwrite'
'set fwrite E:/NCEP/a/E/paper/TC/WS/uwnd.mon.lon_alt.anom.dat'
'set x 1 144'
'set y 1 73'
z0=1
while(z0<=12)
'set z 'z0

*-----求1981-2010年气候月平均-----*
'define jan=ave(uwnd,t=397,t=756,12)'
'define feb=ave(uwnd,t=398,t=756,12)'
'define mar=ave(uwnd,t=399,t=756,12)'
'define apr=ave(uwnd,t=400,t=756,12)'
'define may=ave(uwnd,t=401,t=756,12)'
'define jun=ave(uwnd,t=402,t=756,12)'
'define jul=ave(uwnd,t=403,t=756,12)'
'define aug=ave(uwnd,t=404,t=756,12)'
'define sep=ave(uwnd,t=405,t=756,12)'
'define oct=ave(uwnd,t=406,t=756,12)'
'define nov=ave(uwnd,t=407,t=756,12)'
'define dec=ave(uwnd,t=408,t=756,12)'


*-----求从1948年至最新月份的逐月距平值-----*
t0=1
while(t0<=795)
'set t 't0

if(math_mod(t0,12)=1)
'define u=ave(uwnd-jan,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=2)
'define u=ave(uwnd-feb,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=3)
'define u=ave(uwnd-mar,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=4)
'define u=ave(uwnd-apr,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=5)
'define u=ave(uwnd-may,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=6)
'define u=ave(uwnd-jun,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=7)
'define u=ave(uwnd-jul,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=8)
'define u=ave(uwnd-aug,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=9)
'define u=ave(uwnd-sep,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=10)
'define u=ave(uwnd-oct,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=11)
'define u=ave(uwnd-nov,lat=-5,lat=5)'
endif

if(math_mod(t0,12)=0)
'define u=ave(uwnd-dec,lat=-5,lat=5)'
endif

'd u'
t0=t0+1
endwhile
z0=z0+1
endwhile
'disable fwrite'
'reinit'
我的问题是：我画2014年3月的walker环流对比气候中心的相同月份的walker环流，结果相差很远，请问大虾，我上面的walker环流算法是否正确，如果不正确，该如何计算walker环流呢？
下图：左图为我自己画的50E-60W的经度-高度风场距平剖面图，右图为气候中心网站上的图

姚小娟

补充一句，本想发脚本编辑版的，但是选择那一项发布不了，只能选择图形美化这个版了

wmlydia

也要画这个的飘过

river

d uwnd;-omega.2*1000，您老人家是怎么判断需要在垂直速度上扩大1000倍呢？根据在哪里？
一般人只是扩大100倍，和东西风分量达到一个量级······

sun_shine_Xia

wmlydia 发表于 2014-5-6 18:41
也要画这个的飘过

看好你啊！不要简简单单飘过~

姚小娟

river 发表于 2014-5-7 17:05
d uwnd;-omega.2*1000，您老人家是怎么判断需要在垂直速度上扩大1000倍呢？根据在哪里？
一般人只是扩大10 ...

如果乘以100的话吹着的分量就体现不出来了，我还想请问垂直速度单位为Pa/s，水平速度单位为m/s，两个速度可以直接画风适量的吗？

river

姚小娟 发表于 2014-5-12 09:44
如果乘以100的话吹着的分量就体现不出来了，我还想请问垂直速度单位为Pa/s，水平速度单位为m/s，两个速度 ...

垂直速度乘以100以后和水平速度直接合成但是合适的。你有没有去了解垂直速度Pa/s的量级大概是多少？水平风速的量级是多少？然后在确定垂直速度扩大多少倍吧，而不是想当然的乘以1000吧

更大华

能不能给出最终的图

姚小娟

river 发表于 2014-5-12 09:48
垂直速度乘以100以后和水平速度直接合成但是合适的。你有没有去了解垂直速度Pa/s的量级大概是多少？水平 ...

画出的图如下：
图1是omega*100，图2为omega*1000，如果使用omega*100，基本上是看不出垂直速度的大小的。

river

姚小娟 发表于 2014-5-12 11:27
画出的图如下：
图1是omega*100，图2为omega*1000，如果使用omega*100，基本上是看不出垂直速度的大小的 ...

哦，那您老人自己判断一下吧，反正一般情况下都是扩大100倍