ERA-I的资料是每天四个时次，求日平均的脚本有错吗？两种方法算出来画的图...

jolincai

NC资料里是一天四个时次，我把每天的四个时次求其日平均，出来了GRD，但是当我用同样算法求第N天的日平均，并与同日刚才提取的GRD画图做对比，两个图就不一样了。不知道为什么，我的求平均方法有错吗？资料是2年的，共731天，2924个时次。
'reinit'
'sdfopen E:\data\erai\height\1995_1996erai500h.nc'
'set fwrite E:\data\erai\height\1995_1996erai500hdaily.grd'
'set gxout fwrite'
'set lon -180 180'
'set lat -90 90'
i=1
while (i<=731)
'define rpj=ave(z,t=4*i-3,t=4*i)'
'd rpj'
i=i+1
endwhile
'disable fwrite'
*'define rpj=ave(z,t=5,t=8)'
*'d rpj'
*'disable fwrite'
;

星号部分可以直接求某一天的日平均，我再与之前的做对比，两张图就不一样了，求指点迷津。