求助，grads画站点6小时降水时图形出现大量空白

我喜欢风

rt。我在用站点数据画图时，想画6小时实况降水分布填色图，用了一个循环来实现多时次的操作，但是输出结果中有很多时次的图是空白的，数据均为缺省。但是按照每小时降水的状况，不应该出现这种情况的。。求各位大神救助
下面是我的gs。
'reinit'
'open grid.ctl'
'open 071808_2108.ctl'
'enable print 222r6.gmf'
'set mpdset cnworld'
'set grads off'
'set lon 116 125'
'set lat 22 30'
t=6
****t=26
while(t<=73)
'set t 't
'q time'
res=subwrd(result,3)
'set cmin 3'
'define a=oacres(g.1,rn.2)'
*定义一个i，用来在计算6小时降水时做中间变量
i=t-1
while(i>t-6)
'set t 'i
'define b=oacres(g.1,rn.2)'
'define a=a+b'
i=i-1
endwhile
'set gxout shaded'
'set csmooth on'
'd a'
'cbarn'
'draw title \r6_'res
'print'
'c'
t=t+1
endwhile
'disable print'
'reinit'
;

我喜欢风

另外。这个是站点数据ctl。
DSET 071808_2108.dat
DTYPE station
STNMAP 071808_2108.map
UNDEF -9999.0
TITLE zdz
TDEF 73 LINEAR 00Z18Jul2005 1hr
VARS 6
tmc 0 99 temperature
slp 0 99 sea level pressure
rh  0 99 relative humidity
U 0 99 u componet of wind in 2 min
V 0 99 v componet of wind in 2 min
rn 0 99 rainfall
ENDVARS     
格点背景场的ctl
DSET ^grid.dat
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 131 LINEAR 112 0.1
YDEF 91 LINEAR 21 0.1
ZDEF 1 linear 1 1
TDEF 73 LINEAR 00Z18JUL2005 1hr
VARS 1
g 0 99 Grid Data
ENDVARS

我喜欢风

纠结了好久啊=-=悲了个剧啊

mofangbao

格点场的ctl呢

我喜欢风

格点背景场的ctl
DSET ^grid.dat
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 131 LINEAR 112 0.1
YDEF 91 LINEAR 21 0.1
ZDEF 1 linear 1 1
TDEF 73 LINEAR 00Z18JUL2005 1hr
VARS 1
g 0 99 Grid Data
ENDVARS

我喜欢风

清风大大、土问你作图的时候出现过这种情况么。。

我喜欢风

然后又这么个状况。当我不用循环，注释掉while和endwhile，出来的图形又不是空白的了。。比如我定t=13时图出来了。但是在没注释时。这个时次的图形是空白的。。
土问这是什么状况。。

做个霸气的木头

你这双重循环就是为了求每个时刻开始的6个小时的累计量，其实只要一个循环就行。你这个问题还得自己好好检查，要不楼主将主体求和，试试改成下面这样：'set lon 116 125'
'set lat 22 30'

‘set t 1 73’
''define a=oacres(g.1(t=1),rn.2)''
i=6
while(i<=73)
j=i-5
'set t 1'
'define aa=ave(a,t='j',t='i')'
'define a2=aa*6'
i=i+1
endwhile

凝露灵动

能不能请问一下楼主 你的6小时降水资料是在哪里下载的 不甚感激