#!/bin/ksh


set -A modelc ' ' 'ERA5_00_ComTDry' 'ERA5_01_ComTDry'
set -A varc ' ' 'tp2m' 'mxtp' 'mntp'


contm=1
while (($contm <= 2)) ; do
model=${modelc[${contm}]}

conta=1
while (($conta <= 3)) ; do
var=${varc[${conta}]}

model2=$model

if [ $model = "ERA5_00_ComTDry" ]; then
 model2="ERA5_00"
fi

if [ $model = "ERA5_01_ComTDry" ]; then
 model2="ERA5_01"
fi

dirin="/dados/grpeta/Team/Daniela/temp_Pri/rodadas_cenapad/ciclo_diurno/nc/${model}"
dirout="/dados/grpeta/Team/Daniela/temp_Pri/rodadas_cenapad/ciclo_diurno/txt/${model}"
mkdir -p $dirout

arq=ciclo_diurno_${var}_${model2}_1986_1989


echo ${model} ${model2} ${var}
echo ${dirin}/${arq}.nc

cat <<EOF > txt.gs
'c'
'open ${dirin}/${arq}.ctl'
'sdfopen /dados/grpeta/Team/Priscila/PDoc/FAPESP/rodada/scripts/RMSP.nc'
'sdfopen /dados/grpeta/Team/Priscila/PDoc/FAPESP/rodada/scripts/RMRJ.nc'

'set lat -27 -11'
'set lon -55 -36'

indice=${var}

rc=write('${dirout}/${model2}${var}.txt','${var} - Modelo: ${model}',append)
rc=write('${dirout}/${model2}${var}.txt','RMSP  RMRJ',append)
fdata1='%9.1f'

**************************Indice*******************************************************************
'set dfile 1'

tempo=1
while (tempo<=8)

'define result1=aave(maskout('indice'.1(t='tempo',z=1),mask.2(t=1,z=1)),lon=-55,lon=-36,lat=-27,lat=-11)'
****'define result1=amax('indice'.1(t='tempo',z=1),lon=-47.5,lon=-45.5,lat=-24.2,lat=-23)'
'd result1'
pr1=subwrd(result,4)
prf1=math_format(fdata1,pr1)

'define result2=aave(maskout('indice'.1(t='tempo',z=1),mask.3(t=1,z=1)),lon=-55,lon=-36,lat=-27,lat=-11)'
****'define result2=amax('indice'.1(t='tempo',z=1),lon=-44,lon=-42,lat=-23.2,lat=-22.2)'
'd result2'
pr2=subwrd(result,4)
prf2=math_format(fdata1,pr2)


rc=write('${dirout}/${model2}${var}.txt',prf1''prf2,append)

tempo=tempo+1
endwhile

'close 3'
'close 2'
'close 1'

'quit'

EOF

/usr/local/opengrads-2.2.1.oga.1/Contents//grads -blc "run txt.gs"

rm -f txt.gs

let conta=$conta+1
done
let contm=$contm+1
done