source: tags/ORCHIDEE_1_9_5/ORCHIDEE_OL/FLUXNET/atlas_FLUXNET.cfg @ 8

#!/bin/ksh
#set -vx

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

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 Qh 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 Qle 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 "lwnet : diff with obs" -o  lwnet_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"lwnet : diff with obs\" VVVVVVV RRRRRRR inputfile LWnet LWnet  $FYEAR \"%\" $ssbx "

echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "swnet : diff with obs" -o  swnet_sbx ORCHIDEE_3sbxdif.jnl \
                                                                      "\"swnet : diff with obs\" VVVVVVV RRRRRRR inputfile SWnet SWnet  $FYEAR \"%\" $ssbx "



echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -p -i "netrad : diff with obs" -o  netrad_sbx ORCHIDEE_3sbxdif.jnl \
                                                                       "\"netrad : diff with obs\" VVVVVVV RRRRRRR inputfile netrad (LWnet+SWnet)  $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 GPP*maxvegetfrac*1800) $FYEAR \"gC/m2/day\" $ssbx"
#                                                                     "\"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 "

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 (NEE*maxvegetfrac*1800) $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 Qh 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 Qle 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 Qh 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 Qle 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"
  if [ "${F_FILE}" = "Cabauw" ]; then

      echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "lwnet : climatologie_${ANNEE2}" -o  lwnet_day ORCHIDEE_12day.jnl "\"lwnet\" VVVVVVV RRRRRRR inputfile LWnet lwnet  sloop"

      echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "swnet : climatologie_${ANNEE2}" -o  swnet_day ORCHIDEE_12day.jnl "\"swnet\" VVVVVVV RRRRRRR inputfile SWnet swnet sloop"

  else

      echo $input | ferretloop - ${debug} --ps2gif_option "-g $quality" -s ${ANNEE}:${ANNEE}:1 -p -i "netrad : climatologie_${ANNEE2}" -o  netrad_day ORCHIDEE_12day.jnl "\"netrad\" VVVVVVV RRRRRRR inputfile netrad SWnet+lwnet sloop"

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

      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 (GPP*maxvegetfrac*1800) sloop" 
#                                                                                                   "\"GPP\" VVVVVVV RRRRRRR inputfile (GPP*${ratio}) (GPP*${one_day}) sloop" 

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

      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 NEE*maxvegetfrac*1800 sloop" 

  fi
  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