Intégration de la version LMD de CWRR dans la version 1.9.5.2 de ORCHIDEE

Objectif

Le but de cette modification est de transplanter dans la dernière version stable d'ORCHIDEE la version du modèle CWRR qui a été rendu opérationnelle pendant le projet WATCH. Cela implique en particulier les points suivants :

• Reprendre le module hydrol.f90 du LMD
• Reprendre le module routing.f90
• Introduire le calcule des plaines d'inondation et de l'irrigation dans 1.9.5.2
• Introduire les paramètres de la paramétrisation de la re-infiltration
• Introduire la transpiration potentielle dans diffuco.f90
• Revoir les définitions de veget et frac_bare dans les 2 version d'ORCHIDEE (routing dans 1.9.5.2 reçoit veget_max et veget dans LMD). Il faut faire attention ici de bien conserver l'eau lors de changement de fraction de végétation.

Modules concernés

Variables qui changent dans les interfaces de hydrol et routing ::

• control_in (hydrol.f90) : Structure qui porte les options du modèle.
• frac_bare : variable calculée dans slowproc et qui sert dans hydrol, sechiba, condveg et diffuco.
• floodout (hydrol.f90, routing.f90) : Débit à la sortie des plaines d'inondation. Sortie de routing et doit être INPUT de hydrol et hydrolc.
• vevapflo (hydrol.f90) : Evaporation des plaines d'inondation. Sortie de enerbil et entrée de hydrol et hydrolc.
• reinfiltration (hydrol.f90, routing.f90) : L'eau du fleuve qui se re-infiltre dans l'humidité du sol. Sortie de routing et entrée de hydrol et hydrolc
• flood_frac (hydrol.f90, routing.f90) : fraction de la maille inondée. Sortie de routing et entrée de diffuco, hydrol et hydrolc.
• flood_res (hydrol.f90, routing.f90) : quantité d'eau dans la plaine inondée. Idem que flood_frac.
• k_litt (hydrol.f90, routing.f90) : Litter conductivity : Sortie de hydrol et utilisé dans routing. Fixé pour hydrolc.
• soiltile (hydrol.f90) : distribution des types de sols par type de végétation. Sortie de slowproc et input de hydrol.
• reinf_slope (hydrol.f90, routing.f90) pentes pour le re-infiltration produit par slowproc et utilisé par hydrol.
• transpot (routing.f90) evaporation potentielle estimée dans enerbil et utilisée pour l'irrigation. Nécessite tous le vbeta3pot qui doit être diagnostiqué dans diffuco et puis passé à enerbil.

Modules concernés

• Sechiba.f90 : Modifier les interfaces et les sorties des variables (standard et ALMA)
• slowproc.f90 : Lecture de la carte des pentes et distribution des PFTs sur les types de sol.
• diffuco.f90 : Variables pour le calcule de transpot et évaporation des plaines
• enerbil.f90 : transpot et evaporation des plaines
• hydrol.f90 : tout !
• hydroc.f90 : corrections de l'irrigation et des plaines d'inondation.
• routing.f90 : Tout !

Étapes

Vérifications

• Vérifier la conservation de l'eau … d'abord dans la version Choisnel.
• Re-vérifier la conservation de l'eau en rajoutant progressivement : le routage, l'irrigation et les plaines d'inondation.
• Vérifier le cycle annuel de l'irrigation
• Vérifier l'impact des plaines d'inondation sur le débits de certains grands fleuves : Amazone, Niger et La plata.