垂直水汽通量时间-纬度剖面图

...该...

谢谢，river和jacci_zou的点拨了，图解决了。。。。太感谢。。。有问题还得一起讨论

river

...该... 发表于 2013-3-19 23:50

                               
登录/注册后可看大图

谢谢，river和jacci_zou的点拨了，图解决了。。。。太感谢。。。有问题还得一起讨论

你这个也没点个回复什么的，直接发在你的帖子底下给你自看啊
另外还有一个相关的例子可以学习下http://bbs.06climate.com/forum.p ... p;extra=&page=1

还有就是楼主可以把自己成功的gs贴上来让大家学习嘛

...该...

这里感谢river.下面是我的gs文档。gs 文档中有问题大家一起来交流
#计算5月到九月的时间-纬度的水汽通量散度###
'reinit'
'sdfopen f:/data/201205-09/uwnd.2012.nc'
'sdfopen f:/data/201205-09/vwnd.2012.nc'
'sdfopen f:/data/201205-09/air.2012.nc'
'sdfopen f:/data/201205-09/pres.sfc.2012.nc'
'sdfopen f:/data/201205-09/shum.2012.nc'
'set lon 50 145'
'set lat -30 70'
'set t 122 274'
#当初做不出图来是因为时间放错了位置，现在我个人觉得现在的图还让自己满意#
'set map 1 1 20'
'set annot 1 50'
'set font 4'
'set mpdset cnworld'
'define usq=vint(pres.4(z=1),shum.5*uwnd.1,300)'
'define vsq=vint(pres.4(z=1),shum.5*vwnd.2,300)'
'define qdiv=hdivg(usq,vsq)'

'set lon 1'
'define w1=ave(qdiv,lon=90,lon=120)'
'define w2=ave(uwnd.1,lon=90,lon=120)'
'define w3=ave(vwnd.2,lon=90,lon=120)'
'set lat -10 60'
'set grads off'
'set grid off'
'set xyrev on'
'set timelab off'
'set gxout shaded'
'd w1*1000'
#不知道这里*1000量级对不对？#
'cbarn'
'set xyrev on'
'set gxout vector'
'd w2;w3'

'printim f:\17.gif x4000 y3000 white'

...该...

补充一点，针对set lon 1的设置，我在论坛里面看到有的人提到这个1可任取一个数，结果都不变，我试了，结果真的是不变的

linnamei

{:soso_e100:}

river

...该... 发表于 2013-3-21 23:18

                               
登录/注册后可看大图

这里感谢river.下面是我的gs文档。gs 文档中有问题大家一起来交流
#计算5月到九月的时间-纬度的水汽通量散 ...

'd w1*1000'
#不知道这里*1000量级对不对？#
这个无所谓了，乘与不乘只是扩大的量级不同而已，最主要的是楼主要搞清楚物理量的单位，单位对了就对了，单位搞错了即使算出来的数字是对了也是错的。一般的话都是扩大1e5

...该...

谢谢啦。我知道了！学习了不少！

风荷1

谢谢分享，又学了一招

szx123

不错，学习了。

szx123

不错，学习了。