用日资料求候平均

huangyao

gs文件：
'reinit'
'open d:\bb\14167\ave.ctl'
'set lon 0 360'
'set lat -90 90'
'set lev 850'
'set mpdset cnworld'
'set grid on'
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\tu1\tu3\tu'pen'.'mon'.'year'.png white x1000 y600'
'c'
pen=pen+1
endwhile
mon=mon+1
endwhile
year=year+1
endwhile
;

ctl文件：
dset d:\bb\14167\ave.dat
title bobao
undef -9.99E33
xdef 360 linear 0.0 1.0
ydef 181 linear -90.0 1.0
zdef 1 linear 1 1
tdef 851 linear 01JAN2014 1dy
vars 1
yy 0 99 bobao
endvars

robindct

'set lev 850' 这里是不是有点问题。在ctl里，z并没有850这层呐