关于NECP1°×1°资料的地形使用问题

飞雪残剑

经过了一天的琢磨终于搞明白了如何在grads绘制图形是增加地形，例如天山山区的地形，具体GS如下：
'reinit'
'open E:\20141208sn\fnl_20141208sn.ctl'
'open e:\topo1.ctl'
*设置绘图区域
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.0 10.0 1.5 8.0'
*先绘制地形图
dtop_horzional()
*开始绘制风场
'set dfile 1'
'set lat 34 50'
'set lon 72 98'
'set lev 700'
'set t 12'
'set xlint 5'
'set ylint 100'
'set gxout barb'
'd UGRDprs;VGRDprs'
'set mpdset xinjiang'  
'set map 14 1 7'
'draw map'
*标记ZWWW的坐标
st_lonlat= '43.91 87.47 ZWWW '
'q w2xy '%subwrd(st_lonlat,2)%' '%subwrd(st_lonlat,1)%''
xms1=subwrd(result,3)
yms1=subwrd(result,6)
'draw mark 5 '%xms1%' '%yms1%' 0.14'
'draw string '%xms1%' '%yms1%' '%subwrd(st_lonlat,3)%''

'printim e:\topo1.gif GIF white x1000 y800'
*以下为函数
*horzional位水平地形
function dtop_horzional()
    'set dfile 2'
    'set lat 34 50'
    'set lon 72 98'
    'set lev 850'
    'set t 1'
    'set xlint 5'
    'set ylint 5'
    'set clopts -1 -1 0.25'
    'set xlopts 1 -1 0.25'
    'set ylopts 1 -1 0.25'
    'set cmin 1000'
    'set grads off'  
    'set grid off'  
    'run define_colors.gs'
    'set clevs  700 900 1100 1300 1500 2000 2500 3000'
    'set ccols  0 71 72 73 74 75 76 77 78'
    'set gxout shaded'
    'd smth9(top)'
    'cbar 3 1'
    return
*vertical为垂直地形
function dtop_vertical()
    'set dfile 2'
    'set t 1'
    'set lat 10 60'
    'set lon 100'
    'set lev 1000'
    'set xlint 5'
    'set ylint 100'
    'define L1=smth9(ave(top,lon=100,lon=100))'
    'define L2=top*0.0'
    'set vrange 0 8000'
   'set ylab off'
    'set gxout linefill'
    'set lfcols 1 2'
   'd L1;L2'
  return
绘制图像如下：

WZZ809525579

这售价也太高了吧，望而却步

nunu18

楼主是新疆气象局的吗？我想知道楼主是怎么确定'set vrange 0 8000'中的8000，而不是15000？我也画过地形，我做的的是南疆西部，设置的是16000，显示的地形在500hPa以下，我是自己灵活调的，完了好多人都没搞明白