GrADS剖面加高分辨率率地形，只差一步！求指点。

晚晴山

***************************以下是GS，出的图不对齐，请指出错误之处***************************
'reinit'
*'open F:\lw\WRFcape\0512\wr.ctl'
'open F:\datasets\fnl\200505\fnl_20050512.ctl'
'sdfopen D:\etopo20-0-360.nc'
'enable print D:\2.gmf'
'set lat 33'
'set lon 80 110'
'set t 3'
'set lev 1000 100'
'set vrange 800 200'
'set grid off'
'set grads off'
'set gxout shaded'
'd VVELprs'
'cbarn 0.8 0'
*'set gxout contour'
*'d VVELprs.2'

'set dfile 2'
'set t 1'
'set lat 33'
'set lon 80 110'
'set z 1'
'define pre=1013.25*pow((1-ROSE/44330),5.255)'
'define l1=smth9(ave(pre,lat=29,lat=36))'
'define l2=1000-pre*0.0'
'set vrange 800 200'
'set ylab off'
'set gxout linefill'
'set lfcols 2 15'
'set yflip off'
'd l1;l2'

'print'   
'disable print'
;      

lqouc

你用vpage和parea设定一下

river

把第二个'set vrange 800 200'去掉试试

晚晴山

lqouc 发表于 2013-9-15 21:18

                               
登录/注册后可看大图

你用vpage和parea设定一下

注意你好久了，你是个靠谱的哥们，你是南信大的么？

晚晴山

river 发表于 2013-9-15 22:43

                               
登录/注册后可看大图

把第二个'set vrange 800 200'去掉试试

去掉之后地形就会倒过来

lqouc

晚晴山 发表于 2013-9-16 10:21

                               
登录/注册后可看大图

注意你好久了，你是个靠谱的哥们，你是南信大的么？

过奖啦，我不是南信大的。
其实吧，我说的方法也不一定有用，在grads的使用上面还是river比较在行。

noisy

请教楼主，高分辨率的地形资料是从哪里得到的，能否共享？

晚晴山

lqouc 发表于 2013-9-16 10:28

                               
登录/注册后可看大图

过奖啦，我不是南信大的。
其实吧，我说的方法也不一定有用，在grads的使用上面还是river比较在行。

我搞定了，和你说的差不过。

river

lqouc 发表于 2013-9-16 10:28

                               
登录/注册后可看大图

过奖啦，我不是南信大的。
其实吧，我说的方法也不一定有用，在grads的使用上面还是river比较在行。

过奖了，我也只是稍微懂一些啊

renzi1018

你好请问
‘set z 1'
'define pre=1013.25*pow((1-ROSE/44330),5.255)'
中ROSE是什么变量？