SURFEX v8.1
General documentation of Surfex
modd_co2v_par.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 ! ####################
6  MODULE modd_co2v_par
7 ! ####################
8 !
9 !!*****MODD_CO2V_PAR*
10 !!
11 !! PURPOSE
12 !! -------
13 !! Supply constants and parameter values as
14 !! a function of vegetation class to describe
15 !! the vgegtation cover and for the
16 !! stomatal conductance model of Jacobs.
17 !! (Calvet et al. 1997, Agrig. and For. Met.)
18 !!
19 !! IMPLICIT ARGUMENTS
20 !! ------------------
21 !!
22 !!
23 !! REFERENCE
24 !! ---------
25 !!
26 !! AUTHOR
27 !! ------
28 !! A. Boone * Meteo-France *
29 !!
30 !! MODIFICATIONS
31 !! -------------
32 !! Original 26/10/97
33 !! Modified 12/03/04 by P LeMoigne: XAMMIN in (kg m-2 s-1)
34 !! Modified 12/03/04 by P LeMoigne: XFZERO taken from Calvet 98
35 !! Modified 18/10/04 by P LeMoigne: XFZERO splitted into 2
36 !! P Le Moigne 09/2005 AGS modifs of L. Jarlan
37 ! S. Lafont 03/2009 : change unit of EPSO GAMM ANMAX
38 !! A.L. Gibelin 04/2009 : correction of some AGS parameters and suppress
39 !! useless parameters
40 !! A.L. Gibelin 04/2009 : add parameters for PHOTO='NCB'
41 !! A.L. Gibelin 06/2009 : add parameters for RESPSL='CNT'
42 !! B. Decharme 05/2012 : Add XCC_NITRO and XBIOMASST_LIM (optimization)
43 !! R. Alkama 05/2012 : parameters for 19 vegtype rather than 12
44 !!Seferian & Delire 06/2015 : updating XAMAX peding on TRY database Kattge et al. GCB 2011
45 !!
46 !-------------------------------------------------------------------------------
47 !
48 !* 0. DECLARATIONS
49 ! ------------
50 !
51 ! *** First, CO2 parameters ***
52 ! ---------------------
53 !
54 USE modd_data_cover_par, ONLY : nvegtype_old, nvegtype_ecosg
55 !
56 IMPLICIT NONE
57 !
58 REAL, PARAMETER :: xspin_co2 = 0.1
59 ! fraction of the total spinup period used
60 ! to ramp up CO2 concentration from XCO2_START to XCO2_END
61 !
62 REAL, PARAMETER :: xmco2 = 44.0e-3, xmc = 12.0e-3
63 ! molecular mass of CO2,
64 ! and C (Carbon), respectively
65 ! (used for conversions in kg)
66 !
67 REAL, PARAMETER :: xparcf = 0.48
68 ! coefficient: PAR fraction of incoming solar radiation
69 !
70 REAL, PARAMETER :: xraccf = 1.00
71 ! Factor for aerodynamic resistance for CO2
72 !
73 REAL, PARAMETER :: xpcco2 = 0.40
74 ! proportion of Carbon in dry plant biomass (kgC/kgDM)
75 !
76 REAL, PARAMETER :: xiaopt = 500., xdspopt = 0.0
77 ! optimum/initial values for absorbed global
78 ! radiation, saturation deficit, respectively.
79 
80 REAL, PARAMETER :: xxgt = 0.5
81 ! spherical distribution of leaves
82 !
83 REAL, PARAMETER :: xdifracf = 0.25
84 ! used in computation of fraction of diffusion
85 !
86 REAL, PARAMETER :: xxbomega = 0.9442719
87 ! Factor Assuming the albedo from simple diffusion
88 ! of the leaf (PAR) or 'omega'=0.2
89 !
90 REAL, PARAMETER :: xrdcf = 1./9.
91 ! Dark respiration factor/coefficient
92 !
93 REAL, PARAMETER :: xcondctmin = 0.0002
94 ! minimum canopy conductance
95 ! or (1/5000)
96 !
97 REAL, PARAMETER :: xcondstmin = 1.0e-5
98 ! minimum stomatal conductance for CO2
99 ! or (1/100000)
100 !
101 REAL, PARAMETER :: xanfminit = 1.e-8
102 ! initial maximum leaf assimilation
103 !
104 REAL, PARAMETER :: xairtoh2o = 1.606
105 ! ratio XMD (AIR)/XMV (H2O)
106 !
107 !
108 REAL, PARAMETER :: xco2toh2o = 1.6
109 ! ratio of the binary diffusivities of CO2 and H2O in air
110 !
111 !
112 REAL, PARAMETER :: xaw = 4.7, xasw = 2.8, xbw = 7.0
113 ! coefficient of stress universal relationship
114 ! for Woody species
115 !
116 REAL, PARAMETER :: xdmaxn = 3.0e-2, xdmaxx = 3.0e-1
117 ! minimum and maximum air deficit stress parameters
118 ! used for Herbaceous offensive and defensive strategy
119 !
120 !
121 ! Parameter values [C3,C4] at 25 C:
122 !
123 REAL, PARAMETER, DIMENSION(2) :: xtopt = (/25.0, 35.0/)
124 ! optimum/initial temperatures for evaluating
125 ! compensation points
126 !
127 REAL, PARAMETER, DIMENSION(2) :: xfzero1 = (/0.85, 0.50/) ! AGS LAI
128 REAL, PARAMETER, DIMENSION(2) :: xfzero2 = (/0.95, 0.60/) ! AST, LST, NIT, NCB
129 ! CO2_atm/CO2_int with no photorespiration or saturation deficit,
130 ! used only for crops & herbaceous plants
131 REAL, PARAMETER, DIMENSION(2) :: xfzerotrop = (/0.74, 0.74/) ! Tropical forests with and without
132 ! Carrer et al. radiative transfer
133 !
134 REAL, PARAMETER, DIMENSION(2) :: xepso = (/0.017e-6,0.014e-6/)
135 ! maximum initial quantum use efficiency
136 ! (kgCO2 J-1 PAR )
137 !
138 REAL, PARAMETER, DIMENSION(2) :: xgamm = (/45.0, 2.8/)
139 ! CO2 conpensation concentration (ppmv)
140 !
141 REAL, PARAMETER, DIMENSION(2) :: xqdgamm = (/1.5, 1.5/)
142 ! Q10 function for CO2 conpensation concentration
143 !
144 REAL, PARAMETER, DIMENSION(2) :: xqdgmes = (/2.0, 2.0/)
145 ! Q10 function for mesophyll conductance
146 !
147 REAL, PARAMETER, DIMENSION(2) :: xt1gmes = (/5.0, 13.0/)
148 ! reference temperature for computing compensation
149 ! concentration function for mesophyll conductance:
150 ! minimum temperature
151 !
152 REAL, PARAMETER, DIMENSION(2) :: xt2gmes = (/36.0, 36.0/)
153 ! reference temperature for computing compensation
154 ! concentration function for mesophyll conductance:
155 ! maximum temperature
156 !
157 REAL, PARAMETER, DIMENSION(2) :: xqdamax = (/2.0, 2.0/)
158 ! Q10 function for leaf photosynthetic capacity
159 !
160 REAL, PARAMETER, DIMENSION(2) :: xt1amax = (/8.0, 13.0/)
161 ! reference temperature for computing compensation
162 ! concentration function for leaf photosynthetic
163 ! capacity: minimum temperature
164 !
165 REAL, PARAMETER, DIMENSION(2) :: xt2amax = (/38.0, 38.0/)
166 ! reference temperature for computing compensation
167 ! concentration function for leaf photosynthetic
168 ! capacity: maximum temperature
169 !
170 REAL, PARAMETER, DIMENSION(2) :: xah = (/2.381, 5.323/)
171 !
172 !
173 REAL, PARAMETER, DIMENSION(2) :: xbh = (/-0.6103, -0.8929/)
174 ! coefficient of stress universal relationship
175 ! for herbaceous plants
176 !
177 !
178 ! Nitrogen option (YPHOTO='NIT' OR YPHOTO='NCB') parameters:
179 !
180 REAL, PARAMETER :: xrespfactor_nit = 1.16e-7
181 ! maintenance respiration rate (1% per day)
182 ! of structural biomass (Faurie, 1994) [s-1]
183 !
184 REAL, PARAMETER :: xca_nit = 0.38
185 ! rate of nitrogen dilution of above-ground biomass at all [CO2]
186 ! (Calvet and Soussana 2001, Gibelin et al. 2006)
187 !
188 REAL, PARAMETER :: xcc_nit = 0.753846
189 ! Proportion of active biomass for 1t ha-1
190 ! of total above-ground biomass [-]
191 !
192 REAL, PARAMETER :: xcc_nitro = 0.31425531725
193 ! coef c for biomass in kg/m2
194 ! Old : XCC_NIT/10.**XCA_NIT in nitro_decline.F90
195 !
196 REAL, PARAMETER :: xbiomasst_lim = 4.7540042445e-2
197 ! threshold value for leaf biomass and total
198 ! above ground biomass in nitrogen dilution theory
199 ! Old : XCC_NITRO**(1.0/XCA_NIT) in nitro_decline.F90
200 !
201 !
202 ! Nitrogen option with allocation of carbon (YPHOTO='NCB') parameters:
203 !
204 REAL, PARAMETER :: xcoeff_maint_resp_zero = 1.19e-4/86400.
205 ! value of coefficient for maintenance respiration at 0C (g/g/s)
206 !
207 REAL, PARAMETER :: xslope_maint_resp = 0.16
208 ! slope for maintenance respiration for temperature dependance (1/C)
209 !
210 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xtau_wood = &
211  (/ 0., 0., 0., 40.*365.*86400., 50.*365.*86400., 30.*365.*86400., 0., 0., 0., 0., &
212  0., 0., 40.*365.*86400., 40.*365.*86400., 50.*365.*86400., 40.*365.*86400., &
213  50.*365.*86400., 0., 40.*365.*86400. /)
214 ! Residence time in woody pools (s) (YPHOTO='NCB')
215 !
216 !
217 ! Soil carbon (YPHOTO='NCB' and YRESPSL='CNT') parameters:
218 !
219 REAL, DIMENSION(6) :: xcn
220 ! Biomass Carbon/Nitrogen ratio
221 !
222 REAL, DIMENSION(6) :: xlc
223 ! Biomass Lignin/Carbon ratio
224 !
225 REAL, DIMENSION(6,2) :: xfrac_litter
226 ! Fraction of biomass pools going into litter pools
227 !
228 REAL, DIMENSION(2) :: xtau_litter
229 ! Residence times in litter pools (s)
230 !
231 REAL, DIMENSION(2,3,6) :: xfrac_soilcarb
232 ! Fraction of litter decomposition flux that goes into soil
233 !
234 REAL, DIMENSION(3) :: xtau_soilcarb
235 ! Residence times in carbon pools (s)
236 !
237 ! Radiative transfer parameters
238 !
239 ! single scattering albedo
240 REAL, PARAMETER :: xssa_sup = 0.15 ! single scatering albedo (PAR) for upper layer
241 REAL, PARAMETER :: xssa_inf = 0.15 ! single scatering albedo (PAR) for lower layer
242 !
243 REAL, PARAMETER :: xssa_sup_pir = 0.80 ! single scatering albedo (PAR) for upper layer
244 REAL, PARAMETER :: xssa_inf_pir = 0.80 ! single scatering albedo (PAR) for lower layer
245 
246 ! upper layer (calibration should depend on vegetation type...)
247 REAL, PARAMETER :: xk_sup = 2. ! k_sup = 2*tan(zc); where zc is the crown angle
248 REAL, PARAMETER :: xxsi_sup = 2. ! hight/diameter ratio of the crown
249 !
250 ! lower layer (calibration should depend on vegetation type...)
251 REAL, PARAMETER :: xk_inf = 1. ! k_sup = 2*tan(zc); where zc is the crown angle
252 REAL, PARAMETER :: xxsi_inf = 7. ! height/diameter ratio of the crown
253 !
254 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xlai_shade = & ! above this value shading is activated
255  (/ 0., 0., 0., 2., 3., 6., 2., 2., 2., 2., 2., 2., 2., 2., 3., 2., 3., 2., 2. /)
256 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xxb_sup = & ! b_sup = 1/omega_sup(zs=0) -1 calibration avec omega_sup(zs=0)=0.5
257  (/ 1., 1., 1., 1., 2., 3., 2., 2., 2., 2., 2., 2., 1., 1., 2., 1., 2., 2., 1. /)
258 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xxb_inf = & ! b_sup = 1/omega_sup(zs=0) -1 calibration avec omega_sup(zs=0)=0.5
259  (/ 1., 1., 1., 4., 2., 4., 1., 1.5, 1.5, 1., 1., 1., 4., 4., 2., 4., 2., 1., 4. /)
260 !
261 ! (Calvet et al. 2008) coefs for ratio of biomass to LAI with representation of nitrogen dilution fct of CO2
262 !REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: XPARAM = &
263 !(/ 0., 0., 0., 2.56, 1.81, 1.81, 1.48, 1.48, 1.48, 1.81, 1.81, 1.81, 2.56, 2.56, 1.81, 2.56, &
264 ! 1.81, 1.81, 2.56 /)
265 !
266 ! Kronecker flag for nitrogen dilution hypothesis (based on Yin et al., 2002)
267 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xdiludec = &
268 (/ 0., 0., 0., 1., 0., 0., 1., 1., 1., 1., 1., 0., 1., 0., 0., 1., 1., 0., 1. /)
269 !
270 ! Maximum Leaf photosynthetic capacity (kgCO2 m-2 s-1)
271 ! Modified according to Kattge et al., 2009 median Vcmax at 25C values, except for TRBE : median-std
272 ! For C3 PFTs : Ammax = Vcmax / 2. (Jacobs, p 150)
273 ! For C4 PFTs : Ammax = Vcmax
274 ! Units : [Vcmax]=micromols_CO2 m-2 s-1, [Ammax]=kgCO2 m-2 s-1 --> [Ammax] = [Vcmax] * 44e-3 * 1e-6
275 REAL, PARAMETER, DIMENSION(NVEGTYPE_OLD) :: xamax = &
276 (/ 1., 1., 1., 1.3e-6, 1.4e-6, 0.484e-6, 2.2e-6, 1.7e-6, 1.7e-6, 1.7e-6, 1.7e-6, 1.7e-6, &
277  0.9e-6, 1.3e-6, 1.4e-6, 1.3e-6, 0.9e-6, 1.7e-6, 1.2e-6/)
278 !
279 INTEGER, PARAMETER, DIMENSION(NVEGTYPE_ECOSG) :: itransfert_esg = &
280  (/1,2,3,16,4,13,14,6,5,15,17,19,18,10,11,7,7,8,4,12/)
281 
282 END MODULE modd_co2v_par
283 
284 
285 
286 
287 
288 
Generated on Tue Jan 16 2018 16:23:29 for SURFEX v8.1 by   doxygen 1.8.13