grads为什么不出等值线啊

甜狗狗

这是ctl文件：
dset e:\right.grd
dtype station
stnmap e:\ss.map
undef 99999.9
*title the 160 station Jan t and ind rc
tdef 1 linear jan1982 1mo
vars 1
r 0 99 winter rain anomaly
endvars
这是gs文件：
'reinit'
'open e:\map.ctl'
'open e:\r160.ctl'
'enable print e:\right1.gmf'

'set grads off'
'set grid off'
*'set mpdset lowres'
*'set poli off'
'set mpdset sichuan chongqing guizhou yunnan'
'set lat 20 35'
'set  lon 95 111'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'define aa=oacres(g.1,r.2,4)'
'define rgg=maskout(aa,g.1-0.5)'
'define bb=smth9(rgg)'
*'set csmooth on'
*'set gxout shaded'
*'set clevs    -0.4 -0.3  0  0.3  0.4 '
*'set ccols   4 5 0  0  6  8 '
'set csmooth on'
'd bb'
'set clab forced'
'set gxout contour'
'set cthick 6'
'set cint 0.1'
'd bb'
'print'
'disable print'
;
求大神指导

river

有错误提示吗？

甜狗狗

river 发表于 2013-3-20 21:14

                               
登录/注册后可看大图

有错误提示吗？

我最后搞出来了，是ctl文件写的有问题，谢谢啦