Saturday, January 31, 2009

Northeast monsoon area average (GrADS)

prompt 'Enter south boundary: '
pull s
prompt 'Enter north boundry: '
pull n
prompt 'Enter west boundary: '
pull w
prompt 'Enter east boundary: '
pull e

'!mkdir 'rain'_'s'_'n'_'w'_'e

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='w',lon='e',lat='s',lat='n')'
'd r'
'draw title Rain area ave 'w'e 'e'e 's'n 'n'n'
'printim rain_'s'_'n'_'w'_'e'/rain'yyyy'.jpg white'

yyyy=yyyy+1
endwhile

No comments: