#!/bin/ksh -x
export LANG=US
set -x

mes=oct
ano=2022

grads=/stornext/home/ioper/bin/opengrads-2.2.1.oga.1/Contents/grads

##################### CALCULA DA MEDIA MENSAL DE TMAX TMIN TMED #############################

for VAR in $(echo " TMED TMIN TMAX") ; do
 ctl=/share/ioper/tempo/SAMeT/DAILY/${VAR}/SAMeT_CPTEC_${VAR}.ctl
 d_out=/share/ioper/tempo/SAMeT/CLIMATOLOGY/MONTHLY_AVERAGE_YEARLY/${VAR}
 mkdir -p $d_out
 df=$(cal $mes 2019 |  xargs -n 1 |tail -1) 
 echo $VAR"  "$ano"  "$mes"  "$df
 $grads -blc "run  calc_MONTHLY_AVERAGE_YEARLY.gs $mes $df ${VAR} $d_out/SAMeT_CPTEC_${VAR}_${ano}_${mes}.nc ${ano}"
done

##################### INTERPOLACAO DA TMED MENSAL PARA GERAR O PET ##########################

d_out=/share/ioper/tempo/SAMeT/SCR/MEDIA_MENSAL_for_PET
grads=/stornext/home/ioper/bin/opengrads-2.2.1.oga.1/Contents/grads
cdo=/stornext/home/io/bin/cdo-1.8.2/bin/cdo

mkdir -p $d_out

df=$(cal $mes 2019 |  xargs -n 1 |tail -1)
VAR=TMED
ctl=/share/ioper/tempo/SAMeT/CLIMATOLOGY/MONTHLY_AVERAGE_YEARLY/TMED/SAMeT_CPTEC_TMED_${ano}_${mes}.nc
echo $VAR"  "$mes"  "$df"  "$mm
$grads -lc "run  calc_media_mensal_for_PET.gs $mes $df TMED $ctl $d_out/SAMeT_CPTEC_${VAR}_${ano}_${mes}.nc ${ano}"

exit






for ano in $(seq 2022 2022) ; do
  for mes in $(echo " jul "); do
    mm=$(date +"%m" -d "$mes 01 $ano")
    for VAR in $(echo "TMED") ; do 
      ctl=/share/ioper/tempo/SAMeT/CLIMATOLOGY/MONTHLY_AVERAGE_YEARLY/TMED/SAMeT_CPTEC_TMED_${ano}_${mes}.nc
      df=$(cal $mes 2019 |  xargs -n 1 |tail -1) 
      echo $VAR"  "$mes"  "$df"  "$mm
      $grads -lc "run  calc_media_mensal_for_PET.gs $mes $df TMED $ctl $d_out/SAMeT_CPTEC_${VAR}_${ano}_${mes}.nc ${ano}"
    done
  done
done
exit