24小时热门版块排行榜    

查看: 564  |  回复: 4
当前主题已经存档。

liuzf406

木虫 (著名写手)

[交流] 【求助】grads计算绘图

请各位虫友帮助:我想利用GRADS计算1993年美国及其周边地区冬夏季1000-500百帕的水汽通量,其中夏季的gs文件如下,但运行结果不对,请高人指点,谢谢:
'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen shum.mon.mean.nc'  
'sdfopen uwnd.mon.mean.nc'
'sdfopen vwnd.mon.mean.nc'
'sdfopen pres.mon.mean.nc'
'enable print 10.gmf'
'set lat 10 60'
'set xlint 10'
'set lon -150 -60'
'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 gxout shaded'
'set ylopts 1 4 0.18'
'set xlopts 1 4 0.18'
'set cthick 10'
'define m=shum.1*uwnd.2'
'define n=shum.1*vwnd.3'
'define qu700=vint(pres.4(z=1),m,700)/1000'
'define qv700=vint(pres.4(z=1),n,700)/1000'
'd ave(mag(qv700,qu700),t=546,t=547)'
'set gxout contour'
'set cint 20'
'd ave(mag(qv700,qu700), t=546,t=547)'
'set gxout vector'
'print'
'disable print'
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyong0210

木虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
'sdfopen pres.mon.mean.nc'
是地面气压,是一层吗
pres.4(z=1),z=1是地面层,还是1000hpa层啊???如果是1000hpa层,就错了。好像要地面气压层,不是1000hpa层的气压
王者无敌
2楼2009-08-12 23:51:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄胖胖

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
把运行错误提示贴上来看看
3楼2009-08-13 10:59:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuzf406

木虫 (著名写手)

引用回帖:
Originally posted by wyong0210 at 2009-8-12 23:51:
'sdfopen pres.mon.mean.nc'
是地面气压,是一层吗
pres.4(z=1),z=1是地面层,还是1000hpa层啊???如果是1000hpa层,就错了。好像要地面气压层,不是1000hpa层的气压

谢谢啊,是地面层.
4楼2009-08-13 21:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuzf406

木虫 (著名写手)

引用回帖:
Originally posted by 黄胖胖 at 2009-8-13 10:59:
把运行错误提示贴上来看看

可能是时间有问题,我加了'set t 12', 可以了.
谢谢啊!
5楼2009-08-13 21:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuzf406 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见