问题求助

有个依靠

我下载了1948 至今的slp月平均资料，先提取出1970至2013年12，1,2月的冬季资料
'reinit'
'sdfopen f:\zl\slp.mon.mean.nc'
'set fwrite f:\zl\zliao1.grd'
'set gxout fwrite'
'set lat 20 65'
'set lon 60 240'
tt=264
t=1
while (t<=43)
'set t 'tt+12''
'd slp'
'set t 'tt+13''
'd slp'
'set t 'tt+14''
'd slp'
t=t+1
tt=tt+12
endwhile
'disable fwrite'
;
然后写描述文件

有个依靠

dset f:\zl\zliao1.grd
undef -9.99E+33
title NCEP/NCAR REANALYSIS  PROJECT
xdef 72 linear 60.0 2.5
ydef 18 linear 20.0 2.5
zdef 1 levels 995
tdef 129 linear DEC1970 1mon
vars 1
slp 1 99 presure hpa
endvars

有个依靠

然后想画一个slp多年平均图
'reinit'
'open f:\zl\miaoshu.ctl'
'enable print f:\zl\duonianpingjun1234.gmf'
'set mpdset mres'
'set ccolor rainbow'
'set cint 1'
'set grads off'
'set cterp on'
'set clskip 2'
'set lon 60 240'
'set lat 20 65'
'set t 1 129'
i=1
while(i<=127)
'define have=ave(slp,t='i',t='i+2')'
i=i+3
endwhile
'set t 1 43'
'define hac=ave(have,t=1,t=43)'
'd hac'
'draw title page1 1970-2013 winter slp annual average'
'print'
'disable print'
;

有个依靠

画出来是这样，哪里不对？

有个依靠

想请求大家的帮助

有个依靠

按理说画出的图可以清楚的看出冬季蒙古高压和阿留申低压的分布