'reinit'
'open /cpc/analysis/verif/sat/cld_rad/olr/mon/olrmonth.ctl'
'enable print fig3'
'set vpage 0 8.5 0 11'
'set parea 1.0 7.5 7.5 10.0'
'set grads off'
'set looping off' 
'set annot 1 5'
'set cterp off'
'set csmooth on'
'set map 1 1 5'
'set string 1 tc 5'
'set strsiz 0.1 0.1'
*'draw string 4.25 10.5 FMA 2012'
*'draw string 4.25 10.25 OLR Anomalies'
'draw string 1.6 10.2 (a)'
'set lat -30 30'
'set lon 200 360'
'set t 1 12'
'define olrcli = ave(olr,t+84,t=444,12)'
'modify olrcli seasonal'
'set gxout shaded'
'run /cpc/home/wd52vk/grads/util/rgbset.gs'
'set clevs -20 -15 -10 -5 5 10 15 20'
'set ccols 38 36 34 32 0 72 74 76 78'
'set dfile 1'
'set time mar2012'
'define olra = ave(olr-olrcli,t-1,t+1)'
'set ylint 10'
'd olra'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'open /cpc/analysis/verif/ocean/sst/oi/ctl/monoiv2.ctl'
'open /cpc/analysis/verif/ocean/sst/oi/ctl/aoi.71.00.clim.ctl'
'set vpage 0.0 8.5 0 11.0'
'set parea 1.0 7.5 4.0 6.5'
'set grads off'
'set annot 1 5'
'set dfile 3'
'set t 1 12'
'define sstc=sst'
'modify sstc seasonal'
'set dfile 2'
'set time mar2012'
'set lon 200 360'
'set lat -30 30'
*'set ylopts 1 3 .1'
*'set xlopts 1 3 .1'
'set map 1 1 5'
'set csmooth on'
'set gxout shaded'
'set ccols 59 49 47 45 43 41 21 23 25 27 29 79'
'set clevs -3.0 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 3.0'
'set xlint 20'
'set ylint 10'
'd ave(sst.2-sstc,t-1,t+1)'
'set string 1 tc 5'
'set strsiz 0.1 0.1'
'draw string 1.6 6.7 (b)'
*'draw string 4.25 6.75 SST Anomalies (`ao`nC)'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'open /cpc/analysis/cdas/month/prs/mean/prs.grib.mean.clim.y1981-2010.ctl'
'open /cpc/analysis/cdas/month/prs/mean/prs.grib.mean.y1979-cur.ctl'
'set vpage 0 8.5 0 11'
'set parea 1.5 7.0 0.5 3.0'
'set grads off'
'set annot 1 5'
'set looping off'
'set cterp off'
'set csmooth on'
'set map 1 1 5'
'set dfile 1'
'set lat -15 5'
'set lon 180 365'
'define rad = 3.14159/180.'
'define dx = cdiff(lon,x)*rad*6.37e6'
'set dfile 4'
'set t 1 12'
'set z 1 12'
'define uc = ugrd'
'define vc = vgrd'
'define divc = hdivg(uc,vc)'
'define wc = VVEL*864'
'define chic = VPOT'
'modify uc seasonal'
'modify vc seasonal'
'modify divc seasonal'
'modify wc seasonal'
'modify chic seasonal'
'set lat -5'
'define avedivc = ave(divc,lat=-10,lat=0)'
'define avewc = -ave(wc,lat=-10,lat=0)'
'define avechic = ave(chic,lat=-10,lat=0)'
'define uchic = cdiff(avechic,x)/dx'
'modify avedivc seasonal'
'modify avewc seasonal'
'modify avechic seasonal'
'modify uchic seasonal'
'set t 3'
'define adivc = ave(avedivc,t-1,t+1)'
'define auchic = ave(uchic,t-1,t+1)'
'define awc = ave(avewc,t-1,t+1)'
'set dfile 5'
'set time feb2012 apr2012'
'set lat -15 5'
'set lon 180 365'
'define rad = 3.14159/180.'
'define dx = cdiff(lon,x)*rad*6.37e6'
'define u = ugrd'
'define v = vgrd'
'define div = hdivg(u,v)'
'define w = VVEL*864'
'define chi = VPOT'
'set lat -5'
'define avediv = ave(div,lat=-10,lat=0)'
'define avew = -ave(w,lat=-10,lat=0)'
'define avechi = ave(chi,lat=-10,lat=0)'
'define uchi = cdiff(avechi,x)/dx'
'set time mar2012'
'define adiv = ave(avediv,t-1,t+1)'
'define auchi = ave(uchi,t-1,t+1)'
'define aw = ave(avew,t-1,t+1)'
'set string 1 tc 5'
'set strsiz 0.1 0.1'
*'draw string 4.25 3.25 Anomalous Divergent Circulation EQ-10S'
'draw string 1.6 3.2 (c)'
'set zlog on'
'set gxout shaded'
'set clevs -4 -3 -2 -1 -0.5 0 0.5 1.0 2 3 4'
'set ccols 59 49 47 45 43 41 21 23 25 27 29 79'
'set lon 200 360'
'd smth9((adiv-adivc)*1000000)'
'set gxout contour'
'set ccolor 1'
'set cthick 2'
'set clevs -7 -6 -5 -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 5 6 7'
'd smth9((adiv-adivc)*1000000)'
'set gxout vector'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set cthick 5'
'set ccolor 1'
'set arrlab off'
'd skip(auchi-auchic,2);skip((aw-awc)/10,2)'
'run /cpc/home/wd52vk/grads/util/cbarhov.gs'