候平均

huangyao

用日平均资料画计算候平均，数据是2014年一月到2016年4月的，出来的图没有平均效果，只是每一个候第一天的值，还提示

'reinit'
'open d:\bb\14167\1000\ave.ctl'
'set lon 0 359'
'set lat -90 90'
'set lev 1000'
'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\1000\tu2\tu'pen'.'mon'.'year'.png white x1000 y600'
'c'
   pen=pen+1
endwhile
mon=mon+1
endwhile
year=year+1
endwhile
;