'open /cpc/home/wd52vk/samer/precip/oi-0.5deg-precip-SA-mon-y1979-jun2012.ctl'
'open /cpc/home/wd52vk/samer/precip/oi-0.5deg-precip-SA-mon-clim-y1981-2010.ctl'
'enable print fig1'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 1.0 3.0 7.5 10.0'
'set grads off'
'set mpdset hires'
'set annot 1 5'
'set map 1 1 1'
'set cterp off'
'set csmooth on'
'set string 1 tc 5'
'set strsiz 0.1 0.1'
'draw string 3.25 10.4 Precip. Anoms. (mm)'
'draw string 2.0 10.1 NDJ 2011-2012'
'set lat -35 15'
'set lon 275 330'
'set dfile 2'
'set t 1 12'
'define pclim=pcli.2'
'modify pclim seasonal'
'set dfile 1'
'run /cpc/home/wd52vk/grads/util/rgbset.gs'
'set clevs -400 -300 -200 -100 -50 0 50 100 200 300 400'
'set ccols 21 79 77 75 73 71 31 33 35 37 39 49'
'set xlint 20'
'set ylint 10'
'set gxout shaded'
'set grid off'
'set time dec2011'
'define pseacli=sum(pclim,t-1,t+1)'
'define pseatot=sum(p.1,t-1,t+1)'
'd pseatot-pseacli'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 3.5 5.5 7.5 10.0'
'set grads off'
'set mpdset hires'
'set annot 1 5'
'set map 1 1 1'
'set cterp off'
'set csmooth on'
'set string 1 tc 5'
'set strsiz 0.1 0.1'
'draw string 4.5 10.1 FMA 2012'
'set clevs -400 -300 -200 -100 -50 0 50 100 200 300 400'
'set ccols 21 79 77 75 73 71 31 33 35 37 39 49'
'set xlint 20'
'set ylint 10'
'set gxout shaded'
'set grid off'
'set time mar2012'
'define pseacli=sum(pclim,t-1,t+1)'
'define pseatot=sum(p.1,t-1,t+1)'
'd pseatot-pseacli'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'printim precipa-ndj-fma2011-2012.gif white'
'printim fig1.gif white'