wrf输出求降水的时间-纬度图（求改错）

alxnumb

       首先祝大家周末愉快，感谢看帖人，更感谢回帖好心yin！
       本人想编写一个简单的wrf输出rainc+rainnc在117°E的时间-纬度图，自己编了一个gs出现很多问题，无法画出想要的图，求各位高手指点。以下为我编写的gs

reinit
open Y:\wrfout.ctl
set xyrev on
set lat 10 40
set t 1 80
define aa=0
while(i<80)
aa=rainnc(t='i'+1)+rainc(t='i'+1)-rainnc(t='i')-rainc(t='i')
i=i+1
d aa
endwhile
;

Life

把出现的错误发上来，也许大家能更好的帮助你~~

alxnumb

grads错误为下图

alxnumb

Life 发表于 2012-12-15 15:34

                               
登录/注册后可看大图

把出现的错误发上来，也许大家能更好的帮助你~~

图已经贴下面了

river

你想画时间-纬度图，他和时间--高度、纬度--高度、经度--高度这些剖面图是一个道理的。你现在最基本的剖面画法都不懂呢，先看看基础的东西吧。然后再看看剖面图的画法http://bbs.06climate.com/forum.php?mod=viewthread&tid=10552。最主要还是自己多试一下。

你编写的和你想画的根本不是一回事儿。
'reinit'
'open y:\wrfout.ctl'
'set xyrev on'
'set lon 117'
'set lat 10 40'
'set t 1 80'
'd rainc+rainnc'
;

alxnumb

river 发表于 2012-12-15 19:31

                               
登录/注册后可看大图

你想画时间-纬度图，他和时间--高度、纬度--高度、经度--高度这些剖面图是一个道理的。你现在最基本的剖面画 ...

我已经画出来位势高度的时间-纬度图，但是降水的数据rainc和rainnc是一个累加的量，需要减去前一个时次才行。。所以我就被卡住了

river

alxnumb 发表于 2012-12-15 21:27

                               
登录/注册后可看大图

我已经画出来位势高度的时间-纬度图，但是降水的数据rainc和rainnc是一个累加的量，需要减去前一个时次才 ...

还是对grads不熟悉，你对于什么时候用define和什么时候不用define还不清楚呢，推荐你朱禾老师的GrADS绘图软件实用手册，里面涵盖了grads基础的命令和各种图的画法，只要边看边做一下就会对grads熟悉很多，grads其实还是很强大的。

至于你那个降水量的图这么画就行了：
'reinit'
'open Y:\wrfout.ctl'
'set t 1 80'
'define aa=rainnc(t+1)-rainnc+rainc(t+1)-rainc'
'set xyrev on'
'set lon 117'
'set lat 10 40'
'd aa'
;

alxnumb

river 发表于 2012-12-16 08:46

                               
登录/注册后可看大图

还是对grads不熟悉，你对于什么时候用define和什么时候不用define还不清楚呢，推荐你朱禾老师的GrADS绘图 ...

好的，非常感谢，我一定多加学习！

river

alxnumb 发表于 2012-12-16 20:07

                               
登录/注册后可看大图

好的，非常感谢，我一定多加学习！

互相交流也是一种很好的学习呢，多来论坛逛逛吧