#!/bin/ksh
#set -vx

#**************************************************************
# Author: Martial Mancip, Nicolas Vuichard
# Contact: Nioclas.Vuichard@lsce.jussieu.fr
# $Date:
# $Name: OOL_1_9_5 $
# $Revision: 1.2 $
# History:
# Modification: Nicolas Vuichard
#
#**************************************************************

input=$@
paper=-p
quality=cx1800
#cx800 cx1200
#frequency=1:12:6
echo $input
#reference_file=RRRRRRR
#validation directory=VVVVVVV
debug=
#-d

# Ratio one day over time step
one_day=86400
ratio=$(( one_day/1800 ))

ssbx=481

#----------------------

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "PFT distribution" -o  pft ORCHIDEE_pft.jnl \
                                                                      "\"PFT distribution\" RRRRRRR inputfile MAXVEGETFRAC"


# Smooth over sbx values

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Sensible heat flux : diff with obs" -o  fluxsens_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Sensible heat flux : diff with obs\" VVVVVVV RRRRRRR inputfile Fh Qh $FYEAR \"W/m2\" $ssbx"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Latent heat flux : diff with obs" -o  fluxlat_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Latent heat flux : diff with obs\" VVVVVVV RRRRRRR inputfile Fle Qle $FYEAR \"W/m2\" $ssbx"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Surface Temperature : diff with obs" -o  tsoil_sbx ORCHIDEE_3sbxdif.jnl \
                                                                       "\"Surface Temperature : diff with obs\" VVVVVVV RRRRRRR inputfile TSOIL (AvgSurfT-273.15)  $FYEAR \"degC\" $ssbx "

# slowproc :  veget(ji,jv) = veget_max(ji,jv) * ( 1. - exp( - lai(ji,jv) * ext_coef(jv) ) )
# =>          lai = - ln( 1. - veget/veget_max) / ext_coef
ext_coef=0.5
echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "LAI : lai mean with tsoil/4 obs" -o  lai_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"lai : lai mean with tsoil/4\" VVVVVVV RRRRRRR inputfile (tsoil/4) (0.-maxvegetfrac*ln(1.-vegetfrac/maxvegetfrac)/$ext_coef) $FYEAR \"1\" $ssbx "

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Potential evapotranspiration : diff with obs" -o  epot_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Potential evapotranspiration : diff with obs\" VVVVVVV RRRRRRR inputfile (ET*${ratio}) (PotEvap*${one_day})  $FYEAR \"mm/d\" $ssbx "

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Soil moisture stress : diff with obs" -o  humr_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Soil moisture stress : diff with obs\" VVVVVVV RRRRRRR inputfile RH humrel  $FYEAR \"%\" $ssbx "

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Soil moisture content 1 : diff with obs" -o  humr1_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Soil moisture content 1 : diff with obs\" VVVVVVV RRRRRRR inputfile SWC humrel  $FYEAR \"%\" $ssbx "

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "gpp : diff with obs" -o  gpp_sbx ORCHIDEE_3sbxdif.jnl \
                                                                           "\"gpp : diff with obs\" VVVVVVV RRRRRRR inputfile (GPP*${ratio}) (GPP*${one_day}) $FYEAR \"gC/m2/day\" $ssbx"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "nee : diff with obs" -o  nee_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"nee : diff with obs\" VVVVVVV RRRRRRR inputfile (NEE*${ratio}) (NEE*${one_day}) $FYEAR \"gC/m2/day\" $ssbx"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Ecosystem respiration : diff with obs" -o  reco_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"Ecosystem respiration\" VVVVVVV RRRRRRR inputfile (Reco*${ratio}) ((maint_resp+hetero_resp+growth_resp)*${one_day})  $FYEAR \"gC/m2/day\" $ssbx "


