* This script writes the monthly precipitation anomalies for SRES_A2 runs. 
* anomalies relative to base period 01-60 or 61-90.

 'reinit'
 'xdfopen /large/obs2/cmip5/mo/cnrm/rcp8.5/pr_1850-2100.ctl'

* define clim precc in mm/day
 'set time jan1901 dec1960'
 'define prc=ave(pr,t+0,time=dec1960,1yr)'
 'modify prc seasonal'

 'set time jan1850 dec2100'
 'define pra=(pr-prc)*86400'

* write out jja precip. anom for 21st century
 'set gxout fwrite sequential'
 'set fwrite cnrm_jja_pra_1850-2100_0160.out'

 n=7
 m=1
 while (n<3013)
 'set t 'n''
   'define pajja'm'=ave(pra,t-1,t+1)'
 'set x 1 256'
 'set y 1 128'
 'set z 1'
 'set t 'm''
   'd pajja'm''
  m=m+1
  n=n+12
 endwhile

 'disable fwrite'


**********