MODULE YOS_AGF USE PARKIND1 ,ONLY : JPIM ,JPRB IMPLICIT NONE SAVE ! ----------------------------------------------------------------- !* ** *AGF* - ! ----------------------------------------------------------------- TYPE :: TAGF !REAL(KIND=JPRB), ALLOCATABLE :: temp_growth(:) ! Growth temperature (°C) - Is equal to t2m_month !LOGICAL :: FIRST_CALL !INTEGER(KIND=JPIM) :: KSTEP_SAV INTEGER(KIND=JPIM) :: max_temp ! Maximum temperature for saturated humidity (K) and also used as ! the size of local array to keep saturated humidity (unitless) REAL(KIND=JPRB), DIMENSION(370) :: qsfrict ! Array to keep water vapor pressure at saturation for each temperature level ! (hPa) INTEGER(KIND=JPIM) :: nlai ! Number of LAI levels (unitless) INTEGER(KIND=JPIM) :: NPATH ! Number of photosunthetic pathways, C3 and C4 (unitless) REAL(KIND=JPRB) :: zero ! Numerical constant set to 0 (unitless) REAL(KIND=JPRB) :: one ! Numerical constant set to 1 (unitless) REAL(KIND=JPRB) :: two ! Numerical constant set to 2 (unitless) REAL(KIND=JPRB) :: three ! Numerical constant set to 3 (unitless) REAL(KIND=JPRB) :: four ! Numerical constant set to 4 (unitless) REAL(KIND=JPRB) :: RR ! Ideal gas constant (J.mol^{-1}.K^{-1}) REAL(KIND=JPRB) :: min_wind ! The minimum wind (m.s^{-1}) REAL(KIND=JPRB) :: msmlr_h2o ! Molecular weight of water vapor (kg.mol^{-1}) REAL(KIND=JPRB) :: downregulation_co2_baselevel ! CO2 base level (ppm) REAL(KIND=JPRB) :: downregulation_co2_minimum ! CO2 value above which downregulation is taken into account (ppm) ! Mathematical and numerical constants REAL(KIND=JPRB) :: pi ! pi souce : http://mathworld.wolfram.com/Pi.html (unitless) REAL(KIND=JPRB) :: one_point_five ! Numerical constant set to 1.5 (unitless) REAL(KIND=JPRB) :: nine ! Numerical constant set to 9 (unitless) REAL(KIND=JPRB) :: twenty_seven ! Numerical constant set to 27 (unitless) REAL(KIND=JPRB) :: fifty_four ! Numerical constant set to 54 (unitless) REAL(KIND=JPRB) :: undef_sechiba ! The undef value used in SECHIBA (unitless) REAL(KIND=JPRB) :: undef ! Special value for stomate (unitless) REAL(KIND=JPRB) :: one_month ! One month (s) ! Physics REAL(KIND=JPRB), DIMENSION(2) :: a1 ! Empirical factor involved in the calculation of fvpd (-) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB), DIMENSION(2) :: b1 ! Empirical factor involved in the calculation of fvpd (-) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB), DIMENSION(2) :: arJV ! a coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio (mu mol e- (mu mol CO2)-1) ! See Table 3 of Kattge & Knorr (2007) ! For C4 plants, we assume that there is no ! acclimation and that for a temperature of 25°C, aSV is the same for both C4 and C3 plants (no strong jusitification - need further parametrization) REAL(KIND=JPRB), DIMENSION(2) :: brJV ! b coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio ((mu mol e- (mu mol CO2)-1) (°C)-1) ! See Table 3 of Kattge & Knorr (2007) ! We assume No acclimation term for C4 plants. REAL(KIND=JPRB), DIMENSION(2) :: aSJ ! a coefficient of the linear regression (a+bT) defining the Entropy term for Jmax (J K-1 mol-1) ! See Table 3 of Kattge & Knorr (2007) ! and Table 2 of Yin et al. (2009) for C4 plants REAL(KIND=JPRB), DIMENSION(2) :: bSJ ! b coefficient of the linear regression (a+bT) defining the Entropy term for Jmax (J K-1 mol-1 °C-1) ! See Table 3 of Kattge & Knorr (2007) ! We assume no acclimation term for C4 plants. REAL(KIND=JPRB), DIMENSION(2) :: aSV ! a coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax (J K-1 mol-1) ! See Table 3 of Kattge & Knorr (2007) ! For C4 plants, we assume that there is no ! acclimation and that at for a temperature of 25°C, aSV is the same for both C4 and C3 plants (no strong jusitification - need further parametrization) REAL(KIND=JPRB), DIMENSION(2) :: bSV ! b coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax (J K-1 mol-1 °C-1) ! See Table 3 of Kattge & Knorr (2007) ! We assume No acclimation term for C4 plants. REAL(KIND=JPRB), DIMENSION(2) :: D_Vcmax ! Energy of deactivation for Vcmax (J mol-1) ! Medlyn et al. (2002) also uses 200000. for C3 plants (same value than D_Jmax) ! 'Consequently', we use the value of D_Jmax for C4 plants REAL(KIND=JPRB), DIMENSION(2) :: D_Jmax ! Energy of deactivation for Jmax (J mol-1) ! See Table 2 of Yin et al. (2009) ! Medlyn et al. (2002) also uses 200000. for C3 plants REAL(KIND=JPRB) :: E_gamma_star ! Energy of activation for gamma_star (J mol-1) ! See Medlyn et al. (2002) from Bernacchi al. (2001) ! for C3 plants - We use the same values for C4 plants. REAL(KIND=JPRB), DIMENSION(2) :: E_gm ! Energy of activation for gm (J mol-1) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB), DIMENSION(2) :: S_gm ! Entropy term for gm (J K-1 mol-1) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB), DIMENSION(2) :: D_gm ! Energy of deactivation for gm (J mol-1) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB), DIMENSION(2) :: E_Vcmax ! Energy of activation for Vcmax (J mol-1) ! See Table 2 of Yin et al. (2009) for C4 plants ! and Kattge & Knorr (2007) for C3 plants (table 3) REAL(KIND=JPRB), DIMENSION(2) :: E_Jmax ! Energy of activation for Jmax (J mol-1) ! See Table 2 of Yin et al. (2009) for C4 plants ! and Kattge & Knorr (2007) for C3 plants (table 3) REAL(KIND=JPRB) :: E_KmC ! Energy of activation for KmC (J mol-1) ! See Medlyn et al. (2002) ! from Bernacchi al. (2001) REAL(KIND=JPRB) :: E_KmO ! Energy of activation for KmO (J mol-1) ! See Medlyn et al. (2002) ! from Bernacchi al. (2001) REAL(KIND=JPRB) :: E_Rd ! Energy of activation for Rd (J mol-1) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB) :: E_Sco ! Energy of activation for Sco (J mol-1) ! See Table 2 of Yin et al. (2009) ! Value for C4 plants is not mentioned - We use C3 for all plants. REAL(KIND=JPRB), DIMENSION(2) :: g0 ! Residual stomatal conductance when irradiance approaches zero (mol CO2 m−2 s−1 bar−1) ! Value from ORCHIDEE - No other reference. ! modify to account for the conversion for conductance to H2O to CO2 REAL(KIND=JPRB) :: gamma_star25 ! Ci-based CO2 compensation point in the absence of Rd at 25°C (ubar) ! See Medlyn et al. (2002) for C3 plants - For C4 plants, we use the same value (probably uncorrect) REAL(KIND=JPRB), DIMENSION(2) :: gm25 ! Mesophyll diffusion conductance at 25°C (mol m-2 s-1 bar-1) ! See legend of Figure 6 of Yin et al. (2009) ! and review by Flexas et al. (2008) - gm is not used for C4 plants REAL(KIND=JPRB), DIMENSION(2) :: KmC25 ! Michaelis–Menten constant of Rubisco for CO2 at 25°C (ubar) ! See Table 2 of Yin et al. (2009) for C4 ! and Medlyn et al (2002) for C3 REAL(KIND=JPRB), DIMENSION(2) :: KmO25 ! Michaelis–Menten constant of Rubisco for O2 at 25°C (ubar) ! See Table 2 of Yin et al. (2009) for C4 plants and Medlyn et al. (2002) for C3 REAL(KIND=JPRB), DIMENSION(2) :: Sco25 ! Relative CO2 /O2 specificity factor for Rubisco at 25°C (bar bar-1) ! See Table 2 of Yin et al. (2009) REAL(KIND=JPRB) :: stress_gm ! Water stress on gm REAL(KIND=JPRB) :: stress_gs ! Water stress on gs REAL(KIND=JPRB) :: stress_vcmax ! Water stress on vcmax REAL(KIND=JPRB) :: tphoto_max ! maximum photosynthesis temperature (deg C) REAL(KIND=JPRB) :: tphoto_min ! minimum photosynthesis temperature (deg C) REAL(KIND=JPRB), DIMENSION(20) :: Vcmax25 ! Maximum rate of Rubisco activity-limited carboxylation at 25°C (\mu mol.m^{-2}.s^{-1}) REAL(KIND=JPRB), DIMENSION(20) :: downregulation_co2_coeff ! Coefficient for CO2 downregulation if downregulation_co2 (used for CMIP6 6.1.11) (unitless) REAL(KIND=JPRB), DIMENSION(20) :: rstruct_const ! Structural resistance (s.m^{-1}) REAL(KIND=JPRB), DIMENSION(20) :: CO2TYPE REAL(KIND=JPRB) :: lai_level_depth REAL(KIND=JPRB) :: laimax ! Maximal LAI used for splitting LAI into N layers (m^2.m^{-2}) REAL(KIND=JPRB) :: ext_coeff ! extinction coefficient of the Monsi&Saeki relationship (1953) (unitless) !PFT-dependant REAL(KIND=JPRB) :: gb_ref ! Leaf bulk boundary layer resistance (s m-1) REAL(KIND=JPRB) :: tp_00 ! 0 degree Celsius in degree Kelvin (K) REAL(KIND=JPRB) :: pb_std ! standard pressure (hPa) REAL(KIND=JPRB) :: RG_to_PAR REAL(KIND=JPRB) :: W_to_mol ! W_to_mmol * RG_to_PAR = 2.3 REAL(KIND=JPRB) :: alpha_LL ! Conversion efficiency of absorbed light into J at strictly limiting light (mol e− (mol photon)−1) ! See comment from Yin et al. (2009) after eq. 4 ! alpha value from Medlyn et al. (2002) ! 0.093 mol CO2 fixed per mol absorbed photons ! times 4 mol e- per mol CO2 produced ! PFT-dependant REAL(KIND=JPRB) :: theta ! Convexity factor for response of J to irradiance (-) ! See Table 2 of Yin et al. (2009) ! PFT-dependant !for C4 species REAL(KIND=JPRB) :: fpsir ! Fraction of PSII e− transport rate partitioned to the C4 cycle (-) ! See Table 2 of Yin et al. (2009) - x parameter ! PFT-dependant REAL(KIND=JPRB) :: fQ ! Fraction of electrons at reduced plastoquinone ! that follow the Q-cycle (-) - Values for C3 plants are not used. ! See Table 2 of Yin et al. (2009) ! PFT-dependant REAL(KIND=JPRB) :: fpseudo ! Fraction of electrons at PSI that follow ! pseudocyclic transport (-) - Values for C3 plants are not used. ! See Table 2 of Yin et al. (2009) ! PFT-dependant REAL(KIND=JPRB) :: h_protons ! Number of protons required to produce one ATP (mol mol-1) ! See Table 2 of Yin et al. (2009) - h parameter ! PFT-dependant REAL(KIND=JPRB) :: gbs ! Bundle-sheath conductance (mol m−2 s−1 bar−1) ! See legend of Figure 6 of Yin et al. (2009) ! PFT-dependant REAL(KIND=JPRB) :: Oi ! Intercellular oxygen partial pressure (ubar) REAL(KIND=JPRB) :: alpha ! Fraction of PSII activity in the bundle sheath (-) ! See legend of Figure 6 of Yin et al. (2009) ! PFT-dependant REAL(KIND=JPRB) :: kp ! Initial carboxylation efficiency of the PEP carboxylase (mol m−2 s−1 bar−1) ! See Table 2 of Yin et al. (2009) ! PFT-dependant REAL(KIND=JPRB) :: Tetens_1 ! Ratio between molecular weight of water vapor and molecular weight ! of dry air (unitless) REAL(KIND=JPRB) :: Tetens_2 REAL(KIND=JPRB) :: ref_temp REAL(KIND=JPRB) :: mol_to_m_1 REAL(KIND=JPRB) :: ratio_H2O_to_CO2 ! Ratio of water vapor diffusivity to the CO2 diffusivity (unitless) REAL(KIND=JPRB) :: N_vert_att ! N vertical attenuation factor within the canopy REAL(KIND=JPRB) :: rveg_pft ! Potentiometer to set vegetation resistance (unitless) REAL(KIND=JPRB) :: tau_temp_air_month ! Relaxation constant for computing monthly temperature average (s) REAL(KIND=JPRB) :: cte_grav ! Acceleration of the gravity (m.s^{-2}) REAL(KIND=JPRB) :: cb ! Constant of the Louis scheme (unitless); REAL(KIND=JPRB) :: cc ! Constant of the Louis scheme (unitless); ! reference to Louis (1979) REAL(KIND=JPRB) :: cd ! Constant of the Louis scheme (unitless); ! reference to Louis (1979) REAL(KIND=JPRB) :: cp_air ! Specific heat of dry air (J.kg^{-1}.K^{-1}) REAL(KIND=JPRB) :: cepdu2 ! Squared wind shear (m^2.s^{-2}) REAL(KIND=JPRB) :: ct_karman ! Van Karmann Constant (unitless) REAL(KIND=JPRB) :: min_qc ! The minimum value for qc (qc=drag*wind) used in coupled(enerbil) and forced mode (enerbil and diffuco) END TYPE TAGF END MODULE YOS_AGF