Friday, February 6, 2009

penta-day running mean (GrADS)

'!mkdir rain9'
yyyy=1951
while(yyyy<=2003)

'reinit'
'open rf_5103.ctl'
'set lat 10'
'set lon 75'

'set time 01Oct'yyyy' 31Dec'yyyy
'set grads off'
'define r=aave(rf,lon=70,lon=80,lat=8,lat=13)'
'define pent=ave(r,t-2,t+2)'
'd pent'
'define mean=ave(r,time=01Oct'yyyy',time=31Dec'yyyy')'
'd mean'
'draw title Rain 'yyyy' '
'printim rain9/rain'yyyy'.jpg white'

yyyy=yyyy+1
endwhile

No comments: