1 | !$Header: /home/ssipsl/CVSREP/ORCHIDEE/src_parameters/constantes.f90,v 1.16 2007/08/01 15:19:05 ssipsl Exp $ |
---|
2 | !IPSL (2006) |
---|
3 | ! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC |
---|
4 | !- |
---|
5 | MODULE constantes |
---|
6 | !!-------------------------------------------------------------------- |
---|
7 | !! "constantes" module contains some public technical constants |
---|
8 | !!-------------------------------------------------------------------- |
---|
9 | USE defprec |
---|
10 | USE ioipsl |
---|
11 | !- |
---|
12 | IMPLICIT NONE |
---|
13 | !- |
---|
14 | |
---|
15 | !------------------------- |
---|
16 | ! ORCHIDEE CONSTANTS |
---|
17 | !------------------------ |
---|
18 | |
---|
19 | !---------------- |
---|
20 | ! Global |
---|
21 | !---------------- |
---|
22 | |
---|
23 | ! Unit for output messages |
---|
24 | INTEGER(i_std), SAVE :: numout = 6 |
---|
25 | !- |
---|
26 | ! To set for more printing |
---|
27 | LOGICAL,SAVE :: long_print = .FALSE. |
---|
28 | !- |
---|
29 | ! One of the most frequent problems is a temperature out of range |
---|
30 | ! we provide here a way to catch that in the calling procedure. (JP) |
---|
31 | LOGICAL,PARAMETER :: diag_qsat = .TRUE. |
---|
32 | |
---|
33 | !!$ ! One of the most frequent problems is a temperature out of range |
---|
34 | !!$ ! we provide here a way to catch that in the calling procedure. (JP) |
---|
35 | !!$ LOGICAL,SAVE :: diag_qsat = .TRUE. |
---|
36 | |
---|
37 | !- |
---|
38 | ! Selects the type of output for the model. |
---|
39 | ! Value is read from run.def in intersurf_history. |
---|
40 | LOGICAL :: almaoutput |
---|
41 | |
---|
42 | !- |
---|
43 | ! One day in seconds |
---|
44 | REAL(r_std),SAVE :: one_day |
---|
45 | ! One year in seconds |
---|
46 | REAL(r_std),SAVE :: one_year |
---|
47 | |
---|
48 | ! undef integer for integer arrays |
---|
49 | INTEGER(i_std), PARAMETER :: undef_integer = 999999999 |
---|
50 | |
---|
51 | ! Specific value if no restart value |
---|
52 | REAL(r_std),SAVE :: val_exp = 999999. |
---|
53 | |
---|
54 | ! Special value for stomate |
---|
55 | REAL(r_std),PARAMETER :: undef = -9999. |
---|
56 | |
---|
57 | ! Epsilon to detect a near zero floating point |
---|
58 | REAL(r_std),PARAMETER :: min_sechiba = 1.E-8_r_std |
---|
59 | ! The undef value used in SECHIBA |
---|
60 | REAL(r_std),PARAMETER :: undef_sechiba = 1.E+20_r_std |
---|
61 | |
---|
62 | ! Epsilon to detect a near zero floating point |
---|
63 | REAL(r_std),PARAMETER :: min_stomate = 1.E-8_r_std |
---|
64 | ! some large value (for stomate) |
---|
65 | REAL(r_std),PARAMETER :: large_value = 1.E33_r_std |
---|
66 | |
---|
67 | !- |
---|
68 | TYPE control_type |
---|
69 | LOGICAL :: river_routing |
---|
70 | LOGICAL :: hydrol_cwrr |
---|
71 | LOGICAL :: ok_sechiba |
---|
72 | LOGICAL :: ok_co2 |
---|
73 | LOGICAL :: ok_stomate |
---|
74 | LOGICAL :: ok_dgvm |
---|
75 | LOGICAL :: stomate_watchout |
---|
76 | LOGICAL :: ok_pheno |
---|
77 | END TYPE control_type |
---|
78 | |
---|
79 | ! Flags that (de)activate parts of the model |
---|
80 | TYPE(control_type),SAVE :: control |
---|
81 | !- |
---|
82 | |
---|
83 | !--------------------------------------- |
---|
84 | ! DIMENSIONING AND INDICES PARAMETERS |
---|
85 | !--------------------------------------- |
---|
86 | |
---|
87 | !---------------- |
---|
88 | ! qsat_moisture |
---|
89 | !---------------- |
---|
90 | ! Number of other surface types: land ice (lakes,cities, ...) |
---|
91 | INTEGER(i_std),PARAMETER :: nnobio=1 |
---|
92 | !- |
---|
93 | ! Index for land ice (see nnobio) |
---|
94 | INTEGER(i_std),PARAMETER :: iice = 1 |
---|
95 | |
---|
96 | !------- |
---|
97 | ! Soil |
---|
98 | !------- |
---|
99 | ! Number of soil level |
---|
100 | INTEGER(i_std),PARAMETER :: ngrnd=7 |
---|
101 | !- |
---|
102 | ! Number of diagnostic levels in the soil |
---|
103 | INTEGER(i_std),PARAMETER :: nbdl=11 |
---|
104 | !MM : if you want to compare hydrology variables with old TAG 1.6 and lower, |
---|
105 | ! you must set the Number of diagnostic levels in the soil to 6 : |
---|
106 | ! INTEGER(i_std),PARAMETER :: nbdl=6 |
---|
107 | !- |
---|
108 | ! Number of levels in CWRR |
---|
109 | INTEGER(i_std),PARAMETER :: nslm=11 |
---|
110 | !- |
---|
111 | ! Number of soil types |
---|
112 | INTEGER(i_std),PARAMETER :: nstm = 3 |
---|
113 | !- |
---|
114 | ! Dimensioning parameter for the soil color numbers and their albedo |
---|
115 | INTEGER(i_std), PARAMETER :: classnb = 9 |
---|
116 | |
---|
117 | !- |
---|
118 | ! Diagnostic variables |
---|
119 | !- |
---|
120 | ! The lower limit of the layer on which soil moisture (relative) |
---|
121 | ! and temperature are going to be diagnosed. |
---|
122 | ! These variables are made for transfering the information |
---|
123 | ! to the biogeophyical processes modelled in STOMATE. |
---|
124 | !- |
---|
125 | REAL(r_std),DIMENSION(nbdl),SAVE :: diaglev |
---|
126 | |
---|
127 | !----------------- |
---|
128 | ! STOMATE - LPJ |
---|
129 | !----------------- |
---|
130 | |
---|
131 | ! NV080800 Name of STOMATE forcing file |
---|
132 | CHARACTER(LEN=100) :: stomate_forcing_name='NONE' |
---|
133 | !- |
---|
134 | ! NV080800 Name of soil forcing file |
---|
135 | CHARACTER(LEN=100) :: stomate_Cforcing_name='NONE' |
---|
136 | !- |
---|
137 | INTEGER(i_std),SAVE :: forcing_id |
---|
138 | !- |
---|
139 | ! leaf age discretisation ( 1 = no discretisation ) |
---|
140 | INTEGER(i_std),PARAMETER :: nleafages = 4 |
---|
141 | ! |
---|
142 | !---------------------------- |
---|
143 | ! litter fractions: indices |
---|
144 | !---------------------------- |
---|
145 | INTEGER(i_std),PARAMETER :: ileaf = 1 |
---|
146 | INTEGER(i_std),PARAMETER :: isapabove = 2 |
---|
147 | INTEGER(i_std),PARAMETER :: isapbelow = 3 |
---|
148 | INTEGER(i_std),PARAMETER :: iheartabove = 4 |
---|
149 | INTEGER(i_std),PARAMETER :: iheartbelow = 5 |
---|
150 | INTEGER(i_std),PARAMETER :: iroot = 6 |
---|
151 | INTEGER(i_std),PARAMETER :: ifruit = 7 |
---|
152 | INTEGER(i_std),PARAMETER :: icarbres = 8 |
---|
153 | INTEGER(i_std),PARAMETER :: nparts = 8 |
---|
154 | ! |
---|
155 | !------------------------------------- |
---|
156 | ! indices for assimilation parameters |
---|
157 | !------------------------------------- |
---|
158 | INTEGER(i_std),PARAMETER :: itmin = 1 |
---|
159 | INTEGER(i_std),PARAMETER :: itopt = 2 |
---|
160 | INTEGER(i_std),PARAMETER :: itmax = 3 |
---|
161 | INTEGER(i_std),PARAMETER :: ivcmax = 4 |
---|
162 | INTEGER(i_std),PARAMETER :: ivjmax = 5 |
---|
163 | INTEGER(i_std),PARAMETER :: npco2 = 5 |
---|
164 | !- |
---|
165 | !------------------------------------------ |
---|
166 | ! trees and litter: indices for the parts of heart- and sapwood above |
---|
167 | ! and below the ground |
---|
168 | !----------------------------------------- |
---|
169 | INTEGER(i_std),PARAMETER :: iabove = 1 |
---|
170 | INTEGER(i_std),PARAMETER :: ibelow = 2 |
---|
171 | INTEGER(i_std),PARAMETER :: nlevs = 2 |
---|
172 | !- |
---|
173 | !--------------------------------------------------- |
---|
174 | ! litter: indices for metabolic and structural part |
---|
175 | !-------------------------------------------------- |
---|
176 | INTEGER(i_std),PARAMETER :: imetabolic = 1 |
---|
177 | INTEGER(i_std),PARAMETER :: istructural = 2 |
---|
178 | INTEGER(i_std),PARAMETER :: nlitt = 2 |
---|
179 | ! |
---|
180 | !----------------------- |
---|
181 | ! carbon pools: indices |
---|
182 | !----------------------- |
---|
183 | INTEGER(i_std),PARAMETER :: iactive = 1 |
---|
184 | INTEGER(i_std),PARAMETER :: islow = 2 |
---|
185 | INTEGER(i_std),PARAMETER :: ipassive = 3 |
---|
186 | INTEGER(i_std),PARAMETER :: ncarb = 3 |
---|
187 | |
---|
188 | |
---|
189 | !------------------------------ |
---|
190 | ! MATH AND PHYSICS CONSTANTS |
---|
191 | !------------------------------ |
---|
192 | |
---|
193 | !------------------------------------ |
---|
194 | ! 1 . Maths and numerical constants |
---|
195 | !------------------------------------ |
---|
196 | ! pi |
---|
197 | REAL(r_std), PARAMETER :: pi = 4.*ATAN(1.) |
---|
198 | ! e |
---|
199 | REAL(r_std),PARAMETER :: euler = 2.71828182846 |
---|
200 | !- |
---|
201 | ! Integer constant set to zero |
---|
202 | INTEGER(i_std), PARAMETER :: zero_int = 0 |
---|
203 | !- |
---|
204 | ! Numerical constant set to 0 |
---|
205 | REAL(r_std),PARAMETER :: zero = 0._r_std |
---|
206 | ! Numerical constant set to 1/2 |
---|
207 | REAL(r_std),PARAMETER :: undemi = 0.5_r_std |
---|
208 | ! Numerical constant set to 1 |
---|
209 | REAL(r_std),PARAMETER :: un = 1._r_std |
---|
210 | ! Numerical constant set to -1 |
---|
211 | REAL(r_std),PARAMETER :: moins_un = -1._r_std |
---|
212 | ! Numerical constant set to 2 |
---|
213 | REAL(r_std),PARAMETER :: deux = 2._r_std |
---|
214 | ! Numerical constant set to 3 |
---|
215 | REAL(r_std),PARAMETER :: trois = 3._r_std |
---|
216 | ! Numerical constant set to 4 |
---|
217 | REAL(r_std),PARAMETER :: quatre = 4._r_std |
---|
218 | ! Numerical constant set to 5 |
---|
219 | REAL(r_std),PARAMETER :: cinq = 5._r_std |
---|
220 | ! Numerical constant set to 6 |
---|
221 | REAL(r_std),PARAMETER :: six = 6._r_std |
---|
222 | ! Numerical constant set to 8 |
---|
223 | REAL(r_std),PARAMETER :: huit = 8._r_std |
---|
224 | ! Numerical constant set to 1000 |
---|
225 | REAL(r_std),PARAMETER :: mille = 1000._r_std |
---|
226 | |
---|
227 | !--------------- |
---|
228 | ! 2 . Physics |
---|
229 | !--------------- |
---|
230 | ! |
---|
231 | ! radius of the Earth (m) |
---|
232 | REAL(r_std), PARAMETER :: R_Earth = 6378000. |
---|
233 | ! standard pressure |
---|
234 | REAL(r_std), PARAMETER :: pb_std = 1013. |
---|
235 | !- |
---|
236 | ! Freezing point |
---|
237 | REAL(r_std),PARAMETER :: ZeroCelsius = 273.15 |
---|
238 | !- |
---|
239 | ! 0 degre Celsius in degre Kelvin |
---|
240 | REAL(r_std),PARAMETER :: tp_00=273.15 |
---|
241 | !- |
---|
242 | ! Latent heat of sublimation |
---|
243 | REAL(r_std),PARAMETER :: chalsu0 = 2.8345E06 |
---|
244 | ! Latent heat of evaporation |
---|
245 | REAL(r_std),PARAMETER :: chalev0 = 2.5008E06 |
---|
246 | ! Latent heat of fusion |
---|
247 | REAL(r_std),PARAMETER :: chalfu0 = chalsu0-chalev0 |
---|
248 | !- |
---|
249 | ! Stefan-Boltzman constant |
---|
250 | REAL(r_std),PARAMETER :: c_stefan = 5.6697E-8 |
---|
251 | ! Specific heat of air |
---|
252 | REAL(r_std),PARAMETER :: cp_air = 1004.675 |
---|
253 | ! Constante molere |
---|
254 | REAL(r_std),PARAMETER :: cte_molr = 287.05 |
---|
255 | ! Kappa |
---|
256 | REAL(r_std),PARAMETER :: kappa = cte_molr/cp_air |
---|
257 | ! in -- Kg/mole |
---|
258 | REAL(r_std),PARAMETER :: msmlr_air = 28.964E-03 |
---|
259 | ! in -- Kg/mole |
---|
260 | REAL(r_std),PARAMETER :: msmlr_h2o = 18.02E-03 |
---|
261 | !- |
---|
262 | REAL(r_std),PARAMETER :: cp_h2o = & |
---|
263 | & cp_air*(quatre*msmlr_air)/( 3.5_r_std*msmlr_h2o) |
---|
264 | !- |
---|
265 | REAL(r_std),PARAMETER :: cte_molr_h2o = cte_molr/quatre |
---|
266 | !- |
---|
267 | REAL(r_std),PARAMETER :: retv = msmlr_air/msmlr_h2o-un |
---|
268 | !- |
---|
269 | REAL(r_std),PARAMETER :: rvtmp2 = cp_h2o/cp_air-un |
---|
270 | !- |
---|
271 | REAL(r_std),PARAMETER :: cepdu2 = (0.1_r_std) **2 |
---|
272 | !- |
---|
273 | ! Van Karmann Constante |
---|
274 | REAL(r_std),PARAMETER :: ct_karman = 0.35_r_std |
---|
275 | !- |
---|
276 | ! g acceleration |
---|
277 | REAL(r_std),PARAMETER :: cte_grav = 9.80665_r_std |
---|
278 | !- |
---|
279 | ! Transform pascal into hectopascal |
---|
280 | REAL(r_std),PARAMETER :: pa_par_hpa = 100._r_std |
---|
281 | |
---|
282 | !------------------------------------- |
---|
283 | ! 2.1. Climatic constantes |
---|
284 | !------------------------------------- |
---|
285 | ! |
---|
286 | !$$ To externalise or not ? |
---|
287 | ! |
---|
288 | ! Constantes of the Louis scheme |
---|
289 | REAL(r_std),PARAMETER :: cb = cinq |
---|
290 | REAL(r_std),PARAMETER :: cc = cinq |
---|
291 | REAL(r_std),PARAMETER :: cd = cinq |
---|
292 | !- |
---|
293 | ! Constant in the computation of surface resistance |
---|
294 | REAL(r_std),PARAMETER :: rayt_cste = 125. |
---|
295 | !- |
---|
296 | ! DS :both used in diffuco.f90 |
---|
297 | ! Constant in the computation of surface resistance |
---|
298 | REAL(r_std),PARAMETER :: defc_plus=23.E-3 |
---|
299 | ! Constant in the computation of surface resistance |
---|
300 | REAL(r_std),PARAMETER :: defc_mult=1.5 |
---|
301 | |
---|
302 | !----------------------------------------- |
---|
303 | ! 2.2 Soil thermodynamics constants |
---|
304 | !----------------------------------------- |
---|
305 | ! |
---|
306 | ! Average Thermal Conductivity of soils |
---|
307 | REAL(r_std),PARAMETER :: so_cond = 1.5396 |
---|
308 | ! Average Heat capacity of soils |
---|
309 | REAL(r_std),PARAMETER :: so_capa = 2.0514e+6 |
---|
310 | !- |
---|
311 | ! Values taken from : PIELKE,'MESOSCALE METEOROLOGICAL MODELING',P.384 |
---|
312 | ! Dry soil heat capacity was decreased and conductivity increased. |
---|
313 | !- |
---|
314 | ! To externalise ? |
---|
315 | ! Dry soil Heat capacity of soils |
---|
316 | !*REAL(r_std),PARAMETER :: so_capa_dry = 1.35e+6 |
---|
317 | REAL(r_std),PARAMETER :: so_capa_dry = 1.80e+6 |
---|
318 | ! Dry soil Thermal Conductivity of soils |
---|
319 | !*REAL(r_std),PARAMETER :: so_cond_dry = 0.28 |
---|
320 | REAL(r_std),PARAMETER :: so_cond_dry = 0.40 |
---|
321 | !- |
---|
322 | ! Wet soil Heat capacity of soils |
---|
323 | REAL(r_std),PARAMETER :: so_capa_wet = 3.03e+6 |
---|
324 | ! Wet soil Thermal Conductivity of soils |
---|
325 | REAL(r_std),PARAMETER :: so_cond_wet = 1.89 |
---|
326 | !- |
---|
327 | ! Thermal Conductivity of snow |
---|
328 | REAL(r_std),PARAMETER :: sn_cond = 0.3 |
---|
329 | ! Snow density for the soil thermodynamics |
---|
330 | REAL(r_std),PARAMETER :: sn_dens = 330.0 |
---|
331 | ! Heat capacity for snow |
---|
332 | REAL(r_std),PARAMETER :: sn_capa = 2100.0_r_std*sn_dens |
---|
333 | |
---|
334 | |
---|
335 | !----------------------------------------------- |
---|
336 | !---------------------------------------------- |
---|
337 | ! SCALAR PARAMETERS EXTERNALIZED |
---|
338 | !---------------------------------------------- |
---|
339 | !----------------------------------------------- |
---|
340 | !------------------------------------------ |
---|
341 | ! SECHIBA, SOIL AND VEGETATION parameters |
---|
342 | !----------------------------------------- |
---|
343 | |
---|
344 | !!--------------------------------------- |
---|
345 | !! Parameters for soil type distribution |
---|
346 | !!--------------------------------------- |
---|
347 | ! |
---|
348 | ! Default soil texture distribution in the following order : |
---|
349 | ! sand, loam and clay |
---|
350 | REAL(r_std),SAVE, DIMENSION(nstm) :: soiltype_default = (/ 0.0, 1.0, 0.0 /) |
---|
351 | |
---|
352 | !!---------------------------------------- |
---|
353 | !! Constantes from the Choisnel hydrology |
---|
354 | !!---------------------------------------- |
---|
355 | ! |
---|
356 | ! Wilting point (Has a numerical role for the moment) |
---|
357 | REAL(r_std),SAVE :: qwilt = 5.0 |
---|
358 | ! Total depth of soil reservoir (for hydrolc) |
---|
359 | REAL(r_std),SAVE :: dpu_cste = deux |
---|
360 | ! The minimal size we allow for the upper reservoir (m) |
---|
361 | REAL(r_std),SAVE :: min_resdis = 2.e-5 |
---|
362 | !- |
---|
363 | ! Diffusion constant for the slow regime |
---|
364 | ! (This is for the diffusion between reservoirs) |
---|
365 | REAL(r_std),SAVE :: min_drain = 0.001 |
---|
366 | ! Diffusion constant for the fast regime |
---|
367 | REAL(r_std),SAVE :: max_drain = 0.1 |
---|
368 | ! The exponential in the diffusion law |
---|
369 | REAL(r_std),SAVE :: exp_drain = 1.5 |
---|
370 | !- |
---|
371 | ! Transforms leaf area index into size of interception reservoir |
---|
372 | REAL(r_std),SAVE :: qsintcst = 0.1 |
---|
373 | ! Maximum quantity of water (Kg/M3) |
---|
374 | REAL(r_std),SAVE :: mx_eau_eau = 150. |
---|
375 | !- |
---|
376 | ! Constant in the computation of resistance for bare soil evaporation |
---|
377 | REAL(r_std),SAVE :: rsol_cste = 33.E3 |
---|
378 | ! Scaling depth for litter humidity (m) |
---|
379 | !SZ changed this according to SP from 0.03 to 0.08, 080806 |
---|
380 | REAL(r_std),SAVE :: hcrit_litter=0.08_r_std |
---|
381 | |
---|
382 | !!--------------------------------------------------- |
---|
383 | !! Specific parameters for the CWRR hydrology module |
---|
384 | !!--------------------------------------------------- |
---|
385 | ! |
---|
386 | !!$ DS To externalise ? |
---|
387 | !!$ advice of MM : to put in hydrol |
---|
388 | ! CWRR linearisation |
---|
389 | INTEGER(i_std),PARAMETER :: imin = 1 |
---|
390 | ! number of interval for CWRR |
---|
391 | INTEGER(i_std),PARAMETER :: nbint = 100 |
---|
392 | ! number of points for CWRR |
---|
393 | INTEGER(i_std),PARAMETER :: imax = nbint+1 |
---|
394 | !- |
---|
395 | ! externalise w_time (some bug in hydrol) |
---|
396 | ! Time weighting for discretisation |
---|
397 | REAL(r_std),SAVE :: w_time = un |
---|
398 | !- |
---|
399 | ! Van genuchten coefficient n |
---|
400 | REAL(r_std),SAVE,DIMENSION(nstm) :: nvan = (/ 1.89_r_std, 1.56_r_std, 1.31_r_std /) |
---|
401 | ! Van genuchten coefficient a (mm^{-1}) |
---|
402 | REAL(r_std),SAVE,DIMENSION(nstm) :: avan = (/ 0.0075_r_std, 0.0036_r_std, 0.0019_r_std /) |
---|
403 | !- |
---|
404 | ! Residual soil water content |
---|
405 | REAL(r_std),SAVE,DIMENSION(nstm) :: mcr = (/ 0.065_r_std, 0.078_r_std, 0.095_r_std /) |
---|
406 | ! Saturated soil water content |
---|
407 | REAL(r_std),SAVE,DIMENSION(nstm) :: mcs = (/ 0.41_r_std, 0.43_r_std, 0.41_r_std /) |
---|
408 | ! Total depth of soil reservoir (m) |
---|
409 | REAL(r_std),SAVE,DIMENSION(nstm) :: dpu = (/ 2.0_r_std, 2.0_r_std, 2.0_r_std /) |
---|
410 | !- |
---|
411 | ! dpu must be constant over the different soil types |
---|
412 | ! Hydraulic conductivity Saturation (mm/d) |
---|
413 | REAL(r_std),SAVE,DIMENSION(nstm) :: ks = (/ 1060.8_r_std, 249.6_r_std, 62.4_r_std /) |
---|
414 | ! Soil moisture above which transpir is max |
---|
415 | REAL(r_std),SAVE,DIMENSION(nstm) :: pcent = (/ 0.5_r_std, 0.5_r_std, 0.5_r_std /) |
---|
416 | ! Max value of the permeability coeff at the bottom of the soil |
---|
417 | REAL(r_std),SAVE,DIMENSION(nstm) :: free_drain_max = (/ 1.0_r_std, 1.0_r_std, 1.0_r_std /) |
---|
418 | !- |
---|
419 | ! Volumetric water content field capacity |
---|
420 | REAL(r_std),SAVE,DIMENSION(nstm) :: mcf = (/ 0.32_r_std, 0.32_r_std, 0.32_r_std /) |
---|
421 | ! Volumetric water content Wilting pt |
---|
422 | REAL(r_std),SAVE,DIMENSION(nstm) :: mcw = (/ 0.10_r_std, 0.10_r_std, 0.10_r_std /) |
---|
423 | ! Vol. wat. cont. above which albedo is cst |
---|
424 | REAL(r_std),SAVE,DIMENSION(nstm) :: mc_awet = (/ 0.25_r_std, 0.25_r_std, 0.25_r_std /) |
---|
425 | ! Vol. wat. cont. below which albedo is cst |
---|
426 | REAL(r_std),SAVE,DIMENSION(nstm) :: mc_adry = (/ 0.1_r_std, 0.1_r_std, 0.1_r_std /) |
---|
427 | |
---|
428 | |
---|
429 | !!----------------------------------------------------- |
---|
430 | !! Vegetation parameters (previously in constantes_veg) |
---|
431 | !!----------------------------------------------------- |
---|
432 | ! |
---|
433 | ! Value for frac_nobio for tests in 0-dim simulations |
---|
434 | ! laisser ca tant qu'il n'y a que de la glace (pas de lacs) |
---|
435 | !DS : used in slowproc |
---|
436 | REAL(r_std),SAVE :: frac_nobio_fixed_test_1 = 0.0 |
---|
437 | !- |
---|
438 | ! Is veget_ori array stored in restart file |
---|
439 | !!$ DS: Where is it used ? |
---|
440 | ! LOGICAL,PARAMETER :: ldveget_ori_on_restart = .TRUE. |
---|
441 | !- |
---|
442 | ! Set to .TRUE. if you want q_cdrag coming from GCM |
---|
443 | ! used in diffuco |
---|
444 | LOGICAL,SAVE :: ldq_cdrag_from_gcm = .FALSE. |
---|
445 | !- |
---|
446 | ! allow agricultural PFTs |
---|
447 | LOGICAL,SAVE :: agriculture = .TRUE. |
---|
448 | !- |
---|
449 | ! The maximum mass (kg/m^2) of a glacier. |
---|
450 | REAL(r_std),SAVE :: maxmass_glacier = 3000. |
---|
451 | !- |
---|
452 | ! Minimal fraction of mesh a vegetation type can occupy |
---|
453 | REAL(r_std),SAVE :: min_vegfrac=0.001 |
---|
454 | !- |
---|
455 | !!$ DS not used in the code ? |
---|
456 | ! Limit of air temperature for snow |
---|
457 | REAL(r_std),SAVE :: tsnow=273. |
---|
458 | !- |
---|
459 | ! Sets the amount above which only sublimation occures [Kg/m^2] |
---|
460 | REAL(r_std),SAVE :: snowcri=1.5 |
---|
461 | ! Critical value for computation of snow albedo [Kg/m^2] |
---|
462 | REAL(r_std),SAVE :: snowcri_alb=10. |
---|
463 | ! Lower limit of snow amount |
---|
464 | REAL(r_std),SAVE :: sneige |
---|
465 | !- |
---|
466 | ! The minimum wind |
---|
467 | REAL(r_std),SAVE :: min_wind = 0.1 |
---|
468 | ! bare soil roughness length (m) |
---|
469 | REAL(r_std),SAVE :: z0_bare = 0.01 |
---|
470 | ! ice roughness length (m) |
---|
471 | REAL(r_std),SAVE :: z0_ice = 0.001 |
---|
472 | !- |
---|
473 | ! Time constant of the albedo decay of snow |
---|
474 | REAL(r_std),SAVE :: tcst_snowa = cinq |
---|
475 | ! Maximum period of snow aging |
---|
476 | REAL(r_std),SAVE :: max_snow_age = 50._r_std |
---|
477 | ! Transformation time constant for snow (m) |
---|
478 | REAL(r_std),SAVE :: snow_trans = 0.3_r_std |
---|
479 | !- |
---|
480 | ! albedo of dead leaves, VIS+NIR |
---|
481 | REAL(r_std),DIMENSION(2),SAVE :: alb_deadleaf = (/ .12, .35/) |
---|
482 | ! albedo of ice, VIS+NIR |
---|
483 | REAL(r_std),DIMENSION(2),SAVE :: alb_ice = (/ .60, .20/) |
---|
484 | |
---|
485 | !!-------------------------------- |
---|
486 | !! SECHIBA specific parameters |
---|
487 | !!-------------------------------- |
---|
488 | ! |
---|
489 | !- |
---|
490 | ! condveg |
---|
491 | !- |
---|
492 | ! to get z0 from height |
---|
493 | REAL(r_std), SAVE :: z0_over_height = un/16. |
---|
494 | ! Magic number which relates the height to the displacement height. |
---|
495 | REAL(r_std), SAVE :: height_displacement = 0.75 |
---|
496 | !- |
---|
497 | ! diffuco |
---|
498 | !- |
---|
499 | INTEGER(i_std), SAVE :: nlai = 20 ! dimension de tableau |
---|
500 | ! used in diffuco_trans |
---|
501 | REAL(r_std), SAVE :: laimax = 12. |
---|
502 | REAL(r_std), SAVE :: xc4_1 = .83 |
---|
503 | REAL(r_std), SAVE :: xc4_2 = .93 |
---|
504 | !- |
---|
505 | ! hydrol. |
---|
506 | !- |
---|
507 | ! Allowed moisture above mcs (boundary conditions) |
---|
508 | REAL(r_std), SAVE :: dmcs = 0.002 |
---|
509 | ! Allowed moisture below mcr (boundary conditions) |
---|
510 | REAL(r_std), SAVE :: dmcr = 0.002 |
---|
511 | !- |
---|
512 | ! routing |
---|
513 | !- |
---|
514 | ! Parameter for the Kassel irrigation parametrization linked to the crops |
---|
515 | REAL(r_std), SAVE :: crop_coef = 1.5 |
---|
516 | !- |
---|
517 | ! slowproc |
---|
518 | !- |
---|
519 | REAL(r_std), SAVE :: clayfraction_default = 0.2 |
---|
520 | |
---|
521 | !----------------------------- |
---|
522 | ! STOMATE AND LPJ PARAMETERS |
---|
523 | !----------------------------- |
---|
524 | ! |
---|
525 | !- |
---|
526 | ! lpj_constraints |
---|
527 | !- |
---|
528 | ! longest sustainable time without regeneration (vernalization) |
---|
529 | REAL(r_std), SAVE :: too_long = 5. |
---|
530 | ! |
---|
531 | !- |
---|
532 | ! lpj_fire |
---|
533 | !- |
---|
534 | ! Time scale for memory of the fire index (days). Validated for one year in the DGVM. |
---|
535 | REAL(r_std), SAVE :: tau_fire = 30. |
---|
536 | ! Critical litter quantity for fire |
---|
537 | REAL(r_std), SAVE :: litter_crit = 200. |
---|
538 | ! |
---|
539 | !- |
---|
540 | ! lpj_light |
---|
541 | !- |
---|
542 | ! maximum total number of grass individuals in a closed canopy |
---|
543 | REAL(r_std), SAVE :: grass_mercy = 0.01 |
---|
544 | ! minimum fraction of trees that survive even in a closed canopy |
---|
545 | REAL(r_std), SAVE :: tree_mercy = 0.01 |
---|
546 | ! for diagnosis of fpc increase, compare today's fpc to last year's maximum (T) or |
---|
547 | ! to fpc of last time step (F)? |
---|
548 | LOGICAL, SAVE :: annual_increase = .TRUE. |
---|
549 | ! |
---|
550 | !- |
---|
551 | ! lpj_pftinout |
---|
552 | !- |
---|
553 | ! minimum availability |
---|
554 | REAL(r_std), SAVE :: min_avail = 0.01 |
---|
555 | ! |
---|
556 | !- |
---|
557 | ! stomate_alloc |
---|
558 | !- |
---|
559 | ! Do we try to reach a minimum reservoir even if we are severely stressed? |
---|
560 | LOGICAL, SAVE :: ok_minres = .TRUE. |
---|
561 | ! time (d) to attain the initial foliage using the carbohydrate reserve |
---|
562 | REAL(r_std), SAVE :: tau_leafinit = 10. |
---|
563 | ! maximum time (d) during which reserve is used (trees) |
---|
564 | REAL(r_std), SAVE :: reserve_time_tree = 30. |
---|
565 | ! maximum time (d) during which reserve is used (grasses) |
---|
566 | REAL(r_std), SAVE :: reserve_time_grass = 20. |
---|
567 | ! Standard root allocation |
---|
568 | REAL(r_std), SAVE :: R0 = 0.3 |
---|
569 | ! Standard sapwood allocation |
---|
570 | REAL(r_std), SAVE :: S0 = 0.3 |
---|
571 | ! only used in stomate_alloc |
---|
572 | ! Standard leaf allocation |
---|
573 | REAL(r_std), SAVE :: L0 |
---|
574 | ! Standard fruit allocation |
---|
575 | REAL(r_std), SAVE :: f_fruit = 0.1 |
---|
576 | ! fraction of sapwood allocation above ground (SHOULD BE CALCULATED !!!!) |
---|
577 | REAL(r_std), SAVE :: alloc_sap_above_tree = 0.5 |
---|
578 | REAL(r_std), SAVE :: alloc_sap_above_grass = 1.0 |
---|
579 | ! extrema of leaf allocation fraction |
---|
580 | REAL(r_std), SAVE :: min_LtoLSR = 0.2 |
---|
581 | REAL(r_std), SAVE :: max_LtoLSR = 0.5 |
---|
582 | ! scaling depth for nitrogen limitation (m) |
---|
583 | REAL(r_std), SAVE :: z_nitrogen = 0.2 |
---|
584 | ! |
---|
585 | !- |
---|
586 | ! stomate_data |
---|
587 | !- |
---|
588 | !!------------------------------- |
---|
589 | !! Parameters for the pipe model |
---|
590 | !!------------------------------ |
---|
591 | !- |
---|
592 | ! crown area = pipe_tune1. stem diameter**(1.6) (Reinicke's theory) |
---|
593 | REAL(r_std),SAVE :: pipe_tune1 = 100.0 |
---|
594 | ! height=pipe_tune2 * diameter**pipe_tune3 |
---|
595 | REAL(r_std),SAVE :: pipe_tune2 = 40.0 |
---|
596 | REAL(r_std),SAVE :: pipe_tune3 = 0.5 |
---|
597 | ! needed for stem diameter |
---|
598 | REAL(r_std),SAVE :: pipe_tune4 = 0.3 |
---|
599 | ! Density |
---|
600 | REAL(r_std),SAVE :: pipe_density = 2.e5 |
---|
601 | ! one more SAVE |
---|
602 | REAL(r_std),SAVE :: pipe_k1 = 8.e3 |
---|
603 | ! |
---|
604 | !- |
---|
605 | ! Maximum tree establishment rate |
---|
606 | REAL(r_std),SAVE :: estab_max_tree = 0.12 |
---|
607 | ! Maximum grass establishment rate |
---|
608 | REAL(r_std),SAVE :: estab_max_grass = 0.12 |
---|
609 | ! initial density of individuals |
---|
610 | REAL(r_std),SAVE :: ind_0 = 0.02 |
---|
611 | ! For trees, minimum fraction of crown area occupied |
---|
612 | ! (due to its branches etc.) |
---|
613 | ! This means that only a small fraction of its crown area |
---|
614 | ! can be invaded by other trees. |
---|
615 | REAL(r_std),SAVE :: min_cover = 0.05 |
---|
616 | !- |
---|
617 | ! alpha's : ? |
---|
618 | REAL(r_std),SAVE :: alpha_grass = .5 |
---|
619 | REAL(r_std),SAVE :: alpha_tree = 1. |
---|
620 | !- |
---|
621 | ! maximum reference long term temperature (K) |
---|
622 | REAL(r_std),SAVE :: tlong_ref_max = 303.1 |
---|
623 | ! minimum reference long term temperature (K) |
---|
624 | REAL(r_std),SAVE :: tlong_ref_min = 253.1 |
---|
625 | ! |
---|
626 | !! LOGICAL |
---|
627 | !- |
---|
628 | ! Do we treat PFT expansion across a grid point after introduction? |
---|
629 | ! default = .FALSE. |
---|
630 | LOGICAL,SAVE :: treat_expansion = .FALSE. |
---|
631 | ! |
---|
632 | ! herbivores? |
---|
633 | LOGICAL,SAVE :: ok_herbivores = .FALSE. |
---|
634 | ! |
---|
635 | ! harvesting ? |
---|
636 | LOGICAL,SAVE :: harvest_agri = .TRUE. |
---|
637 | !!---------------------- |
---|
638 | !! climatic parameters |
---|
639 | !!--------------------- |
---|
640 | ! |
---|
641 | ! minimum precip, in mm/year |
---|
642 | REAL(r_std),SAVE :: precip_crit = 100. |
---|
643 | ! minimum gdd for establishment of saplings |
---|
644 | REAL(r_std),SAVE :: gdd_crit_estab = 150. |
---|
645 | ! critical fpc, needed for light competition and establishment |
---|
646 | REAL(r_std),SAVE :: fpc_crit = 0.95 |
---|
647 | !- |
---|
648 | ! fraction of GPP which is lost as growth respiration |
---|
649 | REAL(r_std),SAVE :: frac_growthresp = 0.28 |
---|
650 | ! |
---|
651 | !- |
---|
652 | ! mass ratio (heartwood+sapwood)/sapwood |
---|
653 | REAL(r_std), SAVE :: mass_ratio_heart_sap = 3. |
---|
654 | ! |
---|
655 | !!--------------------------------------------------------- |
---|
656 | ! time scales for phenology and other processes (in days) |
---|
657 | !!--------------------------------------------------------- |
---|
658 | ! |
---|
659 | REAL(r_std), SAVE :: tau_hum_month = 20. |
---|
660 | REAL(r_std), SAVE :: tau_hum_week = 7. |
---|
661 | REAL(r_std), SAVE :: tau_t2m_month = 20. |
---|
662 | REAL(r_std), SAVE :: tau_t2m_week = 7. |
---|
663 | REAL(r_std), SAVE :: tau_tsoil_month = 20. |
---|
664 | REAL(r_std), SAVE :: tau_soilhum_month = 20. |
---|
665 | REAL(r_std), SAVE :: tau_gpp_week = 7. |
---|
666 | REAL(r_std), SAVE :: tau_gdd = 40. |
---|
667 | REAL(r_std), SAVE :: tau_ngd = 50. |
---|
668 | REAL(r_std), SAVE :: coeff_tau_longterm = 3. |
---|
669 | ! used in stomate_data and in stomate_season |
---|
670 | REAL(r_std), SAVE :: tau_longterm |
---|
671 | ! |
---|
672 | !- |
---|
673 | ! stomate_litter |
---|
674 | !- |
---|
675 | ! scaling depth for soil activity (m) |
---|
676 | REAL(r_std), SAVE :: z_decomp = 0.2 |
---|
677 | ! |
---|
678 | !- |
---|
679 | ! stomate_lpj |
---|
680 | !- |
---|
681 | REAL(r_std), SAVE :: frac_turnover_daily = 0.55 |
---|
682 | ! |
---|
683 | !- |
---|
684 | ! stomate_npp |
---|
685 | !- |
---|
686 | ! maximum fraction of allocatable biomass used for maintenance respiration |
---|
687 | REAL(r_std), SAVE :: tax_max = 0.8 |
---|
688 | ! |
---|
689 | !- |
---|
690 | ! stomate_phenology |
---|
691 | !- |
---|
692 | ! take carbon from atmosphere if carbohydrate reserve too small? |
---|
693 | LOGICAL, SAVE :: always_init = .FALSE. |
---|
694 | ! minimum time (d) since last beginning of a growing season |
---|
695 | REAL(r_std), SAVE :: min_growthinit_time = 300. |
---|
696 | ! moisture availability above which moisture tendency doesn't matter |
---|
697 | REAL(r_std), SAVE :: moiavail_always_tree = 1.0 |
---|
698 | REAL(r_std), SAVE :: moiavail_always_grass = 0.6 |
---|
699 | ! monthly temp. above which temp. tendency doesn't matter |
---|
700 | REAL(r_std), SAVE :: t_always |
---|
701 | REAL(r_std), SAVE :: t_always_add = 10. |
---|
702 | ! |
---|
703 | !- |
---|
704 | ! stomate_season |
---|
705 | !- |
---|
706 | ! rapport maximal GPP/GGP_max pour dormance |
---|
707 | REAL(r_std), SAVE :: gppfrac_dormance = 0.2 |
---|
708 | ! minimum gpp considered as not "lowgpp" |
---|
709 | REAL(r_std), SAVE :: min_gpp_allowed = 0.3 |
---|
710 | ! tau (year) for "climatologic variables |
---|
711 | REAL(r_std), SAVE :: tau_climatology = 20 |
---|
712 | ! parameters for herbivore activity |
---|
713 | REAL(r_std), SAVE :: hvc1 = 0.019 |
---|
714 | REAL(r_std), SAVE :: hvc2 = 1.38 |
---|
715 | REAL(r_std), SAVE :: leaf_frac_hvc =.33 |
---|
716 | ! |
---|
717 | !- |
---|
718 | ! stomate_vmax |
---|
719 | !- |
---|
720 | ! offset (minimum relative vcmax) |
---|
721 | REAL(r_std), SAVE :: vmax_offset = 0.3 |
---|
722 | ! leaf age at which vmax attains vcmax_opt (in fraction of critical leaf age) |
---|
723 | REAL(r_std), SAVE :: leafage_firstmax = 0.03 |
---|
724 | ! leaf age at which vmax falls below vcmax_opt (in fraction of critical leaf age) |
---|
725 | REAL(r_std), SAVE :: leafage_lastmax = 0.5 |
---|
726 | ! leaf age at which vmax attains its minimum (in fraction of critical leaf age) |
---|
727 | REAL(r_std), SAVE :: leafage_old = 1. |
---|
728 | |
---|
729 | |
---|
730 | !-------------------------- |
---|
731 | !-------------------------- |
---|
732 | ! ARRAYS-PARAMETERS |
---|
733 | !-------------------------- |
---|
734 | !-------------------------- |
---|
735 | !- |
---|
736 | ! condveg |
---|
737 | !- |
---|
738 | ! The correspondance table for the soil color numbers and their albedo |
---|
739 | ! |
---|
740 | REAL(r_std), DIMENSION(classnb) :: vis_dry = (/0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.27/) |
---|
741 | REAL(r_std), DIMENSION(classnb) :: nir_dry = (/0.48, 0.44, 0.40, 0.36, 0.32, 0.28, 0.24, 0.20, 0.55/) |
---|
742 | REAL(r_std), DIMENSION(classnb) :: vis_wet = (/0.12, 0.11, 0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.15/) |
---|
743 | REAL(r_std), DIMENSION(classnb) :: nir_wet = (/0.24, 0.22, 0.20, 0.18, 0.16, 0.14, 0.12, 0.10, 0.31/) |
---|
744 | ! |
---|
745 | ! Nathalie, introduction d'un albedo moyen, VIS+NIR |
---|
746 | ! Les valeurs suivantes correspondent a la moyenne des valeurs initiales |
---|
747 | ! REAL(stnd), DIMENSION(classnb) :: albsoil_vis = (/0.18, 0.165, 0.15, 0.135, 0.12, 0.105, 0.09, 0.075, 0.21/) |
---|
748 | ! REAL(stnd), DIMENSION(classnb) :: albsoil_nir = (/0.36, 0.33, 0.30, 0.27, 0.24, 0.21, 0.18, 0.15, 0.43/) |
---|
749 | ! les valeurs retenues accentuent le contraste entre equateur et Sahara. |
---|
750 | ! On diminue aussi l'albedo des deserts (tous sauf Sahara) |
---|
751 | REAL(r_std), DIMENSION(classnb) :: albsoil_vis = (/0.18, 0.16, 0.16, 0.15, 0.12, 0.105, 0.09, 0.075, 0.25/) |
---|
752 | REAL(r_std), DIMENSION(classnb) :: albsoil_nir = (/0.36, 0.34, 0.34, 0.33, 0.30, 0.25, 0.20, 0.15, 0.45/) |
---|
753 | |
---|
754 | !- |
---|
755 | ! lpj_fire |
---|
756 | !- |
---|
757 | |
---|
758 | ! What fraction of a burned plant compartment goes into the atmosphere |
---|
759 | ! (rest into litter) |
---|
760 | REAL(r_std), SAVE, DIMENSION(nparts) :: co2frac = (/ .95, .95, 0., 0.3, 0., 0., .95, .95 /) |
---|
761 | |
---|
762 | !- |
---|
763 | ! stomate_litter |
---|
764 | !- |
---|
765 | |
---|
766 | ! C/N ratio |
---|
767 | REAL(r_std), SAVE, DIMENSION(nparts) :: CN = 40.0 |
---|
768 | ! Lignine/C ratio of the different plant parts |
---|
769 | REAL(r_std), SAVE, DIMENSION(nparts) :: LC = (/ 0.22, 0.35, 0.35, 0.35, 0.35, 0.22, 0.22, 0.22 /) |
---|
770 | ! corresponding to frac_soil(istructural,iactive,iabove) |
---|
771 | REAL(r_std), SAVE :: frac_soil_struct_aa = .55 |
---|
772 | ! corresponding to frac_soil(istructural,iactive,ibelow) |
---|
773 | REAL(r_std), SAVE :: frac_soil_struct_ab = .45 |
---|
774 | ! corresponding to frac_soil(istructural,islow,iabove) |
---|
775 | REAL(r_std), SAVE :: frac_soil_struct_sa = .7 |
---|
776 | ! corresponding to frac_soil(istructural,islow,ibelow) |
---|
777 | REAL(r_std), SAVE :: frac_soil_struct_sb = .7 |
---|
778 | ! corresponding to frac_soil(imetabolic,iactive,iabove) |
---|
779 | REAL(r_std), SAVE :: frac_soil_metab_aa = .45 |
---|
780 | ! corresponding to frac_soil(imetabolic,iactive,ibelow) |
---|
781 | REAL(r_std), SAVE :: frac_soil_metab_ab = .45 |
---|
782 | !- |
---|
783 | ! stomate_soilcarbon |
---|
784 | !- |
---|
785 | ! frac_carb_coefficients |
---|
786 | ! from active pool: depends on clay content |
---|
787 | ! correspnding to frac_carb(:,iactive,iactive) |
---|
788 | REAL(r_std), SAVE :: frac_carb_aa = 0.0 |
---|
789 | ! correspnding to frac_carb(:,iactive,ipassive) |
---|
790 | REAL(r_std), SAVE :: frac_carb_ap = 0.004 |
---|
791 | !frac_carb(;;iactive,islow) is computed in stomate_soilcarbon.f90 |
---|
792 | !- |
---|
793 | ! from slow pool |
---|
794 | ! correspnding to frac_carb(:,islow,islow) |
---|
795 | REAL(r_std), SAVE :: frac_carb_ss = 0.0 |
---|
796 | ! correspnding to frac_carb(:,islow,iactive) |
---|
797 | REAL(r_std), SAVE :: frac_carb_sa = .42 |
---|
798 | ! correspnding to frac_carb(:,islow,ipassive) |
---|
799 | REAL(r_std), SAVE :: frac_carb_sp = .03 |
---|
800 | !- |
---|
801 | ! from passive pool |
---|
802 | ! correspnding to frac_carb(:,ipassive,ipassive) |
---|
803 | REAL(r_std), SAVE :: frac_carb_pp = .0 |
---|
804 | ! correspnding to frac_carb(:,ipassive,iactive) |
---|
805 | REAL(r_std), SAVE :: frac_carb_pa = .45 |
---|
806 | ! correspnding to frac_carb(:,ipassive,islow) |
---|
807 | REAL(r_std), SAVE :: frac_carb_ps = .0 |
---|
808 | |
---|
809 | |
---|
810 | !---------------------------------------- |
---|
811 | !--------------------------------------- |
---|
812 | ! COEFFICIENTS OF EQUATIONS |
---|
813 | !------------------------------------- |
---|
814 | !--------------------------------------- |
---|
815 | |
---|
816 | !--------- |
---|
817 | ! SECHIBA |
---|
818 | !--------- |
---|
819 | !- |
---|
820 | ! diffuco |
---|
821 | !- |
---|
822 | REAL(r_std),PARAMETER :: Tetens_1 = 0.622 |
---|
823 | REAL(r_std),PARAMETER :: Tetens_2 = 0.378 |
---|
824 | REAL(r_std),PARAMETER :: std_ci_frac = 0.667 |
---|
825 | REAL(r_std),PARAMETER :: alpha_j = 0.8855 |
---|
826 | REAL(r_std),PARAMETER :: curve_assim = 0.7 |
---|
827 | REAL(r_std),PARAMETER :: WJ_coeff1 = 4.5 |
---|
828 | REAL(r_std),PARAMETER :: WJ_coeff2 = 10.5 |
---|
829 | REAL(r_std),PARAMETER :: Vc_to_Rd_ratio = 0.011 |
---|
830 | REAL(r_std),PARAMETER :: O2toCO2_stoechio = 1.6 |
---|
831 | REAL(r_std),PARAMETER :: mmol_to_m_1 = 0.0244 |
---|
832 | REAL(r_std),PARAMETER :: RG_to_PAR = 0.5 |
---|
833 | REAL(r_std),PARAMETER :: W_to_mmol = 4.6 ! W_to_mmol * RG_to_PAR = 2.3 |
---|
834 | ! |
---|
835 | REAL(r_std), SAVE :: lai_level_depth = .15 |
---|
836 | REAL(r_std), SAVE :: x1_coef = 0.177 |
---|
837 | REAL(r_std), SAVE :: x1_Q10 = 0.069 |
---|
838 | REAL(r_std), SAVE :: quantum_yield = 0.092 |
---|
839 | REAL(r_std), SAVE :: kt_coef = 0.7 |
---|
840 | REAL(r_std), SAVE :: kc_coef = 39.09 |
---|
841 | REAL(r_std), SAVE :: Ko_Q10 = .085 |
---|
842 | REAL(r_std), SAVE :: Oa = 210000. |
---|
843 | REAL(r_std), SAVE :: Ko_coef = 2.412 |
---|
844 | REAL(r_std), SAVE :: CP_0 = 42. |
---|
845 | REAL(r_std), SAVE :: CP_temp_coef = 9.46 |
---|
846 | REAL(r_std), SAVE :: CP_temp_ref = 25. |
---|
847 | ! |
---|
848 | REAL(r_std), SAVE, DIMENSION(2) :: rt_coef = (/ 0.8, 1.3 /) |
---|
849 | REAL(r_std), SAVE, DIMENSION(2) :: vc_coef = (/ 0.39, 0.3 /) |
---|
850 | ! |
---|
851 | ! coefficients of the polynome of degree 5 used inthe equation of coeff_dew_veg |
---|
852 | REAL(r_std), SAVE, DIMENSION(6) :: dew_veg_poly_coeff = & |
---|
853 | & (/ 0.887773, 0.205673, 0.110112, 0.014843, 0.000824, 0.000017 /) |
---|
854 | |
---|
855 | !--------- |
---|
856 | ! LPJ |
---|
857 | !--------- |
---|
858 | !- |
---|
859 | ! lpj_crown |
---|
860 | !- |
---|
861 | REAL(r_std), SAVE :: pipe_tune_exp_coeff = 1.6 |
---|
862 | ! |
---|
863 | !- |
---|
864 | ! lpj_establish |
---|
865 | !- |
---|
866 | REAL(r_std), SAVE :: establish_scal_fact = 15. |
---|
867 | REAL(r_std), SAVE :: fpc_crit_max = .075 |
---|
868 | REAL(r_std), SAVE :: fpc_crit_min= .05 |
---|
869 | ! |
---|
870 | !- |
---|
871 | ! lpj_fire |
---|
872 | !- |
---|
873 | REAL(r_std), SAVE, DIMENSION(3) :: bcfrac_coeff = (/ .3, 1.3, 88.2 /) |
---|
874 | REAL(r_std), SAVE, DIMENSION(4) :: firefrac_coeff = (/ 0.45, 0.8, 0.6, 0.13 /) |
---|
875 | ! |
---|
876 | !- |
---|
877 | ! lpj_gap |
---|
878 | !- |
---|
879 | REAL(r_std), SAVE :: availability_fact = 0.02 |
---|
880 | REAL(r_std), SAVE :: vigour_ref = 0.17 |
---|
881 | REAL(r_std), SAVE :: vigour_coeff = 70. |
---|
882 | !- |
---|
883 | ! lpj_pftinout |
---|
884 | !- |
---|
885 | REAL(r_std), SAVE :: RIP_time_min = 1.25 |
---|
886 | REAL(r_std), SAVE :: npp_longterm_init = 10. |
---|
887 | REAL(r_std), SAVE :: everywhere_init = 0.05 |
---|
888 | ! |
---|
889 | |
---|
890 | !--------- |
---|
891 | ! STOMATE |
---|
892 | !--------- |
---|
893 | !- |
---|
894 | ! stomate_alloc |
---|
895 | !- |
---|
896 | REAL(r_std), PARAMETER :: max_possible_lai = 10. |
---|
897 | REAL(r_std), PARAMETER :: Nlim_Q10 = 10. |
---|
898 | ! |
---|
899 | REAL(r_std), SAVE :: lai_max_to_happy = 0.5 |
---|
900 | REAL(r_std), SAVE :: Nlim_tref = 25. |
---|
901 | ! |
---|
902 | !- |
---|
903 | ! stomate_data |
---|
904 | !- |
---|
905 | REAL(r_std), SAVE :: bm_sapl_carbres = 5. |
---|
906 | REAL(r_std), SAVE :: bm_sapl_sapabove = 0.5 |
---|
907 | REAL(r_std), SAVE :: bm_sapl_heartabove = 2. |
---|
908 | REAL(r_std), SAVE :: bm_sapl_heartbelow = 2. |
---|
909 | REAL(r_std), SAVE :: init_sapl_mass_leaf_nat = 0.1 |
---|
910 | REAL(r_std), SAVE :: init_sapl_mass_leaf_agri = 1. |
---|
911 | REAL(r_std), SAVE :: init_sapl_mass_carbres = 5. |
---|
912 | REAL(r_std), SAVE :: init_sapl_mass_root = 0.1 |
---|
913 | REAL(r_std), SAVE :: init_sapl_mass_fruit = 0.3 |
---|
914 | REAL(r_std), SAVE :: cn_sapl_init = 0.5 |
---|
915 | REAL(r_std), SAVE :: migrate_tree = 10.*1.E3 |
---|
916 | REAL(r_std), SAVE :: migrate_grass = 10.*1.E3 |
---|
917 | REAL(r_std), SAVE :: lai_initmin_tree = 0.3 |
---|
918 | REAL(r_std), SAVE :: lai_initmin_grass = 0.1 |
---|
919 | REAL(r_std), SAVE, DIMENSION(2) :: dia_coeff = (/ 4., 0.5 /) |
---|
920 | REAL(r_std), SAVE, DIMENSION(2) :: maxdia_coeff =(/ 100., 0.01/) |
---|
921 | REAL(r_std), SAVE, DIMENSION(4) :: bm_sapl_leaf = (/ 4., 4., .8, 5./) |
---|
922 | ! |
---|
923 | !- |
---|
924 | ! stomate_litter |
---|
925 | !- |
---|
926 | REAL(r_std), PARAMETER :: Q10 = 10. |
---|
927 | ! |
---|
928 | REAL(r_std), SAVE :: metabolic_ref_frac = 0.85 |
---|
929 | REAL(r_std), SAVE :: metabolic_LN_ratio = 0.018 |
---|
930 | REAL(r_std), SAVE :: tau_metabolic = .066 |
---|
931 | REAL(r_std), SAVE :: tau_struct = .245 |
---|
932 | REAL(r_std), SAVE :: soil_Q10 = .69 != ln 2 |
---|
933 | REAL(r_std), SAVE :: tsoil_ref = 30. |
---|
934 | REAL(r_std), SAVE :: litter_struct_coef = 3. |
---|
935 | REAL(r_std), SAVE, DIMENSION(3) :: moist_coeff = (/ 1.1, 2.4, 0.29 /) |
---|
936 | ! |
---|
937 | !- |
---|
938 | ! stomate_phenology |
---|
939 | !- |
---|
940 | REAL(r_std), SAVE :: gddncd_ref = 603. |
---|
941 | REAL(r_std), SAVE :: gddncd_curve = 0.0091 |
---|
942 | REAL(r_std), SAVE :: gddncd_offset = 64. |
---|
943 | ! |
---|
944 | !- |
---|
945 | ! stomate_prescribe |
---|
946 | !- |
---|
947 | REAL(r_std), SAVE :: cn_tree = 4. |
---|
948 | REAL(r_std), SAVE :: bm_sapl_rescale = 40. |
---|
949 | ! |
---|
950 | !- |
---|
951 | ! stomate_resp |
---|
952 | !- |
---|
953 | REAL(r_std), SAVE :: maint_resp_min_vmax = 0.3 |
---|
954 | REAL(r_std), SAVE :: maint_resp_coeff = 1.4 |
---|
955 | ! |
---|
956 | !- |
---|
957 | ! stomate_season |
---|
958 | !- |
---|
959 | REAL(r_std), SAVE :: ncd_max_year = 3. |
---|
960 | REAL(r_std), SAVE :: gdd_threshold = 5. |
---|
961 | REAL(r_std), SAVE :: green_age_ever = 2. |
---|
962 | REAL(r_std), SAVE :: green_age_dec = 0.5 |
---|
963 | !- |
---|
964 | ! stomate_soilcarbon |
---|
965 | !- |
---|
966 | REAL(r_std), SAVE :: active_to_pass_clay_frac = .68 |
---|
967 | !residence times in carbon pools (days) |
---|
968 | REAL(r_std), SAVE :: carbon_tau_iactive = .149 |
---|
969 | REAL(r_std), SAVE :: carbon_tau_islow = 5.48 |
---|
970 | REAL(r_std), SAVE :: carbon_tau_ipassive = 241. |
---|
971 | ! |
---|
972 | REAL(r_std), SAVE, DIMENSION(3) :: flux_tot_coeff = (/ 1.2, 1.4, .75/) |
---|
973 | ! |
---|
974 | !- |
---|
975 | ! stomate_turnover |
---|
976 | !- |
---|
977 | REAL(r_std), SAVE :: new_turnover_time_ref = 20. |
---|
978 | REAL(r_std), SAVE :: dt_turnover_time = 10. |
---|
979 | REAL(r_std), SAVE :: leaf_age_crit_tref = 20. |
---|
980 | REAL(r_std), SAVE, DIMENSION(3) :: leaf_age_crit_coeff = (/ 1.5, 0.75, 10./) |
---|
981 | |
---|
982 | !************************************************************** |
---|
983 | |
---|
984 | CONTAINS |
---|
985 | |
---|
986 | ! Subroutine called for getin the new parameters values used in sechiba |
---|
987 | ! |
---|
988 | SUBROUTINE getin_sechiba_parameters |
---|
989 | |
---|
990 | IMPLICIT NONE |
---|
991 | ! first call |
---|
992 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
993 | |
---|
994 | IF(first_call) THEN |
---|
995 | |
---|
996 | !!$ CALL getin('DIAG_QSAT',diag_qsat) |
---|
997 | ! |
---|
998 | CALL getin('QWILT',qwilt) |
---|
999 | CALL getin('MIN_RESDIS',min_resdis) |
---|
1000 | CALL getin('MIN_DRAIN',min_drain) |
---|
1001 | CALL getin('MAX_DRAIN',max_drain) |
---|
1002 | CALL getin('EXP_DRAIN',exp_drain) |
---|
1003 | CALL getin('MX_EAU_EAU',mx_eau_eau) |
---|
1004 | CALL getin('RSOL_CSTE',rsol_cste) |
---|
1005 | CALL getin('HCRIT_LITTER',hcrit_litter) |
---|
1006 | !- |
---|
1007 | CALL getin('SOILTYPE_DEFAULT',soiltype_default) |
---|
1008 | !- |
---|
1009 | CALL getin('MAXMASS_GLACIER',maxmass_glacier) |
---|
1010 | CALL getin('MIN_VEGFRAC',min_vegfrac) |
---|
1011 | !- |
---|
1012 | CALL getin('SNOWCRI',snowcri) |
---|
1013 | !- |
---|
1014 | CALL getin('SNOWCRI_ALB',snowcri_alb) |
---|
1015 | CALL getin('MIN_WIND',min_wind) |
---|
1016 | CALL getin('Z0_BARE',z0_bare) |
---|
1017 | CALL getin('Z0_ICE',z0_ice) |
---|
1018 | CALL getin('TCST_SNOWA',tcst_snowa) |
---|
1019 | CALL getin('MAX_SNOW_AGE',max_snow_age) |
---|
1020 | CALL getin('SNOW_TRANS',snow_trans) |
---|
1021 | CALL getin('ALB_DEADLEAF',alb_deadleaf) |
---|
1022 | CALL getin('ALB_ICE',alb_ice) |
---|
1023 | !- |
---|
1024 | CALL getin('Z0_OVER_HEIGHT',z0_over_height) |
---|
1025 | CALL getin('HEIGHT_DISPLACEMENT',height_displacement) |
---|
1026 | !- |
---|
1027 | CALL getin('NLAI',nlai) |
---|
1028 | CALL getin('LAIMAX',laimax) |
---|
1029 | CALL getin('XC4_1',xc4_1) |
---|
1030 | CALL getin('XC4_2',xc4_2) |
---|
1031 | !- |
---|
1032 | CALL getin('DMCS',dmcs) |
---|
1033 | CALL getin('DMCR',dmcr) |
---|
1034 | !- |
---|
1035 | CALL getin('VIS_DRY',vis_dry) |
---|
1036 | CALL getin('NIR_DRY',nir_dry) |
---|
1037 | CALL getin('VIS_WET',vis_wet) |
---|
1038 | CALL getin('NIR_WET',nir_wet) |
---|
1039 | CALL getin('ALBSOIL_VIS',albsoil_vis) |
---|
1040 | CALL getin('ALBSOIL_NIR',albsoil_nir) |
---|
1041 | !- |
---|
1042 | CALL getin('CLAYFRACTION_DEFAULT',clayfraction_default) |
---|
1043 | ! |
---|
1044 | CALL getin('DEW_VEG_POLY_COEFF',dew_veg_poly_coeff) |
---|
1045 | |
---|
1046 | first_call =.FALSE. |
---|
1047 | |
---|
1048 | ENDIF |
---|
1049 | |
---|
1050 | END SUBROUTINE getin_sechiba_parameters |
---|
1051 | |
---|
1052 | !********************************************************* |
---|
1053 | |
---|
1054 | ! Subroutine called only if river_routing is activated |
---|
1055 | |
---|
1056 | SUBROUTINE getin_routing_parameters |
---|
1057 | |
---|
1058 | IMPLICIT NONE |
---|
1059 | |
---|
1060 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
1061 | |
---|
1062 | IF(first_call) THEN |
---|
1063 | |
---|
1064 | CALL getin('CROP_COEF',crop_coef) |
---|
1065 | |
---|
1066 | first_call =.FALSE. |
---|
1067 | |
---|
1068 | ENDIF |
---|
1069 | |
---|
1070 | END SUBROUTINE getin_routing_parameters |
---|
1071 | |
---|
1072 | !******************************************************* |
---|
1073 | |
---|
1074 | ! Subroutine called only if hydrol_cwrr is activated |
---|
1075 | |
---|
1076 | SUBROUTINE getin_hydrol_cwrr_parameters |
---|
1077 | |
---|
1078 | IMPLICIT NONE |
---|
1079 | |
---|
1080 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
1081 | |
---|
1082 | IF(first_call) THEN |
---|
1083 | |
---|
1084 | CALL getin('W_TIME',w_time) |
---|
1085 | CALL getin('NVAN',nvan) |
---|
1086 | CALL getin('AVAN',avan) |
---|
1087 | CALL getin('MCR',mcr) |
---|
1088 | CALL getin('MCS',mcs) |
---|
1089 | CALL getin('KS',ks) |
---|
1090 | CALL getin('PCENT',pcent) |
---|
1091 | CALL getin('FREE_DRAIN_MAX',free_drain_max) |
---|
1092 | CALL getin('MCF',mcf) |
---|
1093 | CALL getin('MCW',mcw) |
---|
1094 | CALL getin('MC_AWET',mc_awet) |
---|
1095 | |
---|
1096 | first_call =.FALSE. |
---|
1097 | |
---|
1098 | ENDIF |
---|
1099 | |
---|
1100 | END SUBROUTINE getin_hydrol_cwrr_parameters |
---|
1101 | !-------------------------------------------- |
---|
1102 | |
---|
1103 | ! Subroutine called only if ok_co2 is activated |
---|
1104 | ! only for diffuco_trans_co2 |
---|
1105 | |
---|
1106 | SUBROUTINE getin_co2_parameters |
---|
1107 | |
---|
1108 | IMPLICIT NONE |
---|
1109 | |
---|
1110 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
1111 | |
---|
1112 | IF(first_call) THEN |
---|
1113 | |
---|
1114 | CALL getin('LAI_LEVEL_DEPTH',lai_level_depth) |
---|
1115 | CALL getin('X1_COEF',x1_coef) |
---|
1116 | CALL getin('X1_Q10',x1_Q10) |
---|
1117 | CALL getin('QUANTUM_YIELD',quantum_yield) |
---|
1118 | CALL getin('KT_COEF',kt_coef) |
---|
1119 | CALL getin('KC_COEF',kc_coef) |
---|
1120 | CALL getin('KO_Q10',Ko_Q10) |
---|
1121 | CALL getin('OA',Oa) |
---|
1122 | CALL getin('KO_COEF',Ko_coef) |
---|
1123 | CALL getin('CP_0',CP_0) |
---|
1124 | CALL getin('CP_TEMP_COEF',cp_temp_coef) |
---|
1125 | CALL getin('CP_TEMP_REF',cp_temp_ref) |
---|
1126 | CALL getin('RT_COEF',rt_coef) |
---|
1127 | CALL getin('VC_COEF',vc_coef) |
---|
1128 | |
---|
1129 | first_call =.FALSE. |
---|
1130 | |
---|
1131 | ENDIF |
---|
1132 | |
---|
1133 | END SUBROUTINE getin_co2_parameters |
---|
1134 | |
---|
1135 | !********************************************************** |
---|
1136 | |
---|
1137 | ! Subroutine called only if stomate is activated |
---|
1138 | |
---|
1139 | SUBROUTINE getin_stomate_parameters |
---|
1140 | |
---|
1141 | IMPLICIT NONE |
---|
1142 | |
---|
1143 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
1144 | |
---|
1145 | IF(first_call) THEN |
---|
1146 | |
---|
1147 | CALL getin('TOO_LONG',too_long) |
---|
1148 | !- |
---|
1149 | CALL getin('TAU_FIRE',tau_fire) |
---|
1150 | CALL getin('LITTER_CRIT',litter_crit) |
---|
1151 | !- |
---|
1152 | CALL getin('OK_MINRES',ok_minres) |
---|
1153 | CALL getin('TAU_LEAFINIT', tau_leafinit) |
---|
1154 | CALL getin('RESERVE_TIME_TREE',reserve_time_tree) |
---|
1155 | CALL getin('RESERVE_TIME_GRASS',reserve_time_grass) |
---|
1156 | CALL getin('R0',R0) |
---|
1157 | CALL getin('S0',S0) |
---|
1158 | CALL getin('F_FRUIT',f_fruit) |
---|
1159 | CALL getin('ALLOC_SAP_ABOVE_TREE',alloc_sap_above_tree) |
---|
1160 | CALL getin('ALLOC_SAP_ABOVE_GRASS',alloc_sap_above_grass) |
---|
1161 | CALL getin('MIN_LTOLSR',min_LtoLSR) |
---|
1162 | CALL getin('MAX_LTOLSR',max_LtoLSR) |
---|
1163 | CALL getin('Z_NITROGEN',z_nitrogen) |
---|
1164 | !- |
---|
1165 | CALL getin('PIPE_TUNE1',pipe_tune1) |
---|
1166 | CALL getin('PIPE_TUNE2',pipe_tune2) |
---|
1167 | CALL getin('PIPE_TUNE3',pipe_tune3) |
---|
1168 | CALL getin('PIPE_TUNE4',pipe_tune4) |
---|
1169 | CALL getin('PIPE_DENSITY',pipe_density) |
---|
1170 | CALL getin('PIPE_K1',pipe_k1) |
---|
1171 | CALL getin('ESTAB_MAX_TREE',estab_max_tree) |
---|
1172 | CALL getin('ESTAB_MAX_GRASS',estab_max_grass) |
---|
1173 | CALL getin('IND_0',ind_0) |
---|
1174 | CALL getin('MIN_COVER',min_cover) |
---|
1175 | CALL getin('PRECIP_CRIT',precip_crit) |
---|
1176 | CALL getin('GDD_CRIT_ESTAB',gdd_crit_estab) |
---|
1177 | CALL getin('FPC_CRIT',fpc_crit) |
---|
1178 | CALL getin('FRAC_GROWTHRESP',frac_growthresp) |
---|
1179 | CALL getin('ALPHA_GRASS',alpha_grass) |
---|
1180 | CALL getin('ALPHA_TREE',alpha_tree) |
---|
1181 | CALL getin('TLONG_REF_MAX',tlong_ref_max) |
---|
1182 | CALL getin('TLONG_REF_MIN',tlong_ref_min) |
---|
1183 | !- |
---|
1184 | CALL getin('MASS_RATIO_HEART_SAP',mass_ratio_heart_sap) |
---|
1185 | CALL getin('TAU_HUM_MONTH',tau_hum_month) |
---|
1186 | CALL getin('TAU_HUM_WEEK',tau_hum_week) |
---|
1187 | CALL getin('TAU_T2M_MONTH',tau_t2m_month) |
---|
1188 | CALL getin('TAU_T2M_WEEK',tau_t2m_week) |
---|
1189 | CALL getin('TAU_TSOIL_MONTH',tau_tsoil_month) |
---|
1190 | CALL getin('TAU_SOILHUM_MONTH',tau_soilhum_month) |
---|
1191 | CALL getin('TAU_GPP_WEEK',tau_gpp_week) |
---|
1192 | CALL getin('TAU_GDD',tau_gdd) |
---|
1193 | CALL getin('TAU_NGD',tau_ngd) |
---|
1194 | CALL getin('COEFF_TAU_LONGTERM',coeff_tau_longterm) |
---|
1195 | ! |
---|
1196 | CALL getin('FRAC_TURNOVER_DAILY',frac_turnover_daily) |
---|
1197 | !- |
---|
1198 | CALL getin('Z_DECOMP',z_decomp) |
---|
1199 | !- |
---|
1200 | CALL getin('TAX_MAX',tax_max) |
---|
1201 | !- |
---|
1202 | CALL getin('ALWAYS_INIT',always_init) |
---|
1203 | CALL getin('MIN_GROWTHINIT_TIME',min_growthinit_time) |
---|
1204 | CALL getin('MOIAVAIL_ALWAYS_TREE',moiavail_always_tree) |
---|
1205 | CALL getin('MOIAVAIL_ALWAYS_GRASS',moiavail_always_grass) |
---|
1206 | CALL getin('T_ALWAYS_ADD',t_always_add) |
---|
1207 | !- |
---|
1208 | CALL getin('VMAX_OFFSET',vmax_offset) |
---|
1209 | CALL getin('LEAFAGE_FIRSTMAX',leafage_firstmax) |
---|
1210 | CALL getin('LEAFAGE_LASTMAX',leafage_lastmax) |
---|
1211 | CALL getin('LEAFAGE_OLD',leafage_old) |
---|
1212 | !- |
---|
1213 | CALL getin('GPPFRAC_DORMANCE',gppfrac_dormance) |
---|
1214 | CALL getin('MIN_GPP_ALLOWED',min_gpp_allowed) |
---|
1215 | CALL getin('TAU_CLIMATOLOGY',tau_climatology) |
---|
1216 | CALL getin('HVC1',hvc1) |
---|
1217 | CALL getin('HVC2',hvc2) |
---|
1218 | CALL getin('LEAF_FRAC_HVC',leaf_frac_hvc) |
---|
1219 | !- |
---|
1220 | CALL getin('CO2FRAC',co2frac) |
---|
1221 | CALL getin('CN',CN) |
---|
1222 | CALL getin('LC',LC) |
---|
1223 | !- |
---|
1224 | CALL getin('FRAC_SOIL_STRUCT_AA',frac_soil_struct_aa) |
---|
1225 | CALL getin('FRAC_SOIL_STRUCT_AB',frac_soil_struct_ab) |
---|
1226 | CALL getin('FRAC_SOIL_STRUCT_SA',frac_soil_struct_sa) |
---|
1227 | CALL getin('FRAC_SOIL_STRUCT_SB',frac_soil_struct_sb) |
---|
1228 | CALL getin('FRAC_SOIL_METAB_AA',frac_soil_metab_aa) |
---|
1229 | CALL getin('FRAC_SOIL_METAB_AB',frac_soil_metab_ab) |
---|
1230 | !- |
---|
1231 | CALL getin('FRAC_CARB_AA',frac_carb_aa) |
---|
1232 | CALL getin('FRAC_CARB_AP',frac_carb_ap) |
---|
1233 | CALL getin('FRAC_CARB_SS',frac_carb_ss) |
---|
1234 | CALL getin('FRAC_CARB_SA',frac_carb_sa) |
---|
1235 | CALL getin('FRAC_CARB_SP',frac_carb_sp) |
---|
1236 | CALL getin('FRAC_CARB_PP',frac_carb_pp) |
---|
1237 | CALL getin('FRAC_CARB_PA',frac_carb_pa) |
---|
1238 | CALL getin('FRAC_CARB_PS',frac_carb_ps) |
---|
1239 | |
---|
1240 | !--------------------------------------- |
---|
1241 | ! COEFFICIENTS OF EQUATIONS |
---|
1242 | !------------------------------------- |
---|
1243 | ! |
---|
1244 | !- |
---|
1245 | CALL getin('BCFRAC_COEFF',bcfrac_coeff) |
---|
1246 | CALL getin('FIREFRAC_COEFF',firefrac_coeff) |
---|
1247 | !- |
---|
1248 | CALL getin('AVAILABILITY_FACT', availability_fact) |
---|
1249 | CALL getin('VIGOUR_REF',vigour_ref) |
---|
1250 | CALL getin('VIGOUR_COEFF',vigour_coeff) |
---|
1251 | !- |
---|
1252 | CALL getin('RIP_TIME_MIN',RIP_time_min) |
---|
1253 | CALL getin('NPP_LONGTERM_INIT',npp_longterm_init) |
---|
1254 | CALL getin('EVERYWHERE_INIT',everywhere_init) |
---|
1255 | ! |
---|
1256 | !- |
---|
1257 | CALL getin('LAI_MAX_TO_HAPPY',lai_max_to_happy) |
---|
1258 | CALL getin('NLIM_TREF',Nlim_tref) |
---|
1259 | !- |
---|
1260 | CALL getin('BM_SAPL_CARBRES',bm_sapl_carbres) |
---|
1261 | CALL getin('BM_SAPL_SAPABOVE',bm_sapl_sapabove) |
---|
1262 | CALL getin('BM_SAPL_HEARTABOVE',bm_sapl_heartabove) |
---|
1263 | CALL getin('BM_SAPL_HEARTBELOW',bm_sapl_heartbelow) |
---|
1264 | CALL getin('INIT_SAPL_MASS_LEAF_NAT',init_sapl_mass_leaf_nat) |
---|
1265 | CALL getin('INIT_SAPL_MASS_LEAF_AGRI',init_sapl_mass_leaf_agri) |
---|
1266 | CALL getin('INIT_SAPL_MASS_CARBRES',init_sapl_mass_carbres) |
---|
1267 | CALL getin('INIT_SAPL_MASS_ROOT',init_sapl_mass_root) |
---|
1268 | CALL getin('INIT_SAPL_MASS_FRUIT',init_sapl_mass_fruit) |
---|
1269 | CALL getin('CN_SAPL_INIT',cn_sapl_init) |
---|
1270 | CALL getin('MIGRATE_TREE',migrate_tree) |
---|
1271 | CALL getin('MIGRATE_GRASS',migrate_grass) |
---|
1272 | CALL getin('MAXDIA_COEFF',maxdia_coeff) |
---|
1273 | CALL getin('LAI_INITMIN_TREE',lai_initmin_tree) |
---|
1274 | CALL getin('LAI_INITMIN_GRASS',lai_initmin_grass) |
---|
1275 | CALL getin('DIA_COEFF',dia_coeff) |
---|
1276 | CALL getin('MAXDIA_COEFF',maxdia_coeff) |
---|
1277 | CALL getin('BM_SAPL_LEAF',bm_sapl_leaf) |
---|
1278 | !- |
---|
1279 | CALL getin('METABOLIC_REF_FRAC',metabolic_ref_frac) |
---|
1280 | CALL getin('METABOLIC_LN_RATIO',metabolic_LN_ratio) |
---|
1281 | CALL getin('TAU_METABOLIC',tau_metabolic) |
---|
1282 | CALL getin('TAU_STRUCT',tau_struct) |
---|
1283 | CALL getin('SOIL_Q10',soil_Q10) |
---|
1284 | CALL getin('TSOIL_REF',tsoil_ref) |
---|
1285 | CALL getin('LITTER_STRUCT_COEF',litter_struct_coef) |
---|
1286 | CALL getin('MOIST_COEFF',moist_coeff) |
---|
1287 | !- |
---|
1288 | CALL getin('GDDNCD_REF',gddncd_ref) |
---|
1289 | CALL getin('GDDNCD_CURVE',gddncd_curve) |
---|
1290 | CALL getin('GDDNCD_OFFSET',gddncd_offset) |
---|
1291 | !- |
---|
1292 | CALL getin('CN_TREE',cn_tree) |
---|
1293 | CALL getin('BM_SAPL_RESCALE',bm_sapl_rescale) |
---|
1294 | !- |
---|
1295 | CALL getin('MAINT_RESP_MIN_VMAX',maint_resp_min_vmax) |
---|
1296 | CALL getin('MAINT_RESP_COEFF',maint_resp_coeff) |
---|
1297 | !- |
---|
1298 | CALL getin('NCD_MAX_YEAR',ncd_max_year) |
---|
1299 | CALL getin('GDD_THRESHOLD',gdd_threshold) |
---|
1300 | CALL getin('GREEN_AGE_EVER',green_age_ever) |
---|
1301 | CALL getin('GREEN_AGE_DEC',green_age_dec) |
---|
1302 | !- |
---|
1303 | CALL getin('ACTIVE_TO_PASS_CLAY_FRAC',active_to_pass_clay_frac) |
---|
1304 | CALL getin('CARBON_TAU_IACTIVE',carbon_tau_iactive) |
---|
1305 | CALL getin('CARBON_TAU_ISLOW',carbon_tau_islow) |
---|
1306 | CALL getin('CARBON_TAU_IPASSIVE',carbon_tau_ipassive) |
---|
1307 | CALL getin('FLUX_TOT_COEFF',flux_tot_coeff) |
---|
1308 | !- |
---|
1309 | CALL getin('NEW_TURNOVER_TIME_REF',new_turnover_time_ref) |
---|
1310 | CALL getin('DT_TURNOVER_TIME',dt_turnover_time) |
---|
1311 | CALL getin('LEAF_AGE_CRIT_TREF',leaf_age_crit_tref) |
---|
1312 | CALL getin('LEAF_AGE_CRIT_COEFF',leaf_age_crit_coeff) |
---|
1313 | |
---|
1314 | first_call = .FALSE. |
---|
1315 | |
---|
1316 | ENDIF |
---|
1317 | |
---|
1318 | END SUBROUTINE getin_stomate_parameters |
---|
1319 | |
---|
1320 | !****************************************** |
---|
1321 | |
---|
1322 | SUBROUTINE getin_dgvm_parameters |
---|
1323 | |
---|
1324 | IMPLICIT NONE |
---|
1325 | |
---|
1326 | LOGICAL, SAVE :: first_call = .TRUE. |
---|
1327 | |
---|
1328 | IF(first_call) THEN |
---|
1329 | |
---|
1330 | CALL getin('PIPE_TUNE_EXP_COEFF',pipe_tune_exp_coeff) |
---|
1331 | ! |
---|
1332 | CALL getin('ESTABLISH_SCAL_FACT',establish_scal_fact) |
---|
1333 | CALL getin('FPC_CRIT_MAX',fpc_crit_max) |
---|
1334 | CALL getin('FPC_CRIT_MIN',fpc_crit_min) |
---|
1335 | ! |
---|
1336 | CALL getin('GRASS_MERCY',grass_mercy) |
---|
1337 | CALL getin('TREE_MERCY',tree_mercy) |
---|
1338 | CALL getin('ANNUAL_INCREASE',annual_increase) |
---|
1339 | ! |
---|
1340 | CALL getin('MIN_AVAIL',min_avail) |
---|
1341 | CALL getin('RIP_TIME_MIN',RIP_time_min) |
---|
1342 | CALL getin('NPP_LONGTERM_INIT',npp_longterm_init) |
---|
1343 | CALL getin('EVERYWHERE_INIT',everywhere_init) |
---|
1344 | |
---|
1345 | first_call = .FALSE. |
---|
1346 | |
---|
1347 | ENDIF |
---|
1348 | |
---|
1349 | |
---|
1350 | END SUBROUTINE getin_dgvm_parameters |
---|
1351 | |
---|
1352 | !-------------------- |
---|
1353 | END MODULE constantes |
---|