水汽通量输送 画图出现问题 求指导

学长吉他

我用grads想画出 1961-2010时间段内指定的几个年份 5-9月份水汽输送通量的距平，编写程序如下;
'reinit'
'sdfopen E:\OpenGrADS\data\shum.mon.mean.nc'
'sdfopen E:\OpenGrADS\data\uwnd.mon.mean.nc'
'sdfopen E:\OpenGrADS\data\vwnd.mon.mean.nc'
'sdfopen E:\OpenGrADS\data\pres.sfc.mon.mean.nc'
'set grads off'
'set grid off'

'set lon 40 160'
'set lat 0 70'
'set map  15'
*'set gxout fwrite'

'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'define suq=vint(pres.4,uq,300)'
'define svq=vint(pres.4,vq,300)'

'set t 1 12'
'have=ave(suq,t+156,t=756,12)'
'modify have seasonal'
'set t 157 756'
'hano=suq-have'
'set t 1'

'define suq1=ave(hano,t=209,t=213)'
'define suq2=ave(hano,t=245,t=249)'
'define suq3=ave(hano,t=413,t=419)'
'define suq4=ave(hano,t=425,t=429)'
'define suq5=ave(hano,t=605,t=609)'
'define AAu=(suq1+suq2+suq3+suq4+suq5)/5'

'set t 1 12'
'have=ave(svq,t+156,t=756,12)'
'modify have seasonal'
'set t 157 756'
'hano=svq-have'
'set t 1'
'define svq1=ave(hano,t=209,t=213)'
'define svq2=ave(hano,t=245,t=249)'
'define svq3=ave(hano,t=413,t=419)'
'define svq4=ave(hano,t=425,t=429)'
'define svq5=ave(hano,t=605,t=609)'
'define AAv=(svq1+svq2+svq3+svq4+svq5)/5'


'd AAu;AAv'
'printim E:\pic\data.png white'
;
图画不出来，

不知道问题出在哪？高手可否解答一下？

lxlxllx89

'set t 157 756'
'hano=suq-have'

have只有1-12 你t=157 不beyond file limits 就怪了

学长吉他

lxlxllx89 发表于 2014-1-20 20:27

                               
登录/注册后可看大图

'set t 157 756'
'hano=suq-have'

那该怎么办呢？