.. _forcing.river_nutriments.config:

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)'