grads

huangyao

我用日平均资料求候平均，ctl和gs文件如下，出来的图老是不对，感觉循环有问题，但是不知道怎么改，求各位大神指点
ctl描述：
dset d:\bb\14167\ave.dat
title bobao
undef -2.56E33
xdef 360 linear 0.0 1.0
ydef 181 linear -90.0 1.0
zdef 1 linear 850 1
tdef 851 linear 01JAN2014 1dy
vars 1
yy 0 99 bobao
endvars

我的gs文件：
'reinit'
'open d:\bb\14167\85h.ctl'
'set lon 0 360'
'set lat -90 90'
'set lev 850'
'set t 1 851'
'd yy'
day=1
year=2014
while(year<=2016)
mon=1
while(mon<=12)
   pen=1
   while(pen<=6)
    if(mon=1|mon=3|mon=5|mon=7|mon=8|mon=10|mon=12)
     if(pen=6)

      'd ave(yy,t='day',t='day'+5)'  
      day=day+6
else
'd ave(yy,t='day',t='day'+4)'
   day=day+5
     endif
endif
    if(mon=2)
     if(pen=6)

      rc=math_fmod(year,4)
      rc1=math_fmod(year,100)
      rc2=math_fmod(year,400)
      if((rc=0&rc1!=0)|(rc2=0))
         'd ave(yy,t='day',t='day'+3)'
          day=day+4
      else
         'd ave(yy,t='day',t='day'+2)'
         day=day+3
      endif
else
'd ave(yy,t='day',t='day'+4)'
   day=day+5
endif
endif
if(mon=4|mon=6|mon=9|mon=11)
   'd ave(yy,t='day',t='day'+4)'
   day=day+5
   endif
'printim d:\bb\14167\qq.png white'
'c'
   pen=pen+1
endwhile
mon=mon+1
endwhile
endwhile
;