#! /bin/ksh -x

DATE=$1 


DIR_WRK=/trab_coroca/meta/worketa_all/posfig
DIR_GRBCTL=/trab_coroca/meta/worketa_all/data/grib
#DIR_GRBCTL=/rede/comum01/sst/E05551
DIR_EXEC=${DIR_WRK}/scripts
DIR_GMP=${DIR_WRK}/gmp
DIR_GIF=${DIR_WRK}/gif
ARQCTL=ll40ganl${DATE}

echo ${ARQCTL}

cat <<EOF> precip.gs

'open ${DIR_GRBCTL}/${ARQCTL}.ctl'
'set mpdset brmap_hires'
'run ${DIR_EXEC}/rgb.gs'
'set rgb 21    0   209    234'
'set rgb 22    0   176    210'
'set rgb 23    0   170    181'
'set rgb 24    0   158    161'
'set rgb 25    0   130    156'
'set rgb 26    0   117    146'
'set rgb 27    0    99    139'
'set rgb 28    0    88    130'
'set rgb 29    0    53    106'
'set rgb 30    0    33     93'
'set rgb 31    0     8     89'
'set rgb 32    0     0     38'

********* + HP precip  **********
'set gxout shaded'
'set cmin 5'
'set cint 5'
'set lat -26.373  -19'
'set lon 310.565 319.973'
*'set clevs 5 10 15 20 30 40 60 80 '

*'set ccols 0 21 22 23 24 25 26 27 28 29 30 31 32'

DTLEUTC.1=`/home/metop/bin/caldate.3.0 ${DATE} + 24h 'hhZddmmmyyyy'`
DTLEUTC.2=`/home/metop/bin/caldate.3.0 ${DATE} + 48h 'hhZddmmmyyyy'`
DTLEUTC.3=`/home/metop/bin/caldate.3.0 ${DATE} + 72h 'hhZddmmmyyyy'`

i=1
ti=2
HP=24
 while (i<4)
 
'enable print precip${DATE}+'HP
'set clevs 0 1 3 5 10 15 20 25 30 40 50 60 80 100'
'set ccols 0 0 9 14 4 11 5 13 3 10 7 12 8 2 6' 
'set grads off' 
*'d (prec(t='ti')+prec(t='ti+1')+prec(t='ti+2')+prec(t='ti+3'))*1000'
'd sum(prec,t='ti',t='ti+3')*1000'
'run cbar.gs'
'set strsiz 0.08'
'set string 1 c 5'
'draw string 0.9 5  CPTEC/INPE/MCT '
'draw string 0.8 4.5  Eta Operacional'
'draw title  Precipitacao Total acumulada em 24h (mm)'
'set strsiz 0.18'
'draw string 5.4 0.3 Previsao ${DATE}+'HP'h, valida para 'DTLEUTC.i
'print'
'disable print precip${DATE}+'HP
'!gxgif -r -x 470 -y 360 -i precip${DATE}+'HP' -o ${DIR_GIF}/precip+${DATE}'HP'.gif'
'!rm precip${DATE}+'HP
*'q pos'
'c'
  i=i+1
  ti=ti+4
  HP=HP+24
 endwhile

'quit'

EOF

/home/metop/GrADS-1.7Beta6/bin/grads -lc "run precip.gs"

exit