'reinit'
'open  MERGE_CPTEC_ciclo_1h.ctl'
'sdfopen Brazil.nc'
'br=lterp(band1.2(t=1,z=1),pmed.1(t=1))'

'open mask_reg.ctl'
'set datawarn off'
'set display color white'
'set xlopts 1 4 0.25'
'set ylopts 1 4 0.25'
'set rgb 22 0 100 0'
'set rgb 23 255 174 66'
'enable print oi.gmt'
'set map 3 1 8'
'set map 1 1 5'
'c'
plota()

dj=25
di=25
i=90
MMAX=0
cont=0
while(i<=482)
  j=240
  while(j<=630) 
    'set dfile 1' 
    'set x ' i
    'set y ' j
    'q dim'
    lat=sublin(result,3)
    lat=subwrd(lat,6)
    lon=sublin(result,2)
    lon=subwrd(lon,6)
*    say lat 
*    say lon
   'q w2xy 'lon ' 'lat
   x=subwrd(result,3)
   y=subwrd(result,6)


*'set dfile 3'
'set lat ' lat
'set lon ' lon
*'d R(t=1)'
co=sublin(result,9)
co=subwrd(co,4)

if (co=1)
 cor=4
endif

if (co=2)
 cor=22
endif

if (co=3)
 cor=1
endif

if (co=4)
 cor=23
endif

if (co=5)
 cor=2
endif

cor=2

'set parea  'x-0.15'  'x+0.15'  'y-0.25'  'y+0.25
'set grid off'
'set xlab off'
'set ylab off'
'set ccolor 22'
'set ccolor  ' cor

*'set baropts outline'
*'set bargap 50' 
'set annot 0'
'set gxout bar'
*if (lat<-20)
*  'set vrange 0.27948410975135425727 2.72737622917779359165  ' * eh dah normalizacao
*  'set vrange 0.2 3 '
*else
  'set vrange  0 0.7'
*endif
'set dfile 1'
'set t 1 24'
'set stat on'
*'d tloop(aave(pmed,x='i'-'di/2',x='i'+'di/2',y='j'-'dj/2',y='j'+'dj/2'))'
* 0.33003844147727272727 eh dah normalizacao
'd (tloop(aave(maskout(pmed,br),x='i'-'di/2',x='i'+'di/2',y='j'-'dj/2',y='j'+'dj/2')))'
a=sublin(result,8)
max=subwrd(a,5)
if(max>=0)
 MMAX=MMAX+max
 cont=cont+1
 say  max"  "MMAX" "cont
endif

*'set strsiz 0.1'
*'draw string  'x-0.21'  'y-0.31'  J  D'
*say result
'set annot 1'
plota()
   j=j+dj
 endwhile
i=i+di
endwhile

'quit
'open BRAZIL.ctl'
'set dfile 4'
'set gxout fgrid'
'set fgvals 0 0 -1 1'
'd img.4(t=1)'
'draw map'
'set strsiz 0.25 0.25'
'set cthick 8'
xi=7.3
yi=9.6 
dx=0.3
dy=0.3
'set line 4' 
'draw recf 'xi'  'yi'  'xi+dx'  'yi+dy
'draw string  'xi+1.5*dx'  'yi' R1' 

yi=yi-dy-0.05
'set line 22' 
'draw recf 'xi'  'yi'  'xi+dx'  'yi+dy
'draw string  'xi+1.5*dx'  'yi' R2' 

yi=yi-dy-0.05
'set line 1' 
'draw recf 'xi'  'yi'  'xi+dx'  'yi+dy
'draw string  'xi+1.5*dx'  'yi' R3' 

yi=yi-dy-0.05
'set line 23' 
'draw recf 'xi'  'yi'  'xi+dx'  'yi+dy
'draw string  'xi+1.5*dx'  'yi' R4' 

yi=yi-dy-0.05
'set line 2' 
'draw recf 'xi'  'yi'  'xi+dx'  'yi+dy
'draw string  'xi+1.5*dx'  'yi' R5' 

'set dfile 1'
i=140
j=170
'set parea  5.55 8.25 1.25 2'    
'set x 140'
'set y 170'
'set t 0 13'
'set vrange  -0.5 10'
'set ccolor 1'
'd tloop(aave(maskout(prec,-br.2(t=1)),x='i'-'di/2',x='i'+'di/2',y='j'-'dj/2',y='j'+'dj/2'))'
'set strsiz 0.11 0.11'
'draw string  5.70439 1.08642  J F M A M J J A S O N D'
'printim clim_prec2.png'
'quit'

************* funcao que apenas plota o mapa do brasil ******************
function plota()
'set dfile 1'
'set mpdset brmap_hires'
'set mpdset hires'
'set grads off'
'set grid off'
'set parea 0.8  8.3  0.5  10.5'
'set t 1'
'set lon  -75  -34'
'set lat  -35    6'
'set clevs -2 -1'
'd pmed'
return