source: branches/publications/ORCHIDEE-GMv3.2/ORCHIDEE/src_parameters/constantes_mtc.f90 @ 6939

Last change on this file since 6939 was 6939, checked in by jinfeng.chang, 20 months ago

Update ORCHIDEE-GMv3.2 for publication

  • Property svn:keywords set to Date Revision
File size: 103.9 KB
Line 
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
43MODULE constantes_mtc
44
45  USE defprec
46  USE constantes_var
47
48  IMPLICIT NONE
49
50  !
51  ! METACLASSES CHARACTERISTICS
52  !
53
54  INTEGER(i_std), PARAMETER :: nvmc = 14                         !! 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  &    '          wheat                   ', &          ! 12
69  &    '          maize                   ', &          ! 13
70  &    '          rice                    '  /)         ! 14
71
72
73  !
74  ! VEGETATION STRUCTURE
75  !
76  INTEGER(i_std),PARAMETER, DIMENSION(nvmc) :: leaf_tab_mtc  =  &                 !! leaf type (1-4, unitless)
77  & (/  4,   1,   1,   2,   1,   1,   2,   &                                      !! 1=broad leaved tree, 2=needle leaved tree
78  &     1,   2,   3,   3,   3,   3,   3   /)                                           !! 3=grass 4=bare ground
79
80  CHARACTER(len=6), PARAMETER, DIMENSION(nvmc) :: pheno_model_mtc  =  &  !! which phenology model is used? (tabulated)
81  & (/  'none  ',   'none  ',   'moi   ',   'none  ',   'none  ',  &
82  &     'ncdgdd',   'none  ',   'ncdgdd',   'ngd   ',   'moigdd',  &
83  &     'moi_C4',   'moigdd',   'moigdd',   'moigdd'  /) 
84
85  LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_tropical_mtc  =  &                       !! Is PFT tropical ? (true/false)
86  & (/ .FALSE.,   .TRUE.,    .TRUE.,    .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,  &
87  &    .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.  /)   
88
89!gmjc
90  ! Is the vegetation a grassland that can be managed ?
91  LOGICAL,PARAMETER, DIMENSION(nvmc) :: is_grassland_manag_mtc   =  &
92  & (/ .FALSE., .FALSE.,  .FALSE., .FALSE., .FALSE.,  &
93  &    .FALSE.,  .FALSE.,  .FALSE., .FALSE., .FALSE., &
94  &    .FALSE., .FALSE., .FALSE.,  .FALSE.  /)
95  ! Is the vegetation a grassland that can be cut ?
96  LOGICAL,PARAMETER, DIMENSION(nvmc) :: is_grassland_cut_mtc   =  &
97  & (/ .FALSE., .FALSE.,  .FALSE., .FALSE., .FALSE.,  &
98  &    .FALSE.,  .FALSE.,  .FALSE., .FALSE., .FALSE., &
99  &    .FALSE., .FALSE., .FALSE. , .FALSE. /)
100  ! Is the vegetation a grassland that can be grazed ?
101  LOGICAL,PARAMETER, DIMENSION(nvmc) :: is_grassland_grazed_mtc   =  &
102  & (/ .FALSE., .FALSE.,  .FALSE., .FALSE., .FALSE.,  &
103  &    .FALSE.,  .FALSE.,  .FALSE., .FALSE., .FALSE., &
104  &    .FALSE., .FALSE., .FALSE. , .FALSE. /)
105  ! Management Intensity reading in MANAGEMENT_MAP
106  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: management_intensity_mtc  =  &
107  & (/ 0,   0,   0,   0,   0,   0,   0,  &
108  &    0,   0,   0,   0,   0,   0,   0  /)
109  ! Start year of management reading in MANAGEMENT_MAP & GRAZING_MAP
110  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: management_start_mtc  =  &
111  & (/ 0,   0,   0,   0,   0,   0,   0,  &
112  &    0,   0,   0,   0,   0,   0,   0  /)
113  ! Start year of N deposition reading in DEPOSITION_MAP
114  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: deposition_start_mtc  =  &
115  & (/ 0,   0,   0,   0,   0,   0,   0,  &
116  &    0,   0,   0,   0,   0,   0,   0  /)
117  ! Number of year that management should be read
118  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: nb_year_management_mtc  =  &
119  & (/ 1,   1,   1,   1,   1,   1,   1,  &
120  &    1,   1,   1,   1,   1,   1,   1  /)
121  ! maximum specific leaf area (m**2/gC)
122   REAL(r_std), PARAMETER, DIMENSION(nvmc)     ::    sla_max_mtc = &  !m2/gC
123(/1.5E-2,1.53E-2, 2.6E-2,9.26E-3, 2E-2, 2.6E-2, 9.26E-3, &
124 &                   2.6E-2, 1.9E-2, 2.6E-2,2.6E-2,2.6E-2,2.6E-2, 2.6E-2 /)
125  ! miniimum specific leaf area (m**2/gC)
126   REAL(r_std), PARAMETER, DIMENSION(nvmc)     ::    sla_min_mtc  = &
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, 2.6E-2 /)
129!end gmjc
130
131!Bioenergy seting
132  ! Is the metaclass bioenergy?
133  LOGICAL,PARAMETER, DIMENSION(nvmc) :: is_bioe1_mtc   =  &
134  & (/ .FALSE., .FALSE.,  .FALSE., .FALSE., .FALSE.,  &
135  &    .FALSE.,  .FALSE.,  .FALSE., .FALSE., .FALSE., &
136  &    .FALSE., .FALSE., .FALSE.,  .FALSE.  /)
137
138  CHARACTER(LEN=5), PARAMETER, DIMENSION(nvmc) :: type_of_lai_mtc  =  &  !! Type of behaviour of the LAI evolution algorithm
139  & (/ 'inter', 'inter', 'inter', 'inter', 'inter',  &                   !! for each vegetation type. (unitless)
140  &    'inter', 'inter', 'inter', 'inter', 'inter',  &                   !! Value of type_of_lai : mean or interp
141  &    'inter', 'inter', 'inter', 'inter' /)
142
143  LOGICAL, PARAMETER, DIMENSION(nvmc) :: natural_mtc =  &                         !! natural?  (true/false)
144  & (/ .TRUE.,   .TRUE.,   .TRUE.,   .TRUE.,   .TRUE.,    .TRUE.,   .TRUE.,  &
145  &    .TRUE.,   .TRUE.,   .TRUE.,   .TRUE.,   .FALSE.,   .FALSE., .FALSE.  /)
146
147! dgvmjc
148  LOGICAL, PARAMETER, DIMENSION(nvmc) :: pasture_mtc =  & !! pasture? (true/false)
149  & (/ .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,    .FALSE.,   .FALSE.,  &
150  &    .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.  /)
151! end dgvmjc
152
153  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: veget_ori_fixed_mtc  =  &  !! Value for veget_ori for tests in
154  & (/ 0.2,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                !! 0-dim simulations (0-1, unitless)
155  &    0.0,   0.0,   0.8,   0.0,   0.0,   0.0,   0.0  /)
156
157  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimax_mtc  =  &          !! laimax for maximum
158  & (/ 0.0,   8.0,   8.0,   4.0,   4.5,   4.5,   4.0,  &                !! See also type of lai interpolation
159  &    4.5,   4.0,   2.0,   2.0,   2.0,   2.0,  2.0  /)                       !! @tex $(m^2.m^{-2})$ @endtex
160
161  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: llaimin_mtc  = &           !! laimin for minimum lai
162  & (/ 0.0,   8.0,   0.0,   4.0,   4.5,   0.0,   4.0,  &                !! See also type of lai interpolation (m^2.m^{-2})
163  &    0.0,   0.0,   0.0,   0.0,   0.0,   0.0, 0.0/)                       !! @tex $(m^2.m^{-2})$ @endtex
164
165  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: height_presc_mtc  =  &     !! prescribed height of vegetation (m)
166  & (/  0.0,   30.0,   30.0,   20.0,   20.0,   20.0,   15.0,  &         !! Value for height_presc : one for each vegetation type
167  &    15.0,   15.0,    0.5,    0.6,    1.0,    1.0,   1.0/)
168
169  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: z0_over_height_mtc = &         !! Factor to calculate roughness height from
170  & (/  0.0, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625,  &         !! vegetation height (unitless)   
171  &  0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625, 0.0625  /)
172
173  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ratio_z0m_z0h_mtc = &      !! Ratio between z0m and z0h values (roughness height for momentum and for heat)
174  & (/  1.0,    1.0,    1.0,    1.0,    1.0,    1.0,    1.0,  &         
175  &     1.0,    1.0,    1.0,    1.0,    1.0,    1.0,    1.0  /)
176
177
178  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rveg_mtc  =  &             !! Potentiometer to set vegetation resistance (unitless)
179  & (/ 1.0,   1.0,   1.0,   1.0,   1.0,   1.0,   1.0,  &                !! Nathalie on March 28th, 2006 - from Fred Hourdin,
180  &    1.0,   1.0,   1.0,   1.0,   1.0,   1.0,   1.0   /)
181
182  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: sla_mtc  =  &                       !! specif leaf area @tex $(m^2.gC^{-1})$ @endtex
183  & (/ 1.5E-2,   1.53E-2,   2.6E-2,   9.26E-3,     2E-2,   2.6E-2,   9.26E-3,  &
184  &    2.6E-2,    1.9E-2,   2.6E-2,    2.6E-2,   2.6E-2,   2.6E-2,   2.6E-2  /) 
185
186  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: availability_fact_mtc  =  &     !! calculate mortality in lpj_gap
187  & (/ undef,   0.14,  0.14,   0.10,   0.10,   0.10,   0.05,  &
188  &     0.05,   0.05,  undef,  undef,  undef,  undef, undef  /)
189
190
191  !
192  ! EVAPOTRANSPIRATION (sechiba)
193  !
194  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: rstruct_const_mtc  =  &  !! Structural resistance.
195  & (/ 0.0,   25.0,   25.0,   25.0,   25.0,   25.0,   25.0,  &        !! @tex $(s.m^{-1})$ @endtex
196  &   25.0,   25.0,    2.5,    2.0,    2.0,    2.0,   2.0   /)
197
198  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kzero_mtc  =  &                  !! A vegetation dependent constant used in the
199  & (/    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.
200  &    25.E-5,   25.E-5,   30.E-5,   30.E-5,   30.E-5,   30.E-5,  30.E-5  /)           !! @tex $(kg.m^2.s^{-1})$ @endtex
201
202
203  !
204  ! WATER (sechiba)
205  !
206  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: wmax_veg_mtc  =  &        !! Volumetric available soil water capacity in each PFT
207  & (/ 150.0,   150.0,   150.0,   150.0,   150.0,   150.0,   150.0,  & !! @tex $(kg.m^{-3} of soil)$ @endtex
208  &    150.0,   150.0,   150.0,   150.0,   150.0,   150.0,   150.0  /)         
209                                                                     
210
211  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref4m  =  &       !! Root profile description for the different
212  & (/ 5.0,   0.4,   0.4,   1.0,   0.8,   0.8,   1.0,  &               !! vegetations types. @tex $(m^{-1})$ @endtex
213  &    1.0,   0.8,   4.0,   1.0,   4.0,   1.0,   4.0  /)                      !! These are the factor in the exponential which gets       
214                                                                       !! the root density as a function of depth
215                                                                       !! Values for zmaxh = 4.0 
216 
217  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: humcste_ref2m  =  &       !! Root profile description for the different
218  & (/ 5.0,   0.8,   0.8,   1.0,   0.8,   0.8,   1.0,  &               !! vegetations types.  @tex $(m^{-1})$ @endtex
219  &    1.0,   0.8,   4.0,   4.0,   4.0,   4.0,   4.0  /)                      !! These are the factor in the exponential which gets       
220                                                                       !! the root density as a function of depth
221                                                                       !! Values for zmaxh = 2.0
222
223  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: throughfall_by_mtc  =  &  !! Fraction of rain intercepted by the canopy
224  & (/ 30.0,   30.0,   30.0,   30.0,   30.0,   30.0,   30.0,  &        !! (0-100, unitless)
225  &    30.0,   30.0,   30.0,   30.0,   30.0,   30.0,   30.0  /)
226
227
228  !
229  ! ALBEDO (sechiba)
230  !
231  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_aged_vis_mtc  =  &    !! Minimum snow albedo value for each vegetation type
232  & (/ 0.50,    0.0,    0.0,   0.15,   0.14,   0.14,   0.15,  &        !! after aging (dirty old snow) (unitless), visible albedo
233  &    0.14,   0.22,   0.35,   0.35,   0.35,   0.35,   0.35   /)       !!  Source : Values are from the Thesis of S. Chalita (1992), optimized on 04/07/2016
234
235  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_aged_nir_mtc  =  &  !! Minimum snow albedo value for each vegetation type
236  & (/ 0.35,    0.0,    0.0,   0.14,   0.14,   0.14,   0.14,  &        !! after aging (dirty old snow) (unitless), near infrared albedo
237  &    0.14,   0.14,   0.18,   0.18,   0.18,   0.18,   0.18  /)                !! Source : Values are from the Thesis of S. Chalita (1992)
238
239  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_dec_vis_mtc  =  &   !! Decay rate of snow albedo value for each vegetation type
240  & (/ 0.45,    0.0,    0.0,   0.10,   0.06,   0.11,   0.10,  &        !! as it will be used in condveg_snow (unitless), visible albedo
241  &    0.11,   0.18,   0.60,   0.60,   0.60,   0.60,   0.60  /)                !! Source : Values are from the Thesis of S. Chalita (1992), optimized on 04/07/2016
242
243  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: snowa_dec_nir_mtc  =  &   !! Decay rate of snow albedo value for each vegetation type
244  & (/ 0.45,    0.0,    0.0,   0.06,   0.06,   0.11,   0.06,  &        !! as it will be used in condveg_snow (unitless), near infrared albedo
245  &    0.11,   0.11,   0.52,   0.52,   0.52,   0.52,   0.52  /)                !! Source : Values are from the Thesis of S. Chalita (1992)
246
247  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_vis_mtc  =  &    !! leaf albedo of vegetation type, visible albedo, optimized on 04/07/2016
248  & (/ 0.00,   0.0397, 0.0474, 0.0386, 0.0484, 0.0411, 0.041,  &       !! (unitless)
249  &    0.0541, 0.0435, 0.0524, 0.0508, 0.0509, 0.0606, 0.0606  /)
250
251  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alb_leaf_nir_mtc  =  &    !! leaf albedo of vegetation type, near infrared albedo, optimized on 04/07/2016
252  & (/ 0.00,   0.227,  0.214,  0.193,  0.208,  0.244,  0.177,  &       !! (unitless)
253  &    0.218,  0.213,  0.252,  0.265,  0.272,  0.244,  0.244  /)
254
255  !
256  ! SOIL - VEGETATION
257  !
258  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pref_soil_veg_mtc  =  &       !! The soil tile number for each vegetation
259  & (/ 1,   2,   2,   2,   2,   2,   2,  &                                   
260  &    2,   2,   3,   3,   3,   3,   3  /)     
261!  LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_crop_soil = &                !! whether it is a crop soil tile
262
263  !
264  ! VEGETATION - AGE CLASSES
265  !
266  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: agec_group_mtc  =  &       !! The age class group that each PFT belongs to.
267       (/ 1,   2,   3,   4,   5,   6,   7,  &                                   
268       8,   9,   10,   11,   12,   13,  14  /)                                         
269
270
271
272  !
273  ! PHOTOSYNTHESIS
274  !
275  !-
276  ! 1 .CO2
277  !-
278  LOGICAL, PARAMETER, DIMENSION(nvmc) :: is_c4_mtc  =  &                            !! flag for C4 vegetation types (true/false)
279  & (/ .FALSE.,  .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,  &
280  &    .FALSE.,  .FALSE.,   .FALSE.,   .TRUE.,    .FALSE.,   .TRUE.,   .FALSE.  /)
281
282  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: vcmax_fix_mtc  =  &     !! values used for vcmax when STOMATE is not
283  & (/  0.0,   40.0,   50.0,   30.0,   35.0,   40.0,   30.0,  &      !! activated @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex
284  &    40.0,   35.0,   60.0,   60.0,   70.0,   70.0,   70.0  /)
285
286! For C4 plant we define a very small downregulation effect as C4 plant are
287! currently saturate with respect to CO2 impact on vcmax
288  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: downregulation_co2_coeff_mtc  =  &  !! coefficient for CO2 downregulation
289  & (/  0.0,   0.38,   0.38,   0.28,   0.28,   0.28,   0.22,  &
290  &     0.22,  0.22,   0.26,   0.03,   0.26,   0.03,   0.26 /)
291
292  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmC_mtc  = &            !! Energy of activation for KmC (J mol-1)
293  & (/undef,  79430.,  79430.,  79430.,  79430.,  79430.,  79430.,  &  !! See Medlyn et al. (2002)
294  &  79430.,  79430.,  79430.,  79430.,  79430.,  79430.,  79430.  /)  !! from Bernacchi al. (2001)
295
296  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_KmO_mtc  = &            !! Energy of activation for KmO (J mol-1)
297  & (/undef,  36380.,  36380.,  36380.,  36380.,  36380.,  36380.,  &  !! See Medlyn et al. (2002)
298  &  36380.,  36380.,  36380.,  36380.,  36380.,  36380.,  36380.  /)           !! from Bernacchi al. (2001)
299
300  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Sco_mtc  = & !! Energy of activation for Sco (J mol-1)
301    & (/undef, -24460., -24460., -24460., -24460., -24460., -24460.,  &  !! See Table 2 of Yin et al. (2009)
302    & -24460., -24460., -24460., -24460., -24460., -24460., -24460.  /) !! Value for C4 plants is not mentioned - We use C3 for all plants
303
304  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_gamma_star_mtc  = &     !! Energy of activation for gamma_star (J mol-1)
305  & (/undef,  37830.,  37830.,  37830.,  37830.,  37830.,  37830.,  &  !! See Medlyn et al. (2002) from Bernacchi al. (2001)
306  &  37830.,  37830.,  37830.,  37830.,  37830.,  37830.,  37830.  /)           !! for C3 plants - We use the same values for C4 plants
307
308  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Vcmax_mtc  = &          !! Energy of activation for Vcmax (J mol-1)
309  & (/undef,  71513.,  71513.,  71513.,  71513.,  71513.,  71513.,  &  !! See Table 2 of Yin et al. (2009) for C4 plants
310  &  71513.,  71513.,  71513.,  67300.,  71513.,  67300.,  71513.  /)           !! and Kattge & Knorr (2007) for C3 plants (table 3)
311
312  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Jmax_mtc  = &            !! Energy of activation for Jmax (J mol-1)
313  & (/undef,  49884.,  49884.,  49884.,  49884.,  49884.,  49884.,  &   !! See Table 2 of Yin et al. (2009) for C4 plants
314  &  49884.,  49884.,  49884.,  77900.,  49884.,  77900.,  49884.  /)            !! and Kattge & Knorr (2007) for C3 plants (table 3)
315
316  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)
317  & (/undef,  668.39,  668.39,  668.39,  668.39,  668.39,  668.39,  &   !! See Table 3 of Kattge & Knorr (2007)
318  &  668.39,  668.39,  668.39,  641.64,  668.39,  641.64,  668.39  /)            !! For C4 plants, we assume that there is no
319                                                                        !! 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)
320
321  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)
322  & (/undef,   -1.07,   -1.07,   -1.07,   -1.07,   -1.07,   -1.07,  &   !! See Table 3 of Kattge & Knorr (2007)
323  &   -1.07,   -1.07,   -1.07,      0.,   -1.07,      0.,   -1.07  /)            !! We assume No acclimation term for C4 plants
324
325  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_min_mtc  =  &  !! minimum photosynthesis temperature (deg C)
326  & (/  undef,   -4.0,    -4.0,   -4.0,   -4.0,   -4.0,   -4.0,  & 
327  &      -4.0,   -4.0,    -4.0,   -4.0,   -4.0,   -4.0,   -4.0  /)
328
329  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tphoto_max_mtc  =  &  !! maximum photosynthesis temperature (deg C)
330  & (/  undef,   55.0,    55.0,   55.0,   55.0,   55.0,   55.0,  & 
331  &      55.0,   55.0,    55.0,   55.0,   55.0,   55.0,   55.0  /)
332
333  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)
334  & (/undef,  659.70,  659.70,  659.70,  659.70,  659.70,  659.70,  &   !! See Table 3 of Kattge & Knorr (2007)
335  &  659.70,  659.70,  659.70,    630.,  659.70,    630.,  659.70  /)            !! and Table 2 of Yin et al. (2009) for C4 plants
336
337  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)
338  & (/undef,   -0.75,   -0.75,   -0.75,   -0.75,   -0.75,   -0.75,  &   !! See Table 3 of Kattge & Knorr (2007)
339  &   -0.75,   -0.75,   -0.75,      0.,   -0.75,      0.,   -0.75  /)            !! We assume no acclimation term for C4 plants
340
341  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Vcmax_mtc  = &           !! Energy of deactivation for Vcmax (J mol-1)
342  & (/undef, 200000., 200000., 200000., 200000., 200000., 200000.,  &   !! Medlyn et al. (2002) also uses 200000. for C3 plants (same value than D_Jmax)
343  & 200000., 200000., 200000., 192000., 200000., 192000., 200000.  /)            !! 'Consequently', we use the value of D_Jmax for C4 plants
344
345  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_Jmax_mtc  = &            !! Energy of deactivation for Jmax (J mol-1)
346  & (/undef, 200000., 200000., 200000., 200000., 200000., 200000.,  &   !! See Table 2 of Yin et al. (2009)
347  & 200000., 200000., 200000., 192000., 200000., 192000., 200000.  /)            !! Medlyn et al. (2002) also uses 200000. for C3 plants
348
349  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_gm_mtc  = & !! Energy of activation for gm (J mol-1) 
350  & (/undef,  49600.,  49600.,  49600.,  49600.,  49600., 49600.,  &   !! See Table 2 of Yin et al. (2009) 
351  &  49600.,  49600.,  49600.,   undef,  49600.,   undef, 49600.  /)             
352                 
353  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: S_gm_mtc  = & !! Entropy term for gm (J K-1 mol-1) 
354  & (/undef,   1400.,   1400.,   1400.,   1400.,   1400., 1400.,  &   !! See Table 2 of Yin et al. (2009) 
355  &   1400.,   1400.,   1400.,   undef,   1400.,   undef, 1400.  /) 
356                 
357  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: D_gm_mtc  = & !! Energy of deactivation for gm (J mol-1) 
358  & (/undef, 437400., 437400., 437400., 437400., 437400., 437400.,  &   !! See Table 2 of Yin et al. (2009) 
359  & 437400., 437400., 437400.,   undef, 437400.,   undef, 437400.  /) 
360
361  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: E_Rd_mtc  = &              !! Energy of activation for Rd (J mol-1)
362  & (/undef,  46390.,  46390.,  46390.,  46390.,  46390.,  46390.,  &   !! See Table 2 of Yin et al. (2009)
363  &  46390.,  46390.,  46390.,  46390.,  46390.,  46390.,  46390.  /)           
364
365  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Vcmax25_mtc  =  &          !! Maximum rate of Rubisco activity-limited carboxylation at 25°C
366  & (/ undef,   50.0,    65.0,    45.0,   45.0,   55.0,   45.0,  &      !! @tex $(\mu mol.m^{-2}.s^{-1})$ @endtex
367  &     45.0,   35.0,    70.0,    70.0,   70.0,   70.0,   70.0  /)
368
369  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)
370  & (/undef,    2.59,    2.59,    2.59,    2.59,    2.59,    2.59,  &   !! See Table 3 of Kattge & Knorr (2007)
371  &    2.59,    2.59,    2.59,   1.715,    2.59,   1.715,   2.59  /)            !! For C4 plants, we assume that there is no
372                                                                        !! 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)
373
374  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)
375  & (/undef,  -0.035,  -0.035,  -0.035,  -0.035,  -0.035,  -0.035,  &   !! See Table 3 of Kattge & Knorr (2007)
376  &  -0.035,  -0.035,  -0.035,      0.,  -0.035,      0.,  -0.035  /)            !! We assume No acclimation term for C4 plants
377
378  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmC25_mtc  = &             !! Michaelis–Menten constant of Rubisco for CO2 at 25°C (ubar)
379  & (/undef,   404.9,   404.9,   404.9,   404.9,  404.9,   404.9,  &    !! See Table 2 of Yin et al. (2009) for C4
380  &   404.9,   404.9,   404.9,    650.,   404.9,   650.,   404.9  /)             !! and Medlyn et al (2002) for C3
381
382  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: KmO25_mtc  = &             !! Michaelis–Menten constant of Rubisco for O2 at 25°C (ubar)
383  & (/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
384  & 278400., 278400., 278400., 450000., 278400., 450000., 278400.  /)           
385
386  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Sco25_mtc  = & !! Relative CO2 /O2 specificity factor for Rubisco at 25???¡ã(bar bar-1)
387    & (/undef,   2800.,   2800.,   2800.,   2800.,   2800., 2800.,  &   !! See Table 2 of Yin et al. (2009)
388    &   2800.,   2800.,   2800.,   2590.,   2800.,   2590., 2800.  /)           
389   
390    REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gm25_mtc  = & !! Mesophyll diffusion conductance at 25???¡ã(mol m-2 s-1 bar-1) 
391    & (/undef,     0.4,     0.4,     0.4,     0.4,    0.4, 0.4,  &   !! See legend of Figure 6 of Yin et al. (2009) 
392    &     0.4,     0.4,     0.4,   undef,     0.4,  undef, 0.4  /) !! and review by Flexas et al. (2008) - gm is not used for C4 plants 
393
394  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gamma_star25_mtc  = &      !! Ci-based CO2 compensation point in the absence of Rd at 25°C (ubar)
395  & (/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)
396  &   42.75,   42.75,   42.75,   42.75,   42.75,   42.75,   42.75  /)   
397
398  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: a1_mtc  = &                !! Empirical factor involved in the calculation of fvpd (-)
399  & (/undef,    0.85,    0.85,    0.85,    0.85,    0.85,  0.85,  &     !! See Table 2 of Yin et al. (2009)
400  &    0.85,    0.85,    0.85,    0.85,    0.72,    0.85,  0.72  /)           
401
402  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: b1_mtc  = &                !! Empirical factor involved in the calculation of fvpd (-)
403  & (/undef,    0.14,    0.14,    0.14,    0.14,    0.14,  0.14,  &     !! See Table 2 of Yin et al. (2009)
404  &    0.14,    0.14,    0.14,    0.20,    0.14,    0.20,  0.14  /)           
405
406  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: g0_mtc  = &                !! Residual stomatal conductance when irradiance approaches zero (mol CO2 m−2 s−1 bar−1)
407  & (/undef, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625, 0.00625,  &   !! Value from ORCHIDEE - No other reference.
408  & 0.00625, 0.00625, 0.00625, 0.01875, 0.00625, 0.01875, 0.00625  /)            !! modofy to account for the conversion for conductance to H2O to CO2
409
410  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: h_protons_mtc  = &         !! Number of protons required to produce one ATP (mol mol-1)
411  & (/undef,      4.,      4.,      4.,      4.,      4.,    4.,  &     !! See Table 2 of Yin et al. (2009) - h parameter
412  &      4.,      4.,      4.,      4.,      4.,      4.,    4.  /)           
413
414  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpsir_mtc = &              !! Fraction of PSII e− transport rate
415  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! partitioned to the C4 cycle (-)
416  &   undef,   undef,   undef,     0.4,   undef,    0.4,  undef  /)             !! See Table 2 of Yin et al. (2009) - x parameter       
417 
418  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fQ_mtc = &                 !! Fraction of electrons at reduced plastoquinone
419  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! that follow the Q-cycle (-) - Values for C3 platns are not used
420  &   undef,   undef,   undef,      1.,   undef,     1.,  undef  /)             !! See Table 2 of Yin et al. (2009)         
421
422  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: fpseudo_mtc = &            !! Fraction of electrons at PSI that follow
423  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! pseudocyclic transport (-) - Values for C3 platns are not used
424  &   undef,   undef,   undef,     0.1,   undef,    0.1,  undef  /)             !! See Table 2 of Yin et al. (2009)   
425
426  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: kp_mtc = &                 !! Initial carboxylation efficiency of the PEP carboxylase (mol m−2 s−1 bar−1)
427  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! See Table 2 of Yin et al. (2009)
428  &   undef,   undef,   undef,     0.7,   undef,    0.7,  undef  /)                 
429
430  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alpha_mtc = &              !! Fraction of PSII activity in the bundle sheath (-)
431  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! See legend of Figure 6 of Yin et al. (2009)
432  &   undef,   undef,   undef,     0.1,   undef,    0.1,  undef  /)                 
433
434  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gbs_mtc = &                !! Bundle-sheath conductance (mol m−2 s−1 bar−1)
435  & (/undef,   undef,   undef,   undef,   undef,  undef,  undef,  &     !! See legend of Figure 6 of Yin et al. (2009)
436  &   undef,   undef,   undef,   0.003,   undef,  0.003,  undef  /)   
437
438  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: theta_mtc = &              !! Convexity factor for response of J to irradiance (-)
439  & (/undef,     0.7,     0.7,     0.7,     0.7,    0.7,    0.7,  &     !! See Table 2 of Yin et al. (2009)
440  &     0.7,     0.7,     0.7,     0.7,     0.7,    0.7,    0.7  /)
441
442  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)
443  & (/undef,     0.3,     0.3,     0.3,     0.3,    0.3,    0.3,  &     !! See comment from Yin et al. (2009) after eq. 4
444  &     0.3,     0.3,     0.3,     0.3,     0.3,    0.3,    0.3  /)             !! alpha value from Medlyn et al. (2002)   
445                                                                        !! 0.093 mol CO2 fixed per mol absorbed photons
446                                                                        !! times 4 mol e- per mol CO2 produced
447
448  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_vcmax_mtc = &       !! Water stress on vcmax
449  & (/    1.,     1.,     1.,       1.,      1.,     1., 1., & 
450  &      1.,     1.,     1.,       1.,      1.,     1. , 1. /) 
451 
452  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_gs_mtc = & !! Water stress on gs
453  & (/    0.,     0.,     0.,       0.,      0.,     0., 0., & 
454  &      0.,     0.,     0.,       0.,      0.,     0. , 0. /) 
455 
456  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: stress_gm_mtc = & !! Water stress on gm
457  & (/    0.,     0.,     0.,       0.,      0.,     0., 0., & 
458  &      0.,     0.,     0.,       0.,      0.,     0. , 0. /) 
459   
460  !-
461  ! 2 .Stomate
462  !-
463  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ext_coeff_mtc  =  &     !! extinction coefficient of the Monsi&Saeki
464  & (/ 0.5,   0.5,   0.5,   0.5,   0.5,   0.5,   0.5,  &             !! relationship (1953) (unitless)
465  &    0.5,   0.5,   0.5,   0.5,   0.5,   0.5,   0.5  /)
466  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ext_coeff_vegetfrac_mtc  =  &     !! extinction coefficient used for defining the fraction
467  & (/ 1.0,   1.0,   1.0,   1.0,   1.0,   1.0,   1.0,  &                       !!  of bare soil (unitless)
468  &    1.0,   1.0,   1.0,   1.0,   1.0,   1.0,   1.0  /)
469
470  !
471  ! ALLOCATION (stomate)
472  !
473  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: R0_mtc = &              !! Default root allocation (0-1, unitless)
474  & (/ undef,   0.30,   0.30,   0.30,   0.30,  0.30,    0.30, &
475  &     0.30,   0.30,   0.30,   0.30,   0.30,  0.30,    0.30 /)                   
476
477  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: S0_mtc = &              !! Default sapwood allocation (0-1, unitless)
478  & (/ undef,   0.25,   0.25,   0.30,   0.30,  0.30,    0.30, &
479  &     0.30,   0.30,   0.30,   0.30,   0.30,  0.30,    0.30 /)                   
480
481  !
482  ! RESPIRATION (stomate)
483  !
484  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: frac_growthresp_mtc  =  &  !! fraction of GPP which is lost as growth respiration
485  & (/  undef,   0.28,   0.28,   0.28,   0.28,   0.28,   0.28,  &
486  &      0.28,   0.28,   0.28,   0.28,   0.28,   0.28,   0.28  /)
487
488  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_c_mtc  =  &  !! slope of maintenance respiration coefficient (1/K),
489  & (/  undef,   0.20,   0.20,   0.16,   0.16,   0.16,   0.16,  &          !! constant c of aT^2+bT+c, tabulated
490  &      0.16,   0.16,   0.16,   0.12,   0.16,   0.12,   0.16  /)
491
492  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_b_mtc  =  &  !! slope of maintenance respiration coefficient (1/K),
493  & (/  undef,   0.0,        0.0,   0.0,        0.0,   0.0,   0.0,  &      !! constant b of aT^2+bT+c, tabulated
494  &       0.0,   0.0,   -0.00133,   0.0,   -0.00133,   0.0,   -0.00133  /)
495
496  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: maint_resp_slope_a_mtc  =  &  !! slope of maintenance respiration coefficient (1/K),
497  & (/  undef,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                !! constant a of aT^2+bT+c, tabulated
498  &       0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0  /)
499
500  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_leaf_mtc  =   &                  !! maintenance respiration coefficient
501  & (/   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,
502  &    2.62E-3,   2.05E-3,   2.62E-3,   2.62E-3,   2.62E-3,   2.62E-3,   2.62E-3  /)  !! @tex $(gC.gC^{-1}.day^{-1})$ @endtex
503
504  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapabove_mtc =  &                !! maintenance respiration coefficient
505  & (/   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,
506  &    1.19E-4,   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
507
508  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_sapbelow_mtc  =  &               !! maintenance respiration coefficient
509  & (/   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,
510  &    1.19E-4,   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
511
512  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartabove_mtc  =  &             !! maintenance respiration coefficient
513  & (/  undef,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                           !! at 0 deg C, for heartwood above,
514  &       0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0  /)                                  !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex
515
516  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_heartbelow_mtc  =  &             !! maintenance respiration coefficient
517  & (/  undef,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                           !! at 0 deg C, for heartwood below,
518  &       0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0  /)                                  !! tabulated, @tex $(gC.gC^{-1}.day^{-1})$ @endtex
519
520  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_root_mtc  =  &                   !! maintenance respiration coefficient
521  & (/   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,
522  &    1.67E-3,   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
523
524  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_fruit_mtc  =  &                  !! maintenance respiration coefficient
525  & (/   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,
526  &    1.19E-4,   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
527
528  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: cm_zero_carbres_mtc  =  &                !! maintenance respiration coefficient
529  & (/   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,
530  &    1.19E-4,   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
531
532
533  !
534  ! FIRE (stomate)
535  !
536  REAL(r_std),PARAMETER, DIMENSION(nvmc) :: flam_mtc  =  &         !! flamability: critical fraction of water
537  & (/  undef,   0.15,   0.25,   0.25,   0.25,   0.25,   0.25,  &  !! holding capacity (0-1, unitless)
538  &      0.25,   0.25,   0.25,   0.25,   0.35,   0.35,   0.35  /)
539
540  REAL(r_std),PARAMETER, DIMENSION(nvmc) :: resist_mtc  =  &       !! fire resistance (0-1, unitless)
541  & (/ undef,   0.95,   0.90,   0.90,   0.90,   0.90,   0.90,  &
542  &    0.90,    0.90,    0.0,    0.0,    0.0,    0.0,   0.0  /) 
543
544!spitfire
545  ! dens_fuel : fuel bulk density
546  REAL(r_std), SAVE, DIMENSION(nvmc)                     ::  dens_fuel_mtc = (/0.0,   25.,      25.,     25.,     10.,      22.,    25.,  &   !!kgm^{-3}
547                                                                       22.,      22.,     2.,      2.,      2.,     2.,     2.   /)
548  ! f_SH : scorch height parameter for crown fire
549  REAL(r_std), SAVE, DIMENSION(nvmc)                     :: f_sh_mtc = (/ 0.0, 0.1487,    0.061,     0.1,   0.371,    0.094,   0.11,  &
550                                                                  0.094,    0.094,      0.,      0.,       0.,     0.,    0.   /)
551  ! crown_length : crown length
552  REAL(r_std), SAVE, DIMENSION(nvmc)                     :: crown_length_mtc = (/ 0.0, 0.3334,     0.10,  0.3334,  0.3334,   0.3334, 0.6667,  &
553                                                                           0.3334,   0.3334,      0.,      0.,       0.,      0.,  0.  /)
554  ! param1 :  Bark thickness parameter
555  REAL(r_std), SAVE, DIMENSION(nvmc)                     :: BTpar1_mtc = (/ 0.0, 0.0301,   0.1085,  0.0670,  0.0451,   0.0347, 0.0292,  &
556                                                                     0.0347,   0.0347,      0.,      0.,       0.,     0.,  0.   /)
557  ! param2 :  Bark thickness parameter
558  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: BTpar2_mtc = (/ 0.0, 0.0281,    0.212,  0.5590,  0.1412,   0.1086, 0.2632,  &
559                                                                      0.1086,   0.1086,      0.,      0.,       0.,     0.,  0.   /)
560  ! r_ck and p_ck parameter for postfire mortality as a result of crown damage
561  ! r_ck
562  REAL(r_std), SAVE, DIMENSION(nvmc)                     :: r_ck_mtc = (/  0.0, 0.95,     0.05,    0.95,    0.95,     0.95,   0.95,  &
563                                                                   0.95,     0.95,      0.,      0.,       0.,     0.,   0.   /)
564  ! p_ck
565  REAL(r_std), SAVE, DIMENSION(nvmc)                     :: p_ck_mtc = (/    0.0, 3.,       3.,     3.75,      3.,       3.,     3.,  &
566                                                                   3.,       3.,       0.,      0.,       0.,     0.,   0.   /)
567  ! emissions factors, unit: g/(kg dry mass burned)
568  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_CO2_mtc = (/ 0.0, 1580.,    1664.,    1568.,   1568.,    1568.,  1568.,  &
569                                                                      1568.,    1568.,    1568.,   1664.,    1568.,  1664.,  1568.   /)
570  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_CO_mtc = (/  0.0, 103.,      63.,     106.,    106.,     106.,   106.,  &
571                                                                     106.,     106.,     106.,     63.,     106.,    63.,  106.   /)
572  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_CH4_mtc = (/  0.0,  6.8,      2.2,      4.8,     4.8,      4.8,    4.8,  &
573                                                                      4.8,      4.8,      4.8,     2.2,      4.8,    2.2, 4.8   /)
574  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_VOC_mtc = (/   0.0, 8.1,      3.4,      5.7,     5.7,      5.7,    5.7,  &
575                                                                       5.7,      5.7,      5.7,     3.4,      5.7,    3.4,  5.7   /)
576  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_TPM_mtc = (/   0.0, 8.5,      8.5,     17.6,    17.6,     17.6,   17.6,  &
577                                                                      17.6,     17.6,     17.6,     8.5,     17.6,    8.5,   17.6   /)
578  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: ef_NOx_mtc = (/ 0.0, 1.999,     2.54,     3.24,    3.24,     3.24,   3.24,  &
579                                                                       3.24,     3.24,     3.24,    2.54,     3.24,   2.54, 3.24   /)
580  ! me : flammability threshold
581  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: me_mtc = (/    1.0, 0.2,     0.3,      0.3,     0.3,      0.3,   0.35,  &
582                                                                   0.35,    0.35,      0.2,     0.2,      0.2,    0.2,  0.2   /)
583
584!  ! Maximum combustion fraction for 100hr and 1000hr fuel for different PFTs.
585!  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: fire_max_cf_mtc = (/    0.0, 0.48,     0.48,      0.45,     0.45,      0.45,   0.41,  &
586!                                                                   0.41,    0.41,      0.85,     0.85,      0.35,    0.35,  0.35   /)
587  !! The 100hr and 1000hr maximum CCs are from Table 2 of Leeuwen et al. (2014),
588  !! with 100hr corresponding to woody debris of 2.55-7.6cm and 1000hr corresponding
589  !! to woody debris of 7.2-20.5cm. The CCs for tropical tress are taken from
590  !! Table 2 (b); The CCs for temperate trees are taken from Table 2 (c). As the
591  !! The equivalent CCs for woody debris were not reported for boreal forests in
592  !! Table 2 (d), we tentatively use those for temperate forests. The CCs for grassland
593  !! are taken from Table 3 by averaging the CCs for savanna and Grassland savanna
594  !! as we don't make distinction between savanna and Grassland savanna in our model
595  !! yet. The CCs for crops are set tentatively as zero to indicate the fact that
596  !! cropland fires are not simulated in the model.
597  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: fire_max_cf_100hr_mtc = (/    0.0, 0.65,     0.65,      0.73,     0.73,      0.73,   0.73,  &
598                                                                     0.73,    0.73,      0.76,     0.76,      0.,    0.,  0.  /)
599  REAL(r_std) , SAVE, DIMENSION(nvmc)                     :: fire_max_cf_1000hr_mtc = (/    0.0, 0.41,     0.41,      0.38,     0.38,      0.38,   0.38,  &
600                                                                   0.38,    0.38,      0.76,     0.76,      0.,    0., 0.  /)
601
602
603!endspit
604
605  !
606  ! FLUX - LUC
607  !
608  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_1_mtc  =  &   !! Coeff of biomass export for the year
609  & (/  undef,   0.897,   0.897,   0.597,   0.597,   0.597,   0.597,  &   !! (0-1, unitless)
610  &     0.597,   0.597,   0.597,   0.597,   0.597,   0.597,   0.597  /)
611
612  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_10_mtc  =  &  !! Coeff of biomass export for the decade
613  & (/  undef,   0.103,   0.103,   0.299,   0.299,   0.299,   0.299,  &   !! (0-1, unitless)
614  &     0.299,   0.299,   0.299,   0.403,   0.299,   0.403,   0.299  /) 
615
616  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_lcchange_100_mtc  =  & !! Coeff of biomass export for the century
617  & (/  undef,     0.0,     0.0,   0.104,   0.104,   0.104,   0.104,  &   !! (0-1, unitless)
618  &     0.104,   0.104,   0.104,     0.0,   0.104,     0.0,   0.104  /)
619
620
621  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_indwood_1_mtc  =  &   !! Coeff of biomass export for the year
622  & (/  undef,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &   !! (0-1, unitless)
623  &     0.0,     0.0,   0.0,   0.0,   0.0,   0.0,   0.0  /)
624
625  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_indwood_10_mtc  =  &  !! Coeff of biomass export for the decade
626  & (/  undef,   0.105,   0.105,   0.19,   0.19,   0.19,   0.19,  &   !! (0-1, unitless)
627  &     0.19,    0.19,    0.0,     0.0,    0.0,    0.0,    0.0  /) 
628
629  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: coeff_indwood_100_mtc  =  & !! Coeff of biomass export for the century
630  & (/  undef,   0.526,   0.526,   0.473,   0.473,   0.473,   0.473,  &   !! (0-1, unitless)
631  &     0.473,   0.473,   0.0,     0.0,     0.0,     0.0,     0.0  /)
632
633  !
634  ! PHENOLOGY
635  !
636  ! The latest modifications regarding leafagecrit, senescence_temp_c, leaffall, hum_min_time and nosenescence_hum are inspired by
637  ! MacBean et al. (2015), following the optimization of phenology parameters using MODIS NDVI (FM/PP).
638  !-
639  ! 1. Stomate
640  !-
641  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: lai_max_to_happy_mtc  =  &  !! threshold of LAI below which plant uses carbohydrate reserves
642  & (/  undef,   0.5,   0.5,   0.5,   0.5,   0.5,   0.5,  &
643  &       0.5,   0.5,   0.5,   0.5,   0.5,   0.5,   0.5  /)
644
645  REAL(r_std), PARAMETER, DIMENSION (nvmc) :: lai_max_mtc  =  &          !! maximum LAI, PFT-specific
646  & (/ undef,   7.0,   7.0,   5.0,   5.0,   5.0,   4.5,  &               !! @tex $(m^2.m^{-2})$ @endtex
647  &      4.5,   3.0,   2.5,   2.5,   5.0,   5.0,   5.0  /)
648
649  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: pheno_type_mtc  =  &     !! type of phenology (0-4, unitless)
650  & (/  0,   1,   3,   1,   1,   2,   1,  &                              !! 0=bare ground 1=evergreen,  2=summergreen,
651  &     2,   2,   4,   4,   2,   3,   2  /)                                   !! 3=raingreen,  4=perennial
652  !-
653  ! 2. Leaf Onset
654  !-
655  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_c_mtc  =  &    !! critical gdd, tabulated (C),
656  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &    !! constant c of aT^2+bT+c
657  &     undef,   undef,   320.0,   400.0,   320.0,   700.0,   320.0  /)
658
659  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_b_mtc  =  &    !! critical gdd, tabulated (C),
660  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &    !! constant b of aT^2+bT+c
661  &     undef,   undef,    6.25,     0.0,    6.25,     0.0,  6.25  /)
662
663  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_gdd_crit_a_mtc  =  &    !! critical gdd, tabulated (C),
664  & (/  undef,   undef,     undef,   undef,   undef,   undef,   undef,  &  !! constant a of aT^2+bT+c
665  &     undef,   undef,   0.03125,     0.0,  0.0315,   0.0,  0.0315  /)
666
667  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: pheno_moigdd_t_crit_mtc  = &  !! temperature threshold for C4 grass(C)
668  & (/  undef,   undef,     undef,   undef,   undef,   undef,   undef,  & 
669  &     undef,   undef,     undef,    22.0,   undef,   undef,   undef  /)
670
671  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ngd_crit_mtc  =  &            !! critical ngd, tabulated.
672  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &    !! Threshold -5 degrees (days)
673  &     undef,    17.0,   undef,   undef,   undef,   undef,   undef  /)
674
675  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ncdgdd_temp_mtc  =  &         !! critical temperature for the ncd vs. gdd
676  & (/  undef,   undef,   undef,   undef,   undef,     5.0,   undef,  &    !! function in phenology (C)
677  &       0.0,   undef,   undef,   undef,   undef,   undef,   undef  /)
678
679  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_frac_mtc  =  &            !! critical humidity (relative to min/max)
680  & (/  undef,   undef,   0.5,   undef,   undef,   undef,   undef, &       !! for phenology (unitless)
681  &     undef,   undef,   0.5,     0.5,     0.5,     0.5,   0.5  /)
682
683  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: hum_min_time_mtc  =  &        !! minimum time elapsed since
684  & (/  undef,   undef,   50.0,   undef,   undef,   undef,   undef,  &     !! moisture minimum (days)
685  &     undef,   undef,   36.0,    35.0,    75.0,    75.0,   75.0  /) 
686
687  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_sap_mtc  =  &             !! time (days) 
688  & (/  undef,   730.0,   730.0,   730.0,   730.0,   730.0,   730.0,  &
689  &     730.0,   730.0,   undef,   undef,   undef,   undef,   undef  /)
690
691  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_leafinit_mtc  =  &  !! time to attain the initial foliage using the carbohydrate reserve
692  & (/  undef,   10.,   10.,   10.,   10.,   10.,   10.,  &
693  &       10.,   10.,   10.,   10.,   10.,   10.,   10.  /)
694
695  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tau_fruit_mtc  =  &           !! fruit lifetime (days)
696  & (/  undef,  90.0,    90.0,    90.0,    90.0,   90.0,   90.0,  &
697  &      90.0,  90.0,   undef,   undef,   undef,   undef,  undef  /)
698
699  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: ecureuil_mtc  =  &            !! fraction of primary leaf and root allocation
700  & (/  undef,   0.0,   1.0,   0.0,   0.0,   1.0,   0.0,  &                !! put into reserve (0-1, unitless)
701  &       1.0,   1.0,   1.0,   1.0,   1.0,   1.0,   1.0  /)
702
703  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_min_mtc  =  &           !! NEW - allocation above/below = f(age)
704  & (/  undef,   0.2,     0.2,     0.2,     0.2,    0.2,   0.2,  &         !! - 30/01/04 NV/JO/PF
705  &       0.2,   0.2,   undef,   undef,   undef,   undef,  undef  /)
706
707  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: alloc_max_mtc  =  &           !! NEW - allocation above/below = f(age)
708  & (/  undef,   0.8,     0.8,     0.8,     0.8,    0.8,   0.8,  &         !! - 30/01/04 NV/JO/PF
709  &       0.8,   0.8,   undef,   undef,   undef,   undef,  undef  /)
710
711  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: demi_alloc_mtc  =  &          !! NEW - allocation above/below = f(age)
712  & (/  undef,   5.0,     5.0,     5.0,     5.0,    5.0,   5.0,  &         !! - 30/01/04 NV/JO/PF
713  &       5.0,   5.0,   undef,   undef,   undef,   undef,  undef  /)
714
715  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaflife_mtc  =  &            !! leaf longevity, tabulated (??units??)
716  & (/  undef,   0.5,   2.0,   0.33,   1.0,   2.0,   0.33,  &
717  &       2.0,   2.0,   2.0,   2.0,    2.0,   2.0,   2.0  /)
718  !-
719  ! 3. Senescence
720  !-
721  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leaffall_mtc  =  &             !! length of death of leaves, tabulated (days)
722  & (/  undef,   undef,   10.0,   undef,   undef,   30.0,   undef,  &
723  &       5.0,    10.0,   10.0,    10.0,    10.0,   10.0,   10.0  /)
724
725  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: leafagecrit_mtc  =  &          !! critical leaf age, tabulated (days)
726  & (/  undef,   730.0,   180.0,   910.0,   730.0,   160.0,   910.0,  &
727  &     220.0,   120.0,   80.0,   120.0,    90.0,    90.0,   90.0  /)
728
729  CHARACTER(LEN=6), PARAMETER, DIMENSION(nvmc) :: senescence_type_mtc  =  & !! type of senescence, tabulated (unitless)
730  & (/  'none  ',  'none  ',   'dry   ',  'none  ',  'none  ',  &
731  &     'cold  ',  'none  ',   'cold  ',  'cold  ',  'mixed ',  &
732  &     'mixed ',  'mixed ',   'mixed ',  'mixed '            /)
733
734  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_hum_mtc  =  &       !! critical relative moisture availability
735  & (/  undef,   undef,   0.3,   undef,   undef,   undef,   undef,  &       !! for senescence (0-1, unitless)
736  &     undef,   undef,   0.2,     0.2,     0.3,     0.2,   0.3  /)
737
738  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: nosenescence_hum_mtc  =  &     !! relative moisture availability above which
739  & (/  undef,   undef,   0.8,   undef,   undef,   undef,   undef,  &       !! there is no humidity-related senescence
740  &     undef,   undef,   0.6,     0.3,     0.3,     0.3,   0.3  /)                !! (0-1, unitless)
741
742  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: max_turnover_time_mtc  =  &    !! maximum turnover time for grasses (days)
743  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &
744  &     undef,   undef,    80.0,    80.0,    80.0,    80.0,   80.0  /)
745
746  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_turnover_time_mtc  =  &    !! minimum turnover time for grasses (days)
747  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &
748  &     undef,   undef,    10.0,    10.0,    10.0,    10.0,   10.0  /)
749 
750  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: min_leaf_age_for_senescence_mtc  =  &  !! minimum leaf age to allow
751  & (/  undef,   undef,   90.0,   undef,   undef,   90.0,   undef,  &               !! senescence g (days)
752  &      60.0,    60.0,   30.0,    30.0,    30.0,   30.0,   30.0  /)
753
754  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_c_mtc  =  &    !! critical temperature for senescence (C)
755  & (/  undef,   undef,    undef,   undef,   undef,  16.0,   undef,  &     !! constant c of aT^2+bT+c, tabulated
756  &      14.0,    10.0,      5.0,     5.0,    5.0,    10.0,   5.0  /)              !! (unitless)
757
758  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_b_mtc  =  &    !! critical temperature for senescence (C),
759  & (/  undef,   undef,   undef,   undef,   undef,   0.0,   undef,  &       !! constant b of aT^2+bT+c, tabulated
760  &       0.0,     0.0,     0.1,     0.0,     0.0,   0.0,   0.0  /)                !! (unitless)
761
762  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: senescence_temp_a_mtc  =  &    !! critical temperature for senescence (C),
763  & (/  undef,   undef,     undef,   undef,   undef,   0.0,   undef,  &     !! constant a of aT^2+bT+c, tabulated
764  &       0.0,     0.0,   0.00375,     0.0,     0.0,   0.0,   0.0  /)              !! (unitless)
765
766  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: gdd_senescence_mtc  =  &       !! minimum gdd to allow senescence of crops (days)
767  & (/  undef,   undef,    undef,   undef,     undef,    undef,    undef,  &
768  &     undef,   undef,    undef,   undef,      950.,    4000.,    950.  /)
769
770
771  !
772  ! DGVM
773  !
774  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: residence_time_mtc  =  &    !! residence time of trees (years)
775  & (/  undef,   30.0,   30.0,   40.0,   40.0,   40.0,   80.0,  &
776  &      80.0,   80.0,    0.0,    0.0,    0.0,    0.0,    0.0  /) 
777
778  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tmin_crit_mtc  =  &
779  & (/  undef,     0.0,     0.0,   -30.0,   -14.0,   -30.0,   -45.0,  &  !! critical tmin, tabulated (C)
780  &     -45.0,   -60.0,   undef,   undef,   undef,   undef,   undef  /)
781
782  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: tcm_crit_mtc  =  &
783  & (/  undef,   undef,   undef,     5.0,    15.5,    15.5,   -8.0,  &   !! critical tcm, tabulated (C)
784  &      -8.0,    -8.0,   undef,   undef,   undef,   undef,  undef  /)
785
786
787
788  !
789  ! Biogenic Volatile Organic Compounds
790  !
791  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_isoprene_mtc = &     !! Isoprene emission factor
792  & (/  0.,    24.,   24.,    8.,   16.,   45.,   8.,  &                    !!
793  &    18.,    0.5,   12.,   18.,    5.,    5.,   5.  /)
794
795  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_monoterpene_mtc = &  !! Monoterpene emission factor
796  & (/   0.,   2.0,    2.0,   1.8,    1.4,    1.6,    1.8,  &               !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
797  &    1.4,    1.8,    0.8,   0.8,    0.22,     0.22,  0.22  /)
798
799  REAL(r_std), PARAMETER :: LDF_mono_mtc = 0.6                                  !! monoterpenes fraction dependancy to light
800  REAL(r_std), PARAMETER :: LDF_sesq_mtc = 0.5                                  !! sesquiterpenes fraction dependancy to light
801  REAL(r_std), PARAMETER :: LDF_meth_mtc = 0.8                                  !! methanol fraction dependancy to light
802  REAL(r_std), PARAMETER :: LDF_acet_mtc = 0.2                                  !! acetone fraction dependancy to light
803
804  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_apinene_mtc = &      !! Alfa pinene emission factor percentage
805  & (/   0.,   0.395,   0.395,   0.354,   0.463,   0.326,   0.354, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
806  &   0.316,   0.662,   0.231,   0.200,   0.277,   0.277,   0.277 /)
807
808  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_bpinene_mtc = &      !! Beta pinene emission factor  percentage     
809  & (/   0.,   0.110,   0.110,   0.146,   0.122,   0.087,   0.146, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
810  &   0.063,   0.150,   0.123,   0.080,   0.154,   0.154,   0.154  /)
811
812  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_limonene_mtc = &     !! Limonene emission factor percentage
813  & (/   0.,   0.092,   0.092,   0.083,   0.122,   0.061,   0.083, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
814  &   0.071,   0.037,   0.146,   0.280,   0.092,   0.092,   0.092  /)
815
816  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_myrcene_mtc = &      !! Myrcene emission factor percentage
817  & (/   0.,   0.073,   0.073,   0.050,   0.054,   0.028,   0.050, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
818  &   0.019,   0.025,   0.062,   0.057,   0.046,   0.046,    0.046  /)
819
820  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sabinene_mtc = &     !! Sabinene emission factor percentage
821  & (/   0.,   0.073,   0.073,   0.050,   0.083,   0.304,   0.050, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
822  &   0.263,   0.030,   0.065,   0.050,   0.062,   0.062,   0.062  /)
823
824  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_camphene_mtc = &     !! Camphene emission factor percentage
825  & (/   0.,   0.055,   0.055,   0.042,   0.049,   0.004,   0.042, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
826  &   0.005,   0.023,   0.054,   0.053,   0.031,   0.031,   0.031  /)
827
828  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_3carene_mtc = &      !! 3-carene emission factor percentage
829  & (/   0.,   0.048,   0.048,   0.175,   0.010,   0.024,   0.175, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
830  &   0.013,   0.042,   0.065,   0.057,   0.200,   0.200,   0.200  /)
831
832  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_tbocimene_mtc = &    !! T-beta-ocimene emission factor percentage
833  & (/   0.,   0.092,   0.092,   0.054,   0.044,   0.113,   0.054, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
834  &   0.105,   0.028,   0.138,   0.120,   0.031,   0.031,   0.031  /)
835
836  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_othermonot_mtc = &   !! Other monoterpenes emission factor percentage
837  & (/   0.,   0.062,   0.062,   0.046,   0.054,   0.052,   0.046, &        !! ATTENTION: for each PFT they are PERCENTAGE of monoterpene EF
838  &   0.144,   0.003,   0.115,   0.103,   0.108,   0.108,   0.108  /)
839
840  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_sesquiterp_mtc = &   !! Sesquiterpene emission factor
841  & (/   0.,  0.45,   0.45,   0.13,   0.30,   0.36,   0.15, &               !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
842  &    0.30,  0.25,   0.60,   0.60,   0.08,   0.08,   0.08  /)
843
844  REAL(r_std), PARAMETER :: beta_mono_mtc = 0.10                            !! Monoterpenes temperature dependency coefficient
845  REAL(r_std), PARAMETER :: beta_sesq_mtc = 0.17                            !! Sesquiterpenes temperature dependency coefficient
846  REAL(r_std), PARAMETER :: beta_meth_mtc = 0.08                            !! Methanol temperature dependency coefficient
847  REAL(r_std), PARAMETER :: beta_acet_mtc = 0.10                            !! Acetone temperature dependency coefficient
848  REAL(r_std), PARAMETER :: beta_oxyVOC_mtc = 0.13                          !! Other oxygenated BVOC temperature dependency coefficient
849
850
851  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_ORVOC_mtc = &        !! ORVOC emissions factor
852  &  (/  0.,    1.5,    1.5,    1.5,    1.5,   1.5,    1.5,  &              !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
853  &     1.5,    1.5,    1.5,    1.5,    1.5,   1.5,    1.5  /) 
854
855  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_OVOC_mtc = &         !! OVOC emissions factor
856  &  (/  0.,    1.5,    1.5,    1.5,    1.5,   1.5,    1.5,  &              !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
857  &     1.5,    1.5,    1.5,    1.5,    1.5,   1.5,    1.5  /)
858 
859  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_MBO_mtc = &          !! MBO emissions factor
860  & (/     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
861  &     2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5, 2.e-5  /) 
862 
863  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_methanol_mtc = &     !! Methanol emissions factor
864  & (/  0.,    0.8,   0.8,   1.8,   0.9,   1.9,   1.8,  &                   !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
865  &    1.8,    1.8,   0.7,   0.9,    2.,     2.,   2.  /) 
866 
867  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetone_mtc = &      !! Acetone emissions factor
868  & (/  0.,   0.25,   0.25,   0.30,   0.20,   0.33,   0.30,  &              !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
869  &   0.25,   0.25,   0.20,   0.20,   0.08,   0.08,   0.08  /)
870 
871  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetal_mtc = &       !! Acetaldehyde emissions factor
872  & (/  0.,   0.2,    0.2,     0.2,   0.2,   0.25,   0.25,   0.16,   &      !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
873  &   0.16,   0.12,   0.12,   0.035,   0.020,  0.020  /) 
874 
875  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formal_mtc = &       !! Formaldehyde emissions factor
876  & (/  0.,   0.04,   0.04,  0.08,    0.04,    0.04,  0.04,  &              !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
877  &   0.04,   0.04,  0.025, 0.025,   0.013,   0.013,  0.013  /) 
878
879  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_acetic_mtc = &       !! Acetic Acid emissions factor
880  & (/   0.,   0.025,   0.025,   0.025,   0.022,   0.08,   0.025,   &      !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
881  &   0.022,   0.013,   0.012,   0.012,   0.008,   0.008,   0.008  /) 
882
883  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: em_factor_formic_mtc = &       !! Formic Acid emissions factor
884  & (/  0.,  0.015,  0.015,   0.02,    0.02,   0.025,  0.025,  &            !! @tex $(\mu gC.g^{-1}.h^{-1})$ @endtex
885  &  0.015,  0.015,  0.010,  0.010,   0.008,   0.008,   0.008  /) 
886
887  REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_wet_mtc = &        !! NOx emissions factor soil emissions and exponential
888  & (/  0.,   2.6,   0.06,   0.03,   0.03,   0.03,   0.03,  &               !! dependancy factor for wet soils
889  &  0.03,   0.03,   0.36,   0.36,   0.36,   0.36,   0.36  /)                       !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex
890
891  REAL(r_std),PARAMETER, DIMENSION(nvmc) :: em_factor_no_dry_mtc = &        !! NOx emissions factor soil emissions and exponential
892  & (/  0.,   8.60,   0.40,   0.22,   0.22,   0.22,   0.22,  &              !! dependancy factor for dry soils
893  &   0.22,   0.22,   2.65,   2.65,   2.65,   2.65,   2.65  /)                      !! @tex $(ngN.m^{-2}.s^{-1})$ @endtex
894
895  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: Larch_mtc = &                  !! Larcher 1991 SAI/LAI ratio (unitless)
896  & (/   0.,   0.015,   0.015,   0.003,   0.005,   0.005,   0.003,  &
897  &   0.005,   0.003,   0.005,   0.005,   0.008,   0.008,   0.008  /) 
898
899
900!
901!WETALND CH4 methane
902!
903!pss+
904  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: sdepth_v_mtc  = &
905  & (/ 0, 129, 129, 129, 129, 129, 129, 129, 129, 79, 79, 162, 162, 162/)
906
907  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: rdepth_v_mtc  = &
908  & (/ 0, 64, 64, 64, 64, 64, 64, 64, 64, 39, 39, 81, 81, 81 /)
909
910  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: tveg_v_mtc  = & 
911  & (/ 0, 1, 1, 1, 1, 1, 1, 1, 1, 10, 10, 15, 15, 15 /)
912!pss-
913
914!chaoyue+
915  LOGICAL,PARAMETER,DIMENSION(nvmc) :: permafrost_veg_exists_mtc = &  !!! set all as true to avoid masking
916   & (/ .TRUE., .TRUE., .TRUE., .TRUE., .TRUE.,  .TRUE., .TRUE., &
917   &    .TRUE., .TRUE.,  .TRUE.,  .TRUE., .TRUE., .TRUE., .TRUE./)
918!chaoyue-
919
920
921!!!!! crop parameters, xuhui
922
923  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: irrig_threshold_mtc  =  &  !! Value for irrig_threshold determine
924  & (/ 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                !! when irrigation started (0-1, vegstress)
925  &    0.0,   0.0,   0.0,   1.0,   1.0,   1.0,   1.0  /)                       !! only useful for non-natural pfts
926
927  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: irrig_fulfill_mtc  =  &  !! Value for irrig_fulfill determine
928  & (/ 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,  &                !! how irrigation demand fulfilled (0-1)
929  &    0.0,   0.0,   0.0,   0.9,   0.9,   0.9,   0.9  /)                       !! only useful for non-natural pfts
930
931
932  !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
933  !-
934  ! 4. STICS
935  !    This part is especially for Stics
936  !-
937  ! 4.0 flag in stomate
938  LOGICAL, PARAMETER, DIMENSION(nvmc) ::  ok_LAIdev_mtc =  &        !! flag for using the LAIdev module for vegetation types (true/false)
939  & (/ .FALSE.,  .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,  &
940  &    .FALSE.,  .FALSE.,   .FALSE.,   .FALSE.,    .FALSE.,    .FALSE.,   .FALSE.  /)
941
942  ! INTEGER(i_std), PARAMETER :: crop_nbjmax = 300                         !! Maximum days of a crop cycle simulation
943 
944  ! 4.1 LAIdev module
945
946  CHARACTER(len = 3), PARAMETER, DIMENSION(nvmc) ::  SP_codeplante_mtc =  &                     !! simulated plante
947  & (/ 'snu',  'snu',   'snu',   'snu',   'snu',   'snu',   'snu',  &
948  &    'snu',  'snu',   'snu',   'snu',    'ble',   'mai',  'ric'  /)
949  CHARACTER(len = 3), PARAMETER, DIMENSION(nvmc) :: SP_stade0_mtc = &
950  & (/ 'snu',  'snu',   'snu',   'snu',   'snu',   'snu',   'snu',  &
951  &    'snu',  'snu',   'snu',   'lev',    'snu',   'snu',  'lev'  /)
952
953
954  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_iplt0_mtc  =  &    !! sowing date
955  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
956  &       undef_int,        undef_int,        undef_int,      120,     284,   100, 120 /)     
957
958  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nbox_mtc  =  &    !! number of agebox for senescence
959  & (/  0,     0,     0,      0,        0,       0,       0,  &   
960  &       0,        0,        0,      3,     3,   3,  3 /)     
961
962  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_iwater_mtc  =  &    !! julian day of the beginning of the simulation
963  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
964  &       undef_int,        undef_int,        undef_int,      1,     1,   1, 1 /)     
965
966  CHARACTER(len = 7), PARAMETER, DIMENSION(nvmc) ::  SP_codesimul_mtc =  &                     !! simulated plante
967  & (/ 'feuille',  'feuille',  'feuille', 'feuille', 'feuille', 'feuille', 'feuille', &
968  &    'feuille',   'feuille', 'feuille', 'culture', 'culture', 'culture',  'culture'/)
969
970  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codelaitr_mtc  =  &    !! option of calculation of LAI: undef_int. LAI, 2. Soil cover
971  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
972  &       undef_int,        undef_int,        undef_int,      1,    1,   1,  1 /)     
973
974  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_slamax_mtc  =  &    !! maximum sla  cm2/g dry mass
975  & (/  undef,   undef,   undef,   undef,   undef,   undef,   undef,  &       !!
976  &       undef,     undef,     undef,     500.0,     350.0,  250.0, 500.0  /)                !!
977
978  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_slamin_mtc  =  &    !! minimum sla
979  & (/  undef,   undef,     undef,   undef,   undef,   undef,   undef,  &     !!
980  &       undef,     undef,    undef,     180.0,     180.0,   180.0,   180.0  /)              !!
981
982  ! STICS:: DEVELOPMENT
983  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeperenne_mtc  =  &    !! annual crop (1) or perennial crop (2)
984  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
985  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
986
987  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codcueille_mtc  =  &    !! harvest option: cutting (1) or picking (2)
988  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
989  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
990
991  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codegdh_mtc  =  &    !!  daily (1) or hourly (2) calculation of development unit
992  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
993  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
994
995  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codetemp_mtc  =  &    !!  air temperature (1) or crop_temperature (2)
996  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
997  &       undef_int,        undef_int,        undef_int,      1,    1,  1,   1 /)     
998 
999  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_coderetflo_mtc  =  &    !! slowing effect of water stress before DRP, YES (1) OR NO (2)
1000  & (/  undef_int,     undef_int,     undef_int,      undef_int,      undef_int,       undef_int,       undef_int,  &   
1001  &       undef_int,        undef_int,        undef_int,      2,     1,   1,   1/)     
1002
1003  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeinnact_mtc  =  &    !! slowing effect of nitrogen stress on crop, active (1) OR NOt active (2)
1004  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1005  &       undef_int,        undef_int,        undef_int,      1,     1,   1,   1 /)     
1006
1007  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeh2oact_mtc  =  &    !! slowing effect of water stress on crop, YES (1) OR NO (2)
1008  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1009  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1/)     
1010
1011  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stressdev_mtc  =  &    !! maximum phasic delay allowed  due to stresses  // SD
1012  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1013  &     undef,        undef,        undef,      0.20,     0.34,        0.20,   0.20 /)     
1014
1015  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_innlai_mtc  =  &    !! nitrogen limitation for leaf growth  // SD
1016  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1017  &     undef,        undef,        undef,      1.0,     0.50,        1.0,   1.0 /)     
1018 
1019  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_innsenes_mtc  =  &    !! nitrogen limitation for leaf senescence  // SD
1020  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1021  &     undef,        undef,        undef,      0.35,     1.0,        1.0,   0.35 /)     
1022
1023  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codebfroid_mtc  =  &    !! PARAMETER option of calculation of chilling requirements,  1 (no need), 2 (vernalization), 3 (development stage)
1024  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &           !!  code undef_int (vernalization), 3 (for perennial)
1025  &       undef_int,        undef_int,        undef_int,      1,    2,   1,   1 /)     
1026
1027
1028  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codephot_mtc  =  &    !!  option of plant photoperiodism: yes (1), no (2)
1029  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1030  &       undef_int,        undef_int,        undef_int,      1,     1,   2,  1 /)     
1031
1032  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codedormance_mtc  =  &    !!option of calculation of dormancy and chilling requirement // code 1/2 
1033  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1034  &       undef_int,        undef_int,        undef_int,      3,    3,   3,  3/)     
1035
1036  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codefauche_mtc  =  &    !! option of cut modes for forage crops: yes (1), no (2)
1037  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1038  &       undef_int,        undef_int,        undef_int,      2,    2,   2,  2 /)     
1039
1040  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codetempfauche_mtc  =  &    !!option of the reference temperature to compute cutting sum of temperatures :
1041  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &    !!  : upvt (1), udevair (2)
1042  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1/)     
1043
1044  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codlainet_mtc  =  &    !!calculation of the LAI (1 : direct LAInet; 2 : LAInet = gross LAI - senescent LAI
1045  & (/  1,     1,     1,      1,        1,       1,       1,  &   
1046  &       1,        1,        1,      3,    3,   3,   3 /)     
1047
1048  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeindetermin_mtc  =  &    !!option of  simulation of the leaf growth and fruit growth 
1049  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &    !!  indeterminate (2) or determinate (1)
1050  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1/)     
1051
1052  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeinitprec_mtc  =  &    !! reinitializing initial status in case of chaining simulations : yes (1), no (2) 
1053  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1054  &       undef_int,        undef_int,        undef_int,      1,    1,   1,  1 /)     
1055
1056  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_culturean_mtc  =  &    !! crop status:  1 = over 1 calendar year
1057  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &    !! other than 1  = on two calendar years (winter crop in northern hemisphere)
1058  &       undef_int,        undef_int,        undef_int,      1,     2,   1,  1/)      !! // code 0/1
1059
1060  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_jvc_mtc  =  &    !! Number of vernalizing days // day
1061  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1062  &       undef,        undef,        undef,      0.0,    38.0,   0.0,  0.0 /)     
1063
1064  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tfroid_mtc  =  &    !!optimal temperature for vernalisation // degree C
1065  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1066  &       undef,        undef,        undef,      6.5,     6.5,   6.5,  6.5 /)     
1067
1068  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_ampfroid_mtc  =  &    !! semi thermal amplitude thermique for vernalising effect // degree C
1069  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1070  &       undef,        undef,        undef,      10.,     10.,   undef, 10. /)     
1071
1072  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_jvcmini_mtc  =  &    !! Minimum number of vernalising days  // day
1073  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1074  &       undef,        undef,        undef,      undef,     7.,   undef,  undef /)     
1075
1076  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgmin_mtc  =  &    !!  Minimum threshold temperature used in emergence stage // degree C
1077  & (/  undef,     undef,     undef,      undef,        undef,       undef,       undef,  &   
1078  &       undef,        undef,        undef,      0.,    0.,   8.,  0. /)     
1079
1080
1081  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stpltger_mtc  =  &            !!Sum of development allowing germination // degree.days
1082  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1083  &       undef,    undef,    undef,    50.,    50.,    35.0,   50.0  /) 
1084
1085  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_profsem_mtc  =  &            !! Sowing depth //cm
1086  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1087  &       undef,    undef,    undef,    3.0,    2.0,   3.0,  3.0  /) 
1088
1089  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_propjgermin_mtc  =  &            !! minimal proportion of the duration P_nbjgerlim
1090  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &  !!  when the temperature is higher than the temperature threshold SP_Tdmax  // %
1091  &       undef,    undef,    undef,    1.0,    1.0,    1.0,  1.0  /) 
1092
1093  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tdmax_mtc  =  &       !! Maximum threshold temperature for development // degree C
1094  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1095  &       undef,    undef,    undef,    30.0,    28.0,    32.0,  30.0  /) 
1096
1097  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nbjgerlim_mtc  =  &    !! Threshold number of day after grain imbibition without germination lack // days
1098  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1099  &       undef_int,        undef_int,        undef_int,      50,     50,   50,  50 /)     
1100 
1101  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_densitesem_mtc  =  &            !! Sowing density  // plants.m-2
1102  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1103  &       undef,    undef,    undef,    66.0,    212.0,    10.0,  66.0  /) 
1104
1105  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_vigueurbat_mtc  =  &            !! indicator of plant vigor allowing to emerge through the crust
1106  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1107  &       undef,    undef,    undef,    1.0,    1.0,    1.0,  1.0  /) 
1108
1109  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codepluiepoquet_mtc  =  &    !!option to replace rainfall by irrigation at poquet depth
1110  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   !! in the case of poquet sowing // code 1/2
1111  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1/)     
1112
1113  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codehypo_mtc  =  &    !! option of simulation of a  phase of hypocotyl growth (1) or planting of plantlets (2)
1114  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1115  &       undef_int,        undef_int,        undef_int,      2,    1,   1, 2 /)     
1116
1117  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_elmax_mtc  =  &            !! Maximum elongation of the coleoptile in darkness condition // cm /
1118  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1119  &       undef,    undef,    undef,    6.5,    8.0,    8.0,  6.5  /) 
1120
1121  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_belong_mtc  =  &            !! parameter of the curve of coleoptile elongation // degree.days -1
1122  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1123  &       undef,    undef,    undef,    0.006,   0.012,    0.022,  0.006  /) 
1124
1125  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_celong_mtc  =  &            !! parameter of the subsoil plantlet elongation curve
1126  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1127  &       undef,    undef,    undef,    2.00,    3.2,    2.04,  2.00  /) 
1128
1129  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nlevlim1_mtc  =  &  !!number of days after germination decreasing the emerged plants if emergence has not occur
1130  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1131  &       undef_int,        undef_int,        undef_int,      10,     10,   10,  10 /)     
1132
1133  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nlevlim2_mtc  =  &    !!number of days after germination after which the emerged plants is 0
1134  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1135  &       undef_int,        undef_int,        undef_int,      50,     50,   50,  50 /)     
1136
1137  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codrecolte_mtc  =  &    !! harvest mode : all the plant (1) or just the fruits (2)
1138  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1139  &       undef_int,        undef_int,        undef_int,      1,    1,  1,  1 /)     
1140
1141  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_variete_mtc  =  &    !!variety number in the technical file
1142  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1143  &       undef_int,        undef_int,        undef_int,      1,    1,  1,  1 /)     
1144
1145  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codegermin_mtc  =  &    !! option of simulation of a germination phase
1146  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &    !! or a delay at the beginning of the crop (undef_int)   or  direct starting (2)
1147  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
1148
1149
1150  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: S_codeulaivernal_mtc  =  &    !! sensitive or not to vernalization
1151  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &    !!
1152  &       undef_int,        undef_int,        undef_int,      0,     1,   0,  0 /)     
1153
1154
1155  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_swfacmin_mtc  =  &    !! minimum of drought stress
1156  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1157  &       undef,    undef,    undef,    0.3,    0.3,    0.3,   0.3  /) 
1158
1159  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_neffmax_mtc  =  &    !! maximum nitrogen effect
1160  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1161  &       undef,    undef,    undef,    0.65,    0.65,    0.65,  0.65  /) 
1162
1163  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_nsatrat_mtc  =  &    !! maximum nitrogen effect
1164  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1165  &       undef,    undef,    undef,    0.91,    0.91,    0.91,   0.91  /) 
1166
1167
1168  ! STICS:: LAI CALCULATION
1169
1170  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_laiplantule_mtc  =  &    !! Plantlet Leaf index at the plantation // m2 leaf  m-2 soil
1171  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1172  &       undef,    undef,    undef,    0.0,    1.0,    1.0,   0.0  /) 
1173
1174  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_vlaimax_mtc  =  &            !!  ULAI at inflection point of the function DELTAI=f(ULAI)
1175  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1176  &       undef,    undef,    undef,    1.5,    2.2,    2.2,  1.5 /) 
1177
1178  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stlevamf_mtc  =  &      !!Sum of development units between the stages LEV and AMF // degree.days
1179  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1180  &       undef,    undef,    undef,    571.0,    185.0,    225.0,  571.0  /) 
1181
1182  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stamflax_mtc  =  &    !! Sum of development units between the stages AMF and LAX // degree.days
1183  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1184  &       undef,    undef,    undef,    409.0,    280.0,    450.0,  409.0  /) 
1185
1186  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_udlaimax_mtc  =  &    !! ulai from which the rate of leaf growth decreases
1187  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1188  &       undef,    undef,    undef,    3.0,    3.0,    3.0,  3.0  /) 
1189
1190  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_laicomp_mtc  =  &     !! LAI from which starts competition inbetween plants // m2 m-2 /
1191  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1192  &       undef,    undef,    undef,    0.305,    0.304,    0.7,  0.305  /) 
1193
1194  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_adens_mtc  =  &     !!  Interplant competition parameter
1195  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1196  &       undef,    undef,    undef,    -0.85,    -0.60,    -0.12,  -0.85  /) 
1197
1198  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_bdens_mtc  =  &     !! minimal density from which interplant competition starts // plants m-2
1199  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1200  &       undef,    undef,    undef,    2.0,    15.0,    5.0,  2.0  /) 
1201
1202  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tcxstop_mtc  =  &     !!threshold temperature beyond which the foliar growth stops // degree C
1203  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1204  &       undef,    undef,    undef,    40.0,    100.0,    35.0,  40.0  /) 
1205
1206  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tcmax_mtc  =  &            !! Maximum temperature of growth
1207  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1208  &       undef,    undef,    undef,    30.0,    40.0,    32.0 ,  30.0 /) 
1209
1210  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tcmin_mtc  =  &            !!Minimum temperature of growth
1211  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1212  &       undef,    undef,    undef,    13.0,    0.0,    6.0,  13.0  /) 
1213
1214  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_dlaimax_mtc  =  &     !! accelerating parameter for the lai growth rate for maize 0.0016
1215  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1216  &       undef,    undef,    undef,    0.0060,    0.00044,    0.0020,  0.0060 /) 
1217
1218  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_dlaimin_mtc  =  &      !!accelerating parameter for the lai growth rate
1219  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1220  &       undef,    undef,    undef,    0.00,    0.00,    0.0012,  0.00  /) 
1221
1222  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_pentlaimax_mtc  =  &    !!parameter of the logistic curve of LAI growth
1223  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1224  &       undef,    undef,    undef,    5.5,    5.5,    5.5,  5.5 /) 
1225
1226
1227  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tigefeuil_mtc  =  &     !! stem (structural part)/leaf proportion // SD
1228  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1229  &       undef,    undef,    undef,    0.5,    0.5,    0.25,  0.5  /) 
1230 
1231  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stlaxsen_mtc  =  &      !!Sum of development units between the stages LEV and AMF // degree.days
1232  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1233  &       undef,    undef,    undef,    390.0,    265.0,    600.0,   390.0  /) 
1234
1235  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stsenlan_mtc  =  &    !! Sum of development units between the stages AMF and LAX // degree.days
1236  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1237  &       undef,    undef,    undef,    89.0,    620.0,    272.0,   89.0  /) 
1238
1239  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stlevdrp_mtc  =  &      !!Sum of development units between the stages LEV and AMF // degree.days
1240  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1241  &       undef,    undef,    undef,    959.0,    690.0,    960.0,   959.0  /) ! 540 for wheat
1242
1243  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stflodrp_mtc  =  &    !! Sum of development units between the stages AMF and LAX // degree.days
1244  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1245  &       undef,    undef,    undef,    50.0,    0.0,    250.0,   50.0  /) 
1246
1247  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stdrpmat_mtc  =  &      !!Sum of development units between the stages LEV and AMF // degree.days
1248  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1249  &       undef,    undef,    undef,    442.0,    940.0,    650.0,   442.0  /)    ! 750 for wheat
1250
1251  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_stdrpdes_mtc  =  &    !! Sum of development units between the stages AMF and LAX // degree.days
1252  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1253  &       undef,    undef,    undef,    350.0,    750.0,    600.0,  350.0  /) 
1254
1255  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_phyllotherme_mtc  =  &    !!
1256  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1257  &       undef,    undef,    undef,    120.0,    120.0,    70.0,   120.0  /) 
1258
1259  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_lai0_mtc  =  &       !!  initial LAI
1260  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1261  &       undef,    undef,    undef,    0.61,    0.0,    0.0,  0.61  /) 
1262
1263  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tustressmin_mtc  =  &       !!rice does not take tustress approach
1264  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1265  &       undef,    undef,    undef,    0.7,    0.7,    0.7,   0.7  /) 
1266 
1267  ! STICS:: LAI SENESCENCE
1268
1269  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nbfgellev_mtc  =  &    !!  leaf number at the end of the juvenile phase (frost sensitivity)  // nb pl-1
1270  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1271  &       undef_int,        undef_int,        undef_int,      2,     2,   50,  2 /)     
1272
1273  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_ratiodurvieI_mtc  =  &      !!  life span of early leaves expressed as a proportion of
1274  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &   !!  the life span of the last leav    es emitted SP_DURVIEF
1275  &       undef,    undef,    undef,    0.8,    0.8,    0.8,   0.8  /) 
1276 
1277  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_durvieF_mtc  =  &       !!  maximal  lifespan of an adult leaf expressed in summation of P_Q10=2 (2**(T-Tbase))
1278  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1279  &       undef,    undef,    undef,    480.0,    480.0,    580.0,  480.0  /) 
1280 
1281
1282  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_ratiosen_mtc  =  &     !!  fraction of senescent biomass (by ratio at the total biomass) // between 0 and 1
1283  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1284  &       undef,    undef,    undef,    0.28,    0.8,    0.8,  0.28  /) 
1285
1286  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tdmin_mtc  =  &       !!Minimum threshold temperature for development // degree C
1287  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1288  &       undef,    undef,    undef,    13.0,    0.0,    6.0,  13.0  /) 
1289 
1290
1291
1292  ! STICS:: f_humerac
1293  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_sensrsec_mtc  =  & 
1294  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1295  &       undef,    undef,    undef,    0.8,    0.5,    0.5,   0.8  /) 
1296
1297  ! STICS:: gel
1298  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codgellev_mtc  =  &   
1299  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1300  &       undef_int,        undef_int,        undef_int,      2,     1,   2,  2 /)     
1301
1302
1303  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codgeljul_mtc  =  &   
1304  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1305  &       undef_int,        undef_int,        undef_int,      2,     1,   2,   2 /)     
1306  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codgelveg_mtc  =  &   
1307  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1308  &       undef_int,        undef_int,        undef_int,      2,     1,   2,  2 /)     
1309
1310
1311
1312  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tletale_mtc  =  & 
1313  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1314  &       undef,    undef,    undef,    -25.0,    -25.0,    -5.0,  -25.0  /) 
1315 
1316  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tdebgel_mtc  =  & 
1317  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1318  &       undef,    undef,    undef,    -4.0,    -4.0,    0.0,   -4.0  /) 
1319 
1320  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgellev10_mtc  =  & 
1321  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1322  &       undef,    undef,    undef,    -4.0,    -4.0,    undef,   -4.0  /) 
1323
1324  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgellev90_mtc  =  & 
1325  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1326  &       undef,    undef,    undef,    -20.0,    -20.0,    undef,  -20.0  /) 
1327
1328
1329  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgeljuv10_mtc  =  & 
1330  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1331  &       undef,    undef,    undef,    -10.0,    -10.0,    undef,   -10.0  /) 
1332
1333  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgeljuv90_mtc  =  & 
1334  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1335  &       undef,    undef,    undef,    -20.0,    -20.0,    undef,   -20.0  /) 
1336  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgelveg10_mtc  =  & 
1337  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1338  &       undef,    undef,    undef,    -4.5,    -4.5,    undef,   -4.5  /) 
1339
1340  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgelveg90_mtc  =  & 
1341  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1342  &       undef,    undef,    undef,    -10.0,    -10.0,    undef,   -10.0  /) 
1343
1344
1345  ! STICS:: Photoperiod
1346
1347  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_sensiphot_mtc  =  & 
1348  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1349  &       undef,    undef,    undef,    1.0,    0.0,    1.0,  1.0 /) 
1350 
1351  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_phosat_mtc  =  & 
1352  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1353  &       undef,    undef,    undef,    20.0,    20.0,    undef,  20.0  /) 
1354
1355  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_phobase_mtc  =  & 
1356  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1357  &       undef,    undef,    undef,    6.3,    6.3,    undef,  6.3  /) 
1358
1359  ! STICS:: CARBON ALLOCATION
1360
1361  CHARACTER(len=3), PARAMETER, DIMENSION(nvmc) :: SP_stoprac_mtc = &  !! description of the MTC (unitless)
1362  & (/ 'GER', &          !  1
1363  &    'GER', &          !  2
1364  &    'GER', &          !  3
1365  &    'GER', &          !  4
1366  &    'GER', &          !  5
1367  &    'GER', &          !  6
1368  &    'GER', &          !  7
1369  &    'GER', &          !  8
1370  &    'GER', &          !  9
1371  &    'GER', &          ! 10
1372  &    'SEN', &          ! 11
1373  &    'SEN', &          ! 12 wheat
1374  &    'LAX', &          ! 13 maize
1375  &    'SEN'  /)         ! 14 rice
1376 
1377  ! initial root depth for planting plant // cm
1378  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_zracplantule_mtc  =  & 
1379  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1380  &       undef,    undef,    undef,    3.0,    3.0,    3.0,  3.0  /) 
1381 
1382  ! trophic effect on root partitioning within the soil // code 1/2/3
1383  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codtrophrac_mtc  =  &   
1384  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1385  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
1386 
1387  !  maximum of root biomass respective to the total biomass (permanent trophic link) // SD /
1388  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_repracpermax_mtc  =  &                      ! not correct
1389  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1390  &       undef,    undef,    undef,    0.24,    0.4,    0.4,   0.24  /) 
1391
1392
1393  ! minimum of root biomass respective to the total biomass (permanent trophic link) // SD /
1394  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_repracpermin_mtc  =  &                      ! not correct
1395  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1396  &       undef,    undef,    undef,    0.047,    0.2,    0.2,    0.047  /) 
1397  !  parameter of biomass root partitioning : evolution of the ratio root/total (permanent trophic link)
1398  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_krepracperm_mtc  =  &                       ! not correct
1399  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1400  &       undef,    undef,    undef,    1.27,    1.27,     1.27,   1.27  /) 
1401
1402  ! maximum of root biomass respective to the total biomass (trophic link by thresholds)
1403  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_repracseumax_mtc  =  &                      ! not correct
1404  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1405  &       undef,    undef,    undef,    0.24,    0.4,    0.4,   0.24  /) 
1406  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_repracseumin_mtc  =  &                      ! not correct
1407  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1408  &       undef,    undef,    undef,    0.047,    0.2,    0.2,   0.047  /)     
1409  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_krepracseu_mtc  =  &                        ! not correct
1410  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1411  &       undef,    undef,    undef,    1.27,    1.27,    1.27,    1.27  /) 
1412 
1413  ! option calculation mode of heat time for the root: with crop temperature (1)  or with soil temperature (2)
1414  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codetemprac_mtc  =  &   
1415  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1416  &       undef_int,        undef_int,        undef_int,      1,     1,   1,   1 /)
1417
1418  ! Activation of the module simulating tiller dynamic: yes (1), no (2) /     
1419  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codedyntalle_mtc  =  &   
1420  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1421  &       undef_int,        undef_int,        undef_int,      1,     1,   2,   1 /)
1422
1423  ! < // PARAMETER // Period to compute NBGRAIN // days // PARPLT     
1424  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_nbjgrain_mtc  =  &   
1425  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1426  &       undef_int,        undef_int,        undef_int,      30,     30,   20,   30 /) 
1427
1428  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_maxgs_mtc  =  &   
1429  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1430  &       undef_int,        undef_int,        undef_int,      300,     300,   235,  300 /) 
1431 
1432  ! activation of frost at anthesis // code 1/2     1, no, 2, yes
1433  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codgelflo_mtc  =  &   
1434  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1435  &       undef_int,        undef_int,        undef_int,      2,     2,   1,  2 /)     
1436  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgelflo10_mtc  =  & 
1437  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1438  &       undef,    undef,    undef,    -4.5,    -4.5,    undef,  -4.5 /) 
1439  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tgelflo90_mtc  =  & 
1440  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1441  &       undef,    undef,    undef,    -6.5,    -6.5,    undef, -6.5  /) 
1442  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_cgrain_mtc  =  & 
1443  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1444  &       undef,    undef,    undef,    0.035,    0.036,    0.05,  0.035  /) 
1445  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_cgrainv0_mtc  =  & 
1446  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1447  &       undef,    undef,    undef,    0.350,    0.0,    0.111,  0.350  /) 
1448
1449  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_nbgrmax_mtc  =  & 
1450  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1451  &       undef,    undef,    undef,    40000.0,    30000.0,    4000.0,  40000.0  /) 
1452  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_nbgrmin_mtc  =  & 
1453  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1454  &       undef,    undef,    undef,    10000.0,    6000.0,    1500.0,  10000.0  /) 
1455
1456  ! option of activation of the direct effect of the nitrogen plant status upon the fruit/grain number , 1: no; 2, yes
1457  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codazofruit_mtc  =  &   
1458  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1459  &       undef_int,        undef_int,        undef_int,      1,     1,   1,   1 /)     
1460  ! option of computing the ratio grain weight/total biomass: proportional to time(1), proportional to sum temperatures (2)
1461  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codeir_mtc  =  &   
1462  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1463  &       undef_int,        undef_int,        undef_int,      2,     1,   1,   2 /)     
1464
1465  ! Rate of increase of the carbon harvest index // g grain g plant -1 day-1
1466  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_vitircarb_mtc  =  & 
1467  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1468  &       undef,    undef,    undef,    0.015,    0.014,    0.020,  0.015  /) ! 0.0107
1469  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_irmax_mtc  =  & 
1470  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1471  &       undef,    undef,    undef,    0.65,    0.55,    0.53,   0.65  /) 
1472  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_vitircarbT_mtc  =  & 
1473  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1474  &       undef,    undef,    undef,    0.0012,    0.0007,    0.00110,  0.0012  /) 
1475 
1476  !  option of heat effect on grain filling: yes (1), no (2) // code 1/2
1477  INTEGER(i_std), PARAMETER, DIMENSION(nvmc) :: SP_codetremp_mtc  =  &   
1478  & (/  undef_int,     undef_int,     undef_int,      undef_int,        undef_int,       undef_int,       undef_int,  &   
1479  &       undef_int,        undef_int,        undef_int,      1,     1,   1,  1 /)     
1480 
1481  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tminremp_mtc  =  & 
1482  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1483  &       undef,    undef,    undef,    14.0,    0.0,    0.0,  14.0  /) 
1484
1485  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_tmaxremp_mtc  =  & 
1486  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1487  &       undef,    undef,    undef,    38.0,    38.0,    42.0,  38.0  /) 
1488 
1489  ! Maximum weight of one grain (at 0% water content) // g // PARPLT // 1, it is depedent on variety (wheat: Thesee, maise: Furio---variety)
1490  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_pgrainmaxi_mtc  =  & 
1491  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1492  &       undef,    undef,    undef,    0.03,    0.052,      0.33,  0.03  /) 
1493
1494  !!===============================SPECIFIC FOR DYNAMIC NITROGEN FOR CROP GROWTH ======================================
1495  LOGICAL, PARAMETER, DIMENSION(nvmc) :: SP_DY_INN_mtc =  &                     !! flag for using the LAIdev module for vegetation types (true/false)
1496  & (/ .FALSE.,  .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,   .FALSE.,  &
1497  &    .FALSE.,  .FALSE.,   .FALSE.,   .FALSE.,    .FALSE.,   .FALSE.,  .FALSE.  /)
1498
1499  ! An average nitrogen fertilization amount for crop // unit is kg N ha-1 year-1
1500  REAL(r_std), PARAMETER, DIMENSION(nvmc) :: SP_avenfert_mtc  =  & 
1501  & (/    undef,    undef,    undef,    undef,    undef,    undef,   undef,  &
1502  &       undef,    undef,    undef,    150.0,   100.0,      100.0,  150.0  /) 
1503
1504  ! >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
1505
1506
1507!!!!! end crop parameters, xuhui
1508
1509END MODULE constantes_mtc
Note: See TracBrowser for help on using the repository browser.