# Annual cumulative values

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Sensible heat flux : Annual diff with obs" -o  fluxsens_year ORCHIDEE_3yeardif.jnl \
                                                                      "\"Sensible heat flux : Annual diff with obs\" VVVVVVV RRRRRRR inputfile Fh Qh $FYEAR $DUREE \"W/m2\" ave"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Latent heat flux : Annual diff with obs" -o  fluxlat_year ORCHIDEE_3yeardif.jnl \
                                                                      "\"Latent heat flux : Annual diff with obs\" VVVVVVV RRRRRRR inputfile Fle Qle $FYEAR $DUREE \"W/m2\" ave"

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Surface Temperature : Annual diff with obs" -o  tsoil_year ORCHIDEE_3yeardif.jnl \
                                                                      "\"Surface Temperature : Annual diff with obs\" VVVVVVV RRRRRRR inputfile TSOIL (AvgSurfT-273.15) $FYEAR $DUREE \"degC\" ave"

if [ "${F_FILE}" != "Cabauw" ] && [ "${FLAG}" != "SECHIBA" ]; then

    echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "gpp : Annual diff with obs" -o  gpp_year ORCHIDEE_3yeardif.jnl \
                                                                          "\"gpp : Annual diff with obs\" VVVVVVV RRRRRRR inputfile (GPP) (GPP*1800) $FYEAR $DUREE \"gC/m2/an\" sum"

fi

if [ "${F_FILE}" != "Cabauw" ] && [ "${FLAG}" = "OK_STOMATE" ]; then

    echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "nee : Annual diff with obs" -o  nee_year ORCHIDEE_3yeardif.jnl \
                                                                          "\"nee : Annual diff with obs\" VVVVVVV RRRRRRR inputfile (NEE) (NEE*1800) $FYEAR $DUREE \"gC/m2/an\" sum"

    echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "Ecosystem respiration : Annual diff with obs" -o  reco_year ORCHIDEE_3yeardif.jnl \
                                                                          "\"Ecosystem respiration : Annual diff with obs\" VVVVVVV RRRRRRR inputfile (Reco) ((maint_resp+hetero_resp+growth_resp)*1800) $FYEAR $DUREE \"gC/m2/an\" sum"
fi


# Diurnal cycle

let ANNEE=0
while [ $ANNEE -lt $DUREE ]; do
  let ANNEE2=$ANNEE+1
  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Sensible heat flux : climatologie_${ANNEE2}" -o  fluxsens_day ORCHIDEE_12day.jnl \
                                                                                               "\"Sensible heat flux\" VVVVVVV RRRRRRR inputfile Fh Qh sloop"

  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Latent heat flux : climatologie_${ANNEE2}" -o  fluxlat_day ORCHIDEE_12day.jnl \
                                                                                               "\"Latent heat flux\" VVVVVVV RRRRRRR inputfile Fle Qle sloop"

  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Surface Temperature : climatologie_${ANNEE2}" -o  tsoil_day ORCHIDEE_12day.jnl \
                                                                                               "\"Surface Temperature\" VVVVVVV RRRRRRR inputfile TSOIL (AvgSurfT-273.15) sloop"


  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "LAI : climatologie_${ANNEE2}" -o  lai_day ORCHIDEE_12day.jnl \
                                                                                               "\"LAI\" VVVVVVV RRRRRRR inputfile (tsoil/4) (0.-maxvegetfrac*ln(1.-vegetfrac/maxvegetfrac)/$ext_coef) sloop"

  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "gpp : climatologie_${ANNEE2}" -o  gpp_day ORCHIDEE_12day.jnl \
                                                                                               "\"GPP\" VVVVVVV RRRRRRR inputfile (GPP*${ratio}) (GPP*${one_day}) sloop" 

  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "nee : climatologie_${ANNEE2}" -o  nee_day ORCHIDEE_12day.jnl \
                                                                                               "\"nee\" VVVVVVV RRRRRRR inputfile (NEE*${ratio}) (NEE*${one_day}) sloop" 

  echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "Ecosystem respiration : climatologie_${ANNEE2}" -o reco_day ORCHIDEE_12day.jnl \
                                                                                               "\"Ecosystem respiration\" VVVVVVV RRRRRRR inputfile (Reco*${ratio}) ((maint_resp+hetero_resp+growth_resp)*${one_day}) sloop"

  let ANNEE=$ANNEE+1
done