1 | ! ================================================================================================================================= |
---|
2 | ! MODULE : constantes_mtc |
---|
3 | ! |
---|
4 | ! CONTACT : orchidee-help _at_ listes.ipsl.fr |
---|
5 | ! |
---|
6 | ! LICENCE : IPSL (2011) |
---|
7 | ! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC |
---|
8 | ! |
---|
9 | !>\BRIEF This module contains the standard values of the parameters for the 13 metaclasses of vegetation used by ORCHIDEE. |
---|
10 | !! |
---|
11 | !!\n DESCRIPTION: None |
---|
12 | !! |
---|
13 | !! RECENT CHANGE(S): |
---|
14 | !! |
---|
15 | !! REFERENCE(S) : |
---|
16 | !! - Kuppel, S. (2012): Doctoral Thesis, Assimilation de mesures de flux turbulents d'eau et de carbone dans un modÚle de la biosphÚre |
---|
17 | !! continentale |
---|
18 | !! - Kuppel, S., Peylin, P., Chevallier, F., Bacour, C., Maignan, F., and Richardson, A. D. (2012). Constraining a global ecosystem |
---|
19 | !! model with multi-site eddy-covariance data, Biogeosciences, 9, 3757-3776, DOI 10.5194/bg-9-3757-2012. |
---|
20 | !! - Wohlfahrt, G., M. Bahn, E. Haubner, I. Horak, W. Michaeler, K.Rottmar, U. Tappeiner, and A. Cemusca, 1999: Inter-specific |
---|
21 | !! variation of the biochemical limitation to photosynthesis and related leaf traits of 30 species from mountain grassland |
---|
22 | !! ecosystems under different land use. Plant Cell Environ., 22, 12811296. |
---|
23 | !! - Malhi, Y., Doughty, C., and Galbraith, D. (2011). The allocation of ecosystem net primary productivity in tropical forests, |
---|
24 | !! Philosophical Transactions of the Royal Society B-Biological Sciences, 366, 3225-3245, DOI 10.1098/rstb.2011.0062. |
---|
25 | !! - Earles, J. M., Yeh, S., and Skog, K. E. (2012). Timing of carbon emissions from global forest clearance, Nature Climate Change, 2, |
---|
26 | !! 682-685, Doi 10.1038/Nclimate1535. |
---|
27 | !! - Piao, S. L., Luyssaert, S., Ciais, P., Janssens, I. A., Chen, A. P., Cao, C., Fang, J. Y., Friedlingstein, P., Luo, Y. Q., and |
---|
28 | !! Wang, S. P. (2010). Forest annual carbon cost: A global-scale analysis of autotrophic respiration, Ecology, 91, 652-661, |
---|
29 | !! Doi 10.1890/08-2176.1. |
---|
30 | !! - Verbeeck, H., Peylin, P., Bacour, C., Bonal, D., Steppe, K., and Ciais, P. (2011). Seasonal patterns of co2 fluxes in amazon |
---|
31 | !! forests: Fusion of eddy covariance data and the orchidee model, Journal of Geophysical Research-Biogeosciences, 116, |
---|
32 | !! Artn G02018, Doi 10.1029/2010jg001544. |
---|
33 | !! - MacBean, N., Maignan, F., Peylin, P., Bacour, C., Breon, F. M., & Ciais, P. (2015). Using satellite data to improve the leaf |
---|
34 | !! phenology of a global terrestrial biosphere model. Biogeosciences, 12(23), 7185-7208. |
---|
35 | !! |
---|
36 | !! SVN : |
---|
37 | !! $HeadURL: $ |
---|
38 | !! $Date$ |
---|
39 | !! $Revision$ |
---|
40 | !! \n |
---|
41 | !_ ================================================================================================================================ |
---|
42 | |
---|
43 | MODULE constantes_mtc |
---|
44 | |
---|
45 | USE defprec |
---|
46 | USE constantes |
---|
47 | |
---|
48 | IMPLICIT NONE |
---|
49 | |
---|
50 | ! |
---|
51 | ! METACLASSES CHARACTERISTICS |
---|
52 | ! |
---|
53 | |
---|
54 | INTEGER(i_std), PARAMETER :: nvmc = 13 !! Number of MTCS fixed in the code (unitless) |
---|
55 | |
---|
56 | CHARACTER(len=34), PARAMETER, DIMENSION(nvmc) :: MTC_name = & !! description of the MTC (unitless) |
---|
57 | & (/ 'bare ground ', & ! 1 |
---|
58 | & 'tropical broad-leaved evergreen ', & ! 2 |
---|
59 | & 'tropical broad-leaved raingreen ', & ! 3 |
---|
60 | & 'temperate needleleaf evergreen ', & ! 4 |
---|
61 | & 'temperate broad-leaved evergreen ', & ! 5 |
---|
62 | & 'temperate broad-leaved summergreen', & ! 6 |
---|
63 | & 'boreal needleleaf evergreen ', & ! 7 |
---|
64 | & 'boreal broad-leaved summergreen', & ! 8 |
---|
65 | & 'boreal needleleaf summergreen', & ! 9 |
---|
66 | & ' C3 grass ', & ! 10 |
---|
67 | & ' C4 grass ', & ! 11 |
---|
68 | & ' C3 agriculture', & ! 12 |
---|
69 | & ' C4 agriculture' /) ! 13 |
---|
70 | |
---|
71 | |
---|
72 | ! |
---|
73 | ! VEGETATION STRUCTURE |
---|
74 | ! |
---|
75 | INTEGER(i_std),PARAMETER, DIMENSION(nvmc) :: leaf_tab_mtc = & !! leaf type (1-4, unitless) |
---|
76 | & (/ 4, 1, 1, 2, 1, 1, 2, & !! 1=broad leaved tree, 2=needle leaved tree |
---|
77 | & 1, 2, 3, 3, 3, 3 /) !! 3=grass 4=bare ground |
---|
78 | |
---|
79 | CHARACTER(len=6), PARAMETER, DIMENSION(nvmc) :: pheno_model_mtc = & !! which phenology model is used? (tabulated) |
---|
80 | & (/ 'none ', 'none ', 'moi ', 'none ', 'none ', & |
---|
81 | & 'ncdgdd', 'none ', 'ncdgdd', 'ngd ', 'moigdd', & |
---|
82 | & 'moigdd', 'moigdd', 'moigdd' /) |
---|
83 | |
---|
84 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_tropical_mtc = & !! Is PFT tropical ? (true/false) |
---|
85 | & (/ .FALSE., .TRUE., .TRUE., .FALSE., .FALSE., .FALSE., .FALSE., & |
---|
86 | & .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE. /) |
---|
87 | |
---|
88 | CHARACTER(LEN=5), PARAMETER, DIMENSION(nvmc) :: type_of_lai_mtc = & !! Type of behaviour of the LAI evolution algorithm |
---|
89 | & (/ 'inter', 'inter', 'inter', 'inter', 'inter', & !! for each vegetation type. (unitless) |
---|
90 | & 'inter', 'inter', 'inter', 'inter', 'inter', & !! Value of type_of_lai : mean or interp |
---|
91 | & 'inter', 'inter', 'inter' /) |
---|
92 | |
---|
93 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: natural_mtc = & !! natural? (true/false) |
---|
94 | & (/ .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., & |
---|
95 | & .TRUE., .TRUE., .TRUE., .TRUE., .FALSE., .FALSE. /) |
---|
96 | |
---|
97 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: veget_ori_fixed_mtc = & !! Value for veget_ori for tests in |
---|
98 | & (/ 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! 0-dim simulations (0-1, unitless) |
---|
99 | & 0.0, 0.0, 0.8, 0.0, 0.0, 0.0 /) |
---|
100 | |
---|
101 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimax_mtc = & !! laimax for maximum |
---|
102 | & (/ 0.0, 8.0, 8.0, 4.0, 4.5, 4.5, 4.0, & !! See also type of lai interpolation |
---|
103 | & 4.5, 4.0, 2.0, 2.0, 2.0, 2.0 /) !! @tex $(m^2.m^{-2})$ @endtex |
---|
104 | |
---|
105 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimin_mtc = & !! laimin for minimum lai |
---|
106 | & (/ 0.0, 8.0, 0.0, 4.0, 4.5, 0.0, 4.0, & !! See also type of lai interpolation (m^2.m^{-2}) |
---|
107 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! @tex $(m^2.m^{-2})$ @endtex |
---|
108 | |
---|
109 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: height_presc_mtc = & !! prescribed height of vegetation (m) |
---|
110 | & (/ 0.0, 30.0, 30.0, 20.0, 20.0, 20.0, 15.0, & !! Value for height_presc : one for each vegetation type |
---|
111 | & 15.0, 15.0, 0.5, 0.6, 1.0, 1.0 /) |
---|
112 | |
---|
113 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: z0_over_height_mtc = & !! Factor to calculate roughness height from |
---|
114 | & (/ 0.0, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, & !! vegetation height (unitless) |
---|
115 | & 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625 /) |
---|
116 | |
---|
117 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ratio_z0m_z0h_mtc = & !! Ratio between z0m and z0h values (roughness height for momentum and for heat) |
---|
118 | & (/ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & |
---|
119 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
120 | |
---|
121 | |
---|
122 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rveg_mtc = & !! Potentiometer to set vegetation resistance (unitless) |
---|
123 | & (/ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & !! Nathalie on March 28th, 2006 - from Fred Hourdin, |
---|
124 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
125 | |
---|
126 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: sla_mtc = & !! specif leaf area @tex $(m^2.gC^{-1})$ @endtex |
---|
127 | & (/ 1.5E-2, 1.53E-2, 2.6E-2, 9.26E-3, 2E-2, 2.6E-2, 9.26E-3, & |
---|
128 | & 2.6E-2, 1.9E-2, 2.6E-2, 2.6E-2, 2.6E-2, 2.6E-2 /) |
---|
129 | |
---|
130 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: availability_fact_mtc = & !! calculate mortality in lpj_gap |
---|
131 | & (/ undef, 0.14, 0.14, 0.10, 0.10, 0.10, 0.05, & |
---|
132 | & 0.05, 0.05, undef, undef, undef, undef /) |
---|
133 | |
---|
134 | ! |
---|
135 | ! EVAPOTRANSPIRATION (sechiba) |
---|
136 | ! |
---|
137 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rstruct_const_mtc = & !! Structural resistance. |
---|
138 | & (/ 0.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, & !! @tex $(s.m^{-1})$ @endtex |
---|
139 | & 25.0, 25.0, 2.5, 2.0, 2.0, 2.0 /) |
---|
140 | |
---|
141 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kzero_mtc = & !! A vegetation dependent constant used in the |
---|
142 | & (/ 0.0, 12.E-5, 12.E-5, 12.E-5, 12.E-5, 25.E-5, 12.E-5, & !! calculation of the surface resistance. |
---|
143 | & 25.E-5, 25.E-5, 30.E-5, 30.E-5, 30.E-5, 30.E-5 /) !! @tex $(kg.m^2.s^{-1})$ @endtex |
---|
144 | |
---|
145 | |
---|
146 | ! |
---|
147 | ! WATER (sechiba) |
---|
148 | ! |
---|
149 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: wmax_veg_mtc = & !! Volumetric available soil water capacity in each PFT |
---|
150 | & (/ 150.0, 150.0, 150.0, 150.0, 150.0, 150.0, 150.0, & !! @tex $(kg.m^{-3} of soil)$ @endtex |
---|
151 | & 150.0, 150.0, 150.0, 150.0, 150.0, 150.0 /) |
---|
152 | |
---|
153 | |
---|
154 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref4m = & !! Root profile description for the different |
---|
155 | & (/ 5.0, 0.4, 0.4, 1.0, 0.8, 0.8, 1.0, & !! vegetations types. @tex $(m^{-1})$ @endtex |
---|
156 | & 1.0, 0.8, 4.0, 1.0, 4.0, 1.0 /) !! These are the factor in the exponential which gets |
---|
157 | !! the root density as a function of depth |
---|
158 | !! Values for zmaxh = 4.0 |
---|
159 | |
---|
160 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref2m = & !! Root profile description for the different |
---|
161 | & (/ 5.0, 0.8, 0.8, 1.0, 0.8, 0.8, 1.0, & !! vegetations types. @tex $(m^{-1})$ @endtex |
---|
162 | & 1.0, 0.8, 4.0, 4.0, 4.0, 4.0 /) !! These are the factor in the exponential which gets |
---|
163 | !! the root density as a function of depth |
---|
164 | !! Values for zmaxh = 2.0 |
---|
165 | |
---|
166 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: throughfall_by_mtc = & !! Percent by PFT of precip that is not intercepted by the canopy |
---|
167 | & (/ 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, & !! (0-100, unitless) |
---|
168 | & 30.0, 30.0, 30.0, 30.0, 30.0, 30.0 /) |
---|
169 | |
---|
170 | |
---|
171 | ! |
---|
172 | ! ALBEDO (sechiba) |
---|
173 | ! |
---|
174 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_aged_vis_mtc = & !! Minimum snow albedo value for each vegetation type |
---|
175 | & (/ 0.74, 0.24, 0.07, 0.08, 0.24, 0.07, 0.18, & !! after aging (dirty old snow) (unitless), visible albedo |
---|
176 | & 0.18, 0.33, 0.57, 0.57, 0.57, 0.57 /) !! Source : Values optimized for ORCHIDEE2.0 |
---|
177 | |
---|
178 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_aged_nir_mtc = & !! Minimum snow albedo value for each vegetation type |
---|
179 | & (/ 0.50, 0.37, 0.08, 0.10, 0.37, 0.08, 0.16, & !! after aging (dirty old snow) (unitless), near infrared albedo |
---|
180 | & 0.17, 0.27, 0.44, 0.44, 0.44, 0.44 /) !! Source : Values optimized for ORCHIDEE2.0 |
---|
181 | |
---|
182 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_dec_vis_mtc = & !! Decay rate of snow albedo value for each vegetation type |
---|
183 | & (/ 0.21, 0.08, 0.17, 0.14, 0.08, 0.17, 0.05, & !! as it will be used in condveg_snow (unitless), visible albedo |
---|
184 | & 0.06, 0.09, 0.15, 0.15, 0.15, 0.15 /) !! Source : Values optimized for ORCHIDEE2.0 |
---|
185 | |
---|
186 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_dec_nir_mtc = & !! Decay rate of snow albedo value for each vegetation type |
---|
187 | & (/ 0.13, 0.10, 0.16, 0.10, 0.10, 0.16, 0.04, & !! as it will be used in condveg_snow (unitless), near infrared albedo |
---|
188 | & 0.07, 0.08, 0.12, 0.12, 0.12, 0.12 /) !! Source : Values optimized for ORCHIDEE2.0 |
---|
189 | |
---|
190 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_vis_mtc = & !! leaf albedo of vegetation type, visible albedo, optimized on 04/07/2016 |
---|
191 | & (/ 0.00, 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, & !! (unitless) |
---|
192 | & 0.03, 0.03, 0.06, 0.06, 0.06, 0.06 /) |
---|
193 | |
---|
194 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_nir_mtc = & !! leaf albedo of vegetation type, near infrared albedo, optimized on 04/07/2016 |
---|
195 | & (/ 0.00, 0.23, 0.18, 0.18, 0.20, 0.24, 0.15, & !! (unitless) |
---|
196 | & 0.26, 0.20, 0.24, 0.27, 0.28, 0.26 /) |
---|
197 | |
---|
198 | ! |
---|
199 | ! SOIL - VEGETATION |
---|
200 | ! |
---|
201 | INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pref_soil_veg_mtc = & !! The soil tile number for each vegetation |
---|
202 | & (/ 1, 2, 2, 2, 2, 2, 2, & |
---|
203 | & 2, 2, 3, 3, 3, 3 /) |
---|
204 | |
---|
205 | |
---|
206 | ! |
---|
207 | ! PHOTOSYNTHESIS |
---|
208 | ! |
---|
209 | !- |
---|
210 | ! 1 .CO2 |
---|
211 | !- |
---|
212 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_c4_mtc = & !! flag for C4 vegetation types (true/false) |
---|
213 | & (/ .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., & |
---|
214 | & .FALSE., .FALSE., .FALSE., .TRUE., .FALSE., .TRUE. /) |
---|
215 | |
---|
216 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: vcmax_fix_mtc = & !! values used for vcmax when STOMATE is not |
---|
217 | & (/ 0.0, 40.0, 50.0, 30.0, 35.0, 40.0, 30.0, & !! activated @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex |
---|
218 | & 40.0, 35.0, 60.0, 60.0, 70.0, 70.0 /) |
---|
219 | |
---|
220 | ! For C4 plant we define a very small downregulation effect as C4 plant are |
---|
221 | ! currently saturate with respect to CO2 impact on vcmax |
---|
222 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: downregulation_co2_coeff_mtc = & !! Coefficient for CO2 downregulation if downregulation_co2 (used for CMIP6 6.1.0-6.1.10) (unitless) |
---|
223 | & (/ 0.0, 0.38, 0.38, 0.28, 0.28, 0.28, 0.22, & |
---|
224 | & 0.22, 0.22, 0.26, 0.03, 0.26, 0.03 /) |
---|
225 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: downregulation_co2_coeff_new_mtc = & !! Coefficient for CO2 downregulation if downregulation_co2_new (used for CMIP6 6.1.11) (unitless) |
---|
226 | & (/ 0.0, 0.35, 0.35, 0.26, 0.26, 0.26, 0.20, & |
---|
227 | & 0.20, 0.20, 0.24, 0.03, 0.24, 0.03 /) |
---|
228 | |
---|
229 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmC_mtc = & !! Energy of activation for KmC (J mol-1) |
---|
230 | & (/undef, 79430., 79430., 79430., 79430., 79430., 79430., & !! See Medlyn et al. (2002) |
---|
231 | & 79430., 79430., 79430., 79430., 79430., 79430. /) !! from Bernacchi al. (2001) |
---|
232 | |
---|
233 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmO_mtc = & !! Energy of activation for KmO (J mol-1) |
---|
234 | & (/undef, 36380., 36380., 36380., 36380., 36380., 36380., & !! See Medlyn et al. (2002) |
---|
235 | & 36380., 36380., 36380., 36380., 36380., 36380. /) !! from Bernacchi al. (2001) |
---|
236 | |
---|
237 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Sco_mtc = & !! Energy of activation for Sco (J mol-1) |
---|
238 | & (/undef, -24460., -24460., -24460., -24460., -24460., -24460., & !! See Table 2 of Yin et al. (2009) |
---|
239 | & -24460., -24460., -24460., -24460., -24460., -24460. /) !! Value for C4 plants is not mentioned - We use C3 for all plants |
---|
240 | |
---|
241 | |
---|
242 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_gamma_star_mtc = & !! Energy of activation for gamma_star (J mol-1) |
---|
243 | & (/undef, 37830., 37830., 37830., 37830., 37830., 37830., & !! See Medlyn et al. (2002) from Bernacchi al. (2001) |
---|
244 | & 37830., 37830., 37830., 37830., 37830., 37830. /) !! for C3 plants - We use the same values for C4 plants |
---|
245 | |
---|
246 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Vcmax_mtc = & !! Energy of activation for Vcmax (J mol-1) |
---|
247 | & (/undef, 71513., 71513., 71513., 71513., 71513., 71513., & !! See Table 2 of Yin et al. (2009) for C4 plants |
---|
248 | & 71513., 71513., 71513., 67300., 71513., 67300. /) !! and Kattge & Knorr (2007) for C3 plants (table 3) |
---|
249 | |
---|
250 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Jmax_mtc = & !! Energy of activation for Jmax (J mol-1) |
---|
251 | & (/undef, 49884., 49884., 49884., 49884., 49884., 49884., & !! See Table 2 of Yin et al. (2009) for C4 plants |
---|
252 | & 49884., 49884., 49884., 77900., 49884., 77900. /) !! and Kattge & Knorr (2007) for C3 plants (table 3) |
---|
253 | |
---|
254 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: aSV_mtc = & !! a coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax (J K-1 mol-1) |
---|
255 | & (/undef, 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, & !! See Table 3 of Kattge & Knorr (2007) |
---|
256 | & 668.39, 668.39, 668.39, 641.64, 668.39, 641.64 /) !! For C4 plants, we assume that there is no |
---|
257 | !! 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) |
---|
258 | |
---|
259 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: bSV_mtc = & !! b coefficient of the linear regression (a+bT) defining the Entropy term for Vcmax (J K-1 mol-1 °C-1) |
---|
260 | & (/undef, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, & !! See Table 3 of Kattge & Knorr (2007) |
---|
261 | & -1.07, -1.07, -1.07, 0., -1.07, 0. /) !! We assume No acclimation term for C4 plants |
---|
262 | |
---|
263 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_min_mtc = & !! minimum photosynthesis temperature (deg C) |
---|
264 | & (/ undef, -4.0, -4.0, -4.0, -4.0, -4.0, -4.0, & |
---|
265 | & -4.0, -4.0, -4.0, -4.0, -4.0, -4.0 /) |
---|
266 | |
---|
267 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_max_mtc = & !! maximum photosynthesis temperature (deg C) |
---|
268 | & (/ undef, 55.0, 55.0, 55.0, 55.0, 55.0, 55.0, & |
---|
269 | & 55.0, 55.0, 55.0, 55.0, 55.0, 55.0 /) |
---|
270 | |
---|
271 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: aSJ_mtc = & !! a coefficient of the linear regression (a+bT) defining the Entropy term for Jmax (J K-1 mol-1) |
---|
272 | & (/undef, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, & !! See Table 3 of Kattge & Knorr (2007) |
---|
273 | & 659.70, 659.70, 659.70, 630., 659.70, 630. /) !! and Table 2 of Yin et al. (2009) for C4 plants |
---|
274 | |
---|
275 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: bSJ_mtc = & !! b coefficient of the linear regression (a+bT) defining the Entropy term for Jmax (J K-1 mol-1 °C-1) |
---|
276 | & (/undef, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, & !! See Table 3 of Kattge & Knorr (2007) |
---|
277 | & -0.75, -0.75, -0.75, 0., -0.75, 0. /) !! We assume no acclimation term for C4 plants |
---|
278 | |
---|
279 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Vcmax_mtc = & !! Energy of deactivation for Vcmax (J mol-1) |
---|
280 | & (/undef, 200000., 200000., 200000., 200000., 200000., 200000., & !! Medlyn et al. (2002) also uses 200000. for C3 plants (same value than D_Jmax) |
---|
281 | & 200000., 200000., 200000., 192000., 200000., 192000. /) !! 'Consequently', we use the value of D_Jmax for C4 plants |
---|
282 | |
---|
283 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Jmax_mtc = & !! Energy of deactivation for Jmax (J mol-1) |
---|
284 | & (/undef, 200000., 200000., 200000., 200000., 200000., 200000., & !! See Table 2 of Yin et al. (2009) |
---|
285 | & 200000., 200000., 200000., 192000., 200000., 192000. /) !! Medlyn et al. (2002) also uses 200000. for C3 plants |
---|
286 | |
---|
287 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_gm_mtc = & !! Energy of activation for gm (J mol-1) |
---|
288 | & (/undef, 49600., 49600., 49600., 49600., 49600., 49600., & !! See Table 2 of Yin et al. (2009) |
---|
289 | & 49600., 49600., 49600., undef, 49600., undef /) |
---|
290 | |
---|
291 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: S_gm_mtc = & !! Entropy term for gm (J K-1 mol-1) |
---|
292 | & (/undef, 1400., 1400., 1400., 1400., 1400., 1400., & !! See Table 2 of Yin et al. (2009) |
---|
293 | & 1400., 1400., 1400., undef, 1400., undef /) |
---|
294 | |
---|
295 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_gm_mtc = & !! Energy of deactivation for gm (J mol-1) |
---|
296 | & (/undef, 437400., 437400., 437400., 437400., 437400., 437400., & !! See Table 2 of Yin et al. (2009) |
---|
297 | & 437400., 437400., 437400., undef, 437400., undef /) |
---|
298 | |
---|
299 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Rd_mtc = & !! Energy of activation for Rd (J mol-1) |
---|
300 | & (/undef, 46390., 46390., 46390., 46390., 46390., 46390., & !! See Table 2 of Yin et al. (2009) |
---|
301 | & 46390., 46390., 46390., 46390., 46390., 46390. /) |
---|
302 | |
---|
303 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Vcmax25_mtc = & !! Maximum rate of Rubisco activity-limited carboxylation at 25°C |
---|
304 | & (/ undef, 45.0, 45.0, 35.0, 40.0, 50.0, 45.0, & !! @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex |
---|
305 | & 35.0, 35.0, 50.0, 50.0, 60.0, 60.0 /) |
---|
306 | |
---|
307 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: arJV_mtc = & !! a coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio (mu mol e- (mu mol CO2)-1) |
---|
308 | & (/undef, 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, & !! See Table 3 of Kattge & Knorr (2007) |
---|
309 | & 2.59, 2.59, 2.59, 1.715, 2.59, 1.715 /) !! For C4 plants, we assume that there is no |
---|
310 | !! 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) |
---|
311 | |
---|
312 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: brJV_mtc = & !! b coefficient of the linear regression (a+bT) defining the Jmax25/Vcmax25 ratio ((mu mol e- (mu mol CO2)-1) (°C)-1) |
---|
313 | & (/undef, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, & !! See Table 3 of Kattge & Knorr (2007) |
---|
314 | & -0.035, -0.035, -0.035, 0., -0.035, 0. /) !! We assume No acclimation term for C4 plants |
---|
315 | |
---|
316 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmC25_mtc = & !! MichaelisâMenten constant of Rubisco for CO2 at 25°C (ubar) |
---|
317 | & (/undef, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, & !! See Table 2 of Yin et al. (2009) for C4 |
---|
318 | & 404.9, 404.9, 404.9, 650., 404.9, 650. /) !! and Medlyn et al (2002) for C3 |
---|
319 | |
---|
320 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmO25_mtc = & !! MichaelisâMenten constant of Rubisco for O2 at 25°C (ubar) |
---|
321 | & (/undef, 278400., 278400., 278400., 278400., 278400., 278400., & !! See Table 2 of Yin et al. (2009) for C4 plants and Medlyn et al. (2002) for C3 |
---|
322 | & 278400., 278400., 278400., 450000., 278400., 450000. /) |
---|
323 | |
---|
324 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Sco25_mtc = & !! Relative CO2 /O2 specificity factor for Rubisco at 25ðC (bar bar-1) |
---|
325 | & (/undef, 2800., 2800., 2800., 2800., 2800., 2800., & !! See Table 2 of Yin et al. (2009) |
---|
326 | & 2800., 2800., 2800., 2590., 2800., 2590. /) |
---|
327 | |
---|
328 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gm25_mtc = & !! Mesophyll diffusion conductance at 25ðC (mol m-2 s-1 bar-1) |
---|
329 | & (/undef, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4, & !! See legend of Figure 6 of Yin et al. (2009) |
---|
330 | & 0.4, 0.4, 0.4, undef, 0.4, undef /) !! and review by Flexas et al. (2008) - gm is not used for C4 plants |
---|
331 | |
---|
332 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gamma_star25_mtc = & !! Ci-based CO2 compensation point in the absence of Rd at 25°C (ubar) |
---|
333 | & (/undef, 42.75, 42.75, 42.75, 42.75, 42.75, 42.75, & !! See Medlyn et al. (2002) for C3 plants - For C4 plants, we use the same value (probably uncorrect) |
---|
334 | & 42.75, 42.75, 42.75, 42.75, 42.75, 42.75 /) |
---|
335 | |
---|
336 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: a1_mtc = & !! Empirical factor involved in the calculation of fvpd (-) |
---|
337 | & (/undef, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, & !! See Table 2 of Yin et al. (2009) |
---|
338 | & 0.85, 0.85, 0.85, 0.72, 0.85, 0.72 /) |
---|
339 | |
---|
340 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: b1_mtc = & !! Empirical factor involved in the calculation of fvpd (-) |
---|
341 | & (/undef, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, & !! See Table 2 of Yin et al. (2009) |
---|
342 | & 0.14, 0.14, 0.14, 0.20, 0.14, 0.20 /) |
---|
343 | |
---|
344 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: g0_mtc = & !! Residual stomatal conductance when irradiance approaches zero (mol CO2 mâ2 sâ1 barâ1) |
---|
345 | & (/undef, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, & !! Value from ORCHIDEE - No other reference. |
---|
346 | & 0.00625, 0.00625, 0.00625, 0.01875, 0.00625, 0.01875 /) !! modofy to account for the conversion for conductance to H2O to CO2 |
---|
347 | |
---|
348 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: h_protons_mtc = & !! Number of protons required to produce one ATP (mol mol-1) |
---|
349 | & (/undef, 4., 4., 4., 4., 4., 4., & !! See Table 2 of Yin et al. (2009) - h parameter |
---|
350 | & 4., 4., 4., 4., 4., 4. /) |
---|
351 | |
---|
352 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpsir_mtc = & !! Fraction of PSII eâ transport rate |
---|
353 | & (/undef, undef, undef, undef, undef, undef, undef, & !! partitioned to the C4 cycle (-) |
---|
354 | & undef, undef, undef, 0.4, undef, 0.4 /) !! See Table 2 of Yin et al. (2009) - x parameter |
---|
355 | |
---|
356 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fQ_mtc = & !! Fraction of electrons at reduced plastoquinone |
---|
357 | & (/undef, undef, undef, undef, undef, undef, undef, & !! that follow the Q-cycle (-) - Values for C3 platns are not used |
---|
358 | & undef, undef, undef, 1., undef, 1. /) !! See Table 2 of Yin et al. (2009) |
---|
359 | |
---|
360 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpseudo_mtc = & !! Fraction of electrons at PSI that follow |
---|
361 | & (/undef, undef, undef, undef, undef, undef, undef, & !! pseudocyclic transport (-) - Values for C3 platns are not used |
---|
362 | & undef, undef, undef, 0.1, undef, 0.1 /) !! See Table 2 of Yin et al. (2009) |
---|
363 | |
---|
364 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kp_mtc = & !! Initial carboxylation efficiency of the PEP carboxylase (mol mâ2 sâ1 barâ1) |
---|
365 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See Table 2 of Yin et al. (2009) |
---|
366 | & undef, undef, undef, 0.7, undef, 0.7 /) |
---|
367 | |
---|
368 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alpha_mtc = & !! Fraction of PSII activity in the bundle sheath (-) |
---|
369 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See legend of Figure 6 of Yin et al. (2009) |
---|
370 | & undef, undef, undef, 0.1, undef, 0.1 /) |
---|
371 | |
---|
372 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gbs_mtc = & !! Bundle-sheath conductance (mol mâ2 sâ1 barâ1) |
---|
373 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See legend of Figure 6 of Yin et al. (2009) |
---|
374 | & undef, undef, undef, 0.003, undef, 0.003 /) |
---|
375 | |
---|
376 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: theta_mtc = & !! Convexity factor for response of J to irradiance (-) |
---|
377 | & (/undef, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, & !! See Table 2 of Yin et al. (2009) |
---|
378 | & 0.7, 0.7, 0.7, 0.7, 0.7, 0.7 /) |
---|
379 | |
---|
380 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alpha_LL_mtc = & !! Conversion efficiency of absorbed light into J at strictly limiting light (mol eâ (mol photon)â1) |
---|
381 | & (/undef, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, & !! See comment from Yin et al. (2009) after eq. 4 |
---|
382 | & 0.3, 0.3, 0.3, 0.3, 0.3, 0.3 /) !! alpha value from Medlyn et al. (2002) |
---|
383 | !! 0.093 mol CO2 fixed per mol absorbed photons |
---|
384 | !! times 4 mol e- per mol CO2 produced |
---|
385 | |
---|
386 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_vcmax_mtc = & !! Water stress on vcmax |
---|
387 | & (/ 1., 1., 1., 1., 1., 1., 1., & |
---|
388 | & 1., 1., 1., 1., 1., 1. /) |
---|
389 | |
---|
390 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_gs_mtc = & !! Water stress on gs |
---|
391 | & (/ 1., 1., 1., 1., 1., 1., 1., & |
---|
392 | & 1., 1., 1., 1., 1., 1. /) |
---|
393 | |
---|
394 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_gm_mtc = & !! Water stress on gm |
---|
395 | & (/ 1., 1., 1., 1., 1., 1., 1., & |
---|
396 | & 1., 1., 1., 1., 1., 1. /) |
---|
397 | |
---|
398 | !- |
---|
399 | ! 2 .Stomate |
---|
400 | !- |
---|
401 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ext_coeff_mtc = & !! extinction coefficient of the Monsi&Saeki |
---|
402 | & (/ 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, & !! relationship (1953) (unitless) |
---|
403 | & 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 /) |
---|
404 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ext_coeff_vegetfrac_mtc = & !! extinction coefficient used for defining the fraction |
---|
405 | & (/ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & !! of bare soil (unitless) |
---|
406 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
407 | |
---|
408 | ! |
---|
409 | ! ALLOCATION (stomate) |
---|
410 | ! |
---|
411 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: R0_mtc = & !! Default root allocation (0-1, unitless) |
---|
412 | & (/ undef, 0.30, 0.30, 0.30, 0.30, 0.30, 0.30, & |
---|
413 | & 0.30, 0.30, 0.30, 0.30, 0.30, 0.30 /) |
---|
414 | |
---|
415 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: S0_mtc = & !! Default sapwood allocation (0-1, unitless) |
---|
416 | & (/ undef, 0.25, 0.25, 0.30, 0.30, 0.30, 0.30, & |
---|
417 | & 0.30, 0.30, 0.30, 0.30, 0.30, 0.30 /) |
---|
418 | |
---|
419 | ! |
---|
420 | ! RESPIRATION (stomate) |
---|
421 | ! |
---|
422 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: frac_growthresp_mtc = & !! fraction of GPP which is lost as growth respiration |
---|
423 | & (/ undef, 0.35, 0.35, 0.28, 0.28, 0.28, 0.35, & |
---|
424 | & 0.35, 0.35, 0.28, 0.28, 0.28, 0.28 /) |
---|
425 | |
---|
426 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_c_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
427 | & (/ undef, 0.12, 0.12, 0.16, 0.16, 0.16, 0.25, & !! constant c of aT^2+bT+c, tabulated |
---|
428 | & 0.25, 0.25, 0.16, 0.12, 0.16, 0.12 /) |
---|
429 | |
---|
430 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_b_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
431 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! constant b of aT^2+bT+c, tabulated |
---|
432 | & 0.0, 0.0, -0.00133, 0.0, -0.00133, 0.0 /) |
---|
433 | |
---|
434 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_a_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
435 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! constant a of aT^2+bT+c, tabulated |
---|
436 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) |
---|
437 | |
---|
438 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_leaf_mtc = & !! maintenance respiration coefficient |
---|
439 | & (/ undef, 2.35E-3, 2.62E-3, 1.01E-3, 2.35E-3, 2.62E-3, 1.01E-3, & !! at 0 deg C,for leaves, tabulated, |
---|
440 | & 2.62E-3, 2.05E-3, 2.62E-3, 2.62E-3, 2.62E-3, 2.62E-3 /) !! @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
441 | |
---|
442 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapabove_mtc = & !! maintenance respiration coefficient |
---|
443 | & (/ undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, & !! at 0 deg C, for sapwood above, |
---|
444 | & 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
445 | |
---|
446 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapbelow_mtc = & !! maintenance respiration coefficient |
---|
447 | & (/ undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, & !! at 0 deg C, for sapwood below, |
---|
448 | & 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
449 | |
---|
450 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartabove_mtc = & !! maintenance respiration coefficient |
---|
451 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! at 0 deg C, for heartwood above, |
---|
452 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
453 | |
---|
454 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartbelow_mtc = & !! maintenance respiration coefficient |
---|
455 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! at 0 deg C, for heartwood below, |
---|
456 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
457 | |
---|
458 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_root_mtc = & !! maintenance respiration coefficient |
---|
459 | & (/ undef, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, & !! at 0 deg C, for roots, tabulated, |
---|
460 | & 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3 /) !! @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
461 | |
---|
462 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_fruit_mtc = & !! maintenance respiration coefficient |
---|
463 | & (/ undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, & !! at 0 deg C, for fruits, tabulated, |
---|
464 | & 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 /) !! @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
465 | |
---|
466 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_carbres_mtc = & !! maintenance respiration coefficient |
---|
467 | & (/ undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, & !! at 0 deg C, for carbohydrate reserve, |
---|
468 | & 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
469 | |
---|
470 | |
---|
471 | ! |
---|
472 | ! FIRE (stomate) |
---|
473 | ! |
---|
474 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: flam_mtc = & !! flamability: critical fraction of water |
---|
475 | & (/ undef, 0.15, 0.25, 0.25, 0.25, 0.25, 0.25, & !! holding capacity (0-1, unitless) |
---|
476 | & 0.25, 0.25, 0.25, 0.25, 0.35, 0.35 /) |
---|
477 | |
---|
478 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: resist_mtc = & !! fire resistance (0-1, unitless) |
---|
479 | & (/ undef, 0.95, 0.90, 0.90, 0.90, 0.90, 0.90, & |
---|
480 | & 0.90, 0.90, 0.0, 0.0, 0.0, 0.0 /) |
---|
481 | |
---|
482 | |
---|
483 | ! |
---|
484 | ! FLUX - LUC |
---|
485 | ! |
---|
486 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_1_mtc = & !! Coeff of biomass export for the year |
---|
487 | & (/ undef, 0.897, 0.897, 0.597, 0.597, 0.597, 0.597, & !! (0-1, unitless) |
---|
488 | & 0.597, 0.597, 0.597, 0.597, 0.597, 0.597 /) |
---|
489 | |
---|
490 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_10_mtc = & !! Coeff of biomass export for the decade |
---|
491 | & (/ undef, 0.103, 0.103, 0.299, 0.299, 0.299, 0.299, & !! (0-1, unitless) |
---|
492 | & 0.299, 0.299, 0.299, 0.403, 0.299, 0.403 /) |
---|
493 | |
---|
494 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_100_mtc = & !! Coeff of biomass export for the century |
---|
495 | & (/ undef, 0.0, 0.0, 0.104, 0.104, 0.104, 0.104, & !! (0-1, unitless) |
---|
496 | & 0.104, 0.104, 0.104, 0.0, 0.104, 0.0 /) |
---|
497 | |
---|
498 | |
---|
499 | ! |
---|
500 | ! PHENOLOGY |
---|
501 | ! |
---|
502 | ! The latest modifications regarding leafagecrit, senescence_temp_c, leaffall, hum_min_time and nosenescence_hum are inspired by |
---|
503 | ! MacBean et al. (2015), following the optimization of phenology parameters using MODIS NDVI (FM/PP). |
---|
504 | !- |
---|
505 | ! 1. Stomate |
---|
506 | !- |
---|
507 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: lai_max_to_happy_mtc = & !! threshold of LAI below which plant uses carbohydrate reserves |
---|
508 | & (/ undef, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, & |
---|
509 | & 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 /) |
---|
510 | |
---|
511 | REAL(r_std), PARAMETER, DIMENSION (nvmc) :: lai_max_mtc = & !! maximum LAI, PFT-specific |
---|
512 | & (/ undef, 7.0, 5.0, 5.0, 4.0, 5.0, 3.5, & !! @tex $(m^2.m^{-2})$ @endtex |
---|
513 | & 4.0, 3.0, 2.5, 2.0, 5.0, 5.0 /) |
---|
514 | |
---|
515 | INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pheno_type_mtc = & !! type of phenology (0-4, unitless) |
---|
516 | & (/ 0, 1, 3, 1, 1, 2, 1, & !! 0=bare ground 1=evergreen, 2=summergreen, |
---|
517 | & 2, 2, 4, 4, 2, 3 /) !! 3=raingreen, 4=perennial |
---|
518 | !- |
---|
519 | ! 2. Leaf Onset |
---|
520 | !- |
---|
521 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_c_mtc = & !! critical gdd, tabulated (C), |
---|
522 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant c of aT^2+bT+c |
---|
523 | & undef, undef, 320.0, 400.0, 320.0, 700.0 /) |
---|
524 | |
---|
525 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_b_mtc = & !! critical gdd, tabulated (C), |
---|
526 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant b of aT^2+bT+c |
---|
527 | & undef, undef, 6.25, 0.0, 6.25, 0.0 /) |
---|
528 | |
---|
529 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_a_mtc = & !! critical gdd, tabulated (C), |
---|
530 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant a of aT^2+bT+c |
---|
531 | & undef, undef, 0.03125, 0.0, 0.0315, 0.0 /) |
---|
532 | |
---|
533 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_moigdd_t_crit_mtc = & !! temperature threshold for C4 grass(C) |
---|
534 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
535 | & undef, undef, undef, 22.0, undef, undef /) |
---|
536 | |
---|
537 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ngd_crit_mtc = & !! critical ngd, tabulated. |
---|
538 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! Threshold -5 degrees (days) |
---|
539 | & undef, 17.0, undef, undef, undef, undef /) |
---|
540 | |
---|
541 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ncdgdd_temp_mtc = & !! critical temperature for the ncd vs. gdd |
---|
542 | & (/ undef, undef, undef, undef, undef, 5.0, undef, & !! function in phenology (C) |
---|
543 | & 0.0, undef, undef, undef, undef, undef /) |
---|
544 | |
---|
545 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_frac_mtc = & !! critical humidity (relative to min/max) |
---|
546 | & (/ undef, undef, 0.5, undef, undef, undef, undef, & !! for phenology (unitless) |
---|
547 | & undef, undef, 0.5, 0.5, 0.5, 0.5 /) |
---|
548 | |
---|
549 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_min_time_mtc = & !! minimum time elapsed since |
---|
550 | & (/ undef, undef, 50.0, undef, undef, undef, undef, & !! moisture minimum (days) |
---|
551 | & undef, undef, 36.0, 35.0, 75.0, 75.0 /) |
---|
552 | |
---|
553 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_sap_mtc = & !! time (days) |
---|
554 | & (/ undef, 730.0, 730.0, 730.0, 730.0, 730.0, 730.0, & |
---|
555 | & 730.0, 730.0, undef, undef, undef, undef /) |
---|
556 | |
---|
557 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_leafinit_mtc = & !! time to attain the initial foliage using the carbohydrate reserve |
---|
558 | & (/ undef, 10., 10., 10., 10., 10., 10., & |
---|
559 | & 10., 10., 10., 10., 10., 10. /) |
---|
560 | |
---|
561 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_fruit_mtc = & !! fruit lifetime (days) |
---|
562 | & (/ undef, 90.0, 90.0, 90.0, 90.0, 90.0, 90.0, & |
---|
563 | & 90.0, 90.0, undef, undef, undef, undef /) |
---|
564 | |
---|
565 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ecureuil_mtc = & !! fraction of primary leaf and root allocation |
---|
566 | & (/ undef, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, & !! put into reserve (0-1, unitless) |
---|
567 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
568 | |
---|
569 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_min_mtc = & !! NEW - allocation above/below = f(age) |
---|
570 | & (/ undef, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, & !! - 30/01/04 NV/JO/PF |
---|
571 | & 0.2, 0.2, undef, undef, undef, undef /) |
---|
572 | |
---|
573 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_max_mtc = & !! NEW - allocation above/below = f(age) |
---|
574 | & (/ undef, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, & !! - 30/01/04 NV/JO/PF |
---|
575 | & 0.8, 0.8, undef, undef, undef, undef /) |
---|
576 | |
---|
577 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: demi_alloc_mtc = & !! NEW - allocation above/below = f(age) |
---|
578 | & (/ undef, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, & !! - 30/01/04 NV/JO/PF |
---|
579 | & 5.0, 5.0, undef, undef, undef, undef /) |
---|
580 | |
---|
581 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaflife_mtc = & !! leaf longevity, tabulated (??units??) |
---|
582 | & (/ undef, 0.5, 2.0, 0.33, 1.0, 2.0, 0.33, & |
---|
583 | & 2.0, 2.0, 2.0, 2.0, 2.0, 2.0 /) |
---|
584 | !- |
---|
585 | ! 3. Senescence |
---|
586 | !- |
---|
587 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaffall_mtc = & !! length of death of leaves, tabulated (days) |
---|
588 | & (/ undef, undef, 10.0, undef, undef, 30.0, undef, & |
---|
589 | & 5.0, 10.0, 10.0, 10.0, 10.0, 10.0 /) |
---|
590 | |
---|
591 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leafagecrit_mtc = & !! critical leaf age, tabulated (days) |
---|
592 | & (/ undef, 730.0, 180.0, 910.0, 730.0, 160.0, 910.0, & |
---|
593 | & 220.0, 120.0, 80.0, 120.0, 90.0, 90.0 /) |
---|
594 | |
---|
595 | CHARACTER(LEN=6), PARAMETER, DIMENSION(nvmc) :: senescence_type_mtc = & !! type of senescence, tabulated (unitless) |
---|
596 | & (/ 'none ', 'none ', 'dry ', 'none ', 'none ', & |
---|
597 | & 'cold ', 'none ', 'cold ', 'cold ', 'mixed ', & |
---|
598 | & 'mixed ', 'mixed ', 'mixed ' /) |
---|
599 | |
---|
600 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_hum_mtc = & !! critical relative moisture availability |
---|
601 | & (/ undef, undef, 0.3, undef, undef, undef, undef, & !! for senescence (0-1, unitless) |
---|
602 | & undef, undef, 0.2, 0.2, 0.3, 0.2 /) |
---|
603 | |
---|
604 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: nosenescence_hum_mtc = & !! relative moisture availability above which |
---|
605 | & (/ undef, undef, 0.8, undef, undef, undef, undef, & !! there is no humidity-related senescence |
---|
606 | & undef, undef, 0.6, 0.3, 0.3, 0.3 /) !! (0-1, unitless) |
---|
607 | |
---|
608 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: max_turnover_time_mtc = & !! maximum turnover time for grasses (days) |
---|
609 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
610 | & undef, undef, 80.0, 80.0, 80.0, 80.0 /) |
---|
611 | |
---|
612 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_turnover_time_mtc = & !! minimum turnover time for grasses (days) |
---|
613 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
614 | & undef, undef, 10.0, 10.0, 10.0, 10.0 /) |
---|
615 | |
---|
616 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_leaf_age_for_senescence_mtc = & !! minimum leaf age to allow |
---|
617 | & (/ undef, undef, 90.0, undef, undef, 90.0, undef, & !! senescence g (days) |
---|
618 | & 60.0, 60.0, 30.0, 30.0, 30.0, 30.0 /) |
---|
619 | |
---|
620 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_c_mtc = & !! critical temperature for senescence (C) |
---|
621 | & (/ undef, undef, undef, undef, undef, 16.0, undef, & !! constant c of aT^2+bT+c, tabulated |
---|
622 | & 14.0, 10.0, 5.0, 5.0, 5.0, 10.0 /) !! (unitless) |
---|
623 | |
---|
624 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_b_mtc = & !! critical temperature for senescence (C), |
---|
625 | & (/ undef, undef, undef, undef, undef, 0.0, undef, & !! constant b of aT^2+bT+c, tabulated |
---|
626 | & 0.0, 0.0, 0.1, 0.0, 0.0, 0.0 /) !! (unitless) |
---|
627 | |
---|
628 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_a_mtc = & !! critical temperature for senescence (C), |
---|
629 | & (/ undef, undef, undef, undef, undef, 0.0, undef, & !! constant a of aT^2+bT+c, tabulated |
---|
630 | & 0.0, 0.0, 0.00375, 0.0, 0.0, 0.0 /) !! (unitless) |
---|
631 | |
---|
632 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gdd_senescence_mtc = & !! minimum gdd to allow senescence of crops (days) |
---|
633 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
634 | & undef, undef, undef, undef, 950., 4000. /) |
---|
635 | |
---|
636 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: always_init_mtc = & !! take carbon from atmosphere if carbohydrate reserve too small (true/false) |
---|
637 | & (/ .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., &!! default is true for all pfts except pft=11 C4 grass |
---|
638 | & .TRUE., .TRUE., .TRUE., .FALSE., .TRUE., .TRUE. /) |
---|
639 | |
---|
640 | |
---|
641 | ! |
---|
642 | ! DGVM |
---|
643 | ! |
---|
644 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: residence_time_mtc = & !! residence time of trees (years) |
---|
645 | & (/ undef, 30.0, 30.0, 40.0, 40.0, 40.0, 80.0, & |
---|
646 | & 80.0, 80.0, 0.0, 0.0, 0.0, 0.0 /) |
---|
647 | |
---|
648 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tmin_crit_mtc = & |
---|
649 | & (/ undef, 0.0, 0.0, -30.0, -14.0, -30.0, -45.0, & !! critical tmin, tabulated (C) |
---|
650 | & -45.0, -60.0, undef, undef, undef, undef /) |
---|
651 | |
---|
652 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tcm_crit_mtc = & |
---|
653 | & (/ undef, undef, undef, 5.0, 15.5, 15.5, -8.0, & !! critical tcm, tabulated (C) |
---|
654 | & -8.0, -8.0, undef, undef, undef, undef /) |
---|
655 | |
---|
656 | |
---|
657 | |
---|
658 | ! |
---|
659 | ! Biogenic Volatile Organic Compounds |
---|
660 | ! |
---|
661 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_isoprene_mtc = & !! Isoprene emission factor |
---|
662 | & (/ 0., 24., 24., 8., 16., 45., 8., & !! |
---|
663 | & 18., 0.5, 12., 18., 5., 5. /) |
---|
664 | |
---|
665 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_monoterpene_mtc = & !! Monoterpene emission factor |
---|
666 | & (/ 0., 2.0, 2.0, 1.8, 1.4, 1.6, 1.8, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
667 | & 1.4, 1.8, 0.8, 0.8, 0.22, 0.22 /) |
---|
668 | |
---|
669 | REAL(r_std), PARAMETER :: LDF_mono_mtc = 0.6 !! monoterpenes fraction dependancy to light |
---|
670 | REAL(r_std), PARAMETER :: LDF_sesq_mtc = 0.5 !! sesquiterpenes fraction dependancy to light |
---|
671 | REAL(r_std), PARAMETER :: LDF_meth_mtc = 0.8 !! methanol fraction dependancy to light |
---|
672 | REAL(r_std), PARAMETER :: LDF_acet_mtc = 0.2 !! acetone fraction dependancy to light |
---|
673 | |
---|
674 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_apinene_mtc = & !! Alfa pinene emission factor percentage |
---|
675 | & (/ 0., 0.395, 0.395, 0.354, 0.463, 0.326, 0.354, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
676 | & 0.316, 0.662, 0.231, 0.200, 0.277, 0.277 /) |
---|
677 | |
---|
678 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_bpinene_mtc = & !! Beta pinene emission factor percentage |
---|
679 | & (/ 0., 0.110, 0.110, 0.146, 0.122, 0.087, 0.146, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
680 | & 0.063, 0.150, 0.123, 0.080, 0.154, 0.154 /) |
---|
681 | |
---|
682 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_limonene_mtc = & !! Limonene emission factor percentage |
---|
683 | & (/ 0., 0.092, 0.092, 0.083, 0.122, 0.061, 0.083, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
684 | & 0.071, 0.037, 0.146, 0.280, 0.092, 0.092 /) |
---|
685 | |
---|
686 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_myrcene_mtc = & !! Myrcene emission factor percentage |
---|
687 | & (/ 0., 0.073, 0.073, 0.050, 0.054, 0.028, 0.050, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
688 | & 0.019, 0.025, 0.062, 0.057, 0.046, 0.046 /) |
---|
689 | |
---|
690 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sabinene_mtc = & !! Sabinene emission factor percentage |
---|
691 | & (/ 0., 0.073, 0.073, 0.050, 0.083, 0.304, 0.050, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
692 | & 0.263, 0.030, 0.065, 0.050, 0.062, 0.062 /) |
---|
693 | |
---|
694 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_camphene_mtc = & !! Camphene emission factor percentage |
---|
695 | & (/ 0., 0.055, 0.055, 0.042, 0.049, 0.004, 0.042, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
696 | & 0.005, 0.023, 0.054, 0.053, 0.031, 0.031 /) |
---|
697 | |
---|
698 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_3carene_mtc = & !! 3-carene emission factor percentage |
---|
699 | & (/ 0., 0.048, 0.048, 0.175, 0.010, 0.024, 0.175, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
700 | & 0.013, 0.042, 0.065, 0.057, 0.200, 0.200 /) |
---|
701 | |
---|
702 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_tbocimene_mtc = & !! T-beta-ocimene emission factor percentage |
---|
703 | & (/ 0., 0.092, 0.092, 0.054, 0.044, 0.113, 0.054, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
704 | & 0.105, 0.028, 0.138, 0.120, 0.031, 0.031 /) |
---|
705 | |
---|
706 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_othermonot_mtc = & !! Other monoterpenes emission factor percentage |
---|
707 | & (/ 0., 0.062, 0.062, 0.046, 0.054, 0.052, 0.046, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
708 | & 0.144, 0.003, 0.115, 0.103, 0.108, 0.108 /) |
---|
709 | |
---|
710 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sesquiterp_mtc = & !! Sesquiterpene emission factor |
---|
711 | & (/ 0., 0.45, 0.45, 0.13, 0.30, 0.36, 0.15, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
712 | & 0.30, 0.25, 0.60, 0.60, 0.08, 0.08 /) |
---|
713 | |
---|
714 | REAL(r_std), PARAMETER :: beta_mono_mtc = 0.10 !! Monoterpenes temperature dependency coefficient |
---|
715 | REAL(r_std), PARAMETER :: beta_sesq_mtc = 0.17 !! Sesquiterpenes temperature dependency coefficient |
---|
716 | REAL(r_std), PARAMETER :: beta_meth_mtc = 0.08 !! Methanol temperature dependency coefficient |
---|
717 | REAL(r_std), PARAMETER :: beta_acet_mtc = 0.10 !! Acetone temperature dependency coefficient |
---|
718 | REAL(r_std), PARAMETER :: beta_oxyVOC_mtc = 0.13 !! Other oxygenated BVOC temperature dependency coefficient |
---|
719 | |
---|
720 | |
---|
721 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_ORVOC_mtc = & !! ORVOC emissions factor |
---|
722 | & (/ 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
723 | & 1.5, 1.5, 1.5, 1.5, 1.5, 1.5 /) |
---|
724 | |
---|
725 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_OVOC_mtc = & !! OVOC emissions factor |
---|
726 | & (/ 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
727 | & 1.5, 1.5, 1.5, 1.5, 1.5, 1.5 /) |
---|
728 | |
---|
729 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_MBO_mtc = & !! MBO emissions factor |
---|
730 | & (/ 0., 2.e-5, 2.e-5, 1.4, 2.e-5, 2.e-5, 0.14, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
731 | & 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5 /) |
---|
732 | |
---|
733 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_methanol_mtc = & !! Methanol emissions factor |
---|
734 | & (/ 0., 0.8, 0.8, 1.8, 0.9, 1.9, 1.8, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
735 | & 1.8, 1.8, 0.7, 0.9, 2., 2. /) |
---|
736 | |
---|
737 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetone_mtc = & !! Acetone emissions factor |
---|
738 | & (/ 0., 0.25, 0.25, 0.30, 0.20, 0.33, 0.30, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
739 | & 0.25, 0.25, 0.20, 0.20, 0.08, 0.08 /) |
---|
740 | |
---|
741 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetal_mtc = & !! Acetaldehyde emissions factor |
---|
742 | & (/ 0., 0.2, 0.2, 0.2, 0.2, 0.25, 0.25, 0.16, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
743 | & 0.16, 0.12, 0.12, 0.035, 0.020 /) |
---|
744 | |
---|
745 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formal_mtc = & !! Formaldehyde emissions factor |
---|
746 | & (/ 0., 0.04, 0.04, 0.08, 0.04, 0.04, 0.04, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
747 | & 0.04, 0.04, 0.025, 0.025, 0.013, 0.013 /) |
---|
748 | |
---|
749 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetic_mtc = & !! Acetic Acid emissions factor |
---|
750 | & (/ 0., 0.025, 0.025, 0.025, 0.022, 0.08, 0.025, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
751 | & 0.022, 0.013, 0.012, 0.012, 0.008, 0.008 /) |
---|
752 | |
---|
753 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formic_mtc = & !! Formic Acid emissions factor |
---|
754 | & (/ 0., 0.015, 0.015, 0.02, 0.02, 0.025, 0.025, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
755 | & 0.015, 0.015, 0.010, 0.010, 0.008, 0.008 /) |
---|
756 | |
---|
757 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_wet_mtc = & !! NOx emissions factor soil emissions and exponential |
---|
758 | & (/ 0., 2.6, 0.06, 0.03, 0.03, 0.03, 0.03, & !! dependancy factor for wet soils |
---|
759 | & 0.03, 0.03, 0.36, 0.36, 0.36, 0.36 /) !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex |
---|
760 | |
---|
761 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_dry_mtc = & !! NOx emissions factor soil emissions and exponential |
---|
762 | & (/ 0., 8.60, 0.40, 0.22, 0.22, 0.22, 0.22, & !! dependancy factor for dry soils |
---|
763 | & 0.22, 0.22, 2.65, 2.65, 2.65, 2.65 /) !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex |
---|
764 | |
---|
765 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Larch_mtc = & !! Larcher 1991 SAI/LAI ratio (unitless) |
---|
766 | & (/ 0., 0.015, 0.015, 0.003, 0.005, 0.005, 0.003, & |
---|
767 | & 0.005, 0.003, 0.005, 0.005, 0.008, 0.008 /) |
---|
768 | |
---|
769 | |
---|
770 | |
---|
771 | END MODULE constantes_mtc |
---|