Changes between Version 2 and Version 3 of DevelopmentActivities/ORCHIDEE-hydraulicArch


Ignore:
Timestamp:
10/09/21 16:24:09 (3 months ago)
Author:
peylin
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-hydraulicArch

    v2 v3  
    33== Meeting the 07 October 2021: comparison of Julien and Yitong schemes =  
    44 
    5 The purpose of the meeting was to compare Yitong Yao’s model with Julien Alléon’s one in order to : 
    6  - Understand the differences between both models ; 
    7  - Prepare a strategy in order to implement only one “mixed” model in ORCHIDEE_trunk. 
     5The purpose of the meeting was to compare Yitong Yao’s model with Julien Alléon’s one in order to i) Understand the differences between both models and ii) prepare a strategy in order to implement only one “mixed” model in ORCHIDEE_trunk. 
    86 
    97The meeting was organised around the presentation of the two models and the discussion on the advantages/drawbacks of each model: 
     
    119 - [attachment:Julien_Présentation.pdf Julien presentation] 
    1210 
    13 Summary of the differences: 
    14  - Starting point 
     11'''Summary of the differences:''' 
     12 - Starting point: 
    1513   - Yitong: Emilie Joetzjer’s previous version of the hydraulic architecture in the module hydraulic_arch.f90. This previous version was more a proof about the necessity of introducing and hydraulic architecture than a clean implemented model as ORCHIDEE would need.  
    1614   - Julien: New model based on ORCHIDEE_trunk version of May 2020. New subroutines introduced in hydrol module. 
    1715 - Overall scheme: Both scheme rely on the description of the water transport from soil/root interface towards leaves thanks to a resistance/capacitance scheme. 
    18  - Root absorption 
     16 - Root absorption: 
    1917  - Yitong: Rely on a development by Emily which, in addition to a soil/root resistance based on soil conductivity and root biomass, adds a weighting of the absorption according to the amount of water present in each layer. The idea is to weight more by the amount of water than by the root biomass.      
    2018  - Julien: Two ways to model root absorption: 
     
    2220   - Classic root absorption scheme relying on a dynamic resistance  
    2321 - Transport through canopy      
    24   - Resistances : 
     22  - Resistances: 
    2523   - Yitong: Dynamic resistances linked to the water potential at each stage 
    2624   - Julien: fixed resistances 
    27   - Capacitances : 
     25  - Capacitances: 
    2826   - Yitong: fixed values 
    2927   - Julien: Dynamic values linked to the water potential of each storage pool. 
    30  - Stomatal conductance T 
    31 wo decoupled stomatal conductances for assimilation and water control. Only the one that controls water is linked to leaf water potential       One stomatal conductance calculated thanks to the leaf water potential. The conductance controls assimilation and water supply. 
    32 Resolution method       Global method: 
    33 The scheme calculates transpiration demand and supply. If the supply does not reach the demand, energy budget is recalculated (mleb has to iterate to converge towards the supply). 
     28 - Stomatal conductance: 
     29  - Yitong: Two decoupled stomatal conductances for C assimilation and water control. Only the one that controls water is linked to leaf water potential         
     30  - Julien: One stomatal conductance calculated thanks to the leaf water potential. The conductance controls assimilation and water supply. 
     31 - Resolution method:    
     32  - Yitong; Global method: the scheme calculates transpiration demand and supply. If the supply does not reach the demand, energy budget is recalculated (enerbil or mleb has to iterate to converge towards the supply). Inside the model: Iteration process in order to converge towards the fluxes at each stages. Uses minpack package for the resolution    
     33  - Julien: Resolution without iterations based on prediction/correction schemes for differential equations and estimation/corrections for non differential equations. ==> Causes instabilities. 
    3434 
    35 Inside the model: 
    36 Iteration process in order to converge towards the fluxes at each stages. Uses minpack package for the resolution       Resolution without iterations based on prediction/correction schemes for differential equations and estimation/corrections for non differential equations. 
    37  
    38 •       Causes instabilities. 
    39  
    40 Main advantage  Dynamic resistances     Resolution method 
    41 Main drawback   Resolution method       Instability 
    42  
    43 Discussion : 
     35'''Discussion:''' 
    4436Both models answer to a precise need which led to the previous choices. 
    4537Both models seem to have good results on the site studied.