1 |
guez |
3 |
module conf_phys_m |
2 |
|
|
|
3 |
|
|
! This module is clean: no C preprocessor directive, no include line. |
4 |
|
|
|
5 |
|
|
implicit none |
6 |
|
|
|
7 |
|
|
integer iflag_pbl |
8 |
|
|
|
9 |
|
|
contains |
10 |
|
|
|
11 |
|
|
subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, & |
12 |
|
|
fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, & |
13 |
|
|
ratqsbas,ratqshaut,if_ebil, & |
14 |
|
|
ok_ade, ok_aie, & |
15 |
|
|
bl95_b0, bl95_b1,& |
16 |
|
|
iflag_thermals,nsplit_thermals) |
17 |
|
|
|
18 |
|
|
! From phylmd/conf_phys.F90,v 1.7 2005/07/05 07:21:23 |
19 |
|
|
|
20 |
guez |
30 |
use getincom, only: getin |
21 |
guez |
3 |
use clesphys, only: solaire, co2_ppm, rco2, ch4_ppb, rch4, n2o_ppb, rn2o, & |
22 |
|
|
cfc11_ppt, rcfc11, cfc12_ppt, rcfc12, top_height, overlap, cdmmax, & |
23 |
|
|
cdhmax, ksta, ksta_ter, ok_kzmin, lev_histhf, lev_histday, & |
24 |
|
|
lev_histmth, type_run, ok_isccp, ok_regdyn, lonmin_ins, lonmax_ins, & |
25 |
|
|
latmin_ins, latmax_ins, ecrit_ins, ecrit_hf, ecrit_hf2mth, & |
26 |
|
|
ecrit_day, ecrit_mth, ecrit_tra, ecrit_reg, bug_ozone |
27 |
guez |
38 |
use yomcst |
28 |
guez |
3 |
use conema3_m |
29 |
|
|
use comfisrtilp |
30 |
|
|
use nuagecom |
31 |
|
|
|
32 |
|
|
!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12 |
33 |
|
|
|
34 |
|
|
! Configuration de la "physique" de LMDZ a l'aide de la fonction |
35 |
|
|
! GETIN de IOIPSL |
36 |
|
|
|
37 |
|
|
! ocean: type d'ocean (force, slab, couple) |
38 |
|
|
! ok_veget: type de modele de vegetation |
39 |
|
|
! ok_journe: sorties journalieres |
40 |
|
|
! ok_mensuel: sorties mensuelles |
41 |
|
|
! ok_instan: sorties instantanees |
42 |
|
|
! ok_ade, ok_aie: apply or not aerosol direct and indirect effects |
43 |
|
|
! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc |
44 |
|
|
|
45 |
|
|
! Sortie: |
46 |
guez |
12 |
character(len=*), intent(out):: ocean |
47 |
guez |
3 |
logical :: ok_veget, ok_newmicro |
48 |
|
|
logical :: ok_journe, ok_mensuel, ok_instan |
49 |
|
|
LOGICAL :: ok_ade, ok_aie |
50 |
|
|
REAL :: bl95_b0, bl95_b1 |
51 |
guez |
17 |
real, intent(out):: fact_cldcon |
52 |
|
|
real, intent(out):: facttemps |
53 |
|
|
real ratqsbas,ratqshaut |
54 |
guez |
3 |
integer :: iflag_cldcon, if_ebil |
55 |
|
|
|
56 |
|
|
! Local |
57 |
|
|
real :: zzz |
58 |
|
|
|
59 |
|
|
integer :: iflag_thermals,nsplit_thermals |
60 |
|
|
|
61 |
|
|
!----------------------------------------------------------- |
62 |
|
|
|
63 |
|
|
print *, "Call sequence information: conf_phys" |
64 |
|
|
|
65 |
|
|
!Config Key = OCEAN |
66 |
|
|
!Config Desc = Type d'ocean |
67 |
|
|
!Config Def = force |
68 |
guez |
12 |
!Config Help = Type d'ocean utilise: force ou slab mais pas couple |
69 |
guez |
3 |
! |
70 |
|
|
ocean = 'force ' |
71 |
|
|
call getin('OCEAN', ocean) |
72 |
|
|
! |
73 |
|
|
!Config Key = VEGET |
74 |
|
|
!Config Desc = Type de modele de vegetation |
75 |
|
|
!Config Def = .false. |
76 |
|
|
!Config Help = Type de modele de vegetation utilise |
77 |
|
|
! |
78 |
|
|
ok_veget = .false. |
79 |
|
|
call getin('VEGET', ok_veget) |
80 |
|
|
! |
81 |
|
|
!Config Key = OK_journe |
82 |
|
|
!Config Desc = Pour des sorties journalieres |
83 |
|
|
!Config Def = .false. |
84 |
|
|
!Config Help = Pour creer le fichier histday contenant les sorties |
85 |
|
|
! journalieres |
86 |
|
|
! |
87 |
|
|
ok_journe = .false. |
88 |
|
|
call getin('OK_journe', ok_journe) |
89 |
|
|
! |
90 |
|
|
!Config Key = OK_mensuel |
91 |
|
|
!Config Desc = Pour des sorties mensuelles |
92 |
|
|
!Config Def = .true. |
93 |
|
|
!Config Help = Pour creer le fichier histmth contenant les sorties |
94 |
|
|
! mensuelles |
95 |
|
|
! |
96 |
|
|
ok_mensuel = .true. |
97 |
|
|
call getin('OK_mensuel', ok_mensuel) |
98 |
|
|
! |
99 |
|
|
!Config Key = OK_instan |
100 |
|
|
!Config Desc = Pour des sorties instantanees |
101 |
|
|
!Config Def = .false. |
102 |
|
|
!Config Help = Pour creer le fichier histins contenant les sorties |
103 |
|
|
! instantanees |
104 |
|
|
! |
105 |
|
|
ok_instan = .false. |
106 |
|
|
call getin('OK_instan', ok_instan) |
107 |
|
|
! |
108 |
|
|
!Config Key = ok_ade |
109 |
|
|
!Config Desc = Aerosol direct effect or not? |
110 |
|
|
!Config Def = .false. |
111 |
|
|
!Config Help = Used in radlwsw.F |
112 |
|
|
! |
113 |
|
|
ok_ade = .false. |
114 |
|
|
call getin('ok_ade', ok_ade) |
115 |
|
|
|
116 |
|
|
! |
117 |
|
|
!Config Key = ok_aie |
118 |
|
|
!Config Desc = Aerosol indirect effect or not? |
119 |
|
|
!Config Def = .false. |
120 |
|
|
!Config Help = Used in nuage.F and radlwsw.F |
121 |
|
|
! |
122 |
|
|
ok_aie = .false. |
123 |
|
|
call getin('ok_aie', ok_aie) |
124 |
|
|
|
125 |
|
|
! |
126 |
|
|
!Config Key = bl95_b0 |
127 |
|
|
!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
128 |
|
|
!Config Def = .false. |
129 |
|
|
!Config Help = Used in nuage.F |
130 |
|
|
! |
131 |
|
|
bl95_b0 = 2. |
132 |
|
|
call getin('bl95_b0', bl95_b0) |
133 |
|
|
|
134 |
|
|
!Config Key = bl95_b1 |
135 |
|
|
!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
136 |
|
|
!Config Def = .false. |
137 |
|
|
!Config Help = Used in nuage.F |
138 |
|
|
! |
139 |
|
|
bl95_b1 = 0.2 |
140 |
|
|
call getin('bl95_b1', bl95_b1) |
141 |
|
|
|
142 |
|
|
! |
143 |
|
|
! |
144 |
|
|
!Config Key = if_ebil |
145 |
|
|
!Config Desc = Niveau de sortie pour les diags bilan d'energie |
146 |
|
|
!Config Def = 0 |
147 |
|
|
!Config Help = |
148 |
|
|
! |
149 |
|
|
! |
150 |
|
|
if_ebil = 0 |
151 |
|
|
call getin('if_ebil', if_ebil) |
152 |
|
|
!! |
153 |
|
|
!! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG |
154 |
|
|
!! |
155 |
|
|
!Config Key = R_ecc |
156 |
|
|
!Config Desc = Excentricite |
157 |
|
|
!Config Def = 0.016715 |
158 |
|
|
!Config Help = |
159 |
|
|
! |
160 |
|
|
!valeur AMIP II |
161 |
|
|
R_ecc = 0.016715 |
162 |
|
|
call getin('R_ecc', R_ecc) |
163 |
|
|
!! |
164 |
|
|
!Config Key = R_peri |
165 |
|
|
!Config Desc = Equinoxe |
166 |
|
|
!Config Def = |
167 |
|
|
!Config Help = |
168 |
|
|
! |
169 |
|
|
! |
170 |
|
|
!valeur AMIP II |
171 |
|
|
R_peri = 102.7 |
172 |
|
|
call getin('R_peri', R_peri) |
173 |
|
|
!! |
174 |
|
|
!Config Key = R_incl |
175 |
|
|
!Config Desc = Inclinaison |
176 |
|
|
!Config Def = |
177 |
|
|
!Config Help = |
178 |
|
|
! |
179 |
|
|
! |
180 |
|
|
!valeur AMIP II |
181 |
|
|
R_incl = 23.441 |
182 |
|
|
call getin('R_incl', R_incl) |
183 |
|
|
!! |
184 |
|
|
!Config Key = solaire |
185 |
|
|
!Config Desc = Constante solaire en W/m2 |
186 |
|
|
!Config Def = 1365. |
187 |
|
|
!Config Help = |
188 |
|
|
! |
189 |
|
|
! |
190 |
|
|
!valeur AMIP II |
191 |
|
|
solaire = 1365. |
192 |
|
|
call getin('solaire', solaire) |
193 |
|
|
!! |
194 |
|
|
!Config Key = co2_ppm |
195 |
|
|
!Config Desc = concentration du gaz carbonique en ppmv |
196 |
|
|
!Config Def = 348. |
197 |
|
|
!Config Help = |
198 |
|
|
! |
199 |
|
|
! |
200 |
|
|
!valeur AMIP II |
201 |
|
|
co2_ppm = 348. |
202 |
|
|
call getin('co2_ppm', co2_ppm) |
203 |
|
|
!! |
204 |
|
|
!Config Key = RCO2 |
205 |
|
|
!Config Desc = Concentration du CO2 |
206 |
|
|
!Config Def = co2_ppm * 1.0e-06 * 44.011/28.97 |
207 |
|
|
!Config Def = 348. * 1.0e-06 * 44.011/28.97 |
208 |
|
|
!Config Help = |
209 |
|
|
! |
210 |
|
|
! RCO2 = 5.286789092164308E-04 |
211 |
|
|
!ancienne valeur |
212 |
|
|
RCO2 = co2_ppm * 1.0e-06 * 44.011/28.97 ! pour co2_ppm=348. |
213 |
|
|
|
214 |
|
|
!! call getin('RCO2', RCO2) |
215 |
|
|
!! |
216 |
|
|
!Config Key = RCH4 |
217 |
|
|
!Config Desc = Concentration du CH4 |
218 |
|
|
!Config Def = 1.65E-06* 16.043/28.97 |
219 |
|
|
!Config Help = |
220 |
|
|
! |
221 |
|
|
! |
222 |
|
|
!valeur AMIP II |
223 |
|
|
!OK RCH4 = 1.65E-06* 16.043/28.97 |
224 |
|
|
! RCH4 = 9.137366240938903E-07 |
225 |
|
|
! |
226 |
|
|
!ancienne valeur |
227 |
|
|
! RCH4 = 1.72E-06* 16.043/28.97 |
228 |
|
|
!OK call getin('RCH4', RCH4) |
229 |
|
|
zzz = 1650. |
230 |
|
|
call getin('CH4_ppb', zzz) |
231 |
|
|
CH4_ppb = zzz |
232 |
|
|
RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97 |
233 |
|
|
!! |
234 |
|
|
!Config Key = RN2O |
235 |
|
|
!Config Desc = Concentration du N2O |
236 |
|
|
!Config Def = 306.E-09* 44.013/28.97 |
237 |
|
|
!Config Help = |
238 |
|
|
! |
239 |
|
|
! |
240 |
|
|
!valeur AMIP II |
241 |
|
|
!OK RN2O = 306.E-09* 44.013/28.97 |
242 |
|
|
! RN2O = 4.648939592682085E-07 |
243 |
|
|
! |
244 |
|
|
!ancienne valeur |
245 |
|
|
! RN2O = 310.E-09* 44.013/28.97 |
246 |
|
|
!OK call getin('RN2O', RN2O) |
247 |
|
|
zzz=306. |
248 |
|
|
call getin('N2O_ppb', zzz) |
249 |
|
|
N2O_ppb = zzz |
250 |
|
|
RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97 |
251 |
|
|
!! |
252 |
|
|
!Config Key = RCFC11 |
253 |
|
|
!Config Desc = Concentration du CFC11 |
254 |
|
|
!Config Def = 280.E-12* 137.3686/28.97 |
255 |
|
|
!Config Help = |
256 |
|
|
! |
257 |
|
|
! |
258 |
|
|
!OK RCFC11 = 280.E-12* 137.3686/28.97 |
259 |
|
|
zzz = 280. |
260 |
|
|
call getin('CFC11_ppt',zzz) |
261 |
|
|
CFC11_ppt = zzz |
262 |
|
|
RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97 |
263 |
|
|
! RCFC11 = 1.327690990680013E-09 |
264 |
|
|
!OK call getin('RCFC11', RCFC11) |
265 |
|
|
!! |
266 |
|
|
!Config Key = RCFC12 |
267 |
|
|
!Config Desc = Concentration du CFC12 |
268 |
|
|
!Config Def = 484.E-12* 120.9140/28.97 |
269 |
|
|
!Config Help = |
270 |
|
|
! |
271 |
|
|
! |
272 |
|
|
!OK RCFC12 = 484.E-12* 120.9140/28.97 |
273 |
|
|
zzz = 484. |
274 |
|
|
call getin('CFC12_ppt',zzz) |
275 |
|
|
CFC12_ppt = zzz |
276 |
|
|
RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97 |
277 |
|
|
! RCFC12 = 2.020102726958923E-09 |
278 |
|
|
!OK call getin('RCFC12', RCFC12) |
279 |
|
|
!! |
280 |
|
|
!! Constante solaire & Parametres orbitaux & taux gaz effet de serre END |
281 |
|
|
!! |
282 |
|
|
!! KE |
283 |
|
|
! |
284 |
|
|
!Config Key = epmax |
285 |
|
|
!Config Desc = Efficacite precip |
286 |
|
|
!Config Def = 0.993 |
287 |
|
|
!Config Help = |
288 |
|
|
! |
289 |
|
|
epmax = .993 |
290 |
|
|
call getin('epmax', epmax) |
291 |
|
|
! |
292 |
|
|
!Config Key = ok_adj_ema |
293 |
|
|
!Config Desc = |
294 |
|
|
!Config Def = false |
295 |
|
|
!Config Help = |
296 |
|
|
! |
297 |
|
|
ok_adj_ema = .false. |
298 |
|
|
call getin('ok_adj_ema',ok_adj_ema) |
299 |
|
|
! |
300 |
|
|
!Config Key = iflag_clw |
301 |
|
|
!Config Desc = |
302 |
|
|
!Config Def = 0 |
303 |
|
|
!Config Help = |
304 |
|
|
! |
305 |
|
|
iflag_clw = 0 |
306 |
|
|
call getin('iflag_clw',iflag_clw) |
307 |
|
|
! |
308 |
|
|
!Config Key = cld_lc_lsc |
309 |
|
|
!Config Desc = |
310 |
|
|
!Config Def = 2.6e-4 |
311 |
|
|
!Config Help = |
312 |
|
|
! |
313 |
|
|
cld_lc_lsc = 2.6e-4 |
314 |
|
|
call getin('cld_lc_lsc',cld_lc_lsc) |
315 |
|
|
! |
316 |
|
|
!Config Key = cld_lc_con |
317 |
|
|
!Config Desc = |
318 |
|
|
!Config Def = 2.6e-4 |
319 |
|
|
!Config Help = |
320 |
|
|
! |
321 |
|
|
cld_lc_con = 2.6e-4 |
322 |
|
|
call getin('cld_lc_con',cld_lc_con) |
323 |
|
|
! |
324 |
|
|
!Config Key = cld_tau_lsc |
325 |
|
|
!Config Desc = |
326 |
|
|
!Config Def = 3600. |
327 |
|
|
!Config Help = |
328 |
|
|
! |
329 |
|
|
cld_tau_lsc = 3600. |
330 |
|
|
call getin('cld_tau_lsc',cld_tau_lsc) |
331 |
|
|
! |
332 |
|
|
!Config Key = cld_tau_con |
333 |
|
|
!Config Desc = |
334 |
|
|
!Config Def = 3600. |
335 |
|
|
!Config Help = |
336 |
|
|
! |
337 |
|
|
cld_tau_con = 3600. |
338 |
|
|
call getin('cld_tau_con',cld_tau_con) |
339 |
|
|
! |
340 |
|
|
!Config Key = ffallv_lsc |
341 |
|
|
!Config Desc = |
342 |
|
|
!Config Def = 1. |
343 |
|
|
!Config Help = |
344 |
|
|
! |
345 |
|
|
ffallv_lsc = 1. |
346 |
|
|
call getin('ffallv_lsc',ffallv_lsc) |
347 |
|
|
! |
348 |
|
|
!Config Key = ffallv_con |
349 |
|
|
!Config Desc = |
350 |
|
|
!Config Def = 1. |
351 |
|
|
!Config Help = |
352 |
|
|
! |
353 |
|
|
ffallv_con = 1. |
354 |
|
|
call getin('ffallv_con',ffallv_con) |
355 |
|
|
! |
356 |
|
|
!Config Key = coef_eva |
357 |
|
|
!Config Desc = |
358 |
|
|
!Config Def = 2.e-5 |
359 |
|
|
!Config Help = |
360 |
|
|
! |
361 |
|
|
coef_eva = 2.e-5 |
362 |
|
|
call getin('coef_eva',coef_eva) |
363 |
|
|
! |
364 |
|
|
!Config Key = reevap_ice |
365 |
|
|
!Config Desc = |
366 |
|
|
!Config Def = .false. |
367 |
|
|
!Config Help = |
368 |
|
|
! |
369 |
|
|
reevap_ice = .false. |
370 |
|
|
call getin('reevap_ice',reevap_ice) |
371 |
|
|
! |
372 |
|
|
!Config Key = iflag_cldcon |
373 |
|
|
!Config Desc = |
374 |
|
|
!Config Def = 1 |
375 |
|
|
!Config Help = |
376 |
|
|
! |
377 |
|
|
iflag_cldcon = 1 |
378 |
|
|
call getin('iflag_cldcon',iflag_cldcon) |
379 |
|
|
|
380 |
|
|
! |
381 |
|
|
!Config Key = iflag_pdf |
382 |
|
|
!Config Desc = |
383 |
|
|
!Config Def = 0 |
384 |
|
|
!Config Help = |
385 |
|
|
! |
386 |
|
|
iflag_pdf = 0 |
387 |
|
|
call getin('iflag_pdf',iflag_pdf) |
388 |
|
|
! |
389 |
|
|
!Config Key = fact_cldcon |
390 |
|
|
!Config Desc = |
391 |
|
|
!Config Def = 0.375 |
392 |
|
|
!Config Help = |
393 |
|
|
! |
394 |
|
|
fact_cldcon = 0.375 |
395 |
|
|
call getin('fact_cldcon',fact_cldcon) |
396 |
|
|
|
397 |
|
|
! |
398 |
|
|
!Config Key = facttemps |
399 |
|
|
!Config Desc = |
400 |
|
|
!Config Def = 1.e-4 |
401 |
|
|
!Config Help = |
402 |
|
|
! |
403 |
|
|
facttemps = 1.e-4 |
404 |
|
|
call getin('facttemps',facttemps) |
405 |
|
|
|
406 |
|
|
! |
407 |
|
|
!Config Key = ok_newmicro |
408 |
|
|
!Config Desc = |
409 |
|
|
!Config Def = .true. |
410 |
|
|
!Config Help = |
411 |
|
|
! |
412 |
|
|
ok_newmicro = .true. |
413 |
|
|
call getin('ok_newmicro',ok_newmicro) |
414 |
|
|
! |
415 |
|
|
!Config Key = ratqsbas |
416 |
|
|
!Config Desc = |
417 |
|
|
!Config Def = 0.01 |
418 |
|
|
!Config Help = |
419 |
|
|
! |
420 |
|
|
ratqsbas = 0.01 |
421 |
|
|
call getin('ratqsbas',ratqsbas) |
422 |
|
|
! |
423 |
|
|
!Config Key = ratqshaut |
424 |
|
|
!Config Desc = |
425 |
|
|
!Config Def = 0.3 |
426 |
|
|
!Config Help = |
427 |
|
|
! |
428 |
|
|
ratqshaut = 0.3 |
429 |
|
|
call getin('ratqshaut',ratqshaut) |
430 |
|
|
|
431 |
|
|
! |
432 |
|
|
!Config Key = rad_froid |
433 |
|
|
!Config Desc = |
434 |
|
|
!Config Def = 35.0 |
435 |
|
|
!Config Help = |
436 |
|
|
! |
437 |
|
|
rad_froid = 35.0 |
438 |
|
|
call getin('rad_froid',rad_froid) |
439 |
|
|
|
440 |
|
|
! |
441 |
|
|
!Config Key = rad_chau1 |
442 |
|
|
!Config Desc = |
443 |
|
|
!Config Def = 13.0 |
444 |
|
|
!Config Help = |
445 |
|
|
! |
446 |
|
|
rad_chau1 = 13.0 |
447 |
|
|
call getin('rad_chau1',rad_chau1) |
448 |
|
|
|
449 |
|
|
! |
450 |
|
|
!Config Key = rad_chau2 |
451 |
|
|
!Config Desc = |
452 |
|
|
!Config Def = 9.0 |
453 |
|
|
!Config Help = |
454 |
|
|
! |
455 |
|
|
rad_chau2 = 9.0 |
456 |
|
|
call getin('rad_chau2',rad_chau2) |
457 |
|
|
|
458 |
|
|
! |
459 |
|
|
!Config Key = top_height |
460 |
|
|
!Config Desc = |
461 |
|
|
!Config Def = 3 |
462 |
|
|
!Config Help = |
463 |
|
|
! |
464 |
|
|
top_height = 3 |
465 |
|
|
call getin('top_height',top_height) |
466 |
|
|
|
467 |
|
|
! |
468 |
|
|
!Config Key = overlap |
469 |
|
|
!Config Desc = |
470 |
|
|
!Config Def = 3 |
471 |
|
|
!Config Help = |
472 |
|
|
! |
473 |
|
|
overlap = 3 |
474 |
|
|
call getin('overlap',overlap) |
475 |
|
|
|
476 |
|
|
|
477 |
|
|
! |
478 |
|
|
! |
479 |
|
|
!Config Key = cdmmax |
480 |
|
|
!Config Desc = |
481 |
|
|
!Config Def = 1.3E-3 |
482 |
|
|
!Config Help = |
483 |
|
|
! |
484 |
|
|
cdmmax = 1.3E-3 |
485 |
|
|
call getin('cdmmax',cdmmax) |
486 |
|
|
|
487 |
|
|
! |
488 |
|
|
!Config Key = cdhmax |
489 |
|
|
!Config Desc = |
490 |
|
|
!Config Def = 1.1E-3 |
491 |
|
|
!Config Help = |
492 |
|
|
! |
493 |
|
|
cdhmax = 1.1E-3 |
494 |
|
|
call getin('cdhmax',cdhmax) |
495 |
|
|
|
496 |
|
|
!261103 |
497 |
|
|
! |
498 |
|
|
!Config Key = ksta |
499 |
|
|
!Config Desc = |
500 |
|
|
!Config Def = 1.0e-10 |
501 |
|
|
!Config Help = |
502 |
|
|
! |
503 |
|
|
ksta = 1.0e-10 |
504 |
|
|
call getin('ksta',ksta) |
505 |
|
|
|
506 |
|
|
! |
507 |
|
|
!Config Key = ksta_ter |
508 |
|
|
!Config Desc = |
509 |
|
|
!Config Def = 1.0e-10 |
510 |
|
|
!Config Help = |
511 |
|
|
! |
512 |
|
|
ksta_ter = 1.0e-10 |
513 |
|
|
call getin('ksta_ter',ksta_ter) |
514 |
|
|
|
515 |
|
|
! |
516 |
|
|
!Config Key = ok_kzmin |
517 |
|
|
!Config Desc = |
518 |
|
|
!Config Def = .true. |
519 |
|
|
!Config Help = |
520 |
|
|
! |
521 |
|
|
ok_kzmin = .true. |
522 |
|
|
call getin('ok_kzmin',ok_kzmin) |
523 |
|
|
|
524 |
|
|
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
525 |
|
|
! PARAMETER FOR THE PLANETARY BOUNDARY LAYER |
526 |
|
|
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
527 |
|
|
|
528 |
|
|
!Config Key = iflag_pbl |
529 |
|
|
!Config Desc = |
530 |
|
|
!Config Def = 1 |
531 |
|
|
!Config Help = |
532 |
|
|
! |
533 |
|
|
iflag_pbl = 1 |
534 |
|
|
call getin('iflag_pbl',iflag_pbl) |
535 |
|
|
! |
536 |
|
|
!Config Key = iflag_thermals |
537 |
|
|
!Config Desc = |
538 |
|
|
!Config Def = 0 |
539 |
|
|
!Config Help = |
540 |
|
|
! |
541 |
|
|
iflag_thermals = 0 |
542 |
|
|
call getin('iflag_thermals',iflag_thermals) |
543 |
|
|
! |
544 |
|
|
! |
545 |
|
|
!Config Key = nsplit_thermals |
546 |
|
|
!Config Desc = |
547 |
|
|
!Config Def = 1 |
548 |
|
|
!Config Help = |
549 |
|
|
! |
550 |
|
|
nsplit_thermals = 1 |
551 |
|
|
call getin('nsplit_thermals',nsplit_thermals) |
552 |
|
|
|
553 |
|
|
|
554 |
|
|
|
555 |
|
|
! |
556 |
|
|
!Config Key = lev_histhf |
557 |
|
|
!Config Desc = |
558 |
|
|
!Config Def = 0 |
559 |
|
|
!Config Help = |
560 |
|
|
! |
561 |
|
|
lev_histhf = 0 |
562 |
|
|
call getin('lev_histhf',lev_histhf) |
563 |
|
|
|
564 |
|
|
! |
565 |
|
|
!Config Key = lev_histday |
566 |
|
|
!Config Desc = |
567 |
|
|
!Config Def = 1 |
568 |
|
|
!Config Help = |
569 |
|
|
! |
570 |
|
|
lev_histday = 1 |
571 |
|
|
call getin('lev_histday',lev_histday) |
572 |
|
|
|
573 |
|
|
! |
574 |
|
|
!Config Key = lev_histmth |
575 |
|
|
!Config Desc = |
576 |
|
|
!Config Def = 2 |
577 |
|
|
!Config Help = |
578 |
|
|
! |
579 |
|
|
lev_histmth = 2 |
580 |
|
|
call getin('lev_histmth',lev_histmth) |
581 |
|
|
|
582 |
|
|
! |
583 |
|
|
!Config Key = type_run |
584 |
|
|
!Config Desc = |
585 |
|
|
!Config Def = 'AMIP' ou 'ENSP' |
586 |
|
|
!Config Help = |
587 |
|
|
! |
588 |
|
|
type_run = 'AMIP' |
589 |
|
|
call getin('type_run',type_run) |
590 |
|
|
|
591 |
|
|
! |
592 |
|
|
!Config Key = ok_isccp |
593 |
|
|
!Config Desc = |
594 |
|
|
!Config Def = .true. |
595 |
|
|
!Config Help = |
596 |
|
|
! |
597 |
|
|
! ok_isccp = .true. |
598 |
|
|
ok_isccp = .false. |
599 |
|
|
call getin('ok_isccp',ok_isccp) |
600 |
|
|
|
601 |
|
|
! |
602 |
|
|
! |
603 |
|
|
!Config Key = ok_regdyn |
604 |
|
|
!Config Desc = |
605 |
|
|
!Config Def = 'AMIP' |
606 |
|
|
!Config Help = |
607 |
|
|
! |
608 |
|
|
! ok_regdyn = .true. |
609 |
|
|
ok_regdyn = .false. |
610 |
|
|
call getin('ok_regdyn',ok_regdyn) |
611 |
|
|
! |
612 |
|
|
! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone |
613 |
|
|
! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc" |
614 |
|
|
! |
615 |
|
|
!Config Key = lonmin_ins |
616 |
|
|
!Config Desc = 100. |
617 |
|
|
!Config Def = longitude minimale sorties "bilKP_ins" |
618 |
|
|
!Config Help = |
619 |
|
|
! |
620 |
|
|
lonmin_ins = 100. |
621 |
|
|
call getin('lonmin_ins',lonmin_ins) |
622 |
|
|
! |
623 |
|
|
!Config Key = lonmax_ins |
624 |
|
|
!Config Desc = 130. |
625 |
|
|
!Config Def = longitude maximale sorties "bilKP_ins" |
626 |
|
|
!Config Help = |
627 |
|
|
! |
628 |
|
|
lonmax_ins = 130. |
629 |
|
|
call getin('lonmax_ins',lonmax_ins) |
630 |
|
|
! |
631 |
|
|
!Config Key = latmin_ins |
632 |
|
|
!Config Desc = -20. |
633 |
|
|
!Config Def = latitude minimale sorties "bilKP_ins" |
634 |
|
|
!Config Help = |
635 |
|
|
! |
636 |
|
|
latmin_ins = -20. |
637 |
|
|
call getin('latmin_ins',latmin_ins) |
638 |
|
|
! |
639 |
|
|
!Config Key = latmax_ins |
640 |
|
|
!Config Desc = 20. |
641 |
|
|
!Config Def = latitude maximale sorties "bilKP_ins" |
642 |
|
|
!Config Help = |
643 |
|
|
! |
644 |
|
|
latmax_ins = 20. |
645 |
|
|
call getin('latmax_ins',latmax_ins) |
646 |
|
|
! |
647 |
|
|
!Config Key = ecrit_ins |
648 |
|
|
!Config Desc = |
649 |
|
|
!Config Def = NINT(86400./dtime/48.) !a chaque pas de temps physique |
650 |
|
|
!Config Help = |
651 |
|
|
! |
652 |
|
|
! ecrit_ins = NINT(86400./dtime/48.) |
653 |
|
|
ecrit_ins = NINT(86400./48.) |
654 |
|
|
call getin('ecrit_ins',ecrit_ins) |
655 |
|
|
! |
656 |
|
|
!Config Key = ecrit_hf |
657 |
|
|
!Config Desc = |
658 |
|
|
!Config Def = NINT(86400./dtime *0.25) !toutes les 6h |
659 |
|
|
!Config Help = |
660 |
|
|
! |
661 |
|
|
! ecrit_hf = NINT(86400./dtime *0.25) |
662 |
|
|
ecrit_hf = NINT(86400. *0.25) |
663 |
|
|
call getin('ecrit_hf',ecrit_hf) |
664 |
|
|
! |
665 |
|
|
!Config Key = ecrit_hf2mth |
666 |
|
|
!Config Desc = |
667 |
|
|
!Config Def = 4*30 !ecriture mens. a partir de val. inst. toutes les 6h |
668 |
|
|
!Config Help = |
669 |
|
|
! |
670 |
|
|
ecrit_hf2mth = 4*30 |
671 |
|
|
call getin('ecrit_hf2mth',ecrit_hf2mth) |
672 |
|
|
! |
673 |
|
|
!Config Key = ecrit_day |
674 |
|
|
ecrit_day = 86400 |
675 |
|
|
call getin('ecrit_day',ecrit_day) |
676 |
|
|
! |
677 |
|
|
ecrit_mth = 86400 |
678 |
|
|
call getin('ecrit_mth',ecrit_mth) |
679 |
|
|
! |
680 |
|
|
ecrit_tra = 1 |
681 |
|
|
call getin('ecrit_tra',ecrit_tra) |
682 |
|
|
! |
683 |
|
|
!Config Key = ecrit_reg |
684 |
|
|
!Config Desc = |
685 |
|
|
!Config Def = NINT(86400./dtime *0.25) !4 fois par jour |
686 |
|
|
!Config Help = |
687 |
|
|
! |
688 |
|
|
! ecrit_reg = NINT(86400./dtime *0.25) !4 fois par jour |
689 |
|
|
ecrit_reg = NINT(86400. *0.25) !4 fois par jour |
690 |
|
|
call getin('ecrit_reg',ecrit_reg) |
691 |
|
|
! |
692 |
|
|
! |
693 |
|
|
! |
694 |
|
|
!Config Key = bug_ozone |
695 |
|
|
!Config Desc = Pour retrouver le bug de l'ozone (IPCC), mettre a true |
696 |
|
|
!Config Def = false |
697 |
|
|
!Config Help = |
698 |
|
|
! |
699 |
|
|
bug_ozone = .false. |
700 |
|
|
call getin('bug_ozone',bug_ozone) |
701 |
|
|
|
702 |
guez |
17 |
print *, ' *********' |
703 |
|
|
print *, ' Configuration des parametres de la physique: ' |
704 |
|
|
print *, ' Config ocean = ', ocean |
705 |
|
|
print *, ' Config veget = ', ok_veget |
706 |
|
|
print *, ' Sortie journaliere = ', ok_journe |
707 |
|
|
print *, ' Sortie mensuelle = ', ok_mensuel |
708 |
|
|
print *, ' Sortie instantanee = ', ok_instan |
709 |
|
|
print *, ' Sortie bilan d''energie, if_ebil =', if_ebil |
710 |
|
|
print *, ' Excentricite = ',R_ecc |
711 |
|
|
print *, ' Equinoxe = ',R_peri |
712 |
|
|
print *, ' Inclinaison =',R_incl |
713 |
|
|
print *, ' Constante solaire =',solaire |
714 |
|
|
print *, ' co2_ppm =',co2_ppm |
715 |
|
|
print *, ' RCO2 = ',RCO2 |
716 |
|
|
print *, ' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4 |
717 |
|
|
print *, ' N2O_ppb =',N2O_ppb,' RN2O = ',RN2O |
718 |
|
|
print *, ' CFC11_ppt=',CFC11_ppt,' RCFC11 = ',RCFC11 |
719 |
|
|
print *, ' CFC12_ppt=',CFC12_ppt,' RCFC12 = ',RCFC12 |
720 |
|
|
print *, ' epmax = ', epmax |
721 |
|
|
print *, ' ok_adj_ema = ', ok_adj_ema |
722 |
|
|
print *, ' iflag_clw = ', iflag_clw |
723 |
|
|
print *, ' cld_lc_lsc = ', cld_lc_lsc |
724 |
|
|
print *, ' cld_lc_con = ', cld_lc_con |
725 |
|
|
print *, ' cld_tau_lsc = ', cld_tau_lsc |
726 |
|
|
print *, ' cld_tau_con = ', cld_tau_con |
727 |
|
|
print *, ' ffallv_lsc = ', ffallv_lsc |
728 |
|
|
print *, ' ffallv_con = ', ffallv_con |
729 |
|
|
print *, ' coef_eva = ', coef_eva |
730 |
|
|
print *, ' reevap_ice = ', reevap_ice |
731 |
|
|
print *, ' iflag_pdf = ', iflag_pdf |
732 |
|
|
print *, ' iflag_cldcon = ', iflag_cldcon |
733 |
|
|
print *, ' fact_cldcon = ', fact_cldcon |
734 |
|
|
print *, ' facttemps = ', facttemps |
735 |
|
|
print *, ' ok_newmicro = ',ok_newmicro |
736 |
|
|
print *, ' ratqsbas = ',ratqsbas |
737 |
|
|
print *, ' ratqshaut = ',ratqshaut |
738 |
|
|
print *, ' top_height = ',top_height |
739 |
|
|
print *, ' overlap = ',overlap |
740 |
|
|
print *, ' cdmmax = ',cdmmax |
741 |
|
|
print *, ' cdhmax = ',cdhmax |
742 |
|
|
print *, ' ksta = ',ksta |
743 |
|
|
print *, ' ksta_ter = ',ksta_ter |
744 |
|
|
print *, ' ok_kzmin = ',ok_kzmin |
745 |
|
|
print *, ' ok_ade = ',ok_ade |
746 |
|
|
print *, ' ok_aie = ',ok_aie |
747 |
|
|
print *, ' bl95_b0 = ',bl95_b0 |
748 |
|
|
print *, ' bl95_b1 = ',bl95_b1 |
749 |
|
|
print *, ' lev_histhf = ',lev_histhf |
750 |
|
|
print *, ' lev_histday = ',lev_histday |
751 |
|
|
print *, ' lev_histmth = ',lev_histmth |
752 |
|
|
print *, ' iflag_pbl = ', iflag_pbl |
753 |
|
|
print *, ' iflag_thermals = ', iflag_thermals |
754 |
|
|
print *, ' type_run = ',type_run |
755 |
|
|
print *, ' ok_isccp = ',ok_isccp |
756 |
|
|
print *, ' ok_regdyn = ',ok_regdyn |
757 |
|
|
print *, ' lonmin lonmax latmin latmax bilKP_ins =',& |
758 |
guez |
3 |
lonmin_ins, lonmax_ins, latmin_ins, latmax_ins |
759 |
guez |
17 |
print *, 'ecrit_ ins, hf, hf2mth, day, mth, reg, tra', ecrit_ins, & |
760 |
guez |
3 |
ecrit_hf, ecrit_hf2mth, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra |
761 |
guez |
17 |
print *, ' bug_ozone = ', bug_ozone |
762 |
guez |
3 |
|
763 |
|
|
end subroutine conf_phys |
764 |
|
|
|
765 |
|
|
end module conf_phys_m |