'open data/camsopi.ctl'
'set display color white'
'clear'
'set lat -40 40'
'set lon -20 55'
'set t 1 12'
'set gxout shaded'
*---------------------------
* 1961-1990 clim
'set t 1 12'
t=25
'define pclim=ave(comb,t+0,time=dec2010,12)'
'modify pclim seasonal'
'define pcltot=365*ave(pclim,t=1,t=12)'
'set t 1'
'set csmooth on'
'set cterp on'
'define tot2015=365*ave(comb,time=jan2015,time=dec2015)'
*
'set t 1'
'define pcncur=(tot2015/pcltot)*100'
'set parea 0 8.5 2 10'
'set mproj scaled'
'set grid off'
'set grads off'
'set mpdset hires'
'set gxout shaded'
'exec colors'
'set csmooth on'
'set cterp on'
'set clevs 50 75 125 150'
'set ccols 27 24 0 34 38'
* Filter regions (with annual precip < 100mm) out for correction 
'define ones = const(const(maskout(pcltot,pcltot-100),1),0,-u)'
'define pcnc= maskout(pcncur,ones-1)'
'd pcnc'
'run cbarhov2.gs'
'printim camsopi_precip_pcnp_2015.png'
*
*
'set fwrite camsopi_precip_pcnp_2015.dat'
'set gxout fwrite'
'set undef -999.0'
'd pcnc'
*
'set mpdraw off'
'set geotiff camsopi_precip_pcnp_2015'
'set gxout geotiff'
'set mproj latlon'
'd pcnc'
*
'quit'