#!/bin/ksh

typeset  -Z2 ldm 


ano=2003
mes=10
dia=15
hora=12

Run_Date=${ano}${mes}${dia}${hora}


############################### Variaveis do namelist_clima ################################
let ldm=${mes}+2
if [ "${ldm}" == "13" ] ; then
ldm=01
fi

if [ "${ldm}" == "14" ] ; then
ldm=02
fi
mesc1=`/scratchin/grupos/grpeta/fixos/util/caldate.3.0 ${ano}${ldm}${dia}${hora} + 0d   "mmm"`


let ldm=${mes}+4
if [ "${ldm}" == "13" ] ; then
ldm=01
fi

if [ "${ldm}" == "14" ] ; then
ldm=02
fi

if [ "${ldm}" == "15" ] ; then
ldm=03
fi

if [ "${ldm}" == "16" ] ; then
ldm=04
fi
mesc2=`/scratchin/grupos/grpeta/fixos/util/caldate.3.0 ${ano}${ldm}${dia}${hora} + 0d   "mmm"`


ldm1=`/scratchin/grupos/grpeta/fixos/util/caldate.3.0 ${ano}${mes}${dia}${hora} + 120d   "ld"`


mes_atual=`/scratchin/grupos/grpeta/fixos/util/caldate.3.0 ${ano}${mes}${dia}${hora} + 31d   "mm"`

trimes=ND


cat << EOF > namelist_clima2
0
/stornext/online7/pnc/Eta
${ano}
${mes}
0
${mesc1}
${ldm1}
${mesc2}
${mes_atual}
${ano}
${trimes}
${ano}
${ano}
${ano}
EOF

############################### Variaveis do namelist_clima ################################


cat << EOF > media_clima.gs

'reinit'
'set display color white'
'c'

dado=read(namelist_clima)
_incano=sublin(dado,2)

say _incano

dado=read(namelist_clima)
_DIRETAOPER=sublin(dado,2)

say _DIRETAOPER

dado=read(namelist_clima)
__ano=sublin(dado,2)

say __ano

dado=read(namelist_clima)
_mes=sublin(dado,2)

say _mes

dado=read(namelist_clima)
_inc=sublin(dado,2)

say _inc

dado=read(namelist_clima)
_mmmtri1=sublin(dado,2)

say _mmmtri1

dado=read(namelist_clima)
_ldaystri3=sublin(dado,2)

say _ldaystri3

dado=read(namelist_clima)
_mmmtri3=sublin(dado,2)

say _mmmtri3

dado=read(namelist_clima)
_mesatual=sublin(dado,2)

say _mesatual

dado=read(namelist_clima)
_anoatual=sublin(dado,2)

say _anoatual

dado=read(namelist_clima)
_TRI=sublin(dado,2)

say _TRI

dado=read(namelist_clima)
_anotitulo=sublin(dado,2)

say _anotitulo

dado=read(namelist_clima)
_yyyytri1=sublin(dado,2)

say _yyyytri1

dado=read(namelist_clima)
_yyyytri3=sublin(dado,2)

say _yyyytri3


_ano=2001
_anod=_ano+_incano
say 'open '_DIRETAOPER'/'__ano''_mes'1312/etasa40km'__ano''_mes'1312.ctl'
'open '_DIRETAOPER'/'__ano''_mes'1312/etasa40km'__ano''_mes'1312.ctl'

'define prc=0.0'
'define tmp=0.0'


