'reinit'
'sdfopen /stornext/online11/eta/Adan_Silva/Observacoes/CRU/0.5x0.5/cru_ts3.22.1901.2013.pre.dat.nc'

'set display color white'
'c'

'set lat -50 27.8'
'set lon -100 -29.2'
'set font 4'

'set mpdset brmap_hires'

* Ano inicial e final da media sazonal
_anoi=1981
_anof=1990

'jan=ave(pre,time=JAN'%_anoi%',time=JAN'%_anof%',12)/31'
say result                                         
'feb=ave(pre,time=FEB'%_anoi%',time=FEB'%_anof%',12)/28'
say result                                         
'mar=ave(pre,time=MAR'%_anoi%',time=MAR'%_anof%',12)/31'
say result                                         
'apr=ave(pre,time=APR'%_anoi%',time=APR'%_anof%',12)/30'
say result                                         
'may=ave(pre,time=MAY'%_anoi%',time=MAY'%_anof%',12)/31'
say result                                         
'jun=ave(pre,time=JUN'%_anoi%',time=JUN'%_anof%',12)/30'
say result                                         
'jul=ave(pre,time=JUL'%_anoi%',time=JUL'%_anof%',12)/31'
say result                           
'aug=ave(pre,time=AUG'%_anoi%',time=AUG'%_anof%',12)/31'
say result                           
'sep=ave(pre,time=SEP'%_anoi%',time=SEP'%_anof%',12)/30'
say result                           
'oct=ave(pre,time=OCT'%_anoi%',time=OCT'%_anof%',12)/31'
say result
'nov=ave(pre,time=NOV'%_anoi%',time=NOV'%_anof%',12)/30'
say result
'dec=ave(pre,time=DEC'%_anoi%',time=DEC'%_anof%',12)/31'
say result

* -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - -- - --


tri=1
while(tri<=4)

if (tri=1)
lag1='DJF'
'define precz'tri'=(dec+jan+feb)/3'
endif

if (tri=2)
lag1='MAM'
'define precz'tri'=(mar+apr+may)/3'
endif

if (tri=3) 
lag1='JJA'
'define precz'tri'=(jun+jul+aug)/3'
endif

if (tri=4) 
lag1='SON'
'define precz'tri'=(sep+oct+nov)/3'
endif

'set rgb 40  30  35 244' 
'set rgb 46  84  98 255'
'set rgb 50 120 140 255'
'set rgb 55 162 189 255'
'set rgb 59 198 231 255'

'set grads off'
'set gxout shaded'
'set clevs 1 3 6 9 12'
'set ccols 0 59 55 50 46 40' 

'set clopts 1 12 0.15'
'set xlopts 1 12 0.15'
'set ylopts 1 12 0.15'
'd precz'tri''

'set gxout contour'
'set cthick 2'
'set clab off'
'set clevs 1 3 6 9 12'
'd precz'tri''
'run cbarn_Left.gs 1 1 1.8 4'

'set string 1 c 10'
'set strsiz .15 .14'

'draw string 5.5 8.4  CRU 50km 'lag1''
'draw string 5.5 8.0  Total Precipitation (mm)'
'gxprint PREC_50km_CRU_'lag1'.png png white x1024 y768'


'c'

tri=tri+1
endwhile

'quit'