Friday, October 2, 2009

Climatology of season JJAS (grads script)

'reinit'
prompt 'Enter the name of new folder : '
pull fld
'!mkdir 'fld''
'open cmip1.ctl'
'q file'
'set lat -20 50'
'set lon 40 130'
'set time 1jun2022 30sep2101'
'define jjas=ave(prate, time=1jun, time=30sep)'
'set gxout shaded'
'd jjas'
'cbarn'
'define clim1=ave(jjas, time=1jun2022, time=30sep2101)'
'set time 1jun2022'
'd clim1*86400'
'set gxout shaded'
'd clim1*86400'
'cbarn'
'draw title cmip1 clomatology of JJAS'
'print climatology.gif'
'print 'fld'/clim1.eps'
'print 'fld'/clim1.gif'
'print 'fld'/clim1.jpg'
'!cd 'fld''
'!ls'