while(_ano <= 2010)

   _anof=_anod+_inc

   say 'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1312/etasa40km'_ano''_mes'1312.ctl'
   say 'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1412/etasa40km'_ano''_mes'1412.ctl'
   say 'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1512/etasa40km'_ano''_mes'1512.ctl'
   say 'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1612/etasa40km'_ano''_mes'1612.ctl'
   say 'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1712/etasa40km'_ano''_mes'1712.ctl'

   'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1312/etasa40km'_ano''_mes'1312.ctl'
   'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1412/etasa40km'_ano''_mes'1412.ctl'
   'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1512/etasa40km'_ano''_mes'1512.ctl'
   'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1612/etasa40km'_ano''_mes'1612.ctl'
   'open /stornext/online3/eta/climatologia_nova/BMJ_ZHAO/'_ano''_mes'1712/etasa40km'_ano''_mes'1712.ctl'
   
   say 'define prc2=sum(prec.2,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define prc3=sum(prec.3,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define prc4=sum(prec.4,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define prc5=sum(prec.5,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define prc6=sum(prec.6,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'  

   'define prc2=sum(prec.2,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define prc3=sum(prec.3,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define prc4=sum(prec.4,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define prc5=sum(prec.5,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define prc6=sum(prec.6,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'     
   
   say 'define tmp2=ave(tp2m.2,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define tmp3=ave(tp2m.3,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define tmp4=ave(tp2m.4,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define tmp5=ave(tp2m.5,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   say 'define tmp6=ave(tp2m.6,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'

   'define tmp2=ave(tp2m.2,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define tmp3=ave(tp2m.3,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define tmp4=ave(tp2m.4,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define tmp5=ave(tp2m.5,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'
   'define tmp6=ave(tp2m.6,time=01'_mmmtri1''_anod',time='_ldaystri3''_mmmtri3''_anof')'

   
   say 'define prc=prc+((prc2+prc3+prc4+prc5+prc6)/5)'
   say 'define tmp=tmp+((tmp2+tmp3+tmp4+tmp5+tmp6)/5)'
   
   'define prc=prc+((prc2+prc3+prc4+prc5+prc6)/5)'
   'define tmp=tmp+((tmp2+tmp3+tmp4+tmp5+tmp6)/5)'
   
   _ano=_ano+1
   _anod=_anod+1

   'q files'
   say result


   'close 6'
   'close 5'
   'close 4'
   'close 3'
   'close 2'
   
endwhile

say 'define prc=(prc/10)*1000'
say 'define tmp=(tmp/10)-273.15'

'define prc=(prc/10)*1000'
'define tmp=(tmp/10)-273.15'


#*********America do Sul**************
'set lon -90.2 -30.2'
'set lat -58 15.2'
#*************************************

# Perfumaria:
'set mpdset brmap_hires'
'set grads off'
'set grid off'
'run rgbset.gs'

'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'
#'set clevs  10 50 100 150 200 250 300 400 500'
#'set ccols  0 21 22 23 41 42 44 45 34 36 '
'set clevs  30 100 200 300 400 500 600 700 800 '
'set ccols  0 21 22 23 41 42 44 45 34 36 '

say 'Plotando a Climatologia PREC'

'set gxout shaded'

*'d maskout(prc,0.9999-smav.1(t=1))'
'd prc'

'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' CLIMATOLOGIA DE PRECIPITACAO (mm/2 meses)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: Precipitacao '  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/climatologia_prec_eta_'_TRI'.gif'
'c'

'set gxout fwrite'
'set fwrite dataout/climatologia_prec_eta_'_TRI'.bin'
'd prc'
'disable fwrite'

# Perfumaria:
'set mpdset brmap_hires'
'set grads off'
'set grid off'
'run rgbset.gs'

'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'

'set clevs   13 15 17 19 21 23 25 27 29 31 33'
***GSM'set clevs   22 23 24 25 26 27 28 29 30 31 32 33'
'set ccols  45 44 43 42 41 21 22 23 24 25 27 28 29'

say 'Plotando a Climatologia TEMP'

'set gxout shaded'

*'d maskout(tmp,0.9999-smav.1(t=1))'
'd tmp'


'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' CLIMATOLOGIA DE TEMPERATURA a 2m (C)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: Temperatura a 2m'  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/climatologia_tp2m_eta_'_TRI'.gif'

'c'

'set gxout fwrite'
'set fwrite dataout/climatologia_tp2m_eta_'_TRI'.bin'
'd tmp'
'disable fwrite'




# Calculando a anomalia:
say 'open '_DIRETAOPER'/'__ano''_mes'1312/etasa40km'__ano''_mes'1312.ctl'
say 'open '_DIRETAOPER'/'__ano''_mes'1412/etasa40km'__ano''_mes'1412.ctl'
say 'open '_DIRETAOPER'/'__ano''_mes'1512/etasa40km'__ano''_mes'1512.ctl'
say 'open '_DIRETAOPER'/'__ano''_mes'1612/etasa40km'__ano''_mes'1612.ctl'
say 'open '_DIRETAOPER'/'__ano''_mes'1712/etasa40km'__ano''_mes'1712.ctl'

'open '_DIRETAOPER'/'__ano''_mes'1312/etasa40km'__ano''_mes'1312.ctl'
'open '_DIRETAOPER'/'__ano''_mes'1412/etasa40km'__ano''_mes'1412.ctl'
'open '_DIRETAOPER'/'__ano''_mes'1512/etasa40km'__ano''_mes'1512.ctl'
'open '_DIRETAOPER'/'__ano''_mes'1612/etasa40km'__ano''_mes'1612.ctl'
'open '_DIRETAOPER'/'__ano''_mes'1712/etasa40km'__ano''_mes'1712.ctl'

'q files'
say result

say 'define anomprc2=sum(prec.2,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomprc3=sum(prec.3,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomprc4=sum(prec.4,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomprc5=sum(prec.5,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomprc6=sum(prec.6,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'

'define anomprc2=sum(prec.2,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomprc3=sum(prec.3,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomprc4=sum(prec.4,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomprc5=sum(prec.5,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomprc6=sum(prec.6,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'

say 'define anomtmp2=ave(tp2m.2,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomtmp3=ave(tp2m.3,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomtmp4=ave(tp2m.4,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomtmp5=ave(tp2m.5,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
say 'define anomtmp6=ave(tp2m.6,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'

'define anomtmp2=ave(tp2m.2,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomtmp3=ave(tp2m.3,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomtmp4=ave(tp2m.4,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomtmp5=ave(tp2m.5,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'
'define anomtmp6=ave(tp2m.6,time=01'_mmmtri1''_yyyytri1',time='_ldaystri3''_mmmtri3''_yyyytri3')'

say 'define prevprc=((anomprc2+anomprc3+anomprc4+anomprc5+anomprc6)/5)*1000'
say 'define prevtmp=((anomtmp2+anomtmp3+anomtmp4+anomtmp5+anomtmp6)/5)-273.15'

'define prevprc=((anomprc2+anomprc3+anomprc4+anomprc5+anomprc6)/5)*1000'
'define prevtmp=((anomtmp2+anomtmp3+anomtmp4+anomtmp5+anomtmp6)/5)-273.15'

say 'define anomprc=prevprc-prc'
say 'define anomtmp=prevtmp-tmp'

'define anomprc=prevprc-prc'
'define anomtmp=prevtmp-tmp'

'set mpdset brmap_hires'
'set grads off'
'set grid off'
'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'

'set clevs  -400  -300    -200   -100   -50   -10    10    50     100     200   300   400'
'set ccols 26  25      24      23     22    21     0    41    42      43      44      45         35'

say 'Plotando a anomalia prevista PREC'

'set gxout shaded'

*'d maskout(anomprc,0.9999-smav.1(t=1))'
'd anomprc'

'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' PREVISAO: ANOMALIA DE PRECIPITACAO TOTAL (mm/2 meses)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: Precipitacao '  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/anomalia_prec_eta_'_TRI'.gif'
'c'

'set mpdset brmap_hires'
'set grads off'
'set grid off'
'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'

***GSM'set clevs -4 -3 -2 -1 -0.5 0.50 1 2 3 4'

'set clevs     -1.5     -1.25     -1     -0.75     -0.5    -0.25    0.25       0.5     0.75     1     1.25       1.5  '
'set ccols  46       45       44       43      42        41      0        21       22       23    24        25        26'

*** ok'set ccols 45   44    43      42     41      0      21   22   23  24 25  '
say 'Plotando a anomalia prevista TEMP'

'set gxout shaded'

*'d maskout(anomtmp,0.9999-smav.1(t=1))'
'd anomtmp'

'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' PREVISAO: ANOMALIA DE TEMPERATURA a 2m (C)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: TEMPERATURA a 2m'  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/anomalia_tp2m_eta_'_TRI'.gif'
'c'



# plotando a previsao mensal:
'set mpdset brmap_hires'
'set grads off'
'set grid off'
'run rgbset.gs'

'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'

'set clevs    13 15 17 19 21 23 25 27 29 31 33'
***GSM'set clevs    22 23 24 25 26 27 28 29 30 31 32 33'
'set ccols  45 44 43 42 41 21 22 23 24 25 27 28'
***GSM'set ccols  45 44 43 42 41 21 22 23 24 25 27 28 29'

say 'Plotando a previsao TEMP'

'set gxout shaded'

*'d maskout(prevtmp,0.9999-smav.1(t=1))'
'd prevtmp'

'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' PREVISAO: TEMPERATURA a 2m (C)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: Temperatura a 2m'  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/previsao_tp2m_eta_'_TRI'.gif'

'c'
'set gxout fwrite'
'set fwrite dataout/previsao_tp2m_eta_'_TRI'.bin'
'd prevtmp'
'disable fwrite'


'set mpdset brmap_hires'
'set grads off'
'set grid off'
'run rgbset.gs'

'set xlopts 1 7 .15'
'set ylopts 1 7 .15'
'set xlint 10'
'set ylint 10'

'set clevs  30 100 200 300 400 500 600 700 800 '
'set ccols  0 21 22 23 41 42 44 45 34 36 '

say 'Plotando a previsao PREC'

'set gxout shaded'

*'d maskout(prevprc,0.9999-smav.1(t=1))'
'd prevprc'

'q gxinfo'
lin=sublin(result,4)
ymax=subwrd(lin,6)
   
a.2=ymax+0.64
b.2=ymax+0.45
c.2=ymax+0.28
d.2=ymax+0.11
  
'set string 1 c 10 0' 
'set strsiz 0.14 0.16' 
'draw string 5.6 'a.2' PREVISAO: PRECIPITACAO TOTAL (mm/2 meses)'
'draw string 5.4 'd.2' America do Sul'
  
  
'set strsiz 0.10 0.12'
'set string 1 c 5 0'
'draw string 5.6 'b.2' ETA-CPTEC: Precipitacao'  
'draw string 5.4 'c.2' Produzida em '_mesatual' '_anoatual', valido para '_TRI' '_anotitulo
'run cbarn.gs 0.9 0 5.5 0.2'

'printim dataout/previsao_prec_eta_'_TRI'.gif'

'q dims'
say result

'c'
'set gxout fwrite'
'set fwrite dataout/previsao_prec_eta_'_TRI'.bin'
'd prevprc'
'disable fwrite'


'quit'

EOF


grads -lc "run media_clima.gs"
rm -f media_clima.gs


cat << EOF > dataout/climatologia_tp2m_eta_${trimes}.ctl
dset ^climatologia_tp2m_eta_${trimes}.bin
undef 1e+20
xdef 64 linear -35.0  0.400000
ydef 64 linear   0.0  0.400000
zdef 1 levels 1000
tdef 1 linear 12Z13${mesc1}${ano}  1hr
vars 1
tp2m    0  99  2 METRE TEMPERATURE [C]
endvars
EOF


cat << EOF > dataout/climatologia_prec_eta_${trimes}.ctl
dset ^climatologia_prec_eta_${trimes}.bin
undef 1e+20
xdef 64 linear -35.0  0.400000
ydef 64 linear   0.0  0.400000
zdef 1 levels 1000
tdef 1 linear 12Z13${mesc1}${ano}  1hr
vars 1
prec    0  99  precipitation
endvars
EOF


cat << EOF > dataout/previsao_prec_eta_${trimes}.ctl
dset ^previsao_prec_eta_${trimes}.bin
undef 1e+20
xdef 64 linear -35.0  0.400000
ydef 64 linear   0.0  0.400000
zdef 1 levels 1000
tdef 1 linear 12Z13${mesc1}${ano}  1hr
vars 1
prec    0  99  precipitation
endvars
EOF


cat << EOF > dataout/previsao_tp2m_eta_${trimes}.ctl
dset ^previsao_tp2m_eta_${trimes}.bin
undef 1e+20
xdef 64 linear -35.0  0.400000
ydef 64 linear   0.0  0.400000
zdef 1 levels 1000
tdef 1 linear 12Z13${mesc1}${ano}  1hr
vars 1
tp2m    0  99  2 METRE TEMPERATURE [C]
endvars
EOF

chmod -R 2755 dataout


exit

