地形叠加显示出错

酒国

想将这两张图合在一起，但总是出错。单独画是可以的。不知道gs里面有什么问题？
下面是gs:
'reinit'
'sdfopen d:/uwnd.mon.mean.nc'
'sdfopen d:/omega.mon.mean.nc'
'set grads off'
'set grid off'
'set lon 90 180'
'set lat 30'
'set z 1 12'
'set t 37'
'set csmooth on'
'set xlopts 1 4 0.14'
'set ylopts 1 4 0.14'
'set gxout stream'
'set cthick 6'
'set ccolor 1'
'set clopts 0 4 0.10'
'define u1=ave(uwnd,t=37,t=625,12)'
'define u2=ave(uwnd,t=38,t=626,12)'
'define u3=ave(uwnd,t=39,t=627,12)'
'define u4=ave(uwnd,t=40,t=628,12)'
'define u5=ave(uwnd,t=41,t=629,12)'
'define u=(u1+u2+u3+u4)/4'
'define o1=ave(omega.2,t=37,t=625,12)'
'define o2=ave(omega.2,t=38,t=626,12)'
'define o3=ave(omega.2,t=39,t=627,12)'
'define o4=ave(omega.2,t=40,t=628,12)'
'define o5=ave(omega.2,t=41,t=629,12)'
'define o=(o1+o2+o3+o4)/4'
'd u;-o*100'
'sdfopen d:/topo.nc'
'define pressure=1000-0.76*htopo/8'
'set lon 90 180'
'set gxout bar'
'set lat 30'
'set axlim 100 1000'
'set yflip on'
'set ccolor 1'
'd smth9(pressure)'
pull '     '
'!gxps -cr -i vor_lon_vertical.gmf -o vor_lon_vertical.ps'

'draw title 1000-100 u,w'
'print'
'disable print'
;

numa

看报错情况似乎是指向htopo这个变量。你确定这个变量存在么？

酒国

numa 发表于 2017-9-25 12:09
看报错情况似乎是指向htopo这个变量。你确定这个变量存在么？

存在的。。。。。单独画都是可以的。。。。。和在一块就说定义的名字有问题

numa

酒国 发表于 2017-9-25 12:12
存在的。。。。。单独画都是可以的。。。。。和在一块就说定义的名字有问题

不妨把脚本中涉及的文件上传一下让大家测试一下。

柚子皮

'sdfopen d:/uwnd.mon.mean.nc'
'sdfopen d:/omega.mon.mean.nc'

看文件有两个，u，v属于第一个，htopo属于第二个是吗？那么使用的时候htopo需要表明为“htopo.2”表示来自第二个文件

酒国

柚子皮 发表于 2017-9-26 12:17
'sdfopen d:/uwnd.mon.mean.nc'
'sdfopen d:/omega.mon.mean.nc'

这个问题我之前试过了，还是不对的。

酒国

numa 发表于 2017-9-25 12:33
不妨把脚本中涉及的文件上传一下让大家测试一下。

好啊。月平均的可以在网站下。地形的就是这个。