|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 ! #################### 00002 MODULE MODD_CO2V_PAR 00003 ! #################### 00004 ! 00005 !!*****MODD_CO2V_PAR* 00006 !! 00007 !! PURPOSE 00008 !! ------- 00009 !! Supply constants and parameter values as 00010 !! a function of vegetation class to describe 00011 !! the vgegtation cover and for the 00012 !! stomatal conductance model of Jacobs. 00013 !! (Calvet et al. 1997, Agrig. and For. Met.) 00014 !! 00015 !! IMPLICIT ARGUMENTS 00016 !! ------------------ 00017 !! 00018 !! 00019 !! REFERENCE 00020 !! --------- 00021 !! 00022 !! AUTHOR 00023 !! ------ 00024 !! A. Boone * Meteo-France * 00025 !! 00026 !! MODIFICATIONS 00027 !! ------------- 00028 !! Original 26/10/97 00029 !! Modified 12/03/04 by P LeMoigne: XAMMIN in (kg m-2 s-1) 00030 !! Modified 12/03/04 by P LeMoigne: XFZERO taken from Calvet 98 00031 !! Modified 18/10/04 by P LeMoigne: XFZERO splitted into 2 00032 !! P Le Moigne 09/2005 AGS modifs of L. Jarlan 00033 ! S. Lafont 03/2009 : change unit of EPSO GAMM ANMAX 00034 !! A.L. Gibelin 04/2009 : correction of some AGS parameters and suppress 00035 !! useless parameters 00036 !! A.L. Gibelin 04/2009 : add parameters for PHOTO='NCB' 00037 !! A.L. Gibelin 06/2009 : add parameters for RESPSL='CNT' 00038 !! B. Decahrme 05/2012 : Add XCC_NITRO and XBIOMASST_LIM (optimization) 00039 !! 00040 !------------------------------------------------------------------------------- 00041 ! 00042 !* 0. DECLARATIONS 00043 ! ------------ 00044 ! 00045 ! *** First, CO2 parameters *** 00046 ! --------------------- 00047 ! 00048 IMPLICIT NONE 00049 ! 00050 INTEGER :: NGROUND_CO2_GAUS_LAYERS 00051 ! number of Gaussian intervals for integration over canopy 00052 ! 00053 ! 00054 REAL, PARAMETER :: XMCO2 = 44.0E-3, XMC = 12.0E-3 00055 ! molecular mass of CO2, 00056 ! and C (Carbon), respectively 00057 ! (used for conversions in kg) 00058 ! 00059 REAL, PARAMETER :: XDMAX_AGS = 0.045 00060 ! maximum specific humidity deficit (kg kg-1) 00061 ! for standard ISBA AGS (CPHOTO = 'AGS' or 'LAI') 00062 ! 00063 REAL, PARAMETER :: XPARCF = 0.48 00064 ! coefficient: PAR fraction of incoming solar radiation 00065 ! 00066 REAL, PARAMETER :: XRACCF = 1.00 00067 ! Factor for aerodynamic resistance for CO2 00068 ! 00069 REAL, PARAMETER :: XPCCO2 = 0.40 00070 ! proportion of Carbon in dry plant biomass (kgC/kgDM) 00071 ! 00072 REAL, PARAMETER :: XIAOPT = 500., XDSPOPT = 0.0 00073 ! optimum/initial values for absorbed global 00074 ! radiation, saturation deficit, respectively. 00075 00076 REAL, PARAMETER :: XXGT = 0.5 00077 ! spherical distribution of leaves 00078 ! 00079 REAL, PARAMETER :: XDIFRACF = 0.25 00080 ! used in computation of fraction of diffusion 00081 ! 00082 REAL, PARAMETER :: XXBOMEGA = 0.9442719 00083 ! Factor Assuming the albedo from simple diffusion 00084 ! of the leaf (PAR) or 'omega'=0.2 00085 ! 00086 REAL, PARAMETER :: XRDCF = 1./9. 00087 ! Dark respiration factor/coefficient 00088 ! 00089 REAL, PARAMETER :: XCONDCTMIN = 0.0002 00090 ! minimum canopy conductance 00091 ! or (1/5000) 00092 ! 00093 REAL, PARAMETER :: XCONDSTMIN = 1.0E-5 00094 ! minimum stomatal conductance for CO2 00095 ! or (1/100000) 00096 ! 00097 REAL, PARAMETER :: XANFMINIT = 1.E-8 00098 ! initial maximum leaf assimilation 00099 ! 00100 REAL, PARAMETER :: XAIRTOH2O = 1.606 00101 ! ratio XMD (AIR)/XMV (H2O) 00102 ! 00103 ! 00104 REAL, PARAMETER :: XCO2TOH2O = 1.6 00105 ! ratio of the binary diffusivities of CO2 and H2O in air 00106 ! 00107 ! 00108 REAL, PARAMETER :: XAW = 4.7, XASW = 2.8, XBW = 7.0 00109 ! coefficient of stress universal relationship 00110 ! for Woody species 00111 ! 00112 REAL, PARAMETER :: XDMAXN = 3.0E-2, XDMAXX = 3.0E-1 00113 ! minimum and maximum air deficit stress parameters 00114 ! used for Herbaceous offensive and defensive strategy 00115 ! 00116 ! 00117 ! Parameter values [C3,C4] at 25 C: 00118 ! 00119 REAL, PARAMETER, DIMENSION(2) :: XTOPT = (/25.0, 35.0/) 00120 ! optimum/initial temperatures for evaluating 00121 ! compensation points 00122 ! 00123 REAL, PARAMETER, DIMENSION(2) :: XFZERO1 = (/0.85, 0.50/) ! AGS LAI 00124 REAL, PARAMETER, DIMENSION(2) :: XFZERO2 = (/0.95, 0.60/) ! AST, LST, NIT, NCB 00125 ! ideal value of F, no photorespiration or 00126 ! saturation deficit 00127 ! 00128 REAL, PARAMETER, DIMENSION(2) :: XEPSO = (/0.017E-6,0.014E-6/) 00129 ! maximum initial quantum use efficiency 00130 ! (kgCO2 J-1 PAR ) 00131 ! 00132 REAL, PARAMETER, DIMENSION(2) :: XGAMM = (/45.0, 2.8/) 00133 ! CO2 conpensation concentration (ppmv) 00134 ! 00135 REAL, PARAMETER, DIMENSION(2) :: XQDGAMM = (/1.5, 1.5/) 00136 ! Q10 function for CO2 conpensation concentration 00137 ! 00138 REAL, PARAMETER, DIMENSION(2) :: XQDGMES = (/2.0, 2.0/) 00139 ! Q10 function for mesophyll conductance 00140 ! 00141 REAL, PARAMETER, DIMENSION(2) :: XT1GMES = (/5.0, 13.0/) 00142 ! reference temperature for computing compensation 00143 ! concentration function for mesophyll conductance: 00144 ! minimum temperature 00145 ! 00146 REAL, PARAMETER, DIMENSION(2) :: XT2GMES = (/36.0, 36.0/) 00147 ! reference temperature for computing compensation 00148 ! concentration function for mesophyll conductance: 00149 ! maximum temperature 00150 ! 00151 00152 REAL, PARAMETER, DIMENSION(2) :: XAMAX = (/2.2E-6, 1.7E-6/) 00153 ! leaf photosynthetic capacity (kgCO2 m-2 s-1) 00154 ! 00155 REAL, PARAMETER, DIMENSION(2) :: XQDAMAX = (/2.0, 2.0/) 00156 ! Q10 function for leaf photosynthetic capacity 00157 ! 00158 REAL, PARAMETER, DIMENSION(2) :: XT1AMAX = (/8.0, 13.0/) 00159 ! reference temperature for computing compensation 00160 ! concentration function for leaf photosynthetic 00161 ! capacity: minimum temperature 00162 ! 00163 REAL, PARAMETER, DIMENSION(2) :: XT2AMAX = (/38.0, 38.0/) 00164 ! reference temperature for computing compensation 00165 ! concentration function for leaf photosynthetic 00166 ! capacity: maximum temperature 00167 ! 00168 REAL, PARAMETER, DIMENSION(2) :: XAH = (/2.381, 5.323/) 00169 ! 00170 ! 00171 REAL, PARAMETER, DIMENSION(2) :: XBH = (/-0.6103, -0.8929/) 00172 ! coefficient of stress universal relationship 00173 ! for herbaceous plants 00174 ! 00175 ! 00176 ! Nitrogen option (YPHOTO='NIT' OR YPHOTO='NCB') parameters: 00177 ! 00178 REAL, PARAMETER :: XRESPFACTOR_NIT = 1.16E-7 00179 ! maintenance respiration rate (1% per day) 00180 ! of structural biomass (Faurie, 1994) [s-1] 00181 ! 00182 REAL, PARAMETER :: XCA_NIT = 0.38 00183 ! rate of nitrogen dilution of above-ground biomass at all [CO2] 00184 ! (Calvet and Soussana 2001, Gibelin et al. 2006) 00185 ! 00186 REAL, PARAMETER :: XCC_NIT = 0.753846 00187 ! Proportion of active biomass for 1t ha-1 00188 ! of total above-ground biomass [-] 00189 ! 00190 REAL, PARAMETER :: XCC_NITRO = 0.31425531725 00191 ! coef c for biomass in kg/m2 00192 ! Old : XCC_NIT/10.**XCA_NIT in nitro_decline.F90 00193 ! 00194 REAL, PARAMETER :: XBIOMASST_LIM = 4.7540042445E-2 00195 ! threshold value for leaf biomass and total 00196 ! above ground biomass in nitrogen dilution theory 00197 ! Old : XCC_NITRO**(1.0/XCA_NIT) in nitro_decline.F90 00198 ! 00199 ! 00200 ! Nitrogen option with allocation of carbon (YPHOTO='NCB') parameters: 00201 ! 00202 REAL, PARAMETER :: XCOEFF_MAINT_RESP_ZERO = 1.19E-4/86400. 00203 ! value of coefficient for maintenance respiration at 0°C (g/g/s) 00204 ! 00205 REAL, PARAMETER :: XSLOPE_MAINT_RESP = 0.16 00206 ! slope for maintenance respiration for temperature dependance (1/°C) 00207 ! 00208 REAL, PARAMETER, DIMENSION(12) :: XTAU_WOOD = 00209 (/ 0., 0., 0., 40.*365.*86400., 50.*365.*86400., 30.*365.*86400., 00210 0., 0., 0., 0., 0., 0. /) 00211 ! Residence time in woody pools (s) (YPHOTO='NCB') 00212 ! 00213 ! 00214 ! Soil carbon (YPHOTO='NCB' and YRESPSL='CNT') parameters: 00215 ! 00216 REAL, ALLOCATABLE, DIMENSION(:) :: XCN 00217 ! Biomass Carbon/Nitrogen ratio 00218 ! 00219 REAL, ALLOCATABLE, DIMENSION(:) :: XLC 00220 ! Biomass Lignin/Carbon ratio 00221 ! 00222 REAL, ALLOCATABLE, DIMENSION(:,:) :: XFRAC_LITTER 00223 ! Fraction of biomass pools going into litter pools 00224 ! 00225 REAL, ALLOCATABLE, DIMENSION(:) :: XTAU_LITTER 00226 ! Residence times in litter pools (s) 00227 ! 00228 REAL, ALLOCATABLE, DIMENSION(:,:,:) :: XFRAC_SOILCARB 00229 ! Fraction of litter decomposition flux that goes into soil 00230 ! 00231 REAL, ALLOCATABLE, DIMENSION(:) :: XTAU_SOILCARB 00232 ! Residence times in carbon pools (s) 00233 ! 00234 ! Radiative transfer parameters 00235 ! 00236 ! single scattering albedo 00237 REAL, PARAMETER :: XSSA_SUP = 0.15 ! single scatering albedo (PAR) for upper layer 00238 REAL, PARAMETER :: XSSA_INF = 0.15 ! single scatering albedo (PAR) for lower layer 00239 ! 00240 REAL, PARAMETER :: XSSA_SUP_PIR = 0.80 ! single scatering albedo (PAR) for upper layer 00241 REAL, PARAMETER :: XSSA_INF_PIR = 0.80 ! single scatering albedo (PAR) for lower layer 00242 00243 ! upper layer (calibration should depend on vegetation type...) 00244 REAL, PARAMETER :: XK_SUP = 2. ! k_sup = 2*tan(zc); where zc is the crown angle 00245 REAL, PARAMETER :: XXSI_SUP = 2. ! hight/diameter ratio of the crown 00246 ! 00247 ! lower layer (calibration should depend on vegetation type...) 00248 REAL, PARAMETER :: XK_INF = 1. ! k_sup = 2*tan(zc); where zc is the crown angle 00249 REAL, PARAMETER :: XXSI_INF = 7. ! hight/diameter ratio of the crown 00250 ! 00251 REAL, PARAMETER, DIMENSION(12) :: XLAI_SHADE = ! above this value sahding is activated 00252 (/ 0., 0., 0., 2., 3., 6., 2., 2., 2., 2., 2., 2. /) 00253 REAL, PARAMETER, DIMENSION(12) :: XXB_SUP = ! b_sup = 1/omega_sup(zs=0) -1 calibration avec omega_sup(zs=0)=0.5 00254 (/ 1., 1., 1., 1., 2., 3., 2., 2., 2., 2., 2., 2. /) 00255 REAL, PARAMETER, DIMENSION(12) :: XXB_INF = ! b_sup = 1/omega_sup(zs=0) -1 calibration avec omega_sup(zs=0)=0.5 00256 (/ 1., 1., 1., 4., 2., 4., 1., 1.5, 1.5, 1., 1., 1. /) 00257 ! 00258 ! (Calvet et al. 2008) coefs for ratio of biomass to LAI with representation of nitrogen dilution fct of CO2 00259 REAL, PARAMETER, DIMENSION(12) :: XCNAREF = 00260 (/ 0., 0., 0., 2.0, 2.8, 2.5, 1.3, 1.9, 1.9, 1.3, 1.3, 1.3 /) 00261 00262 00263 REAL, PARAMETER, DIMENSION(12) :: XPARAM = 00264 (/ 0., 0., 0., 2.56, 1.81, 1.81, 1.48, 1.48, 1.48, 1.81, 1.81, 1.81 /) 00265 00266 END MODULE MODD_CO2V_PAR 00267 00268 00269 00270 00271 00272
1.8.0
