'!mkdir rain7'
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)'
'd r'
'define mean=ave(r,time=01Oct'yyyy',time=31Dec'yyyy')'
'define dev=((r-mean)*(r-mean))'
'define rmsdev=sqrt(ave(dev,time=01Oct'yyyy',time=31Dec'yyyy'))'
'd mean-rmsdev'
'd mean'
'd mean+rmsdev'
'draw title Rain 'yyyy' '
'printim rain7/rain'yyyy'.jpg white'
yyyy=yyyy+1
endwhile
Thursday, February 5, 2009
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment