'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 marengo'
'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 7'
'set strsiz 0.1 0.1'
'draw string 3.25 10.5 Precip. Anoms. (mm)'
'draw string 2.0 10.2 DJF 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 jan2012'
'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 7'
'set strsiz 0.1 0.1'
'draw string 4.5 10.2  MAM 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 apr2012'
'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-djf-mam2011-2012.gif white'