Friday, January 23, 2009

Indian ocean Dipole using grads

yyyy=1960
while(yyyy<=1979)

'reinit'
'open sst6079.ctl'

'set lat 10'
'set lon 75'

'set time 01Oct'yyyy' 31Dec'yyyy
'set grads off'
'define sste=aave(sst,lon=90,lon=110,lat=-10,lat=0)'
'd sste'
'define sstw=aave(sst,lon=50,lon=70,lat=-10,lat=10)'
'd sstw'
'define sstd=sstw-sste'
'd sstd'
'draw title sst control run'
'printim sst1/sst'yyyy'.jpg white'

yyyy=yyyy+1
endwhile


yyyy=1980
while(yyyy<=1998)

'reinit'
'open sst8098.ctl'

'set lat 10'
'set lon 75'

'set time 01Oct'yyyy' 31Dec'yyyy
'set grads off'
'define sste=aave(sst.2,lon=90,lon=110,lat=-10,lat=0)'
'd sste'
'define sstw=aave(sst.2,lon=50,lon=70,lat=-10,lat=10)'
'd sstw'
'define sstd=sstw-sste'
'd sstd'
'draw title sst control run'
'printim sst1/sst'yyyy'.jpg white'

yyyy=yyyy+1
endwhile

No comments: