River nutriments configuration

• To generate river forcing from 01/01/2014 to 05/01/2014, run treatment river:

launch_mangae config_mangae.cfg -s 2014-01-01,00:00 -e 2014-01-05,00:00 --only river

• If nutriments option activated, nutriments output files are in inputs/MANGAE2500-r1630/river/concentrations/:

MES-adour_st-vincent-paul.dat    N-org-adour_st-vincent-paul.dat   Si-adour_st-vincent-paul.dat
N-NH4-adour_st-vincent-paul.dat  P-PO4-adour_st-vincent-paul.dat
N-NO3-adour_st-vincent-paul.dat  P-part-adour_st-vincent-paul.dat

• Configuration is precised in configuration file config_mangae.cfg, Section river.

• To activate nutriments, set ‘nutriments = True’ and precise output_data_format. Find an example here for Adour:

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# Section river
# Préparation des forçages hydrologiques
# Entrées disponibles :
#     depend       : liste contenant les noms des traitements dont dépend le traitement courant, liste vide si aucune dépendance
#     description  : description du traitement (récupéré via l'option --list du script launcher.py)
#     name         : nom du job soumis sur le calculateur (!!!! ne pas dépasser 15 caractères !!!)
#     launch       : procédure de lancement associée
#     nutriments   : logique pour activation global du calcul des concentrations en nutriments
#     function     : fonction de préparation des données hydrologiques (�|  choisir parmi 'prepare' ,'combine' ou 'fusion')
#     startshift   : optionnel, réél indiquant le décalage en jours du début du traitement par rapport �|  la date de démarrage de référence
#     flow_dirname : optionnel, nom du dossier de stockage des fichiers de débits
#     concentration_dirname : optionnel, nom du dossier de stockage des fichiers de concentrations en nutriments
# Sous-sections :
#     rivers       : sous-section contenant les paramètres de traitements des rivières
#       nameriver  : sous-section concernant la rivière nameriver, contient les entrées suivantes
#         formula  : chaîne de caractères indiquant la méthode de calcul du débit
#         nutriments : logique pour activation du calcul des concentrations en nutriments
#     output_data_format : sous-section décrivant le format des fichiers de sortie
#       flow             : sous-section décrivant le format  des fichiers de débits
#         date_format : format des dates
#         file_format : agencement des données
#         header      : description du �~@~\header�~@~]
#         filename_format : nomenclature du fichier de sortie
#       nutriments    : sous-section décrivant le format des fichiers de débits, vide si aucun calcul de nutriments
#         date_format : format des dates
#         file_format : agencement des données
#         header      : description du �~@~\header�~@~]
#         filename_format : nomenclature du fichier de sortie
#         filename_nutriment_aliases : sous-section donnant les alias pour chaque type de nutriments
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
[[river]]
  depend = []
  description = 'Preparation of river forcings'
  name = 'mang25_river'
  launch = 'tools.river.launch_river'
  nutriments = True
  function = 'treatRivers'
  concentration_dirname='concentrations'
  [[[rivers]]]
    [[[[adour_st-vincent-paul]]]]
      formula = 'adour_saint_vincent_de_paul'
      nutriments = True
  [[[output_data_format]]]
    [[[[flow]]]]
      date_format = '%d/%m/%Y %H:%M:%S'
      file_format = 'DATE VALUE'
      header = '''["DATE VALUE (RIVER)","","",""]'''
      filename_format = 'RIVER.dat'
    [[[[nutriments]]]]
      date_format = '%d/%m/%Y %H:%M:%S'
      file_format = 'DATE VALUE'
      header = ['NUTRIMENT-RIVER.dat', '"source: relation debits/temps J.F. Guillaud - Septembre 2008"', '"unite: micromol/l sauf pour MES en g/l"', '---------------------------------------']
     filename_format = 'NUTRIMENT-RIVER.dat'

• Formula for nutriments computation are defined in file /home13/caparmor/previmer/op/lib_operational_chain/data_oco/config_nutriments.cfg, an extraction of this file:

[config]
  [[adour]]
    NNO3 = '53.2 * (Q**0.28)'
    NH4 = '207 * (Q**-0.49)'
    Si = '125'
    PPO4 = '14.7 * (Q**-0.32)'
    Ppart = '0.068 * MES*1000. + 2.6'
    Norg = '63.6'
    MES = '(1/1000.)*(0.97 * (Q**0.58))'
  [[garonne]]
    NNO3 = '26.8 * log(Q) - 15.7'
    NH4 = '4.67'
    Si = '30.74 * (Q**0.163)'
    PPO4 = '1.76'
    Ppart = '2.44'
    Norg = '33.8'
    MES = '(1/1000.)*(8.12 * exp(0.0012 * Q))'
  [[dordogne]]
    NNO3 = '-17 * log(Q) + 219'
    NH4 = '3.19'
    Si = '146'
    PPO4 = '0.73'
    Ppart = '0.42 * log(MES*1000.) + 0.45'
    Norg = '33.8'
    MES = '(1/1000.)*numpy.where(0.031 * Q - 2.09 >1, 0.031 * Q - 2.09 , 1)'
 [[charente]]
    NNO3 = '20.4 + 104.7 * log(Q) - 0.70 * Q'
    NH4 = '7.75 * (Q**-0.18)'
    Si = '104.7 * (Q**0.103)'
    PPO4 = '1.3'
    Ppart = '2.0'
    Norg = '22.3 * (Q**0.10)'
    MES = '(1/1000.)*(7.9)'