1 | |
---|
2 | ! |
---|
3 | ! $Id: conf_phys.F90 1539 2011-06-08 22:13:33Z fhourdin $ |
---|
4 | ! |
---|
5 | ! |
---|
6 | ! |
---|
7 | module conf_phys_m |
---|
8 | |
---|
9 | implicit none |
---|
10 | |
---|
11 | contains |
---|
12 | |
---|
13 | subroutine conf_phys(ok_journe, ok_mensuel, ok_instan, ok_hf, & |
---|
14 | ok_LES,& |
---|
15 | callstats,& |
---|
16 | solarlong0,seuil_inversion, & |
---|
17 | fact_cldcon, facttemps,ok_newmicro,iflag_radia,& |
---|
18 | iflag_cldcon, & |
---|
19 | iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, & |
---|
20 | ok_ade, ok_aie, aerosol_couple, & |
---|
21 | flag_aerosol, new_aod, & |
---|
22 | bl95_b0, bl95_b1,& |
---|
23 | read_climoz, & |
---|
24 | alp_offset) |
---|
25 | |
---|
26 | use IOIPSL |
---|
27 | USE surface_data |
---|
28 | USE phys_cal_mod |
---|
29 | USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl |
---|
30 | use control_mod |
---|
31 | |
---|
32 | include "conema3.h" |
---|
33 | include "fisrtilp.h" |
---|
34 | include "nuage.h" |
---|
35 | include "YOMCST.h" |
---|
36 | include "YOMCST2.h" |
---|
37 | |
---|
38 | include "thermcell.h" |
---|
39 | |
---|
40 | !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12 |
---|
41 | include "clesphys.h" |
---|
42 | include "compbl.h" |
---|
43 | include "comsoil.h" |
---|
44 | ! |
---|
45 | ! Configuration de la "physique" de LMDZ a l'aide de la fonction |
---|
46 | ! GETIN de IOIPSL |
---|
47 | ! |
---|
48 | ! LF 05/2001 |
---|
49 | ! |
---|
50 | |
---|
51 | ! |
---|
52 | ! type_ocean: type d'ocean (force, slab, couple) |
---|
53 | ! version_ocean: version d'ocean (opa8/nemo pour type_ocean=couple ou |
---|
54 | ! sicOBS pour type_ocean=slab) |
---|
55 | ! ok_veget: type de modele de vegetation |
---|
56 | ! ok_journe: sorties journalieres |
---|
57 | ! ok_hf: sorties haute frequence |
---|
58 | ! ok_mensuel: sorties mensuelles |
---|
59 | ! ok_instan: sorties instantanees |
---|
60 | ! ok_ade, ok_aie: apply or not aerosol direct and indirect effects |
---|
61 | ! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc |
---|
62 | ! |
---|
63 | |
---|
64 | |
---|
65 | ! Sortie: |
---|
66 | logical :: ok_newmicro |
---|
67 | integer :: iflag_radia |
---|
68 | logical :: ok_journe, ok_mensuel, ok_instan, ok_hf |
---|
69 | logical :: ok_LES |
---|
70 | LOGICAL :: callstats |
---|
71 | LOGICAL :: ok_ade, ok_aie, aerosol_couple |
---|
72 | INTEGER :: flag_aerosol |
---|
73 | LOGICAL :: new_aod |
---|
74 | REAL :: bl95_b0, bl95_b1 |
---|
75 | real :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs |
---|
76 | integer :: iflag_cldcon |
---|
77 | integer :: iflag_ratqs |
---|
78 | |
---|
79 | character (len = 6),SAVE :: type_ocean_omp, version_ocean_omp, ocean_omp |
---|
80 | CHARACTER(len = 8),SAVE :: aer_type_omp |
---|
81 | logical,SAVE :: ok_veget_omp, ok_newmicro_omp |
---|
82 | logical,SAVE :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp |
---|
83 | logical,SAVE :: ok_LES_omp |
---|
84 | LOGICAL,SAVE :: callstats_omp |
---|
85 | LOGICAL,SAVE :: ok_ade_omp, ok_aie_omp, aerosol_couple_omp |
---|
86 | INTEGER, SAVE :: flag_aerosol_omp |
---|
87 | LOGICAL, SAVE :: new_aod_omp |
---|
88 | REAL,SAVE :: bl95_b0_omp, bl95_b1_omp |
---|
89 | REAL,SAVE :: freq_ISCCP_omp, ecrit_ISCCP_omp |
---|
90 | REAL,SAVE :: freq_COSP_omp |
---|
91 | real,SAVE :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp |
---|
92 | real,SAVE :: ratqshaut_omp |
---|
93 | real,SAVE :: tau_ratqs_omp |
---|
94 | integer,SAVE :: iflag_radia_omp |
---|
95 | integer,SAVE :: iflag_rrtm_omp |
---|
96 | integer,SAVE :: iflag_cldcon_omp, ip_ebil_phy_omp |
---|
97 | integer,SAVE :: iflag_ratqs_omp |
---|
98 | |
---|
99 | Real,SAVE :: f_cdrag_ter_omp,f_cdrag_oce_omp |
---|
100 | Real,SAVE :: f_rugoro_omp |
---|
101 | |
---|
102 | ! Local |
---|
103 | integer :: numout = 6 |
---|
104 | real :: zzz |
---|
105 | |
---|
106 | real :: seuil_inversion |
---|
107 | real,save :: seuil_inversion_omp |
---|
108 | |
---|
109 | integer,SAVE :: iflag_thermals_ed_omp,iflag_thermals_optflux_omp |
---|
110 | integer,SAVE :: iflag_thermals_omp,nsplit_thermals_omp |
---|
111 | real,save :: tau_thermals_omp,alp_bl_k_omp |
---|
112 | real :: alp_offset |
---|
113 | REAL, SAVE :: alp_offset_omp |
---|
114 | integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp |
---|
115 | integer,SAVE :: iflag_cvl_sigd_omp |
---|
116 | REAL, SAVE :: supcrit1_omp, supcrit2_omp |
---|
117 | INTEGER, SAVE :: iflag_mix_omp |
---|
118 | real, save :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp |
---|
119 | |
---|
120 | REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp |
---|
121 | REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp |
---|
122 | REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp |
---|
123 | REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp |
---|
124 | REAL,SAVE :: CFC11_ppt_omp,RCFC11_omp,CFC11_ppt_per_omp,RCFC11_per_omp |
---|
125 | REAL,SAVE :: CFC12_ppt_omp,RCFC12_omp,CFC12_ppt_per_omp,RCFC12_per_omp |
---|
126 | REAL,SAVE :: epmax_omp |
---|
127 | LOGICAL,SAVE :: ok_adj_ema_omp |
---|
128 | INTEGER,SAVE :: iflag_clw_omp |
---|
129 | REAL,SAVE :: cld_lc_lsc_omp,cld_lc_con_omp,cld_tau_lsc_omp,cld_tau_con_omp |
---|
130 | REAL,SAVE :: ffallv_lsc_omp, ffallv_con_omp,coef_eva_omp |
---|
131 | LOGICAL,SAVE :: reevap_ice_omp |
---|
132 | INTEGER,SAVE :: iflag_pdf_omp |
---|
133 | REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp |
---|
134 | REAL,SAVE :: t_glace_min_omp, t_glace_max_omp |
---|
135 | REAL,SAVE :: rei_min_omp, rei_max_omp |
---|
136 | REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp |
---|
137 | REAL,SAVE :: qsol0_omp |
---|
138 | REAL :: solarlong0 |
---|
139 | REAL,SAVE :: solarlong0_omp |
---|
140 | INTEGER,SAVE :: top_height_omp,overlap_omp |
---|
141 | REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp |
---|
142 | LOGICAL,SAVE :: ok_kzmin_omp |
---|
143 | REAL, SAVE :: fmagic_omp, pmagic_omp |
---|
144 | INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp |
---|
145 | Integer, save :: lev_histins_omp, lev_histLES_omp |
---|
146 | INTEGER, SAVE :: lev_histdayNMC_omp |
---|
147 | LOGICAL, SAVE :: ok_histNMC_omp(3) |
---|
148 | REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3) |
---|
149 | CHARACTER*4, SAVE :: type_run_omp |
---|
150 | LOGICAL,SAVE :: ok_isccp_omp |
---|
151 | LOGICAL,SAVE :: ok_cosp_omp |
---|
152 | LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp |
---|
153 | REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp |
---|
154 | REAL,SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp |
---|
155 | REAL,SAVE :: ecrit_ins_omp |
---|
156 | REAL,SAVE :: ecrit_LES_omp |
---|
157 | REAL,SAVE :: ecrit_tra_omp |
---|
158 | REAL,SAVE :: cvl_corr_omp |
---|
159 | LOGICAL,SAVE :: ok_lic_melt_omp |
---|
160 | ! |
---|
161 | LOGICAL,SAVE :: cycle_diurne_omp,soil_model_omp,new_oliq_omp |
---|
162 | LOGICAL,SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp |
---|
163 | INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp |
---|
164 | LOGICAL,SAVE :: ok_strato_omp |
---|
165 | LOGICAL,SAVE :: ok_hines_omp |
---|
166 | LOGICAL,SAVE :: carbon_cycle_tr_omp |
---|
167 | LOGICAL,SAVE :: carbon_cycle_cpl_omp |
---|
168 | |
---|
169 | integer, intent(out):: read_climoz ! read ozone climatology, OpenMP shared |
---|
170 | ! Allowed values are 0, 1 and 2 |
---|
171 | ! 0: do not read an ozone climatology |
---|
172 | ! 1: read a single ozone climatology that will be used day and night |
---|
173 | ! 2: read two ozone climatologies, the average day and night |
---|
174 | ! climatology and the daylight climatology |
---|
175 | |
---|
176 | !$OMP MASTER |
---|
177 | !Config Key = type_ocean |
---|
178 | !Config Desc = Type d'ocean |
---|
179 | !Config Def = force |
---|
180 | !Config Help = Type d'ocean utilise: force, slab,couple |
---|
181 | ! |
---|
182 | type_ocean_omp = 'force ' |
---|
183 | call getin('type_ocean', type_ocean_omp) |
---|
184 | ! |
---|
185 | !Config Key = version_ocean |
---|
186 | !Config Desc = Version d'ocean |
---|
187 | !Config Def = xxxxxx |
---|
188 | !Config Help = Version d'ocean utilise: opa8/nemo/sicOBS/xxxxxx |
---|
189 | ! |
---|
190 | version_ocean_omp = 'xxxxxx' |
---|
191 | call getin('version_ocean', version_ocean_omp) |
---|
192 | |
---|
193 | !Config Key = OCEAN |
---|
194 | !Config Desc = Old parameter name for type_ocean |
---|
195 | !Config Def = yyyyyy |
---|
196 | !Config Help = This is only for testing purpose |
---|
197 | ! |
---|
198 | ocean_omp = 'yyyyyy' |
---|
199 | call getin('OCEAN', ocean_omp) |
---|
200 | IF (ocean_omp /= 'yyyyyy') THEN |
---|
201 | WRITE(numout,*)'ERROR!! Old variable name OCEAN used in parmeter file.' |
---|
202 | WRITE(numout,*)'Variable OCEAN has been replaced by the variable type_ocean.' |
---|
203 | WRITE(numout,*)'You have to update your parameter file physiq.def to succed running' |
---|
204 | CALL abort_gcm('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1) |
---|
205 | END IF |
---|
206 | |
---|
207 | ! |
---|
208 | !Config Key = VEGET |
---|
209 | !Config Desc = Type de modele de vegetation |
---|
210 | !Config Def = .false. |
---|
211 | !Config Help = Type de modele de vegetation utilise |
---|
212 | ! |
---|
213 | ok_veget_omp = .false. |
---|
214 | call getin('VEGET', ok_veget_omp) |
---|
215 | ! |
---|
216 | !Config Key = OK_journe |
---|
217 | !Config Desc = Pour des sorties journalieres |
---|
218 | !Config Def = .false. |
---|
219 | !Config Help = Pour creer le fichier histday contenant les sorties |
---|
220 | ! journalieres |
---|
221 | ! |
---|
222 | ok_journe_omp = .false. |
---|
223 | call getin('OK_journe', ok_journe_omp) |
---|
224 | ! |
---|
225 | !Config Key = ok_hf |
---|
226 | !Config Desc = Pour des sorties haute frequence |
---|
227 | !Config Def = .false. |
---|
228 | !Config Help = Pour creer le fichier histhf contenant les sorties |
---|
229 | ! haute frequence ( 3h ou 6h) |
---|
230 | ! |
---|
231 | ok_hf_omp = .false. |
---|
232 | call getin('ok_hf', ok_hf_omp) |
---|
233 | ! |
---|
234 | !Config Key = OK_mensuel |
---|
235 | !Config Desc = Pour des sorties mensuelles |
---|
236 | !Config Def = .true. |
---|
237 | !Config Help = Pour creer le fichier histmth contenant les sorties |
---|
238 | ! mensuelles |
---|
239 | ! |
---|
240 | ok_mensuel_omp = .true. |
---|
241 | call getin('OK_mensuel', ok_mensuel_omp) |
---|
242 | ! |
---|
243 | !Config Key = OK_instan |
---|
244 | !Config Desc = Pour des sorties instantanees |
---|
245 | !Config Def = .false. |
---|
246 | !Config Help = Pour creer le fichier histins contenant les sorties |
---|
247 | ! instantanees |
---|
248 | ! |
---|
249 | ok_instan_omp = .false. |
---|
250 | call getin('OK_instan', ok_instan_omp) |
---|
251 | ! |
---|
252 | !Config Key = ok_ade |
---|
253 | !Config Desc = Aerosol direct effect or not? |
---|
254 | !Config Def = .false. |
---|
255 | !Config Help = Used in radlwsw.F |
---|
256 | ! |
---|
257 | ok_ade_omp = .false. |
---|
258 | call getin('ok_ade', ok_ade_omp) |
---|
259 | |
---|
260 | ! |
---|
261 | !Config Key = ok_aie |
---|
262 | !Config Desc = Aerosol indirect effect or not? |
---|
263 | !Config Def = .false. |
---|
264 | !Config Help = Used in nuage.F and radlwsw.F |
---|
265 | ! |
---|
266 | ok_aie_omp = .false. |
---|
267 | call getin('ok_aie', ok_aie_omp) |
---|
268 | |
---|
269 | ! |
---|
270 | !Config Key = aerosol_couple |
---|
271 | !Config Desc = read aerosol in file or calcul by inca |
---|
272 | !Config Def = .false. |
---|
273 | !Config Help = Used in physiq.F |
---|
274 | ! |
---|
275 | aerosol_couple_omp = .false. |
---|
276 | CALL getin('aerosol_couple',aerosol_couple_omp) |
---|
277 | |
---|
278 | ! |
---|
279 | !Config Key = flag_aerosol |
---|
280 | !Config Desc = which aerosol is use for coupled model |
---|
281 | !Config Def = 1 |
---|
282 | !Config Help = Used in physiq.F |
---|
283 | ! |
---|
284 | ! - flag_aerosol=1 => so4 only (defaut) |
---|
285 | ! - flag_aerosol=2 => bc only |
---|
286 | ! - flag_aerosol=3 => pom only |
---|
287 | ! - flag_aerosol=4 => seasalt only |
---|
288 | ! - flag_aerosol=5 => dust only |
---|
289 | ! - flag_aerosol=6 => all aerosol |
---|
290 | |
---|
291 | flag_aerosol_omp = 1 |
---|
292 | CALL getin('flag_aerosol',flag_aerosol_omp) |
---|
293 | |
---|
294 | ! Temporary variable for testing purpose!! |
---|
295 | !Config Key = new_aod |
---|
296 | !Config Desc = which calcul of aeropt |
---|
297 | !Config Def = false |
---|
298 | !Config Help = Used in physiq.F |
---|
299 | ! |
---|
300 | new_aod_omp = .true. |
---|
301 | CALL getin('new_aod',new_aod_omp) |
---|
302 | |
---|
303 | ! |
---|
304 | !Config Key = aer_type |
---|
305 | !Config Desc = Use a constant field for the aerosols |
---|
306 | !Config Def = scenario |
---|
307 | !Config Help = Used in readaerosol.F90 |
---|
308 | ! |
---|
309 | aer_type_omp = 'scenario' |
---|
310 | call getin('aer_type', aer_type_omp) |
---|
311 | |
---|
312 | ! |
---|
313 | !Config Key = bl95_b0 |
---|
314 | !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
315 | !Config Def = .false. |
---|
316 | !Config Help = Used in nuage.F |
---|
317 | ! |
---|
318 | bl95_b0_omp = 2. |
---|
319 | call getin('bl95_b0', bl95_b0_omp) |
---|
320 | |
---|
321 | !Config Key = bl95_b1 |
---|
322 | !Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995) |
---|
323 | !Config Def = .false. |
---|
324 | !Config Help = Used in nuage.F |
---|
325 | ! |
---|
326 | bl95_b1_omp = 0.2 |
---|
327 | call getin('bl95_b1', bl95_b1_omp) |
---|
328 | |
---|
329 | !Config Key = freq_ISCCP |
---|
330 | !Config Desc = Frequence d'appel du simulateur ISCCP en secondes; |
---|
331 | ! par defaut 10800, i.e. 3 heures |
---|
332 | !Config Def = 10800. |
---|
333 | !Config Help = Used in ini_histISCCP.h |
---|
334 | ! |
---|
335 | freq_ISCCP_omp = 10800. |
---|
336 | call getin('freq_ISCCP', freq_ISCCP_omp) |
---|
337 | ! |
---|
338 | !Config Key = ecrit_ISCCP |
---|
339 | !Config Desc = Frequence d'ecriture des resultats du simulateur ISCCP en nombre de jours; |
---|
340 | ! par defaut 1., i.e. 1 jour |
---|
341 | !Config Def = 1. |
---|
342 | !Config Help = Used in ini_histISCCP.h |
---|
343 | ! |
---|
344 | ! |
---|
345 | ecrit_ISCCP_omp = 1. |
---|
346 | call getin('ecrit_ISCCP', ecrit_ISCCP_omp) |
---|
347 | |
---|
348 | !Config Key = freq_COSP |
---|
349 | !Config Desc = Frequence d'appel du simulateur COSP en secondes; |
---|
350 | ! par defaut 10800, i.e. 3 heures |
---|
351 | !Config Def = 10800. |
---|
352 | !Config Help = Used in ini_histdayCOSP.h |
---|
353 | ! |
---|
354 | freq_COSP_omp = 10800. |
---|
355 | call getin('freq_COSP', freq_COSP_omp) |
---|
356 | |
---|
357 | ! |
---|
358 | !Config Key = ip_ebil_phy |
---|
359 | !Config Desc = Niveau de sortie pour les diags bilan d'energie |
---|
360 | !Config Def = 0 |
---|
361 | !Config Help = |
---|
362 | ! |
---|
363 | ip_ebil_phy_omp = 0 |
---|
364 | call getin('ip_ebil_phy', ip_ebil_phy_omp) |
---|
365 | ! |
---|
366 | !Config Key = seuil_inversion |
---|
367 | !Config Desc = Seuil ur dTh pour le choix entre les schemas de CL |
---|
368 | !Config Def = -0.1 |
---|
369 | !Config Help = |
---|
370 | ! |
---|
371 | seuil_inversion_omp = -0.1 |
---|
372 | call getin('seuil_inversion', seuil_inversion_omp) |
---|
373 | |
---|
374 | !! |
---|
375 | !! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG |
---|
376 | !! |
---|
377 | !Config Key = R_ecc |
---|
378 | !Config Desc = Excentricite |
---|
379 | !Config Def = 0.016715 |
---|
380 | !Config Help = |
---|
381 | ! |
---|
382 | !valeur AMIP II |
---|
383 | R_ecc_omp = 0.016715 |
---|
384 | call getin('R_ecc', R_ecc_omp) |
---|
385 | !! |
---|
386 | !Config Key = R_peri |
---|
387 | !Config Desc = Equinoxe |
---|
388 | !Config Def = |
---|
389 | !Config Help = |
---|
390 | ! |
---|
391 | ! |
---|
392 | !valeur AMIP II |
---|
393 | R_peri_omp = 102.7 |
---|
394 | call getin('R_peri', R_peri_omp) |
---|
395 | !! |
---|
396 | !Config Key = R_incl |
---|
397 | !Config Desc = Inclinaison |
---|
398 | !Config Def = |
---|
399 | !Config Help = |
---|
400 | ! |
---|
401 | ! |
---|
402 | !valeur AMIP II |
---|
403 | R_incl_omp = 23.441 |
---|
404 | call getin('R_incl', R_incl_omp) |
---|
405 | !! |
---|
406 | !Config Key = solaire |
---|
407 | !Config Desc = Constante solaire en W/m2 |
---|
408 | !Config Def = 1365. |
---|
409 | !Config Help = |
---|
410 | ! |
---|
411 | ! |
---|
412 | !valeur AMIP II |
---|
413 | solaire_omp = 1365. |
---|
414 | call getin('solaire', solaire_omp) |
---|
415 | !! |
---|
416 | !Config Key = co2_ppm |
---|
417 | !Config Desc = concentration du gaz carbonique en ppmv |
---|
418 | !Config Def = 348. |
---|
419 | !Config Help = |
---|
420 | ! |
---|
421 | ! |
---|
422 | !valeur AMIP II |
---|
423 | co2_ppm_omp = 348. |
---|
424 | call getin('co2_ppm', co2_ppm_omp) |
---|
425 | !! |
---|
426 | !Config Key = RCO2 |
---|
427 | !Config Desc = Concentration du CO2 |
---|
428 | !Config Def = co2_ppm * 1.0e-06 * 44.011/28.97 |
---|
429 | !Config Def = 348. * 1.0e-06 * 44.011/28.97 |
---|
430 | !Config Help = |
---|
431 | ! |
---|
432 | ! RCO2 = 5.286789092164308E-04 |
---|
433 | !ancienne valeur |
---|
434 | RCO2_omp = co2_ppm_omp * 1.0e-06 * 44.011/28.97 ! pour co2_ppm=348. |
---|
435 | |
---|
436 | !! call getin('RCO2', RCO2) |
---|
437 | !! |
---|
438 | !Config Key = RCH4 |
---|
439 | !Config Desc = Concentration du CH4 |
---|
440 | !Config Def = 1.65E-06* 16.043/28.97 |
---|
441 | !Config Help = |
---|
442 | ! |
---|
443 | ! |
---|
444 | !valeur AMIP II |
---|
445 | !OK RCH4 = 1.65E-06* 16.043/28.97 |
---|
446 | ! RCH4 = 9.137366240938903E-07 |
---|
447 | ! |
---|
448 | !ancienne valeur |
---|
449 | ! RCH4 = 1.72E-06* 16.043/28.97 |
---|
450 | !OK call getin('RCH4', RCH4) |
---|
451 | zzz = 1650. |
---|
452 | call getin('CH4_ppb', zzz) |
---|
453 | CH4_ppb_omp = zzz |
---|
454 | RCH4_omp = CH4_ppb_omp * 1.0E-09 * 16.043/28.97 |
---|
455 | !! |
---|
456 | !Config Key = RN2O |
---|
457 | !Config Desc = Concentration du N2O |
---|
458 | !Config Def = 306.E-09* 44.013/28.97 |
---|
459 | !Config Help = |
---|
460 | ! |
---|
461 | ! |
---|
462 | !valeur AMIP II |
---|
463 | !OK RN2O = 306.E-09* 44.013/28.97 |
---|
464 | ! RN2O = 4.648939592682085E-07 |
---|
465 | ! |
---|
466 | !ancienne valeur |
---|
467 | ! RN2O = 310.E-09* 44.013/28.97 |
---|
468 | !OK call getin('RN2O', RN2O) |
---|
469 | zzz=306. |
---|
470 | call getin('N2O_ppb', zzz) |
---|
471 | N2O_ppb_omp = zzz |
---|
472 | RN2O_omp = N2O_ppb_omp * 1.0E-09 * 44.013/28.97 |
---|
473 | !! |
---|
474 | !Config Key = RCFC11 |
---|
475 | !Config Desc = Concentration du CFC11 |
---|
476 | !Config Def = 280.E-12* 137.3686/28.97 |
---|
477 | !Config Help = |
---|
478 | ! |
---|
479 | ! |
---|
480 | !OK RCFC11 = 280.E-12* 137.3686/28.97 |
---|
481 | zzz = 280. |
---|
482 | call getin('CFC11_ppt',zzz) |
---|
483 | CFC11_ppt_omp = zzz |
---|
484 | RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * 137.3686/28.97 |
---|
485 | ! RCFC11 = 1.327690990680013E-09 |
---|
486 | !OK call getin('RCFC11', RCFC11) |
---|
487 | !! |
---|
488 | !Config Key = RCFC12 |
---|
489 | !Config Desc = Concentration du CFC12 |
---|
490 | !Config Def = 484.E-12* 120.9140/28.97 |
---|
491 | !Config Help = |
---|
492 | ! |
---|
493 | ! |
---|
494 | !OK RCFC12 = 484.E-12* 120.9140/28.97 |
---|
495 | zzz = 484. |
---|
496 | call getin('CFC12_ppt',zzz) |
---|
497 | CFC12_ppt_omp = zzz |
---|
498 | RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * 120.9140/28.97 |
---|
499 | ! RCFC12 = 2.020102726958923E-09 |
---|
500 | !OK call getin('RCFC12', RCFC12) |
---|
501 | |
---|
502 | !ajout CFMIP begin |
---|
503 | !! |
---|
504 | !Config Key = co2_ppm_per |
---|
505 | !Config Desc = concentration du co2_ppm_per |
---|
506 | !Config Def = 348. |
---|
507 | !Config Help = |
---|
508 | ! |
---|
509 | co2_ppm_per_omp = co2_ppm_omp |
---|
510 | call getin('co2_ppm_per', co2_ppm_per_omp) |
---|
511 | !! |
---|
512 | !Config Key = RCO2_per |
---|
513 | !Config Desc = Concentration du CO2_per |
---|
514 | !Config Def = co2_ppm_per * 1.0e-06 * 44.011/28.97 |
---|
515 | !Config Def = 348. * 1.0e-06 * 44.011/28.97 |
---|
516 | !Config Help = |
---|
517 | ! |
---|
518 | RCO2_per_omp = co2_ppm_per_omp * 1.0e-06 * 44.011/28.97 |
---|
519 | |
---|
520 | !Config Key = RCH4_per |
---|
521 | !Config Desc = Concentration du CH4_per |
---|
522 | !Config Def = 1.65E-06* 16.043/28.97 |
---|
523 | !Config Help = |
---|
524 | ! |
---|
525 | zzz = CH4_ppb_omp |
---|
526 | call getin('CH4_ppb_per', zzz) |
---|
527 | CH4_ppb_per_omp = zzz |
---|
528 | RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * 16.043/28.97 |
---|
529 | !! |
---|
530 | !Config Key = RN2O_per |
---|
531 | !Config Desc = Concentration du N2O_per |
---|
532 | !Config Def = 306.E-09* 44.013/28.97 |
---|
533 | !Config Help = |
---|
534 | ! |
---|
535 | zzz = N2O_ppb_omp |
---|
536 | call getin('N2O_ppb_per', zzz) |
---|
537 | N2O_ppb_per_omp = zzz |
---|
538 | RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * 44.013/28.97 |
---|
539 | !! |
---|
540 | !Config Key = RCFC11_per |
---|
541 | !Config Desc = Concentration du CFC11_per |
---|
542 | !Config Def = 280.E-12* 137.3686/28.97 |
---|
543 | !Config Help = |
---|
544 | ! |
---|
545 | zzz = CFC11_ppt_omp |
---|
546 | call getin('CFC11_ppt_per',zzz) |
---|
547 | CFC11_ppt_per_omp = zzz |
---|
548 | RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * 137.3686/28.97 |
---|
549 | !! |
---|
550 | !Config Key = RCFC12_per |
---|
551 | !Config Desc = Concentration du CFC12_per |
---|
552 | !Config Def = 484.E-12* 120.9140/28.97 |
---|
553 | !Config Help = |
---|
554 | ! |
---|
555 | zzz = CFC12_ppt_omp |
---|
556 | call getin('CFC12_ppt_per',zzz) |
---|
557 | CFC12_ppt_per_omp = zzz |
---|
558 | RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * 120.9140/28.97 |
---|
559 | !ajout CFMIP end |
---|
560 | |
---|
561 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
562 | ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique |
---|
563 | ! Constantes precedemment dans dyn3d/conf_gcm |
---|
564 | |
---|
565 | !Config Key = cycle_diurne |
---|
566 | !Config Desc = Cycle ddiurne |
---|
567 | !Config Def = y |
---|
568 | !Config Help = Cette option permet d'eteidre le cycle diurne. |
---|
569 | !Config Peut etre util pour accelerer le code ! |
---|
570 | cycle_diurne_omp = .TRUE. |
---|
571 | CALL getin('cycle_diurne',cycle_diurne_omp) |
---|
572 | |
---|
573 | !Config Key = soil_model |
---|
574 | !Config Desc = Modele de sol |
---|
575 | !Config Def = y |
---|
576 | !Config Help = Choix du modele de sol (Thermique ?) |
---|
577 | !Config Option qui pourait un string afin de pouvoir |
---|
578 | !Config plus de choix ! Ou meme une liste d'options ! |
---|
579 | soil_model_omp = .TRUE. |
---|
580 | CALL getin('soil_model',soil_model_omp) |
---|
581 | |
---|
582 | !Config Key = new_oliq |
---|
583 | !Config Desc = Nouvelle eau liquide |
---|
584 | !Config Def = y |
---|
585 | !Config Help = Permet de mettre en route la |
---|
586 | !Config nouvelle parametrisation de l'eau liquide ! |
---|
587 | new_oliq_omp = .TRUE. |
---|
588 | CALL getin('new_oliq',new_oliq_omp) |
---|
589 | |
---|
590 | !Config Key = ok_orodr |
---|
591 | !Config Desc = Orodr ??? |
---|
592 | !Config Def = y |
---|
593 | !Config Help = Y en a pas comprendre ! |
---|
594 | !Config |
---|
595 | ok_orodr_omp = .TRUE. |
---|
596 | CALL getin('ok_orodr',ok_orodr_omp) |
---|
597 | |
---|
598 | !Config Key = ok_orolf |
---|
599 | !Config Desc = Orolf ?? |
---|
600 | !Config Def = y |
---|
601 | !Config Help = Connais pas ! |
---|
602 | ok_orolf_omp = .TRUE. |
---|
603 | CALL getin('ok_orolf', ok_orolf_omp) |
---|
604 | |
---|
605 | !Config Key = ok_limitvrai |
---|
606 | !Config Desc = Force la lecture de la bonne annee |
---|
607 | !Config Def = n |
---|
608 | !Config Help = On peut forcer le modele a lire le |
---|
609 | !Config fichier SST de la bonne annee. C'est une tres bonne |
---|
610 | !Config idee, pourquoi ne pas mettre toujours a y ??? |
---|
611 | ok_limitvrai_omp = .FALSE. |
---|
612 | CALL getin('ok_limitvrai',ok_limitvrai_omp) |
---|
613 | |
---|
614 | !Config Key = nbapp_rad |
---|
615 | !Config Desc = Frequence d'appel au rayonnement |
---|
616 | !Config Def = 12 |
---|
617 | !Config Help = Nombre d'appels des routines de rayonnements |
---|
618 | !Config par jour. |
---|
619 | nbapp_rad_omp = 12 |
---|
620 | CALL getin('nbapp_rad',nbapp_rad_omp) |
---|
621 | |
---|
622 | !Config Key = iflag_con |
---|
623 | !Config Desc = Flag de convection |
---|
624 | !Config Def = 2 |
---|
625 | !Config Help = Flag pour la convection les options suivantes existent : |
---|
626 | !Config 1 pour LMD, |
---|
627 | !Config 2 pour Tiedtke, |
---|
628 | !Config 3 pour CCM(NCAR) |
---|
629 | iflag_con_omp = 2 |
---|
630 | CALL getin('iflag_con',iflag_con_omp) |
---|
631 | |
---|
632 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
633 | !! |
---|
634 | !! Constante solaire & Parametres orbitaux & taux gaz effet de serre END |
---|
635 | !! |
---|
636 | !! KE |
---|
637 | ! |
---|
638 | |
---|
639 | !Config key = cvl_corr |
---|
640 | !Config Desc = Facteur multiplication des precip convectives dans KE |
---|
641 | !Config Def = 1.00 |
---|
642 | !Config Help = 1.02 pour un moderne ou un pre-ind. A ajuster pour un glaciaire |
---|
643 | cvl_corr_omp = 1.00 |
---|
644 | CALL getin('cvl_corr', cvl_corr_omp) |
---|
645 | |
---|
646 | |
---|
647 | !Config Key = epmax |
---|
648 | !Config Desc = Efficacite precip |
---|
649 | !Config Def = 0.993 |
---|
650 | !Config Help = |
---|
651 | ! |
---|
652 | epmax_omp = .993 |
---|
653 | call getin('epmax', epmax_omp) |
---|
654 | ! |
---|
655 | !Config Key = ok_adj_ema |
---|
656 | !Config Desc = |
---|
657 | !Config Def = false |
---|
658 | !Config Help = |
---|
659 | ! |
---|
660 | ok_adj_ema_omp = .false. |
---|
661 | call getin('ok_adj_ema',ok_adj_ema_omp) |
---|
662 | ! |
---|
663 | !Config Key = iflag_clw |
---|
664 | !Config Desc = |
---|
665 | !Config Def = 0 |
---|
666 | !Config Help = |
---|
667 | ! |
---|
668 | iflag_clw_omp = 0 |
---|
669 | call getin('iflag_clw',iflag_clw_omp) |
---|
670 | ! |
---|
671 | !Config Key = cld_lc_lsc |
---|
672 | !Config Desc = |
---|
673 | !Config Def = 2.6e-4 |
---|
674 | !Config Help = |
---|
675 | ! |
---|
676 | cld_lc_lsc_omp = 2.6e-4 |
---|
677 | call getin('cld_lc_lsc',cld_lc_lsc_omp) |
---|
678 | ! |
---|
679 | !Config Key = cld_lc_con |
---|
680 | !Config Desc = |
---|
681 | !Config Def = 2.6e-4 |
---|
682 | !Config Help = |
---|
683 | ! |
---|
684 | cld_lc_con_omp = 2.6e-4 |
---|
685 | call getin('cld_lc_con',cld_lc_con_omp) |
---|
686 | ! |
---|
687 | !Config Key = cld_tau_lsc |
---|
688 | !Config Desc = |
---|
689 | !Config Def = 3600. |
---|
690 | !Config Help = |
---|
691 | ! |
---|
692 | cld_tau_lsc_omp = 3600. |
---|
693 | call getin('cld_tau_lsc',cld_tau_lsc_omp) |
---|
694 | ! |
---|
695 | !Config Key = cld_tau_con |
---|
696 | !Config Desc = |
---|
697 | !Config Def = 3600. |
---|
698 | !Config Help = |
---|
699 | ! |
---|
700 | cld_tau_con_omp = 3600. |
---|
701 | call getin('cld_tau_con',cld_tau_con_omp) |
---|
702 | ! |
---|
703 | !Config Key = ffallv_lsc |
---|
704 | !Config Desc = |
---|
705 | !Config Def = 1. |
---|
706 | !Config Help = |
---|
707 | ! |
---|
708 | ffallv_lsc_omp = 1. |
---|
709 | call getin('ffallv_lsc',ffallv_lsc_omp) |
---|
710 | ! |
---|
711 | !Config Key = ffallv_con |
---|
712 | !Config Desc = |
---|
713 | !Config Def = 1. |
---|
714 | !Config Help = |
---|
715 | ! |
---|
716 | ffallv_con_omp = 1. |
---|
717 | call getin('ffallv_con',ffallv_con_omp) |
---|
718 | ! |
---|
719 | !Config Key = coef_eva |
---|
720 | !Config Desc = |
---|
721 | !Config Def = 2.e-5 |
---|
722 | !Config Help = |
---|
723 | ! |
---|
724 | coef_eva_omp = 2.e-5 |
---|
725 | call getin('coef_eva',coef_eva_omp) |
---|
726 | ! |
---|
727 | !Config Key = reevap_ice |
---|
728 | !Config Desc = |
---|
729 | !Config Def = .false. |
---|
730 | !Config Help = |
---|
731 | ! |
---|
732 | reevap_ice_omp = .false. |
---|
733 | call getin('reevap_ice',reevap_ice_omp) |
---|
734 | |
---|
735 | !Config Key = iflag_ratqs |
---|
736 | !Config Desc = |
---|
737 | !Config Def = 1 |
---|
738 | !Config Help = |
---|
739 | ! |
---|
740 | iflag_ratqs_omp = 1 |
---|
741 | call getin('iflag_ratqs',iflag_ratqs_omp) |
---|
742 | |
---|
743 | ! |
---|
744 | !Config Key = iflag_radia |
---|
745 | !Config Desc = |
---|
746 | !Config Def = 1 |
---|
747 | !Config Help = |
---|
748 | ! |
---|
749 | iflag_radia_omp = 1 |
---|
750 | call getin('iflag_radia',iflag_radia_omp) |
---|
751 | |
---|
752 | ! |
---|
753 | !Config Key = iflag_rrtm |
---|
754 | !Config Desc = |
---|
755 | !Config Def = 0 |
---|
756 | !Config Help = |
---|
757 | ! |
---|
758 | iflag_rrtm_omp = 0 |
---|
759 | call getin('iflag_rrtm',iflag_rrtm_omp) |
---|
760 | |
---|
761 | ! |
---|
762 | !Config Key = iflag_cldcon |
---|
763 | !Config Desc = |
---|
764 | !Config Def = 1 |
---|
765 | !Config Help = |
---|
766 | ! |
---|
767 | iflag_cldcon_omp = 1 |
---|
768 | call getin('iflag_cldcon',iflag_cldcon_omp) |
---|
769 | |
---|
770 | ! |
---|
771 | !Config Key = iflag_pdf |
---|
772 | !Config Desc = |
---|
773 | !Config Def = 0 |
---|
774 | !Config Help = |
---|
775 | ! |
---|
776 | iflag_pdf_omp = 0 |
---|
777 | call getin('iflag_pdf',iflag_pdf_omp) |
---|
778 | ! |
---|
779 | !Config Key = fact_cldcon |
---|
780 | !Config Desc = |
---|
781 | !Config Def = 0.375 |
---|
782 | !Config Help = |
---|
783 | ! |
---|
784 | fact_cldcon_omp = 0.375 |
---|
785 | call getin('fact_cldcon',fact_cldcon_omp) |
---|
786 | |
---|
787 | ! |
---|
788 | !Config Key = facttemps |
---|
789 | !Config Desc = |
---|
790 | !Config Def = 1.e-4 |
---|
791 | !Config Help = |
---|
792 | ! |
---|
793 | facttemps_omp = 1.e-4 |
---|
794 | call getin('facttemps',facttemps_omp) |
---|
795 | |
---|
796 | ! |
---|
797 | !Config Key = ok_newmicro |
---|
798 | !Config Desc = |
---|
799 | !Config Def = .true. |
---|
800 | !Config Help = |
---|
801 | ! |
---|
802 | ok_newmicro_omp = .true. |
---|
803 | call getin('ok_newmicro',ok_newmicro_omp) |
---|
804 | ! |
---|
805 | !Config Key = ratqsbas |
---|
806 | !Config Desc = |
---|
807 | !Config Def = 0.01 |
---|
808 | !Config Help = |
---|
809 | ! |
---|
810 | ratqsbas_omp = 0.01 |
---|
811 | call getin('ratqsbas',ratqsbas_omp) |
---|
812 | ! |
---|
813 | !Config Key = ratqshaut |
---|
814 | !Config Desc = |
---|
815 | !Config Def = 0.3 |
---|
816 | !Config Help = |
---|
817 | ! |
---|
818 | ratqshaut_omp = 0.3 |
---|
819 | call getin('ratqshaut',ratqshaut_omp) |
---|
820 | |
---|
821 | !Config Key = tau_ratqs |
---|
822 | !Config Desc = |
---|
823 | !Config Def = 1800. |
---|
824 | !Config Help = |
---|
825 | ! |
---|
826 | tau_ratqs_omp = 1800. |
---|
827 | call getin('tau_ratqs',tau_ratqs_omp) |
---|
828 | |
---|
829 | ! |
---|
830 | !----------------------------------------------------------------------- |
---|
831 | ! Longitude solaire pour le calcul de l'ensoleillement en degre |
---|
832 | ! si on veut imposer la saison. Sinon, solarlong0=-999.999 |
---|
833 | !Config Key = solarlong0 |
---|
834 | !Config Desc = |
---|
835 | !Config Def = -999.999 |
---|
836 | !Config Help = |
---|
837 | ! |
---|
838 | solarlong0_omp = -999.999 |
---|
839 | call getin('solarlong0',solarlong0_omp) |
---|
840 | ! |
---|
841 | !----------------------------------------------------------------------- |
---|
842 | ! Valeur imposee de l'humidite du sol pour le modele bucket. |
---|
843 | !Config Key = qsol0 |
---|
844 | !Config Desc = |
---|
845 | !Config Def = -1. |
---|
846 | !Config Help = |
---|
847 | ! |
---|
848 | qsol0_omp = -1. |
---|
849 | call getin('qsol0',qsol0_omp) |
---|
850 | ! |
---|
851 | !----------------------------------------------------------------------- |
---|
852 | ! |
---|
853 | !Config Key = inertie_ice |
---|
854 | !Config Desc = |
---|
855 | !Config Def = 2000. |
---|
856 | !Config Help = |
---|
857 | ! |
---|
858 | inertie_ice_omp = 2000. |
---|
859 | call getin('inertie_ice',inertie_ice_omp) |
---|
860 | ! |
---|
861 | !Config Key = inertie_sno |
---|
862 | !Config Desc = |
---|
863 | !Config Def = 2000. |
---|
864 | !Config Help = |
---|
865 | ! |
---|
866 | inertie_sno_omp = 2000. |
---|
867 | call getin('inertie_sno',inertie_sno_omp) |
---|
868 | ! |
---|
869 | !Config Key = inertie_sol |
---|
870 | !Config Desc = |
---|
871 | !Config Def = 2000. |
---|
872 | !Config Help = |
---|
873 | ! |
---|
874 | inertie_sol_omp = 2000. |
---|
875 | call getin('inertie_sol',inertie_sol_omp) |
---|
876 | |
---|
877 | ! |
---|
878 | !Config Key = rad_froid |
---|
879 | !Config Desc = |
---|
880 | !Config Def = 35.0 |
---|
881 | !Config Help = |
---|
882 | ! |
---|
883 | rad_froid_omp = 35.0 |
---|
884 | call getin('rad_froid',rad_froid_omp) |
---|
885 | |
---|
886 | ! |
---|
887 | !Config Key = rad_chau1 |
---|
888 | !Config Desc = |
---|
889 | !Config Def = 13.0 |
---|
890 | !Config Help = |
---|
891 | ! |
---|
892 | rad_chau1_omp = 13.0 |
---|
893 | call getin('rad_chau1',rad_chau1_omp) |
---|
894 | |
---|
895 | ! |
---|
896 | !Config Key = rad_chau2 |
---|
897 | !Config Desc = |
---|
898 | !Config Def = 9.0 |
---|
899 | !Config Help = |
---|
900 | ! |
---|
901 | rad_chau2_omp = 9.0 |
---|
902 | call getin('rad_chau2',rad_chau2_omp) |
---|
903 | |
---|
904 | ! |
---|
905 | !Config Key = t_glace_min |
---|
906 | !Config Desc = |
---|
907 | !Config Def = 258. |
---|
908 | !Config Help = |
---|
909 | ! |
---|
910 | t_glace_min_omp = 258. |
---|
911 | call getin('t_glace_min',t_glace_min_omp) |
---|
912 | |
---|
913 | ! |
---|
914 | !Config Key = t_glace_max |
---|
915 | !Config Desc = |
---|
916 | !Config Def = 273.13 |
---|
917 | !Config Help = |
---|
918 | ! |
---|
919 | t_glace_max_omp = 273.13 |
---|
920 | call getin('t_glace_max',t_glace_max_omp) |
---|
921 | |
---|
922 | !Config Key = rei_min |
---|
923 | !Config Desc = |
---|
924 | !Config Def = 3.5 |
---|
925 | !Config Help = |
---|
926 | ! |
---|
927 | rei_min_omp = 3.5 |
---|
928 | call getin('rei_min',rei_min_omp) |
---|
929 | |
---|
930 | ! |
---|
931 | !Config Key = rei_max |
---|
932 | !Config Desc = |
---|
933 | !Config Def = 61.29 |
---|
934 | !Config Help = |
---|
935 | ! |
---|
936 | rei_max_omp = 61.29 |
---|
937 | call getin('rei_max',rei_max_omp) |
---|
938 | |
---|
939 | ! |
---|
940 | !Config Key = top_height |
---|
941 | !Config Desc = |
---|
942 | !Config Def = 3 |
---|
943 | !Config Help = |
---|
944 | ! |
---|
945 | top_height_omp = 3 |
---|
946 | call getin('top_height',top_height_omp) |
---|
947 | |
---|
948 | ! |
---|
949 | !Config Key = overlap |
---|
950 | !Config Desc = |
---|
951 | !Config Def = 3 |
---|
952 | !Config Help = |
---|
953 | ! |
---|
954 | overlap_omp = 3 |
---|
955 | call getin('overlap',overlap_omp) |
---|
956 | |
---|
957 | |
---|
958 | ! |
---|
959 | ! |
---|
960 | !Config Key = cdmmax |
---|
961 | !Config Desc = |
---|
962 | !Config Def = 1.3E-3 |
---|
963 | !Config Help = |
---|
964 | ! |
---|
965 | cdmmax_omp = 1.3E-3 |
---|
966 | call getin('cdmmax',cdmmax_omp) |
---|
967 | |
---|
968 | ! |
---|
969 | !Config Key = cdhmax |
---|
970 | !Config Desc = |
---|
971 | !Config Def = 1.1E-3 |
---|
972 | !Config Help = |
---|
973 | ! |
---|
974 | cdhmax_omp = 1.1E-3 |
---|
975 | call getin('cdhmax',cdhmax_omp) |
---|
976 | |
---|
977 | !261103 |
---|
978 | ! |
---|
979 | !Config Key = ksta |
---|
980 | !Config Desc = |
---|
981 | !Config Def = 1.0e-10 |
---|
982 | !Config Help = |
---|
983 | ! |
---|
984 | ksta_omp = 1.0e-10 |
---|
985 | call getin('ksta',ksta_omp) |
---|
986 | |
---|
987 | ! |
---|
988 | !Config Key = ksta_ter |
---|
989 | !Config Desc = |
---|
990 | !Config Def = 1.0e-10 |
---|
991 | !Config Help = |
---|
992 | ! |
---|
993 | ksta_ter_omp = 1.0e-10 |
---|
994 | call getin('ksta_ter',ksta_ter_omp) |
---|
995 | |
---|
996 | ! |
---|
997 | !Config Key = ok_kzmin |
---|
998 | !Config Desc = |
---|
999 | !Config Def = .true. |
---|
1000 | !Config Help = |
---|
1001 | ! |
---|
1002 | ok_kzmin_omp = .true. |
---|
1003 | call getin('ok_kzmin',ok_kzmin_omp) |
---|
1004 | |
---|
1005 | ! |
---|
1006 | !Config Key = fmagic |
---|
1007 | !Config Desc = additionnal multiplicator factor used for albedo |
---|
1008 | !Config Def = 1. |
---|
1009 | !Config Help = additionnal multiplicator factor used in albedo.F |
---|
1010 | ! |
---|
1011 | fmagic_omp = 1. |
---|
1012 | call getin('fmagic',fmagic_omp) |
---|
1013 | |
---|
1014 | ! |
---|
1015 | !Config Key = pmagic |
---|
1016 | !Config Desc = additional factor used for albedo |
---|
1017 | !Config Def = 0. |
---|
1018 | !Config Help = additional factor used in albedo.F |
---|
1019 | ! |
---|
1020 | pmagic_omp = 0. |
---|
1021 | call getin('pmagic',pmagic_omp) |
---|
1022 | |
---|
1023 | |
---|
1024 | !Config Key = ok_lic_melt |
---|
1025 | !Config Desc = Prise en compte de la fonte de la calotte dans le bilan d'eau |
---|
1026 | !Config Def = .false. |
---|
1027 | !Config Help = mettre a .false. pour assurer la conservation en eau |
---|
1028 | ok_lic_melt_omp = .false. |
---|
1029 | call getin('ok_lic_melt', ok_lic_melt_omp) |
---|
1030 | |
---|
1031 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
1032 | ! PARAMETER FOR THE PLANETARY BOUNDARY LAYER |
---|
1033 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
1034 | |
---|
1035 | !Config Key = iflag_pbl |
---|
1036 | !Config Desc = |
---|
1037 | !Config Def = 1 |
---|
1038 | !Config Help = |
---|
1039 | ! |
---|
1040 | iflag_pbl_omp = 1 |
---|
1041 | call getin('iflag_pbl',iflag_pbl_omp) |
---|
1042 | ! |
---|
1043 | !Config Key = iflag_thermals |
---|
1044 | !Config Desc = |
---|
1045 | !Config Def = 0 |
---|
1046 | !Config Help = |
---|
1047 | ! |
---|
1048 | iflag_thermals_omp = 0 |
---|
1049 | call getin('iflag_thermals',iflag_thermals_omp) |
---|
1050 | ! |
---|
1051 | ! |
---|
1052 | !Config Key = iflag_thermals_ed |
---|
1053 | !Config Desc = |
---|
1054 | !Config Def = 0 |
---|
1055 | !Config Help = |
---|
1056 | ! |
---|
1057 | iflag_thermals_ed_omp = 0 |
---|
1058 | call getin('iflag_thermals_ed',iflag_thermals_ed_omp) |
---|
1059 | ! |
---|
1060 | ! |
---|
1061 | !Config Key = iflag_thermals_optflux |
---|
1062 | !Config Desc = |
---|
1063 | !Config Def = 0 |
---|
1064 | !Config Help = |
---|
1065 | ! |
---|
1066 | iflag_thermals_optflux_omp = 0 |
---|
1067 | call getin('iflag_thermals_optflux',iflag_thermals_optflux_omp) |
---|
1068 | ! |
---|
1069 | ! |
---|
1070 | !Config Key = nsplit_thermals |
---|
1071 | !Config Desc = |
---|
1072 | !Config Def = 1 |
---|
1073 | !Config Help = |
---|
1074 | ! |
---|
1075 | nsplit_thermals_omp = 1 |
---|
1076 | call getin('nsplit_thermals',nsplit_thermals_omp) |
---|
1077 | |
---|
1078 | !Config Key = alp_bl_k |
---|
1079 | !Config Desc = |
---|
1080 | !Config Def = 0. |
---|
1081 | !Config Help = |
---|
1082 | ! |
---|
1083 | alp_bl_k_omp = 1. |
---|
1084 | call getin('alp_bl_k',alp_bl_k_omp) |
---|
1085 | |
---|
1086 | ! |
---|
1087 | !Config Key = tau_thermals |
---|
1088 | !Config Desc = |
---|
1089 | !Config Def = 0. |
---|
1090 | !Config Help = |
---|
1091 | ! |
---|
1092 | tau_thermals_omp = 0. |
---|
1093 | call getin('tau_thermals',tau_thermals_omp) |
---|
1094 | |
---|
1095 | ! |
---|
1096 | !Config Key = iflag_coupl |
---|
1097 | !Config Desc = |
---|
1098 | !Config Def = 0 |
---|
1099 | !Config Help = |
---|
1100 | ! |
---|
1101 | iflag_coupl_omp = 0 |
---|
1102 | call getin('iflag_coupl',iflag_coupl_omp) |
---|
1103 | |
---|
1104 | ! |
---|
1105 | !Config Key = iflag_clos |
---|
1106 | !Config Desc = |
---|
1107 | !Config Def = 0 |
---|
1108 | !Config Help = |
---|
1109 | ! |
---|
1110 | iflag_clos_omp = 1 |
---|
1111 | call getin('iflag_clos',iflag_clos_omp) |
---|
1112 | ! |
---|
1113 | !Config Key = iflag_cvl_sigd |
---|
1114 | !Config Desc = |
---|
1115 | !Config Def = 0 |
---|
1116 | !Config Help = |
---|
1117 | ! |
---|
1118 | iflag_cvl_sigd_omp = 0 |
---|
1119 | call getin('iflag_cvl_sigd',iflag_cvl_sigd_omp) |
---|
1120 | |
---|
1121 | !Config Key = iflag_wake |
---|
1122 | !Config Desc = |
---|
1123 | !Config Def = 0 |
---|
1124 | !Config Help = |
---|
1125 | ! |
---|
1126 | iflag_wake_omp = 0 |
---|
1127 | call getin('iflag_wake',iflag_wake_omp) |
---|
1128 | |
---|
1129 | !Config Key = alp_offset |
---|
1130 | !Config Desc = |
---|
1131 | !Config Def = 0 |
---|
1132 | !Config Help = |
---|
1133 | ! |
---|
1134 | alp_offset_omp = 0. |
---|
1135 | call getin('alp_offset',alp_offset_omp) |
---|
1136 | |
---|
1137 | ! |
---|
1138 | !Config Key = lev_histhf |
---|
1139 | !Config Desc = |
---|
1140 | !Config Def = 1 |
---|
1141 | !Config Help = |
---|
1142 | ! |
---|
1143 | lev_histhf_omp = 1 |
---|
1144 | call getin('lev_histhf',lev_histhf_omp) |
---|
1145 | |
---|
1146 | ! |
---|
1147 | !Config Key = lev_histday |
---|
1148 | !Config Desc = |
---|
1149 | !Config Def = 1 |
---|
1150 | !Config Help = |
---|
1151 | ! |
---|
1152 | lev_histday_omp = 1 |
---|
1153 | call getin('lev_histday',lev_histday_omp) |
---|
1154 | |
---|
1155 | ! |
---|
1156 | !Config Key = lev_histmth |
---|
1157 | !Config Desc = |
---|
1158 | !Config Def = 2 |
---|
1159 | !Config Help = |
---|
1160 | ! |
---|
1161 | lev_histmth_omp = 2 |
---|
1162 | call getin('lev_histmth',lev_histmth_omp) |
---|
1163 | ! |
---|
1164 | !Config Key = lev_histins |
---|
1165 | !Config Desc = |
---|
1166 | !Config Def = 1 |
---|
1167 | !Config Help = |
---|
1168 | ! |
---|
1169 | lev_histins_omp = 1 |
---|
1170 | call getin('lev_histins',lev_histins_omp) |
---|
1171 | ! |
---|
1172 | !Config Key = lev_histLES |
---|
1173 | !Config Desc = |
---|
1174 | !Config Def = 1 |
---|
1175 | !Config Help = |
---|
1176 | ! |
---|
1177 | lev_histLES_omp = 1 |
---|
1178 | call getin('lev_histLES',lev_histLES_omp) |
---|
1179 | ! |
---|
1180 | !Config Key = lev_histdayNMC |
---|
1181 | !Config Desc = |
---|
1182 | !Config Def = 8 |
---|
1183 | !Config Help = |
---|
1184 | ! |
---|
1185 | lev_histdayNMC_omp = 8 |
---|
1186 | call getin('lev_histdayNMC',lev_histdayNMC_omp) |
---|
1187 | ! |
---|
1188 | !histNMC BEG |
---|
1189 | !Config Key = ok_histNMC |
---|
1190 | !Config Desc = ok_histNMC(1) = frequence de sortie fichiers histmthNMC |
---|
1191 | !Config Desc = ok_histNMC(2) = frequence de sortie fichiers histdayNMC |
---|
1192 | !Config Desc = ok_histNMC(3) = frequence de sortie fichiers histhfNMC |
---|
1193 | !Config Def = n, n, n |
---|
1194 | !Config Help = |
---|
1195 | ! |
---|
1196 | ok_histNMC_omp(1) = .false. |
---|
1197 | ok_histNMC_omp(2) = .false. |
---|
1198 | ok_histNMC_omp(3) = .false. |
---|
1199 | call getin('ok_histNMC',ok_histNMC_omp) |
---|
1200 | ! |
---|
1201 | !Config Key = freq_outNMC |
---|
1202 | !Config Desc = freq_outNMC(1) = frequence de sortie fichiers histmthNMC |
---|
1203 | !Config Desc = freq_outNMC(2) = frequence de sortie fichiers histdayNMC |
---|
1204 | !Config Desc = freq_outNMC(3) = frequence de sortie fichiers histhfNMC |
---|
1205 | !Config Def = 2592000., 86400., 21600. |
---|
1206 | !Config Help = |
---|
1207 | ! |
---|
1208 | ! freq_outNMC_omp(1) = 2592000. |
---|
1209 | freq_outNMC_omp(1) = mth_len*86400. |
---|
1210 | freq_outNMC_omp(2) = 86400. |
---|
1211 | freq_outNMC_omp(3) = 21600. |
---|
1212 | call getin('freq_outNMC',freq_outNMC_omp) |
---|
1213 | ! |
---|
1214 | !Config Key = freq_calNMC |
---|
1215 | !Config Desc = freq_calNMC(1) = frequence de calcul fichiers histmthNMC |
---|
1216 | !Config Desc = freq_calNMC(2) = frequence de calcul fichiers histdayNMC |
---|
1217 | !Config Desc = freq_calNMC(3) = frequence de calcul fichiers histhfNMC |
---|
1218 | !Config Def = pasphys |
---|
1219 | !Config Help = |
---|
1220 | ! |
---|
1221 | freq_calNMC_omp(1) = pasphys |
---|
1222 | freq_calNMC_omp(2) = pasphys |
---|
1223 | freq_calNMC_omp(3) = pasphys |
---|
1224 | call getin('freq_calNMC',freq_calNMC_omp) |
---|
1225 | ! |
---|
1226 | !Config Key = type_run |
---|
1227 | !Config Desc = |
---|
1228 | !Config Def = 'AMIP'/'CFMIP' ou 'CLIM'/'ENSP' |
---|
1229 | !Config Help = |
---|
1230 | ! |
---|
1231 | type_run_omp = 'AMIP' |
---|
1232 | call getin('type_run',type_run_omp) |
---|
1233 | |
---|
1234 | ! |
---|
1235 | !Config Key = ok_isccp |
---|
1236 | !Config Desc = |
---|
1237 | !Config Def = .true. |
---|
1238 | !Config Help = |
---|
1239 | ! |
---|
1240 | ! ok_isccp = .true. |
---|
1241 | ok_isccp_omp = .false. |
---|
1242 | call getin('ok_isccp',ok_isccp_omp) |
---|
1243 | |
---|
1244 | ! |
---|
1245 | !Config Key = ok_cosp |
---|
1246 | !Config Desc = |
---|
1247 | !Config Def = .false. |
---|
1248 | !Config Help = |
---|
1249 | ! |
---|
1250 | ok_cosp_omp = .false. |
---|
1251 | call getin('ok_cosp',ok_cosp_omp) |
---|
1252 | |
---|
1253 | ! |
---|
1254 | !Config Key = ok_mensuelCOSP |
---|
1255 | !Config Desc = |
---|
1256 | !Config Def = .true. |
---|
1257 | !Config Help = |
---|
1258 | ! |
---|
1259 | ok_mensuelCOSP_omp = .true. |
---|
1260 | call getin('ok_mensuelCOSP',ok_mensuelCOSP_omp) |
---|
1261 | |
---|
1262 | ! |
---|
1263 | !Config Key = ok_journeCOSP |
---|
1264 | !Config Desc = |
---|
1265 | !Config Def = .true. |
---|
1266 | !Config Help = |
---|
1267 | ! |
---|
1268 | ok_journeCOSP_omp = .true. |
---|
1269 | call getin('ok_journeCOSP',ok_journeCOSP_omp) |
---|
1270 | |
---|
1271 | ! |
---|
1272 | !Config Key = ok_hfCOSP |
---|
1273 | !Config Desc = |
---|
1274 | !Config Def = .false. |
---|
1275 | !Config Help = |
---|
1276 | ! |
---|
1277 | ok_hfCOSP_omp = .false. |
---|
1278 | call getin('ok_hfCOSP',ok_hfCOSP_omp) |
---|
1279 | |
---|
1280 | ! |
---|
1281 | ! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone |
---|
1282 | ! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc" |
---|
1283 | ! |
---|
1284 | !Config Key = lonmin_ins |
---|
1285 | !Config Desc = 100. |
---|
1286 | !Config Def = longitude minimale sorties "bilKP_ins" |
---|
1287 | !Config Help = |
---|
1288 | ! |
---|
1289 | lonmin_ins_omp = 100. |
---|
1290 | call getin('lonmin_ins',lonmin_ins_omp) |
---|
1291 | ! |
---|
1292 | !Config Key = lonmax_ins |
---|
1293 | !Config Desc = 130. |
---|
1294 | !Config Def = longitude maximale sorties "bilKP_ins" |
---|
1295 | !Config Help = |
---|
1296 | ! |
---|
1297 | lonmax_ins_omp = 130. |
---|
1298 | call getin('lonmax_ins',lonmax_ins_omp) |
---|
1299 | ! |
---|
1300 | !Config Key = latmin_ins |
---|
1301 | !Config Desc = -20. |
---|
1302 | !Config Def = latitude minimale sorties "bilKP_ins" |
---|
1303 | !Config Help = |
---|
1304 | ! |
---|
1305 | latmin_ins_omp = -20. |
---|
1306 | call getin('latmin_ins',latmin_ins_omp) |
---|
1307 | ! |
---|
1308 | !Config Key = latmax_ins |
---|
1309 | !Config Desc = 20. |
---|
1310 | !Config Def = latitude maximale sorties "bilKP_ins" |
---|
1311 | !Config Help = |
---|
1312 | ! |
---|
1313 | latmax_ins_omp = 20. |
---|
1314 | call getin('latmax_ins',latmax_ins_omp) |
---|
1315 | ! |
---|
1316 | !Config Key = ecrit_hf |
---|
1317 | !Config Desc = |
---|
1318 | !Config Def = 1./8. !toutes les 3h |
---|
1319 | !Config Help = |
---|
1320 | ! |
---|
1321 | ecrit_hf_omp = 1./8. |
---|
1322 | call getin('ecrit_hf',ecrit_hf_omp) |
---|
1323 | ! |
---|
1324 | !Config Key = ecrit_ins |
---|
1325 | !Config Desc = |
---|
1326 | !Config Def = 1./48. ! toutes les 1/2 h |
---|
1327 | !Config Help = |
---|
1328 | ! |
---|
1329 | ecrit_ins_omp = 1./48. |
---|
1330 | call getin('ecrit_ins',ecrit_ins_omp) |
---|
1331 | ! |
---|
1332 | !Config Key = ecrit_day |
---|
1333 | !Config Desc = |
---|
1334 | !Config Def = 1.0 !tous les jours |
---|
1335 | !Config Help = nombre de jours pour ecriture fichier histday.nc |
---|
1336 | ! |
---|
1337 | ecrit_day_omp = 1.0 |
---|
1338 | call getin('ecrit_day',ecrit_day_omp) |
---|
1339 | ! |
---|
1340 | !Config Key = ecrit_mth |
---|
1341 | !Config Desc = |
---|
1342 | !Config Def = 30. !tous les 30jours (1 fois par mois) |
---|
1343 | !Config Help = |
---|
1344 | ! |
---|
1345 | ecrit_mth_omp = 30. |
---|
1346 | call getin('ecrit_mth',ecrit_mth_omp) |
---|
1347 | ! |
---|
1348 | !Config Key = ecrit_tra |
---|
1349 | !Config Desc = |
---|
1350 | !Config Def = 30. !tous les 30jours (1 fois par mois) |
---|
1351 | !Config Help = |
---|
1352 | ! |
---|
1353 | ecrit_tra_omp = 30. |
---|
1354 | call getin('ecrit_tra',ecrit_tra_omp) |
---|
1355 | ! |
---|
1356 | !Config Key = ecrit_reg |
---|
1357 | !Config Desc = |
---|
1358 | !Config Def = 0.25 !4 fois par jour |
---|
1359 | !Config Help = |
---|
1360 | ! |
---|
1361 | ecrit_reg_omp = 0.25 !4 fois par jour |
---|
1362 | call getin('ecrit_reg',ecrit_reg_omp) |
---|
1363 | ! |
---|
1364 | ! |
---|
1365 | ! |
---|
1366 | ! PARAMETRES CDRAG |
---|
1367 | ! |
---|
1368 | !Config Key = f_cdrag_ter |
---|
1369 | !Config Desc = |
---|
1370 | !Config Def = 0.8 |
---|
1371 | !Config Help = |
---|
1372 | ! |
---|
1373 | f_cdrag_ter_omp = 0.8 |
---|
1374 | call getin('f_cdrag_ter',f_cdrag_ter_omp) |
---|
1375 | ! |
---|
1376 | !Config Key = f_cdrag_oce |
---|
1377 | !Config Desc = |
---|
1378 | !Config Def = 0.8 |
---|
1379 | !Config Help = |
---|
1380 | ! |
---|
1381 | f_cdrag_oce_omp = 0.8 |
---|
1382 | call getin('f_cdrag_oce',f_cdrag_oce_omp) |
---|
1383 | ! |
---|
1384 | ! RUGORO |
---|
1385 | !Config Key = f_rugoro |
---|
1386 | !Config Desc = |
---|
1387 | !Config Def = 0. |
---|
1388 | !Config Help = |
---|
1389 | ! |
---|
1390 | f_rugoro_omp = 0. |
---|
1391 | call getin('f_rugoro',f_rugoro_omp) |
---|
1392 | |
---|
1393 | ! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS |
---|
1394 | ! |
---|
1395 | !Config Key = supcrit1 |
---|
1396 | !Config Desc = |
---|
1397 | !Config Def = .540 |
---|
1398 | !Config Help = |
---|
1399 | ! |
---|
1400 | supcrit1_omp = .540 |
---|
1401 | call getin('supcrit1',supcrit1_omp) |
---|
1402 | |
---|
1403 | ! |
---|
1404 | !Config Key = supcrit2 |
---|
1405 | !Config Desc = |
---|
1406 | !Config Def = .600 |
---|
1407 | !Config Help = |
---|
1408 | ! |
---|
1409 | supcrit2_omp = .600 |
---|
1410 | call getin('supcrit2',supcrit2_omp) |
---|
1411 | |
---|
1412 | ! |
---|
1413 | ! PARAMETERS FOR THE MIXING DISTRIBUTION |
---|
1414 | ! iflag_mix: 0=OLD, |
---|
1415 | ! 1=NEW (JYG), |
---|
1416 | ! 2=NEW + conv. depth inhib. by tropos. dryness |
---|
1417 | ! '2' is NOT operationnal and should not be used. |
---|
1418 | ! |
---|
1419 | !Config Key = iflag_mix |
---|
1420 | !Config Desc = |
---|
1421 | !Config Def = 1 |
---|
1422 | !Config Help = |
---|
1423 | ! |
---|
1424 | iflag_mix_omp = 1 |
---|
1425 | call getin('iflag_mix',iflag_mix_omp) |
---|
1426 | |
---|
1427 | ! |
---|
1428 | !Config Key = scut |
---|
1429 | !Config Desc = |
---|
1430 | !Config Def = 0.95 |
---|
1431 | !Config Help = |
---|
1432 | ! |
---|
1433 | scut_omp = 0.95 |
---|
1434 | call getin('scut',scut_omp) |
---|
1435 | |
---|
1436 | ! |
---|
1437 | !Config Key = qqa1 |
---|
1438 | !Config Desc = |
---|
1439 | !Config Def = 1.0 |
---|
1440 | !Config Help = |
---|
1441 | ! |
---|
1442 | qqa1_omp = 1.0 |
---|
1443 | call getin('qqa1',qqa1_omp) |
---|
1444 | |
---|
1445 | ! |
---|
1446 | !Config Key = qqa2 |
---|
1447 | !Config Desc = |
---|
1448 | !Config Def = 0.0 |
---|
1449 | !Config Help = |
---|
1450 | ! |
---|
1451 | qqa2_omp = 0.0 |
---|
1452 | call getin('qqa2',qqa2_omp) |
---|
1453 | |
---|
1454 | ! |
---|
1455 | !Config Key = gammas |
---|
1456 | !Config Desc = |
---|
1457 | !Config Def = 0.05 |
---|
1458 | !Config Help = |
---|
1459 | ! |
---|
1460 | gammas_omp = 0.05 |
---|
1461 | call getin('gammas',gammas_omp) |
---|
1462 | |
---|
1463 | ! |
---|
1464 | !Config Key = Fmax |
---|
1465 | !Config Desc = |
---|
1466 | !Config Def = 0.65 |
---|
1467 | !Config Help = |
---|
1468 | ! |
---|
1469 | Fmax_omp = 0.65 |
---|
1470 | call getin('Fmax',Fmax_omp) |
---|
1471 | |
---|
1472 | ! |
---|
1473 | !Config Key = alphas |
---|
1474 | !Config Desc = |
---|
1475 | !Config Def = -5. |
---|
1476 | !Config Help = |
---|
1477 | ! |
---|
1478 | alphas_omp = -5. |
---|
1479 | call getin('alphas',alphas_omp) |
---|
1480 | |
---|
1481 | !Config key = ok_strato |
---|
1482 | !Config Desc = activation de la version strato |
---|
1483 | !Config Def = .FALSE. |
---|
1484 | !Config Help = active la version stratosphérique de LMDZ de F. Lott |
---|
1485 | |
---|
1486 | ok_strato_omp=.FALSE. |
---|
1487 | CALL getin('ok_strato',ok_strato_omp) |
---|
1488 | |
---|
1489 | !Config key = ok_hines |
---|
1490 | !Config Desc = activation de la parametrisation de hines |
---|
1491 | !Config Def = .FALSE. |
---|
1492 | !Config Help = Clefs controlant la parametrization de Hines |
---|
1493 | ! Et la sponge layer (Runs Stratospheriques) |
---|
1494 | |
---|
1495 | ok_hines_omp=.FALSE. |
---|
1496 | CALL getin('ok_hines',ok_hines_omp) |
---|
1497 | |
---|
1498 | !Config Key = OK_LES |
---|
1499 | !Config Desc = Pour des sorties LES |
---|
1500 | !Config Def = .false. |
---|
1501 | !Config Help = Pour creer le fichier histLES contenant les sorties |
---|
1502 | ! LES |
---|
1503 | ! |
---|
1504 | ok_LES_omp = .false. |
---|
1505 | call getin('OK_LES', ok_LES_omp) |
---|
1506 | |
---|
1507 | !Config Key = callstats |
---|
1508 | !Config Desc = Pour des sorties callstats |
---|
1509 | !Config Def = .false. |
---|
1510 | !Config Help = Pour creer le fichier stats contenant les sorties |
---|
1511 | ! stats |
---|
1512 | ! |
---|
1513 | callstats_omp = .false. |
---|
1514 | call getin('callstats', callstats_omp) |
---|
1515 | ! |
---|
1516 | !Config Key = ecrit_LES |
---|
1517 | !Config Desc = Frequence d'ecriture des resultats du LES en nombre de jours; |
---|
1518 | ! par defaut 1., i.e. 1 jour |
---|
1519 | !Config Def = 1./8. |
---|
1520 | !Config Help = ... |
---|
1521 | ! |
---|
1522 | ! |
---|
1523 | ecrit_LES_omp = 1./8. |
---|
1524 | call getin('ecrit_LES', ecrit_LES_omp) |
---|
1525 | ! |
---|
1526 | read_climoz = 0 ! default value |
---|
1527 | call getin('read_climoz', read_climoz) |
---|
1528 | |
---|
1529 | carbon_cycle_tr_omp=.FALSE. |
---|
1530 | CALL getin('carbon_cycle_tr',carbon_cycle_tr_omp) |
---|
1531 | |
---|
1532 | carbon_cycle_cpl_omp=.FALSE. |
---|
1533 | CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp) |
---|
1534 | |
---|
1535 | !$OMP END MASTER |
---|
1536 | !$OMP BARRIER |
---|
1537 | |
---|
1538 | R_ecc = R_ecc_omp |
---|
1539 | R_peri = R_peri_omp |
---|
1540 | R_incl = R_incl_omp |
---|
1541 | solaire = solaire_omp |
---|
1542 | co2_ppm = co2_ppm_omp |
---|
1543 | RCO2 = RCO2_omp |
---|
1544 | CH4_ppb = CH4_ppb_omp |
---|
1545 | RCH4 = RCH4_omp |
---|
1546 | N2O_ppb = N2O_ppb_omp |
---|
1547 | RN2O = RN2O_omp |
---|
1548 | CFC11_ppt = CFC11_ppt_omp |
---|
1549 | RCFC11 = RCFC11_omp |
---|
1550 | CFC12_ppt = CFC12_ppt_omp |
---|
1551 | RCFC12 = RCFC12_omp |
---|
1552 | RCO2_act = RCO2 |
---|
1553 | RCH4_act = RCH4 |
---|
1554 | RN2O_act = RN2O |
---|
1555 | RCFC11_act = RCFC11 |
---|
1556 | RCFC12_act = RCFC12 |
---|
1557 | RCO2_per = RCO2_per_omp |
---|
1558 | RCH4_per = RCH4_per_omp |
---|
1559 | RN2O_per = RN2O_per_omp |
---|
1560 | RCFC11_per = RCFC11_per_omp |
---|
1561 | RCFC12_per = RCFC12_per_omp |
---|
1562 | |
---|
1563 | cycle_diurne = cycle_diurne_omp |
---|
1564 | soil_model = soil_model_omp |
---|
1565 | new_oliq = new_oliq_omp |
---|
1566 | ok_orodr = ok_orodr_omp |
---|
1567 | ok_orolf = ok_orolf_omp |
---|
1568 | ok_limitvrai = ok_limitvrai_omp |
---|
1569 | nbapp_rad = nbapp_rad_omp |
---|
1570 | iflag_con = iflag_con_omp |
---|
1571 | |
---|
1572 | epmax = epmax_omp |
---|
1573 | ok_adj_ema = ok_adj_ema_omp |
---|
1574 | iflag_clw = iflag_clw_omp |
---|
1575 | cld_lc_lsc = cld_lc_lsc_omp |
---|
1576 | cld_lc_con = cld_lc_con_omp |
---|
1577 | cld_tau_lsc = cld_tau_lsc_omp |
---|
1578 | cld_tau_con = cld_tau_con_omp |
---|
1579 | ffallv_lsc = ffallv_lsc_omp |
---|
1580 | ffallv_con = ffallv_con_omp |
---|
1581 | coef_eva = coef_eva_omp |
---|
1582 | reevap_ice = reevap_ice_omp |
---|
1583 | iflag_pdf = iflag_pdf_omp |
---|
1584 | solarlong0 = solarlong0_omp |
---|
1585 | qsol0 = qsol0_omp |
---|
1586 | inertie_sol = inertie_sol_omp |
---|
1587 | inertie_ice = inertie_ice_omp |
---|
1588 | inertie_sno = inertie_sno_omp |
---|
1589 | rad_froid = rad_froid_omp |
---|
1590 | rad_chau1 = rad_chau1_omp |
---|
1591 | rad_chau2 = rad_chau2_omp |
---|
1592 | t_glace_min = t_glace_min_omp |
---|
1593 | t_glace_max = t_glace_max_omp |
---|
1594 | rei_min = rei_min_omp |
---|
1595 | rei_max = rei_max_omp |
---|
1596 | top_height = top_height_omp |
---|
1597 | overlap = overlap_omp |
---|
1598 | cdmmax = cdmmax_omp |
---|
1599 | cdhmax = cdhmax_omp |
---|
1600 | ksta = ksta_omp |
---|
1601 | ksta_ter = ksta_ter_omp |
---|
1602 | ok_kzmin = ok_kzmin_omp |
---|
1603 | fmagic = fmagic_omp |
---|
1604 | pmagic = pmagic_omp |
---|
1605 | iflag_pbl = iflag_pbl_omp |
---|
1606 | lev_histhf = lev_histhf_omp |
---|
1607 | lev_histday = lev_histday_omp |
---|
1608 | lev_histmth = lev_histmth_omp |
---|
1609 | lev_histins = lev_histins_omp |
---|
1610 | lev_histLES = lev_histLES_omp |
---|
1611 | lev_histdayNMC = lev_histdayNMC_omp |
---|
1612 | ok_histNMC(:) = ok_histNMC_omp(:) |
---|
1613 | freq_outNMC(:) = freq_outNMC_omp(:) |
---|
1614 | freq_calNMC(:) = freq_calNMC_omp(:) |
---|
1615 | |
---|
1616 | type_ocean = type_ocean_omp |
---|
1617 | version_ocean = version_ocean_omp |
---|
1618 | ok_veget = ok_veget_omp |
---|
1619 | ok_newmicro = ok_newmicro_omp |
---|
1620 | ok_journe = ok_journe_omp |
---|
1621 | ok_hf = ok_hf_omp |
---|
1622 | ok_mensuel = ok_mensuel_omp |
---|
1623 | ok_instan = ok_instan_omp |
---|
1624 | freq_ISCCP = freq_ISCCP_omp |
---|
1625 | ecrit_ISCCP = ecrit_ISCCP_omp |
---|
1626 | freq_COSP = freq_COSP_omp |
---|
1627 | ok_ade = ok_ade_omp |
---|
1628 | ok_aie = ok_aie_omp |
---|
1629 | aerosol_couple = aerosol_couple_omp |
---|
1630 | flag_aerosol=flag_aerosol_omp |
---|
1631 | new_aod=new_aod_omp |
---|
1632 | aer_type = aer_type_omp |
---|
1633 | bl95_b0 = bl95_b0_omp |
---|
1634 | bl95_b1 = bl95_b1_omp |
---|
1635 | fact_cldcon = fact_cldcon_omp |
---|
1636 | facttemps = facttemps_omp |
---|
1637 | ratqsbas = ratqsbas_omp |
---|
1638 | ratqshaut = ratqshaut_omp |
---|
1639 | tau_ratqs = tau_ratqs_omp |
---|
1640 | |
---|
1641 | iflag_radia = iflag_radia_omp |
---|
1642 | iflag_rrtm = iflag_rrtm_omp |
---|
1643 | iflag_cldcon = iflag_cldcon_omp |
---|
1644 | iflag_ratqs = iflag_ratqs_omp |
---|
1645 | ip_ebil_phy = ip_ebil_phy_omp |
---|
1646 | iflag_thermals = iflag_thermals_omp |
---|
1647 | iflag_thermals_ed = iflag_thermals_ed_omp |
---|
1648 | iflag_thermals_optflux = iflag_thermals_optflux_omp |
---|
1649 | nsplit_thermals = nsplit_thermals_omp |
---|
1650 | tau_thermals = tau_thermals_omp |
---|
1651 | alp_bl_k = alp_bl_k_omp |
---|
1652 | iflag_coupl = iflag_coupl_omp |
---|
1653 | iflag_clos = iflag_clos_omp |
---|
1654 | iflag_wake = iflag_wake_omp |
---|
1655 | alp_offset = alp_offset_omp |
---|
1656 | iflag_cvl_sigd = iflag_cvl_sigd_omp |
---|
1657 | type_run = type_run_omp |
---|
1658 | ok_isccp = ok_isccp_omp |
---|
1659 | ok_cosp = ok_cosp_omp |
---|
1660 | ok_mensuelCOSP = ok_mensuelCOSP_omp |
---|
1661 | ok_journeCOSP = ok_journeCOSP_omp |
---|
1662 | ok_hfCOSP = ok_hfCOSP_omp |
---|
1663 | seuil_inversion=seuil_inversion_omp |
---|
1664 | lonmin_ins = lonmin_ins_omp |
---|
1665 | lonmax_ins = lonmax_ins_omp |
---|
1666 | latmin_ins = latmin_ins_omp |
---|
1667 | latmax_ins = latmax_ins_omp |
---|
1668 | ecrit_hf = ecrit_hf_omp |
---|
1669 | ecrit_ins = ecrit_ins_omp |
---|
1670 | ecrit_day = ecrit_day_omp |
---|
1671 | ecrit_mth = ecrit_mth_omp |
---|
1672 | ecrit_tra = ecrit_tra_omp |
---|
1673 | ecrit_reg = ecrit_reg_omp |
---|
1674 | cvl_corr = cvl_corr_omp |
---|
1675 | ok_lic_melt = ok_lic_melt_omp |
---|
1676 | f_cdrag_ter=f_cdrag_ter_omp |
---|
1677 | f_cdrag_oce=f_cdrag_oce_omp |
---|
1678 | f_rugoro=f_rugoro_omp |
---|
1679 | supcrit1 = supcrit1_omp |
---|
1680 | supcrit2 = supcrit2_omp |
---|
1681 | iflag_mix = iflag_mix_omp |
---|
1682 | scut = scut_omp |
---|
1683 | qqa1 = qqa1_omp |
---|
1684 | qqa2 = qqa2_omp |
---|
1685 | gammas = gammas_omp |
---|
1686 | Fmax = Fmax_omp |
---|
1687 | alphas = alphas_omp |
---|
1688 | ok_strato = ok_strato_omp |
---|
1689 | ok_hines = ok_hines_omp |
---|
1690 | ok_LES = ok_LES_omp |
---|
1691 | callstats = callstats_omp |
---|
1692 | ecrit_LES = ecrit_LES_omp |
---|
1693 | carbon_cycle_tr = carbon_cycle_tr_omp |
---|
1694 | carbon_cycle_cpl = carbon_cycle_cpl_omp |
---|
1695 | |
---|
1696 | ! Test of coherence between type_ocean and version_ocean |
---|
1697 | IF (type_ocean=='couple' .AND. (version_ocean/='opa8' .AND. version_ocean/='nemo') ) THEN |
---|
1698 | WRITE(numout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration' |
---|
1699 | CALL abort_gcm('conf_phys','version_ocean not valid',1) |
---|
1700 | END IF |
---|
1701 | |
---|
1702 | IF (type_ocean=='slab' .AND. version_ocean=='xxxxxx') THEN |
---|
1703 | version_ocean='sicOBS' |
---|
1704 | ELSE IF (type_ocean=='slab' .AND. version_ocean/='sicOBS') THEN |
---|
1705 | WRITE(numout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean' |
---|
1706 | CALL abort_gcm('conf_phys','version_ocean not valid',1) |
---|
1707 | END IF |
---|
1708 | |
---|
1709 | ! Test sur new_aod. Ce flag permet de retrouver les resultats de l'AR4 |
---|
1710 | ! il n'est utilisable que lors du couplage avec le SO4 seul |
---|
1711 | IF (ok_ade .OR. ok_aie) THEN |
---|
1712 | IF ( .NOT. new_aod .AND. flag_aerosol .NE. 1) THEN |
---|
1713 | CALL abort_gcm('conf_phys','new_aod=.FALSE. not compatible avec flag_aerosol=1',1) |
---|
1714 | END IF |
---|
1715 | END IF |
---|
1716 | |
---|
1717 | !$OMP MASTER |
---|
1718 | |
---|
1719 | write(numout,*)' ##############################################' |
---|
1720 | write(numout,*)' Configuration des parametres de la physique: ' |
---|
1721 | write(numout,*)' Type ocean = ', type_ocean |
---|
1722 | write(numout,*)' Version ocean = ', version_ocean |
---|
1723 | write(numout,*)' Config veget = ', ok_veget |
---|
1724 | write(numout,*)' Sortie journaliere = ', ok_journe |
---|
1725 | write(numout,*)' Sortie haute frequence = ', ok_hf |
---|
1726 | write(numout,*)' Sortie mensuelle = ', ok_mensuel |
---|
1727 | write(numout,*)' Sortie instantanee = ', ok_instan |
---|
1728 | write(numout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP |
---|
1729 | write(numout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP |
---|
1730 | write(numout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP |
---|
1731 | write(numout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy |
---|
1732 | write(numout,*)' Excentricite = ',R_ecc |
---|
1733 | write(numout,*)' Equinoxe = ',R_peri |
---|
1734 | write(numout,*)' Inclinaison =',R_incl |
---|
1735 | write(numout,*)' Constante solaire =',solaire |
---|
1736 | write(numout,*)' co2_ppm =',co2_ppm |
---|
1737 | write(numout,*)' RCO2_act = ',RCO2_act |
---|
1738 | write(numout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act |
---|
1739 | write(numout,*)' N2O_ppb =',N2O_ppb,' RN2O_act= ',RN2O_act |
---|
1740 | write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act= ',RCFC11_act |
---|
1741 | write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act= ',RCFC12_act |
---|
1742 | write(numout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per |
---|
1743 | write(numout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per |
---|
1744 | write(numout,*)' RCFC12_per = ',RCFC12_per |
---|
1745 | write(numout,*)' cvl_corr=', cvl_corr |
---|
1746 | write(numout,*)'ok_lic_melt=', ok_lic_melt |
---|
1747 | write(numout,*)'cycle_diurne=',cycle_diurne |
---|
1748 | write(numout,*)'soil_model=',soil_model |
---|
1749 | write(numout,*)'new_oliq=',new_oliq |
---|
1750 | write(numout,*)'ok_orodr=',ok_orodr |
---|
1751 | write(numout,*)'ok_orolf=',ok_orolf |
---|
1752 | write(numout,*)'ok_limitvrai=',ok_limitvrai |
---|
1753 | write(numout,*)'nbapp_rad=',nbapp_rad |
---|
1754 | write(numout,*)'iflag_con=',iflag_con |
---|
1755 | write(numout,*)' epmax = ', epmax |
---|
1756 | write(numout,*)' ok_adj_ema = ', ok_adj_ema |
---|
1757 | write(numout,*)' iflag_clw = ', iflag_clw |
---|
1758 | write(numout,*)' cld_lc_lsc = ', cld_lc_lsc |
---|
1759 | write(numout,*)' cld_lc_con = ', cld_lc_con |
---|
1760 | write(numout,*)' cld_tau_lsc = ', cld_tau_lsc |
---|
1761 | write(numout,*)' cld_tau_con = ', cld_tau_con |
---|
1762 | write(numout,*)' ffallv_lsc = ', ffallv_lsc |
---|
1763 | write(numout,*)' ffallv_con = ', ffallv_con |
---|
1764 | write(numout,*)' coef_eva = ', coef_eva |
---|
1765 | write(numout,*)' reevap_ice = ', reevap_ice |
---|
1766 | write(numout,*)' iflag_pdf = ', iflag_pdf |
---|
1767 | write(numout,*)' iflag_cldcon = ', iflag_cldcon |
---|
1768 | write(numout,*)' iflag_radia = ', iflag_radia |
---|
1769 | write(numout,*)' iflag_rrtm = ', iflag_rrtm |
---|
1770 | write(numout,*)' iflag_ratqs = ', iflag_ratqs |
---|
1771 | write(numout,*)' seuil_inversion = ', seuil_inversion |
---|
1772 | write(numout,*)' fact_cldcon = ', fact_cldcon |
---|
1773 | write(numout,*)' facttemps = ', facttemps |
---|
1774 | write(numout,*)' ok_newmicro = ',ok_newmicro |
---|
1775 | write(numout,*)' ratqsbas = ',ratqsbas |
---|
1776 | write(numout,*)' ratqshaut = ',ratqshaut |
---|
1777 | write(numout,*)' tau_ratqs = ',tau_ratqs |
---|
1778 | write(numout,*)' top_height = ',top_height |
---|
1779 | write(numout,*)' rad_froid = ',rad_froid |
---|
1780 | write(numout,*)' rad_chau1 = ',rad_chau1 |
---|
1781 | write(numout,*)' rad_chau2 = ',rad_chau2 |
---|
1782 | write(numout,*)' t_glace_min = ',t_glace_min |
---|
1783 | write(numout,*)' t_glace_max = ',t_glace_max |
---|
1784 | write(numout,*)' rei_min = ',rei_min |
---|
1785 | write(numout,*)' rei_max = ',rei_max |
---|
1786 | write(numout,*)' overlap = ',overlap |
---|
1787 | write(numout,*)' cdmmax = ',cdmmax |
---|
1788 | write(numout,*)' cdhmax = ',cdhmax |
---|
1789 | write(numout,*)' ksta = ',ksta |
---|
1790 | write(numout,*)' ksta_ter = ',ksta_ter |
---|
1791 | write(numout,*)' ok_kzmin = ',ok_kzmin |
---|
1792 | write(numout,*)' fmagic = ',fmagic |
---|
1793 | write(numout,*)' pmagic = ',pmagic |
---|
1794 | write(numout,*)' ok_ade = ',ok_ade |
---|
1795 | write(numout,*)' ok_aie = ',ok_aie |
---|
1796 | write(numout,*)' aerosol_couple = ', aerosol_couple |
---|
1797 | write(numout,*)' flag_aerosol = ', flag_aerosol |
---|
1798 | write(numout,*)' new_aod = ', new_aod |
---|
1799 | write(numout,*)' aer_type = ',aer_type |
---|
1800 | write(numout,*)' bl95_b0 = ',bl95_b0 |
---|
1801 | write(numout,*)' bl95_b1 = ',bl95_b1 |
---|
1802 | write(numout,*)' lev_histhf = ',lev_histhf |
---|
1803 | write(numout,*)' lev_histday = ',lev_histday |
---|
1804 | write(numout,*)' lev_histmth = ',lev_histmth |
---|
1805 | write(numout,*)' lev_histins = ',lev_histins |
---|
1806 | write(numout,*)' lev_histLES = ',lev_histLES |
---|
1807 | write(numout,*)' lev_histdayNMC = ',lev_histdayNMC |
---|
1808 | write(numout,*)' ok_histNMC = ',ok_histNMC |
---|
1809 | write(numout,*)' freq_outNMC = ',freq_outNMC |
---|
1810 | write(numout,*)' freq_calNMC = ',freq_calNMC |
---|
1811 | write(numout,*)' iflag_pbl = ', iflag_pbl |
---|
1812 | write(numout,*)' iflag_thermals = ', iflag_thermals |
---|
1813 | write(numout,*)' iflag_thermals_ed = ', iflag_thermals_ed |
---|
1814 | write(numout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux |
---|
1815 | write(numout,*)' iflag_clos = ', iflag_clos |
---|
1816 | write(numout,*)' type_run = ',type_run |
---|
1817 | write(numout,*)' ok_isccp = ',ok_isccp |
---|
1818 | write(numout,*)' ok_cosp = ',ok_cosp |
---|
1819 | write(numout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP |
---|
1820 | write(numout,*)' ok_journeCOSP = ',ok_journeCOSP |
---|
1821 | write(numout,*)' ok_hfCOSP =',ok_hfCOSP |
---|
1822 | write(numout,*)' solarlong0 = ', solarlong0 |
---|
1823 | write(numout,*)' qsol0 = ', qsol0 |
---|
1824 | write(numout,*)' inertie_sol = ', inertie_sol |
---|
1825 | write(numout,*)' inertie_ice = ', inertie_ice |
---|
1826 | write(numout,*)' inertie_sno = ', inertie_sno |
---|
1827 | write(numout,*)' f_cdrag_ter = ',f_cdrag_ter |
---|
1828 | write(numout,*)' f_cdrag_oce = ',f_cdrag_oce |
---|
1829 | write(numout,*)' f_rugoro = ',f_rugoro |
---|
1830 | write(numout,*)' supcrit1 = ', supcrit1 |
---|
1831 | write(numout,*)' supcrit2 = ', supcrit2 |
---|
1832 | write(numout,*)' iflag_mix = ', iflag_mix |
---|
1833 | write(numout,*)' scut = ', scut |
---|
1834 | write(numout,*)' qqa1 = ', qqa1 |
---|
1835 | write(numout,*)' qqa2 = ', qqa2 |
---|
1836 | write(numout,*)' gammas = ', gammas |
---|
1837 | write(numout,*)' Fmax = ', Fmax |
---|
1838 | write(numout,*)' alphas = ', alphas |
---|
1839 | write(numout,*)' iflag_wake = ', iflag_wake |
---|
1840 | write(numout,*)' alp_offset = ', alp_offset |
---|
1841 | |
---|
1842 | write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',& |
---|
1843 | lonmin_ins, lonmax_ins, latmin_ins, latmax_ins |
---|
1844 | write(numout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',& |
---|
1845 | ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES |
---|
1846 | |
---|
1847 | write(numout,*) 'ok_strato = ', ok_strato |
---|
1848 | write(numout,*) 'ok_hines = ', ok_hines |
---|
1849 | write(numout,*) 'read_climoz = ', read_climoz |
---|
1850 | write(numout,*) 'carbon_cycle_tr = ', carbon_cycle_tr |
---|
1851 | write(numout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl |
---|
1852 | |
---|
1853 | !$OMP END MASTER |
---|
1854 | |
---|
1855 | return |
---|
1856 | |
---|
1857 | end subroutine conf_phys |
---|
1858 | |
---|
1859 | end module conf_phys_m |
---|
1860 | ! |
---|
1861 | !################################################################# |
---|
1862 | ! |
---|
1863 | |
---|
1864 | subroutine conf_interface(tau_calv) |
---|
1865 | |
---|
1866 | use IOIPSL |
---|
1867 | implicit none |
---|
1868 | |
---|
1869 | ! Configuration de l'interace atm/surf |
---|
1870 | ! |
---|
1871 | ! tau_calv: temps de relaxation pour la fonte des glaciers |
---|
1872 | |
---|
1873 | REAL :: tau_calv |
---|
1874 | REAL,SAVE :: tau_calv_omp |
---|
1875 | |
---|
1876 | ! Local |
---|
1877 | integer :: numout = 6 |
---|
1878 | ! |
---|
1879 | !Config Key = tau_calv |
---|
1880 | !Config Desc = temps de relaxation pour fonte des glaciers en jours |
---|
1881 | !Config Def = 1 an |
---|
1882 | !Config Help = |
---|
1883 | ! |
---|
1884 | tau_calv_omp = 360.*10. |
---|
1885 | !$OMP MASTER |
---|
1886 | call getin('tau_calv',tau_calv_omp) |
---|
1887 | !$OMP END MASTER |
---|
1888 | !$OMP BARRIER |
---|
1889 | |
---|
1890 | tau_calv=tau_calv_omp |
---|
1891 | |
---|
1892 | !$OMP MASTER |
---|
1893 | write(numout,*)' ##############################################' |
---|
1894 | WRITE(numout,*)' Configuration de l''interface atm/surfaces : ' |
---|
1895 | WRITE(numout,*)' tau_calv = ',tau_calv |
---|
1896 | !$OMP END MASTER |
---|
1897 | |
---|
1898 | return |
---|
1899 | |
---|
1900 | end subroutine conf_interface |
---|