1 | ! ================================================================================================================================= |
---|
2 | ! MODULE : constantes_mtc |
---|
3 | ! |
---|
4 | ! CONTACT : orchidee-help _at_ ipsl.jussieu.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 | !! |
---|
34 | !! SVN : |
---|
35 | !! $HeadURL: $ |
---|
36 | !! $Date$ |
---|
37 | !! $Revision$ |
---|
38 | !! \n |
---|
39 | !_ ================================================================================================================================ |
---|
40 | |
---|
41 | MODULE constantes_mtc |
---|
42 | |
---|
43 | USE defprec |
---|
44 | USE constantes |
---|
45 | |
---|
46 | IMPLICIT NONE |
---|
47 | |
---|
48 | ! |
---|
49 | ! METACLASSES CHARACTERISTICS |
---|
50 | ! |
---|
51 | |
---|
52 | INTEGER(i_std), PARAMETER :: nvmc = 13 !! Number of MTCS fixed in the code (unitless) |
---|
53 | |
---|
54 | CHARACTER(len=34), PARAMETER, DIMENSION(nvmc) :: MTC_name = & !! description of the MTC (unitless) |
---|
55 | & (/ 'bare ground ', & ! 1 |
---|
56 | & 'tropical broad-leaved evergreen ', & ! 2 |
---|
57 | & 'tropical broad-leaved raingreen ', & ! 3 |
---|
58 | & 'temperate needleleaf evergreen ', & ! 4 |
---|
59 | & 'temperate broad-leaved evergreen ', & ! 5 |
---|
60 | & 'temperate broad-leaved summergreen', & ! 6 |
---|
61 | & 'boreal needleleaf evergreen ', & ! 7 |
---|
62 | & 'boreal broad-leaved summergreen', & ! 8 |
---|
63 | & 'boreal needleleaf summergreen', & ! 9 |
---|
64 | & ' C3 grass ', & ! 10 |
---|
65 | & ' C4 grass ', & ! 11 |
---|
66 | & ' C3 agriculture', & ! 12 |
---|
67 | & ' C4 agriculture' /) ! 13 |
---|
68 | |
---|
69 | |
---|
70 | ! |
---|
71 | ! VEGETATION STRUCTURE |
---|
72 | ! |
---|
73 | INTEGER(i_std),PARAMETER, DIMENSION(nvmc) :: leaf_tab_mtc = & !! leaf type (1-4, unitless) |
---|
74 | & (/ 4, 1, 1, 2, 1, 1, 2, & !! 1=broad leaved tree, 2=needle leaved tree |
---|
75 | & 1, 2, 3, 3, 3, 3 /) !! 3=grass 4=bare ground |
---|
76 | |
---|
77 | CHARACTER(len=6), PARAMETER, DIMENSION(nvmc) :: pheno_model_mtc = & !! which phenology model is used? (tabulated) |
---|
78 | & (/ 'none ', 'none ', 'moi ', 'none ', 'none ', & |
---|
79 | & 'ncdgdd', 'none ', 'ncdgdd', 'ngd ', 'moigdd', & |
---|
80 | & 'moigdd', 'moigdd', 'moigdd' /) |
---|
81 | |
---|
82 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_tropical_mtc = & !! Is PFT tropical ? (true/false) |
---|
83 | & (/ .FALSE., .TRUE., .TRUE., .FALSE., .FALSE., .FALSE., .FALSE., & |
---|
84 | & .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE. /) |
---|
85 | |
---|
86 | CHARACTER(LEN=5), PARAMETER, DIMENSION(nvmc) :: type_of_lai_mtc = & !! Type of behaviour of the LAI evolution algorithm |
---|
87 | & (/ 'inter', 'inter', 'inter', 'inter', 'inter', & !! for each vegetation type. (unitless) |
---|
88 | & 'inter', 'inter', 'inter', 'inter', 'inter', & !! Value of type_of_lai : mean or interp |
---|
89 | & 'inter', 'inter', 'inter' /) |
---|
90 | |
---|
91 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: natural_mtc = & !! natural? (true/false) |
---|
92 | & (/ .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., .TRUE., & |
---|
93 | & .TRUE., .TRUE., .TRUE., .TRUE., .FALSE., .FALSE. /) |
---|
94 | |
---|
95 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: veget_ori_fixed_mtc = & !! Value for veget_ori for tests in |
---|
96 | & (/ 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! 0-dim simulations (0-1, unitless) |
---|
97 | & 0.0, 0.0, 0.8, 0.0, 0.0, 0.0 /) |
---|
98 | |
---|
99 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimax_mtc = & !! laimax for maximum |
---|
100 | & (/ 0.0, 8.0, 8.0, 4.0, 4.5, 4.5, 4.0, & !! See also type of lai interpolation |
---|
101 | & 4.5, 4.0, 2.0, 2.0, 2.0, 2.0 /) !! @tex $(m^2.m^{-2})$ @endtex |
---|
102 | |
---|
103 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimin_mtc = & !! laimin for minimum lai |
---|
104 | & (/ 0.0, 8.0, 0.0, 4.0, 4.5, 0.0, 4.0, & !! See also type of lai interpolation (m^2.m^{-2}) |
---|
105 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! @tex $(m^2.m^{-2})$ @endtex |
---|
106 | |
---|
107 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: height_presc_mtc = & !! prescribed height of vegetation (m) |
---|
108 | & (/ 0.0, 30.0, 30.0, 20.0, 20.0, 20.0, 15.0, & !! Value for height_presc : one for each vegetation type |
---|
109 | & 15.0, 15.0, 0.5, 0.6, 1.0, 1.0 /) |
---|
110 | |
---|
111 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rveg_mtc = & !! Potentiometer to set vegetation resistance (unitless) |
---|
112 | & (/ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & !! Nathalie on March 28th, 2006 - from Fred Hourdin, |
---|
113 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
114 | |
---|
115 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: sla_mtc = & !! specif leaf area @tex $(m^2.gC^{-1})$ @endtex |
---|
116 | & (/ 1.5E-2, 1.53E-2, 2.6E-2, 9.26E-3, 2E-2, 2.6E-2, 9.26E-3, & |
---|
117 | & 2.6E-2, 1.9E-2, 2.6E-2, 2.6E-2, 2.6E-2, 2.6E-2 /) |
---|
118 | |
---|
119 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: availability_fact_mtc = & !! calculate mortality in lpj_gap |
---|
120 | & (/ undef, 0.14, 0.14, 0.10, 0.10, 0.10, 0.05, & |
---|
121 | & 0.05, 0.05, undef, undef, undef, undef /) |
---|
122 | |
---|
123 | ! |
---|
124 | ! EVAPOTRANSPIRATION (sechiba) |
---|
125 | ! |
---|
126 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rstruct_const_mtc = & !! Structural resistance. |
---|
127 | & (/ 0.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, & !! @tex $(s.m^{-1})$ @endtex |
---|
128 | & 25.0, 25.0, 2.5, 2.0, 2.0, 2.0 /) |
---|
129 | |
---|
130 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kzero_mtc = & !! A vegetation dependent constant used in the |
---|
131 | & (/ 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. |
---|
132 | & 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 |
---|
133 | |
---|
134 | |
---|
135 | ! |
---|
136 | ! WATER (sechiba) |
---|
137 | ! |
---|
138 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: wmax_veg_mtc = & !! Volumetric available soil water capacity in each PFT |
---|
139 | & (/ 150.0, 150.0, 150.0, 150.0, 150.0, 150.0, 150.0, & !! @tex $(kg.m^{-3} of soil)$ @endtex |
---|
140 | & 150.0, 150.0, 150.0, 150.0, 150.0, 150.0 /) |
---|
141 | |
---|
142 | |
---|
143 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref4m = & !! Root profile description for the different |
---|
144 | & (/ 5.0, 0.4, 0.4, 1.0, 0.8, 0.8, 1.0, & !! vegetations types. @tex $(m^{-1})$ @endtex |
---|
145 | & 1.0, 0.8, 4.0, 1.0, 4.0, 1.0 /) !! These are the factor in the exponential which gets |
---|
146 | !! the root density as a function of depth |
---|
147 | !! Values for zmaxh = 4.0 |
---|
148 | |
---|
149 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref2m = & !! Root profile description for the different |
---|
150 | & (/ 5.0, 0.8, 0.8, 1.0, 0.8, 0.8, 1.0, & !! vegetations types. @tex $(m^{-1})$ @endtex |
---|
151 | & 1.0, 0.8, 4.0, 4.0, 4.0, 4.0 /) !! These are the factor in the exponential which gets |
---|
152 | !! the root density as a function of depth |
---|
153 | !! Values for zmaxh = 2.0 |
---|
154 | |
---|
155 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: throughfall_by_mtc = & !! Fraction of rain intercepted by the canopy |
---|
156 | & (/ 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, & !! (0-100, unitless) |
---|
157 | & 30.0, 30.0, 30.0, 30.0, 30.0, 30.0 /) |
---|
158 | |
---|
159 | |
---|
160 | ! |
---|
161 | ! ALBEDO (sechiba) |
---|
162 | ! |
---|
163 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_aged_mtc = & !! Minimum snow albedo value for each vegetation type |
---|
164 | & (/ 0.35, 0.0, 0.0, 0.14, 0.14, 0.14, 0.14, & !! after aging (dirty old snow) (unitless) |
---|
165 | & 0.14, 0.14, 0.18, 0.18, 0.18, 0.18 /) !! Source : Values are from the Thesis of S. Chalita (1992) |
---|
166 | |
---|
167 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_dec_mtc = & !! Decay rate of snow albedo value for each vegetation type |
---|
168 | & (/ 0.45, 0.0, 0.0, 0.06, 0.06, 0.11, 0.06, & !! as it will be used in condveg_snow (unitless) |
---|
169 | & 0.11, 0.11, 0.52, 0.52, 0.52, 0.52 /) !! Source : Values are from the Thesis of S. Chalita (1992) |
---|
170 | |
---|
171 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_vis_mtc = & !! leaf albedo of vegetation type, visible albedo |
---|
172 | & (/ 0.00, 0.04, 0.06, 0.06, 0.06, 0.06, 0.06, & !! (unitless) |
---|
173 | & 0.06, 0.06, 0.10, 0.10, 0.10, 0.10 /) |
---|
174 | ! REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_vis_mtc = & !! optimized set of parameter values, to be switched to after validation |
---|
175 | ! & (/ 0.00, 0.037, 0.047, 0.048, 0.06, 0.06, 0.06, & !! (values are different only for PFTs 2,3,4) |
---|
176 | ! & 0.06, 0.06, 0.10, 0.10, 0.10, 0.10 /) |
---|
177 | |
---|
178 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_nir_mtc = & !! leaf albedo of vegetation type, near infrared albedo |
---|
179 | & (/ 0.00, 0.20, 0.22, 0.22, 0.22, 0.22, 0.22, & !! (unitless) |
---|
180 | & 0.22, 0.22, 0.30, 0.30, 0.30, 0.30 /) |
---|
181 | ! REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_nir_mtc = & !! optimized set of parameter values, to be switched to after validation |
---|
182 | ! & (/ 0.00, 0.218, 0.198, 0.208, 0.197, 0.232, 0.2, & !! (for PFTs 2-13) |
---|
183 | ! & 0.217, 0.244, 0.246, 0.28, 0.291, 0.245 /) |
---|
184 | |
---|
185 | |
---|
186 | ! |
---|
187 | ! SOIL - VEGETATION |
---|
188 | ! |
---|
189 | INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pref_soil_veg_mtc = & !! The soil tile number for each vegetation |
---|
190 | & (/ 1, 2, 2, 2, 2, 2, 2, & |
---|
191 | & 2, 2, 3, 3, 3, 3 /) |
---|
192 | |
---|
193 | |
---|
194 | ! |
---|
195 | ! PHOTOSYNTHESIS |
---|
196 | ! |
---|
197 | !- |
---|
198 | ! 1 .CO2 |
---|
199 | !- |
---|
200 | LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_c4_mtc = & !! flag for C4 vegetation types (true/false) |
---|
201 | & (/ .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., & |
---|
202 | & .FALSE., .FALSE., .FALSE., .TRUE., .FALSE., .TRUE. /) |
---|
203 | |
---|
204 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: vcmax_fix_mtc = & !! values used for vcmax when STOMATE is not |
---|
205 | & (/ 0.0, 40.0, 50.0, 30.0, 35.0, 40.0, 30.0, & !! activated @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex |
---|
206 | & 40.0, 35.0, 60.0, 60.0, 70.0, 70.0 /) |
---|
207 | |
---|
208 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: downregulation_co2_coeff_mtc = & !! coefficient for CO2 downregulation |
---|
209 | & (/ 0.0, 0.38, 0.38, 0.28, 0.28, 0.28, 0.22, & |
---|
210 | & 0.22, 0.22, 0.26, 0.26, 0.26, 0.26 /) |
---|
211 | |
---|
212 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmC_mtc = & !! Energy of activation for KmC (J mol-1) |
---|
213 | & (/undef, 79430., 79430., 79430., 79430., 79430., 79430., & !! See Medlyn et al. (2002) |
---|
214 | & 79430., 79430., 79430., 79430., 79430., 79430. /) !! from Bernacchi al. (2001) |
---|
215 | |
---|
216 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmO_mtc = & !! Energy of activation for KmO (J mol-1) |
---|
217 | & (/undef, 36380., 36380., 36380., 36380., 36380., 36380., & !! See Medlyn et al. (2002) |
---|
218 | & 36380., 36380., 36380., 36380., 36380., 36380. /) !! from Bernacchi al. (2001) |
---|
219 | |
---|
220 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_gamma_star_mtc = & !! Energy of activation for gamma_star (J mol-1) |
---|
221 | & (/undef, 37830., 37830., 37830., 37830., 37830., 37830., & !! See Medlyn et al. (2002) from Bernacchi al. (2001) |
---|
222 | & 37830., 37830., 37830., 37830., 37830., 37830. /) !! for C3 plants - We use the same values for C4 plants |
---|
223 | |
---|
224 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Vcmax_mtc = & !! Energy of activation for Vcmax (J mol-1) |
---|
225 | & (/undef, 71513., 71513., 71513., 71513., 71513., 71513., & !! See Table 2 of Yin et al. (2009) for C4 plants |
---|
226 | & 71513., 71513., 71513., 67300., 71513., 67300. /) !! and Kattge & Knorr (2007) for C3 plants (table 3) |
---|
227 | |
---|
228 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Jmax_mtc = & !! Energy of activation for Jmax (J mol-1) |
---|
229 | & (/undef, 49884., 49884., 49884., 49884., 49884., 49884., & !! See Table 2 of Yin et al. (2009) for C4 plants |
---|
230 | & 49884., 49884., 49884., 77900., 49884., 77900. /) !! and Kattge & Knorr (2007) for C3 plants (table 3) |
---|
231 | |
---|
232 | 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) |
---|
233 | & (/undef, 668.39, 668.39, 668.39, 668.39, 668.39, 668.39, & !! See Table 3 of Kattge & Knorr (2007) |
---|
234 | & 668.39, 668.39, 668.39, 641.64, 668.39, 641.64 /) !! For C4 plants, we assume that there is no |
---|
235 | !! 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) |
---|
236 | |
---|
237 | 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) |
---|
238 | & (/undef, -1.07, -1.07, -1.07, -1.07, -1.07, -1.07, & !! See Table 3 of Kattge & Knorr (2007) |
---|
239 | & -1.07, -1.07, -1.07, 0., -1.07, 0. /) !! We assume No acclimation term for C4 plants |
---|
240 | |
---|
241 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_min_mtc = & !! minimum photosynthesis temperature (deg C) |
---|
242 | & (/ undef, -4.0, -4.0, -4.0, -4.0, -4.0, -4.0, & |
---|
243 | & -4.0, -4.0, -4.0, -4.0, -4.0, -4.0 /) |
---|
244 | |
---|
245 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_max_mtc = & !! maximum photosynthesis temperature (deg C) |
---|
246 | & (/ undef, 55.0, 55.0, 55.0, 55.0, 55.0, 55.0, & |
---|
247 | & 55.0, 55.0, 55.0, 55.0, 55.0, 55.0 /) |
---|
248 | |
---|
249 | 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) |
---|
250 | & (/undef, 659.70, 659.70, 659.70, 659.70, 659.70, 659.70, & !! See Table 3 of Kattge & Knorr (2007) |
---|
251 | & 659.70, 659.70, 659.70, 630., 659.70, 630. /) !! and Table 2 of Yin et al. (2009) for C4 plants |
---|
252 | |
---|
253 | 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) |
---|
254 | & (/undef, -0.75, -0.75, -0.75, -0.75, -0.75, -0.75, & !! See Table 3 of Kattge & Knorr (2007) |
---|
255 | & -0.75, -0.75, -0.75, 0., -0.75, 0. /) !! We assume no acclimation term for C4 plants |
---|
256 | |
---|
257 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Vcmax_mtc = & !! Energy of deactivation for Vcmax (J mol-1) |
---|
258 | & (/undef, 200000., 200000., 200000., 200000., 200000., 200000., & !! Medlyn et al. (2002) also uses 200000. for C3 plants (same value than D_Jmax) |
---|
259 | & 200000., 200000., 200000., 192000., 200000., 192000. /) !! 'Consequently', we use the value of D_Jmax for C4 plants |
---|
260 | |
---|
261 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Jmax_mtc = & !! Energy of deactivation for Jmax (J mol-1) |
---|
262 | & (/undef, 200000., 200000., 200000., 200000., 200000., 200000., & !! See Table 2 of Yin et al. (2009) |
---|
263 | & 200000., 200000., 200000., 192000., 200000., 192000. /) !! Medlyn et al. (2002) also uses 200000. for C3 plants |
---|
264 | |
---|
265 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Rd_mtc = & !! Energy of activation for Rd (J mol-1) |
---|
266 | & (/undef, 46390., 46390., 46390., 46390., 46390., 46390., & !! See Table 2 of Yin et al. (2009) |
---|
267 | & 46390., 46390., 46390., 46390., 46390., 46390. /) |
---|
268 | |
---|
269 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Vcmax25_mtc = & !! Maximum rate of Rubisco activity-limited carboxylation at 25°C |
---|
270 | & (/ undef, 65.0, 65.0, 45.0, 45.0, 55.0, 45.0, & !! @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex |
---|
271 | & 45.0, 35.0, 70.0, 70.0, 70.0, 70.0 /) |
---|
272 | |
---|
273 | 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) |
---|
274 | & (/undef, 2.59, 2.59, 2.59, 2.59, 2.59, 2.59, & !! See Table 3 of Kattge & Knorr (2007) |
---|
275 | & 2.59, 2.59, 2.59, 1.715, 2.59, 1.715 /) !! For C4 plants, we assume that there is no |
---|
276 | !! 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) |
---|
277 | |
---|
278 | 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) |
---|
279 | & (/undef, -0.035, -0.035, -0.035, -0.035, -0.035, -0.035, & !! See Table 3 of Kattge & Knorr (2007) |
---|
280 | & -0.035, -0.035, -0.035, 0., -0.035, 0. /) !! We assume No acclimation term for C4 plants |
---|
281 | |
---|
282 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmC25_mtc = & !! MichaelisâMenten constant of Rubisco for CO2 at 25°C (ubar) |
---|
283 | & (/undef, 404.9, 404.9, 404.9, 404.9, 404.9, 404.9, & !! See Table 2 of Yin et al. (2009) for C4 |
---|
284 | & 404.9, 404.9, 404.9, 650., 404.9, 650. /) !! and Medlyn et al (2002) for C3 |
---|
285 | |
---|
286 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmO25_mtc = & !! MichaelisâMenten constant of Rubisco for O2 at 25°C (ubar) |
---|
287 | & (/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 |
---|
288 | & 278400., 278400., 278400., 450000., 278400., 450000. /) |
---|
289 | |
---|
290 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gamma_star25_mtc = & !! Ci-based CO2 compensation point in the absence of Rd at 25°C (ubar) |
---|
291 | & (/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) |
---|
292 | & 42.75, 42.75, 42.75, 42.75, 42.75, 42.75 /) |
---|
293 | |
---|
294 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: a1_mtc = & !! Empirical factor involved in the calculation of fvpd (-) |
---|
295 | & (/undef, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, & !! See Table 2 of Yin et al. (2009) |
---|
296 | & 0.85, 0.85, 0.85, 0.85, 0.85, 0.85 /) |
---|
297 | |
---|
298 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: b1_mtc = & !! Empirical factor involved in the calculation of fvpd (-) |
---|
299 | & (/undef, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, & !! See Table 2 of Yin et al. (2009) |
---|
300 | & 0.14, 0.14, 0.14, 0.20, 0.14, 0.20 /) |
---|
301 | |
---|
302 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: g0_mtc = & !! Residual stomatal conductance when irradiance approaches zero (mol CO2 mâ2 sâ1 barâ1) |
---|
303 | & (/undef, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, & !! Value from ORCHIDEE - No other reference. |
---|
304 | & 0.00625, 0.00625, 0.00625, 0.01875, 0.00625, 0.01875 /) !! modofy to account for the conversion for conductance to H2O to CO2 |
---|
305 | |
---|
306 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: h_protons_mtc = & !! Number of protons required to produce one ATP (mol mol-1) |
---|
307 | & (/undef, 4., 4., 4., 4., 4., 4., & !! See Table 2 of Yin et al. (2009) - h parameter |
---|
308 | & 4., 4., 4., 4., 4., 4. /) |
---|
309 | |
---|
310 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpsir_mtc = & !! Fraction of PSII eâ transport rate |
---|
311 | & (/undef, undef, undef, undef, undef, undef, undef, & !! partitioned to the C4 cycle (-) |
---|
312 | & undef, undef, undef, 0.4, undef, 0.4 /) !! See Table 2 of Yin et al. (2009) - x parameter |
---|
313 | |
---|
314 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fQ_mtc = & !! Fraction of electrons at reduced plastoquinone |
---|
315 | & (/undef, undef, undef, undef, undef, undef, undef, & !! that follow the Q-cycle (-) - Values for C3 platns are not used |
---|
316 | & undef, undef, undef, 1., undef, 1. /) !! See Table 2 of Yin et al. (2009) |
---|
317 | |
---|
318 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpseudo_mtc = & !! Fraction of electrons at PSI that follow |
---|
319 | & (/undef, undef, undef, undef, undef, undef, undef, & !! pseudocyclic transport (-) - Values for C3 platns are not used |
---|
320 | & undef, undef, undef, 0.1, undef, 0.1 /) !! See Table 2 of Yin et al. (2009) |
---|
321 | |
---|
322 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kp_mtc = & !! Initial carboxylation efficiency of the PEP carboxylase (mol mâ2 sâ1 barâ1) |
---|
323 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See Table 2 of Yin et al. (2009) |
---|
324 | & undef, undef, undef, 0.7, undef, 0.7 /) |
---|
325 | |
---|
326 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alpha_mtc = & !! Fraction of PSII activity in the bundle sheath (-) |
---|
327 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See legend of Figure 6 of Yin et al. (2009) |
---|
328 | & undef, undef, undef, 0.1, undef, 0.1 /) |
---|
329 | |
---|
330 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gbs_mtc = & !! Bundle-sheath conductance (mol mâ2 sâ1 barâ1) |
---|
331 | & (/undef, undef, undef, undef, undef, undef, undef, & !! See legend of Figure 6 of Yin et al. (2009) |
---|
332 | & undef, undef, undef, 0.003, undef, 0.003 /) |
---|
333 | |
---|
334 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: theta_mtc = & !! Convexity factor for response of J to irradiance (-) |
---|
335 | & (/undef, 0.7, 0.7, 0.7, 0.7, 0.7, 0.7, & !! See Table 2 of Yin et al. (2009) |
---|
336 | & 0.7, 0.7, 0.7, 0.7, 0.7, 0.7 /) |
---|
337 | |
---|
338 | 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) |
---|
339 | & (/undef, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, & !! See comment from Yin et al. (2009) after eq. 4 |
---|
340 | & 0.3, 0.3, 0.3, 0.3, 0.3, 0.3 /) !! alpha value from Medlyn et al. (2002) |
---|
341 | !! 0.093 mol CO2 fixed per mol absorbed photons |
---|
342 | !! times 4 mol e- per mol CO2 produced |
---|
343 | |
---|
344 | !- |
---|
345 | ! 2 .Stomate |
---|
346 | !- |
---|
347 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ext_coeff_mtc = & !! extinction coefficient of the Monsi&Saeki |
---|
348 | & (/ 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, & !! relationship (1953) (unitless) |
---|
349 | & 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 /) |
---|
350 | |
---|
351 | ! |
---|
352 | ! ALLOCATION (stomate) |
---|
353 | ! |
---|
354 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: R0_mtc = & !! Default root allocation (0-1, unitless) |
---|
355 | & (/ undef, 0.30, 0.30, 0.30, 0.30, 0.30, 0.30, & |
---|
356 | & 0.30, 0.30, 0.30, 0.30, 0.30, 0.30 /) |
---|
357 | |
---|
358 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: S0_mtc = & !! Default sapwood allocation (0-1, unitless) |
---|
359 | & (/ undef, 0.25, 0.25, 0.30, 0.30, 0.30, 0.30, & |
---|
360 | & 0.30, 0.30, 0.30, 0.30, 0.30, 0.30 /) |
---|
361 | |
---|
362 | ! |
---|
363 | ! RESPIRATION (stomate) |
---|
364 | ! |
---|
365 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: frac_growthresp_mtc = & !! fraction of GPP which is lost as growth respiration |
---|
366 | & (/ undef, 0.28, 0.28, 0.28, 0.28, 0.28, 0.28, & |
---|
367 | & 0.28, 0.28, 0.28, 0.28, 0.28, 0.28 /) |
---|
368 | |
---|
369 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_c_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
370 | & (/ undef, 0.20, 0.20, 0.16, 0.16, 0.16, 0.16, & !! constant c of aT^2+bT+c, tabulated |
---|
371 | & 0.16, 0.16, 0.16, 0.12, 0.16, 0.12 /) |
---|
372 | |
---|
373 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_b_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
374 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! constant b of aT^2+bT+c, tabulated |
---|
375 | & 0.0, 0.0, -0.00133, 0.0, -0.00133, 0.0 /) |
---|
376 | |
---|
377 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_a_mtc = & !! slope of maintenance respiration coefficient (1/K), |
---|
378 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! constant a of aT^2+bT+c, tabulated |
---|
379 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) |
---|
380 | |
---|
381 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_leaf_mtc = & !! maintenance respiration coefficient |
---|
382 | & (/ 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, |
---|
383 | & 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 |
---|
384 | |
---|
385 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapabove_mtc = & !! maintenance respiration coefficient |
---|
386 | & (/ 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, |
---|
387 | & 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 |
---|
388 | |
---|
389 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapbelow_mtc = & !! maintenance respiration coefficient |
---|
390 | & (/ 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, |
---|
391 | & 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 |
---|
392 | |
---|
393 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartabove_mtc = & !! maintenance respiration coefficient |
---|
394 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! at 0 deg C, for heartwood above, |
---|
395 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
396 | |
---|
397 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartbelow_mtc = & !! maintenance respiration coefficient |
---|
398 | & (/ undef, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & !! at 0 deg C, for heartwood below, |
---|
399 | & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /) !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex |
---|
400 | |
---|
401 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_root_mtc = & !! maintenance respiration coefficient |
---|
402 | & (/ 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, |
---|
403 | & 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 |
---|
404 | |
---|
405 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_fruit_mtc = & !! maintenance respiration coefficient |
---|
406 | & (/ 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, |
---|
407 | & 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 |
---|
408 | |
---|
409 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_carbres_mtc = & !! maintenance respiration coefficient |
---|
410 | & (/ 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, |
---|
411 | & 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 |
---|
412 | |
---|
413 | |
---|
414 | ! |
---|
415 | ! FIRE (stomate) |
---|
416 | ! |
---|
417 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: flam_mtc = & !! flamability: critical fraction of water |
---|
418 | & (/ undef, 0.15, 0.25, 0.25, 0.25, 0.25, 0.25, & !! holding capacity (0-1, unitless) |
---|
419 | & 0.25, 0.25, 0.25, 0.25, 0.35, 0.35 /) |
---|
420 | |
---|
421 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: resist_mtc = & !! fire resistance (0-1, unitless) |
---|
422 | & (/ undef, 0.95, 0.90, 0.90, 0.90, 0.90, 0.90, & |
---|
423 | & 0.90, 0.90, 0.0, 0.0, 0.0, 0.0 /) |
---|
424 | |
---|
425 | |
---|
426 | ! |
---|
427 | ! FLUX - LUC |
---|
428 | ! |
---|
429 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_1_mtc = & !! Coeff of biomass export for the year |
---|
430 | & (/ undef, 0.897, 0.897, 0.597, 0.597, 0.597, 0.597, & !! (0-1, unitless) |
---|
431 | & 0.597, 0.597, 0.597, 0.597, 0.597, 0.597 /) |
---|
432 | |
---|
433 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_10_mtc = & !! Coeff of biomass export for the decade |
---|
434 | & (/ undef, 0.103, 0.103, 0.299, 0.299, 0.299, 0.299, & !! (0-1, unitless) |
---|
435 | & 0.299, 0.299, 0.299, 0.403, 0.299, 0.403 /) |
---|
436 | |
---|
437 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_100_mtc = & !! Coeff of biomass export for the century |
---|
438 | & (/ undef, 0.0, 0.0, 0.104, 0.104, 0.104, 0.104, & !! (0-1, unitless) |
---|
439 | & 0.104, 0.104, 0.104, 0.0, 0.104, 0.0 /) |
---|
440 | |
---|
441 | |
---|
442 | ! |
---|
443 | ! PHENOLOGY |
---|
444 | ! |
---|
445 | !- |
---|
446 | ! 1. Stomate |
---|
447 | !- |
---|
448 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: lai_max_to_happy_mtc = & !! threshold of LAI below which plant uses carbohydrate reserves |
---|
449 | & (/ undef, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, & |
---|
450 | & 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 /) |
---|
451 | |
---|
452 | REAL(r_std), PARAMETER, DIMENSION (nvmc) :: lai_max_mtc = & !! maximum LAI, PFT-specific |
---|
453 | & (/ undef, 7.0, 7.0, 5.0, 5.0, 5.0, 4.5, & !! @tex $(m^2.m^{-2})$ @endtex |
---|
454 | & 4.5, 3.0, 2.5, 2.5, 5.0, 5.0 /) |
---|
455 | |
---|
456 | INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pheno_type_mtc = & !! type of phenology (0-4, unitless) |
---|
457 | & (/ 0, 1, 3, 1, 1, 2, 1, & !! 0=bare ground 1=evergreen, 2=summergreen, |
---|
458 | & 2, 2, 4, 4, 2, 3 /) !! 3=raingreen, 4=perennial |
---|
459 | !- |
---|
460 | ! 2. Leaf Onset |
---|
461 | !- |
---|
462 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_c_mtc = & !! critical gdd, tabulated (C), |
---|
463 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant c of aT^2+bT+c |
---|
464 | & undef, undef, 320.0, 400.0, 320.0, 700.0 /) |
---|
465 | |
---|
466 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_b_mtc = & !! critical gdd, tabulated (C), |
---|
467 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant b of aT^2+bT+c |
---|
468 | & undef, undef, 6.25, 0.0, 6.25, 0.0 /) |
---|
469 | |
---|
470 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_a_mtc = & !! critical gdd, tabulated (C), |
---|
471 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! constant a of aT^2+bT+c |
---|
472 | & undef, undef, 0.03125, 0.0, 0.0315, 0.0 /) |
---|
473 | |
---|
474 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_moigdd_t_crit_mtc = & !! temperature threshold for C4 grass(C) |
---|
475 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
476 | & undef, undef, undef, 22.0, undef, undef /) |
---|
477 | |
---|
478 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ngd_crit_mtc = & !! critical ngd, tabulated. |
---|
479 | & (/ undef, undef, undef, undef, undef, undef, undef, & !! Threshold -5 degrees (days) |
---|
480 | & undef, 17.0, undef, undef, undef, undef /) |
---|
481 | |
---|
482 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ncdgdd_temp_mtc = & !! critical temperature for the ncd vs. gdd |
---|
483 | & (/ undef, undef, undef, undef, undef, 5.0, undef, & !! function in phenology (C) |
---|
484 | & 0.0, undef, undef, undef, undef, undef /) |
---|
485 | |
---|
486 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_frac_mtc = & !! critical humidity (relative to min/max) |
---|
487 | & (/ undef, undef, 0.5, undef, undef, undef, undef, & !! for phenology (unitless) |
---|
488 | & undef, undef, 0.5, 0.5, 0.5, 0.5 /) |
---|
489 | |
---|
490 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_min_time_mtc = & !! minimum time elapsed since |
---|
491 | & (/ undef, undef, 50.0, undef, undef, undef, undef, & !! moisture minimum (days) |
---|
492 | & undef, undef, 35.0, 35.0, 75.0, 75.0 /) |
---|
493 | |
---|
494 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_sap_mtc = & !! time (days) |
---|
495 | & (/ undef, 730.0, 730.0, 730.0, 730.0, 730.0, 730.0, & |
---|
496 | & 730.0, 730.0, undef, undef, undef, undef /) |
---|
497 | |
---|
498 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_leafinit_mtc = & !! time to attain the initial foliage using the carbohydrate reserve |
---|
499 | & (/ undef, 10., 10., 10., 10., 10., 10., & |
---|
500 | & 10., 10., 10., 10., 10., 10. /) |
---|
501 | |
---|
502 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_fruit_mtc = & !! fruit lifetime (days) |
---|
503 | & (/ undef, 90.0, 90.0, 90.0, 90.0, 90.0, 90.0, & |
---|
504 | & 90.0, 90.0, undef, undef, undef, undef /) |
---|
505 | |
---|
506 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ecureuil_mtc = & !! fraction of primary leaf and root allocation |
---|
507 | & (/ undef, 0.0, 1.0, 0.0, 0.0, 1.0, 0.0, & !! put into reserve (0-1, unitless) |
---|
508 | & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /) |
---|
509 | |
---|
510 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_min_mtc = & !! NEW - allocation above/below = f(age) |
---|
511 | & (/ undef, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, & !! - 30/01/04 NV/JO/PF |
---|
512 | & 0.2, 0.2, undef, undef, undef, undef /) |
---|
513 | |
---|
514 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_max_mtc = & !! NEW - allocation above/below = f(age) |
---|
515 | & (/ undef, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, & !! - 30/01/04 NV/JO/PF |
---|
516 | & 0.8, 0.8, undef, undef, undef, undef /) |
---|
517 | |
---|
518 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: demi_alloc_mtc = & !! NEW - allocation above/below = f(age) |
---|
519 | & (/ undef, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, & !! - 30/01/04 NV/JO/PF |
---|
520 | & 5.0, 5.0, undef, undef, undef, undef /) |
---|
521 | |
---|
522 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaflife_mtc = & !! leaf longevity, tabulated (??units??) |
---|
523 | & (/ undef, 0.5, 2.0, 0.33, 1.0, 2.0, 0.33, & |
---|
524 | & 2.0, 2.0, 2.0, 2.0, 2.0, 2.0 /) |
---|
525 | !- |
---|
526 | ! 3. Senescence |
---|
527 | !- |
---|
528 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaffall_mtc = & !! length of death of leaves, tabulated (days) |
---|
529 | & (/ undef, undef, 10.0, undef, undef, 10.0, undef, & |
---|
530 | & 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 /) |
---|
531 | |
---|
532 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leafagecrit_mtc = & !! critical leaf age, tabulated (days) |
---|
533 | & (/ undef, 730.0, 180.0, 910.0, 730.0, 180.0, 910.0, & |
---|
534 | & 180.0, 180.0, 120.0, 120.0, 90.0, 90.0 /) |
---|
535 | |
---|
536 | CHARACTER(LEN=6), PARAMETER, DIMENSION(nvmc) :: senescence_type_mtc = & !! type of senescence, tabulated (unitless) |
---|
537 | & (/ 'none ', 'none ', 'dry ', 'none ', 'none ', & |
---|
538 | & 'cold ', 'none ', 'cold ', 'cold ', 'mixed ', & |
---|
539 | & 'mixed ', 'mixed ', 'mixed ' /) |
---|
540 | |
---|
541 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_hum_mtc = & !! critical relative moisture availability |
---|
542 | & (/ undef, undef, 0.3, undef, undef, undef, undef, & !! for senescence (0-1, unitless) |
---|
543 | & undef, undef, 0.2, 0.2, 0.3, 0.2 /) |
---|
544 | |
---|
545 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: nosenescence_hum_mtc = & !! relative moisture availability above which |
---|
546 | & (/ undef, undef, 0.8, undef, undef, undef, undef, & !! there is no humidity-related senescence |
---|
547 | & undef, undef, 0.3, 0.3, 0.3, 0.3 /) !! (0-1, unitless) |
---|
548 | |
---|
549 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: max_turnover_time_mtc = & !! maximum turnover time for grasses (days) |
---|
550 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
551 | & undef, undef, 80.0, 80.0, 80.0, 80.0 /) |
---|
552 | |
---|
553 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_turnover_time_mtc = & !! minimum turnover time for grasses (days) |
---|
554 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
555 | & undef, undef, 10.0, 10.0, 10.0, 10.0 /) |
---|
556 | |
---|
557 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_leaf_age_for_senescence_mtc = & !! minimum leaf age to allow |
---|
558 | & (/ undef, undef, 90.0, undef, undef, 90.0, undef, & !! senescence g (days) |
---|
559 | & 60.0, 60.0, 30.0, 30.0, 30.0, 30.0 /) |
---|
560 | |
---|
561 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_c_mtc = & !! critical temperature for senescence (C) |
---|
562 | & (/ undef, undef, undef, undef, undef, 12.0, undef, & !! constant c of aT^2+bT+c, tabulated |
---|
563 | & 7.0, 2.0, -1.375, 5.0, 5.0, 10.0 /) !! (unitless) |
---|
564 | |
---|
565 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_b_mtc = & !! critical temperature for senescence (C), |
---|
566 | & (/ undef, undef, undef, undef, undef, 0.0, undef, & !! constant b of aT^2+bT+c, tabulated |
---|
567 | & 0.0, 0.0, 0.1, 0.0, 0.0, 0.0 /) !! (unitless) |
---|
568 | |
---|
569 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_a_mtc = & !! critical temperature for senescence (C), |
---|
570 | & (/ undef, undef, undef, undef, undef, 0.0, undef, & !! constant a of aT^2+bT+c, tabulated |
---|
571 | & 0.0, 0.0, 0.00375, 0.0, 0.0, 0.0 /) !! (unitless) |
---|
572 | |
---|
573 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gdd_senescence_mtc = & !! minimum gdd to allow senescence of crops (days) |
---|
574 | & (/ undef, undef, undef, undef, undef, undef, undef, & |
---|
575 | & undef, undef, undef, undef, 950., 4000. /) |
---|
576 | |
---|
577 | |
---|
578 | ! |
---|
579 | ! DGVM |
---|
580 | ! |
---|
581 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: residence_time_mtc = & !! residence time of trees (years) |
---|
582 | & (/ undef, 30.0, 30.0, 40.0, 40.0, 40.0, 80.0, & |
---|
583 | & 80.0, 80.0, 0.0, 0.0, 0.0, 0.0 /) |
---|
584 | |
---|
585 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tmin_crit_mtc = & |
---|
586 | & (/ undef, 0.0, 0.0, -30.0, -14.0, -30.0, -45.0, & !! critical tmin, tabulated (C) |
---|
587 | & -45.0, undef, undef, undef, undef, undef /) |
---|
588 | |
---|
589 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tcm_crit_mtc = & |
---|
590 | & (/ undef, undef, undef, 5.0, 15.5, 15.5, -8.0, & !! critical tcm, tabulated (C) |
---|
591 | & -8.0, -8.0, undef, undef, undef, undef /) |
---|
592 | |
---|
593 | |
---|
594 | |
---|
595 | ! |
---|
596 | ! Biogenic Volatile Organic Compounds |
---|
597 | ! |
---|
598 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_isoprene_mtc = & !! Isoprene emission factor |
---|
599 | & (/ 0., 24., 24., 8., 16., 45., 8., & !! |
---|
600 | & 18., 0.5, 12., 18., 5., 5. /) |
---|
601 | |
---|
602 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_monoterpene_mtc = & !! Monoterpene emission factor |
---|
603 | & (/ 0., 2.0, 2.0, 1.8, 1.4, 1.6, 1.8, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
604 | & 1.4, 1.8, 0.8, 0.8, 0.22, 0.22 /) |
---|
605 | |
---|
606 | REAL(r_std), PARAMETER :: LDF_mono_mtc = 0.6 !! monoterpenes fraction dependancy to light |
---|
607 | REAL(r_std), PARAMETER :: LDF_sesq_mtc = 0.5 !! sesquiterpenes fraction dependancy to light |
---|
608 | REAL(r_std), PARAMETER :: LDF_meth_mtc = 0.8 !! methanol fraction dependancy to light |
---|
609 | REAL(r_std), PARAMETER :: LDF_acet_mtc = 0.2 !! acetone fraction dependancy to light |
---|
610 | |
---|
611 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_apinene_mtc = & !! Alfa pinene emission factor percentage |
---|
612 | & (/ 0., 0.395, 0.395, 0.354, 0.463, 0.326, 0.354, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
613 | & 0.316, 0.662, 0.231, 0.200, 0.277, 0.277 /) |
---|
614 | |
---|
615 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_bpinene_mtc = & !! Beta pinene emission factor percentage |
---|
616 | & (/ 0., 0.110, 0.110, 0.146, 0.122, 0.087, 0.146, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
617 | & 0.063, 0.150, 0.123, 0.080, 0.154, 0.154 /) |
---|
618 | |
---|
619 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_limonene_mtc = & !! Limonene emission factor percentage |
---|
620 | & (/ 0., 0.092, 0.092, 0.083, 0.122, 0.061, 0.083, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
621 | & 0.071, 0.037, 0.146, 0.280, 0.092, 0.092 /) |
---|
622 | |
---|
623 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_myrcene_mtc = & !! Myrcene emission factor percentage |
---|
624 | & (/ 0., 0.073, 0.073, 0.050, 0.054, 0.028, 0.050, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
625 | & 0.019, 0.025, 0.062, 0.057, 0.046, 0.046 /) |
---|
626 | |
---|
627 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sabinene_mtc = & !! Sabinene emission factor percentage |
---|
628 | & (/ 0., 0.073, 0.073, 0.050, 0.083, 0.304, 0.050, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
629 | & 0.263, 0.030, 0.065, 0.050, 0.062, 0.062 /) |
---|
630 | |
---|
631 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_camphene_mtc = & !! Camphene emission factor percentage |
---|
632 | & (/ 0., 0.055, 0.055, 0.042, 0.049, 0.004, 0.042, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
633 | & 0.005, 0.023, 0.054, 0.053, 0.031, 0.031 /) |
---|
634 | |
---|
635 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_3carene_mtc = & !! 3-carene emission factor percentage |
---|
636 | & (/ 0., 0.048, 0.048, 0.175, 0.010, 0.024, 0.175, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
637 | & 0.013, 0.042, 0.065, 0.057, 0.200, 0.200 /) |
---|
638 | |
---|
639 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_tbocimene_mtc = & !! T-beta-ocimene emission factor percentage |
---|
640 | & (/ 0., 0.092, 0.092, 0.054, 0.044, 0.113, 0.054, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
641 | & 0.105, 0.028, 0.138, 0.120, 0.031, 0.031 /) |
---|
642 | |
---|
643 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_othermonot_mtc = & !! Other monoterpenes emission factor percentage |
---|
644 | & (/ 0., 0.062, 0.062, 0.046, 0.054, 0.052, 0.046, & !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF |
---|
645 | & 0.144, 0.003, 0.115, 0.103, 0.108, 0.108 /) |
---|
646 | |
---|
647 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sesquiterp_mtc = & !! Sesquiterpene emission factor |
---|
648 | & (/ 0., 0.45, 0.45, 0.13, 0.30, 0.36, 0.15, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
649 | & 0.30, 0.25, 0.60, 0.60, 0.08, 0.08 /) |
---|
650 | |
---|
651 | REAL(r_std), PARAMETER :: beta_mono_mtc = 0.10 !! Monoterpenes temperature dependency coefficient |
---|
652 | REAL(r_std), PARAMETER :: beta_sesq_mtc = 0.17 !! Sesquiterpenes temperature dependency coefficient |
---|
653 | REAL(r_std), PARAMETER :: beta_meth_mtc = 0.08 !! Methanol temperature dependency coefficient |
---|
654 | REAL(r_std), PARAMETER :: beta_acet_mtc = 0.10 !! Acetone temperature dependency coefficient |
---|
655 | REAL(r_std), PARAMETER :: beta_oxyVOC_mtc = 0.13 !! Other oxygenated BVOC temperature dependency coefficient |
---|
656 | |
---|
657 | |
---|
658 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_ORVOC_mtc = & !! ORVOC emissions factor |
---|
659 | & (/ 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
660 | & 1.5, 1.5, 1.5, 1.5, 1.5, 1.5 /) |
---|
661 | |
---|
662 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_OVOC_mtc = & !! OVOC emissions factor |
---|
663 | & (/ 0., 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
664 | & 1.5, 1.5, 1.5, 1.5, 1.5, 1.5 /) |
---|
665 | |
---|
666 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_MBO_mtc = & !! MBO emissions factor |
---|
667 | & (/ 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 |
---|
668 | & 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5 /) |
---|
669 | |
---|
670 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_methanol_mtc = & !! Methanol emissions factor |
---|
671 | & (/ 0., 0.8, 0.8, 1.8, 0.9, 1.9, 1.8, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
672 | & 1.8, 1.8, 0.7, 0.9, 2., 2. /) |
---|
673 | |
---|
674 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetone_mtc = & !! Acetone emissions factor |
---|
675 | & (/ 0., 0.25, 0.25, 0.30, 0.20, 0.33, 0.30, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
676 | & 0.25, 0.25, 0.20, 0.20, 0.08, 0.08 /) |
---|
677 | |
---|
678 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetal_mtc = & !! Acetaldehyde emissions factor |
---|
679 | & (/ 0., 0.2, 0.2, 0.2, 0.2, 0.25, 0.25, 0.16, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
680 | & 0.16, 0.12, 0.12, 0.035, 0.020 /) |
---|
681 | |
---|
682 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formal_mtc = & !! Formaldehyde emissions factor |
---|
683 | & (/ 0., 0.04, 0.04, 0.08, 0.04, 0.04, 0.04, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
684 | & 0.04, 0.04, 0.025, 0.025, 0.013, 0.013 /) |
---|
685 | |
---|
686 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetic_mtc = & !! Acetic Acid emissions factor |
---|
687 | & (/ 0., 0.025, 0.025, 0.025, 0.022, 0.08, 0.025, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
688 | & 0.022, 0.013, 0.012, 0.012, 0.008, 0.008 /) |
---|
689 | |
---|
690 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formic_mtc = & !! Formic Acid emissions factor |
---|
691 | & (/ 0., 0.015, 0.015, 0.02, 0.02, 0.025, 0.025, & !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex |
---|
692 | & 0.015, 0.015, 0.010, 0.010, 0.008, 0.008 /) |
---|
693 | |
---|
694 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_wet_mtc = & !! NOx emissions factor soil emissions and exponential |
---|
695 | & (/ 0., 2.6, 0.06, 0.03, 0.03, 0.03, 0.03, & !! dependancy factor for wet soils |
---|
696 | & 0.03, 0.03, 0.36, 0.36, 0.36, 0.36 /) !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex |
---|
697 | |
---|
698 | REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_dry_mtc = & !! NOx emissions factor soil emissions and exponential |
---|
699 | & (/ 0., 8.60, 0.40, 0.22, 0.22, 0.22, 0.22, & !! dependancy factor for dry soils |
---|
700 | & 0.22, 0.22, 2.65, 2.65, 2.65, 2.65 /) !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex |
---|
701 | |
---|
702 | REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Larch_mtc = & !! Larcher 1991 SAI/LAI ratio (unitless) |
---|
703 | & (/ 0., 0.015, 0.015, 0.003, 0.005, 0.005, 0.003, & |
---|
704 | & 0.005, 0.003, 0.005, 0.005, 0.008, 0.008 /) |
---|
705 | |
---|
706 | |
---|
707 | |
---|
708 | END MODULE constantes_mtc |
---|