求1948-2011年的每月平均水汽通量，每月平均水汽通量图一样，不知道是哪儿错了？

冰点

hustlh 发表于 2012-8-15 15:11

                               
登录/注册后可看大图

也就是说求时间平均的话，还是得用while循环，是吧？？？  每循环一次，改一次循环时间，就每得一张图？？ ...

http://bbs.06climate.com/forum.p ... d=82&extra=page%3D1 可以参考下这个帖子，挺不错的，版主还给予了注释。或者你不进行时间循环，假定先计算1月的，一个GS文件只输出一张图，完了修改下时间计算2、3月……，再看得出的图是否还一样

hustlh

冰点 发表于 2012-8-15 15:22

                               
登录/注册后可看大图

http://bbs.06climate.com/forum.p ... d=82&extra=page%3D1 可以参考下这个帖子，挺不错的，版主还给予 ...

好的

我试试看

hustlh

冰点 发表于 2012-8-15 15:22

                               
登录/注册后可看大图

http://bbs.06climate.com/forum.p ... d=82&extra=page%3D1 可以参考下这个帖子，挺不错的，版主还给予 ...

这个能够得出来。  按照你的方法我先不管64年的月平均，直接设置t从1变到12.即先设置set t 1   ，画1948年1月份的，之后再将t改成2，得出1948年2月份的，就这样一直改到set t 12  。这个到没有什么问题，画出的图都不一样，也应该是正确的。    就是，这个帖子里面求的水汽通量是700hpa的水汽通量，那个求整层的的话，z怎样设置啊？？？？     拜托你了！！！    单层（也就是某层）的水汽通量我会做，但是整层的话，我就有点晕了！


你是搞气象的吗？？能不能留个QQ，我单独请教啊，拜托了！ 因为我是学水的，最近才接触这个，什么都不懂，老板只交代让我做，具体什么都没说，做的过程中遇到了很多问题啊！

做个霸气的木头

hustlh 发表于 2012-8-15 15:11

                               
登录/注册后可看大图

也就是说求时间平均的话，还是得用while循环，是吧？？？  每循环一次，改一次循环时间，就每得一张图？？ ...

你的gs我修改了下 ，有断开的地反是修改了的，设定了相应的时间和层次，根据我上面说的你理解理解，应该没问题了。如果你想12个月的一起画出来就得用循环了，只需要在最后那段加循环。
'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen f:\data\monthly\shum.mon.mean.nc'  
'sdfopen f:\data\monthly\uwnd.mon.mean.nc'
'sdfopen f:\data\monthly\vwnd.mon.mean.nc'
'sdfopen f:\data\monthly\pres.mon.mean.nc'
'set lat 4 53'
'set xlint 10'
'set lon 73 135'
'set ylint 10'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld mres'
'set map 1 1 1'
'draw map'
'set ylopts 1 4 0.12'
'set xlopts 1 4 0.12'
'set cthick 4'
'set t 1 757'
'set z 1 8'
'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'set z 1'
'set t 1 757'
'define uqs=vint(pres.4(z=1),uq,300)/1000'
'define vqs=vint(pres.4(z=1),vq,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 1 12'
'define w1=ave(sq,t+0,t=757,12)'
'define w2=ave(uqs,t+0,t=757,12)'
'define w3=ave(vqs,t+0,t=757,12)'

'set t 1'
'set gxout shaded'
'd w1'
'set grads off'
'set gxout contour'
'set cint 40'
'd w1'
'set grads off'
'set gxout vector'
'd w2;w3'
'run cbarn'
'draw title 1948-2011 Jan average vopor flux china'
'printim e:\1.tif x1000 y800 white'
;

hustlh

做个霸气的木头 发表于 2012-8-15 16:33

                               
登录/注册后可看大图

你的gs我修改了下 ，有断开的地反是修改了的，设定了相应的时间和层次，根据我上面说的你理解理解，应该没 ...

哦，我大致懂了。前2个变量和中间3个变量是时间维相同，层次维z不同。后3个变量和中间3个变量是层次维z相同，时间维不同。  所以每类都得重新设置时间t和层次z。   最后一段开头的‘set t 1’  ，就说明最后得出的图是这64年所有一月的平均，改动相应t从2到12，就可以分别得出其他2月到12月的。   是不是我理解的这样的？？？

还有你刚刚说最后一段用t的while循环就可以一次性得出12张图？？？

冰点

hustlh 发表于 2012-8-15 16:32

                               
登录/注册后可看大图

这个能够得出来。  按照你的方法我先不管64年的月平均，直接设置t从1变到12.即先设置set t 1   ，画1948年 ...

遇到了问题，才能进步啊！版主帮你修改的GS文件，你试下。

hustlh

冰点 发表于 2012-8-15 16:51

                               
登录/注册后可看大图

遇到了问题，才能进步啊！版主帮你修改的GS文件，你试下。

恩

正在试啊，真是感谢大家的帮助啊。

谢谢了！

做个霸气的木头

hustlh 发表于 2012-8-15 16:48

                               
登录/注册后可看大图

哦，我大致懂了。前2个变量和中间3个变量是时间维相同，层次维z不同。后3个变量和中间3个变量是层次维z相 ...

总算理解对了，循环的东西你试试吧，不难的！

冰点

做个霸气的木头 发表于 2012-8-15 17:00

                               
登录/注册后可看大图

总算理解对了，循环的东西你试试吧，不难的！

最近版主运气不怎么好啊！

hustlh

真的可以呢，每个月都不一样。

太感谢木头版主和冰点了，