Context Navigation

conf_phys_m.F90 @ 2311

Last change on this file since 2311 was 2311, checked in by Ehouarn Millour, 9 years ago

Further modifications to enforce physics/dynamics separation:

moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to physics (should be used rather than iniprint.h)
created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead when in physics.
reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and print_control_mod flags.

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

File size: 68.3 KB

Line
1	!
2	! $Id: conf_phys.F90 1668 2012-10-12 10:47:37Z idelkadi $
3	!
4	!
5	!
6	module conf_phys_m
7
8	implicit none
9
10	contains
11
12	subroutine conf_phys(ok_journe, ok_mensuel, ok_instan, ok_hf, &
13	ok_LES,&
14	callstats,&
15	solarlong0,seuil_inversion, &
16	fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
17	iflag_cld_th, &
18	iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
19	ok_ade, ok_aie, ok_cdnc, aerosol_couple, &
20	flag_aerosol, flag_aerosol_strat, new_aod, &
21	bl95_b0, bl95_b1,&
22	read_climoz, &
23	alp_offset)
24
25	use IOIPSL
26	USE surface_data
27	USE phys_cal_mod
28	USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl
29	USE control_mod
30	USE mod_grid_phy_lmdz, only: klon_glo
31	USE print_control_mod, ONLY: lunout
32
33
34	include "conema3.h"
35	include "fisrtilp.h"
36	include "nuage.h"
37	include "YOMCST.h"
38	include "YOMCST2.h"
39
40	include "thermcell.h"
41
42
43	!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
44	include "clesphys.h"
45	include "compbl.h"
46	include "comsoil.h"
47	include "YOEGWD.h"
48	!
49	! Configuration de la "physique" de LMDZ a l'aide de la fonction
50	! GETIN de IOIPSL
51	!
52	! LF 05/2001
53	!
54
55	!
56	! type_ocean: type d'ocean (force, slab, couple)
57	! version_ocean: version d'ocean (opa8/nemo pour type_ocean=couple ou
58	! sicOBS,sicINT,sicNO pour type_ocean=slab)
59	! ok_veget: type de modele de vegetation
60	! ok_journe: sorties journalieres
61	! ok_hf: sorties haute frequence
62	! ok_mensuel: sorties mensuelles
63	! ok_instan: sorties instantanees
64	! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
65	! ok_cdnc, ok cloud droplet number concentration
66	! flag_aerosol_strat : flag pour les aerosols stratos
67	! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
68	!
69
70
71	! Sortie:
72	logical :: ok_newmicro
73	integer :: iflag_radia
74	logical :: ok_journe, ok_mensuel, ok_instan, ok_hf
75	logical :: ok_LES
76	LOGICAL :: callstats
77	LOGICAL :: ok_ade, ok_aie, ok_cdnc, aerosol_couple
78	INTEGER :: flag_aerosol
79	LOGICAL :: flag_aerosol_strat
80	LOGICAL :: new_aod
81	REAL :: bl95_b0, bl95_b1
82	real :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
83	integer :: iflag_cld_th
84	integer :: iflag_ratqs
85
86	character (len = 6),SAVE :: type_ocean_omp, version_ocean_omp, ocean_omp
87	character (len = 10),SAVE :: type_veget_omp
88	CHARACTER(len = 8),SAVE :: aer_type_omp
89	logical,SAVE :: ok_snow_omp
90	logical,SAVE :: ok_newmicro_omp
91	logical,SAVE :: ok_all_xml_omp
92	logical,SAVE :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp
93	logical,SAVE :: ok_LES_omp
94	LOGICAL,SAVE :: callstats_omp
95	LOGICAL,SAVE :: ok_ade_omp, ok_aie_omp, ok_cdnc_omp, aerosol_couple_omp
96	INTEGER, SAVE :: flag_aerosol_omp
97	LOGICAL, SAVE :: flag_aerosol_strat_omp
98	LOGICAL, SAVE :: new_aod_omp
99	REAL,SAVE :: bl95_b0_omp, bl95_b1_omp
100	REAL,SAVE :: freq_ISCCP_omp, ecrit_ISCCP_omp
101	REAL,SAVE :: freq_COSP_omp
102	real,SAVE :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
103	real,SAVE :: tau_cld_cv_omp, coefw_cld_cv_omp
104	integer,SAVE :: iflag_cld_cv_omp
105
106
107	real,SAVE :: ratqshaut_omp
108	real,SAVE :: tau_ratqs_omp
109	REAL, SAVE :: t_coupl_omp
110	integer,SAVE :: iflag_radia_omp
111	integer,SAVE :: iflag_rrtm_omp
112	integer,SAVE :: iflag_albedo_omp !albedo SB
113	logical,save :: ok_chlorophyll_omp ! albedo SB
114	integer,SAVE :: NSW_omp
115	integer,SAVE :: iflag_cld_th_omp, ip_ebil_phy_omp
116	integer,SAVE :: iflag_ratqs_omp
117
118	Real,SAVE :: f_cdrag_ter_omp,f_cdrag_oce_omp
119	Real,SAVE :: f_rugoro_omp , z0min_omp
120	Real,SAVE :: z0m_seaice_omp,z0h_seaice_omp
121	REAL,SAVE :: min_wind_speed_omp,f_gust_wk_omp,f_gust_bl_omp,f_qsat_oce_omp, f_z0qh_oce_omp
122	INTEGER,SAVE :: iflag_gusts_omp,iflag_z0_oce_omp
123
124	! Local
125	real :: zzz
126
127	real :: seuil_inversion
128	real,save :: seuil_inversion_omp
129
130	integer,SAVE :: iflag_thermals_ed_omp,iflag_thermals_optflux_omp,iflag_thermals_closure_omp
131	real, SAVE :: fact_thermals_ed_dz_omp
132	integer,SAVE :: iflag_thermals_omp,nsplit_thermals_omp
133	real,save :: tau_thermals_omp,alp_bl_k_omp
134	! nrlmd le 10/04/2012
135	integer,SAVE :: iflag_trig_bl_omp,iflag_clos_bl_omp
136	integer,SAVE :: tau_trig_shallow_omp,tau_trig_deep_omp
137	real,SAVE :: s_trig_omp
138	! fin nrlmd le 10/04/2012
139	real :: alp_offset
140	REAL, SAVE :: alp_offset_omp
141	integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
142	integer,SAVE :: iflag_cvl_sigd_omp
143	REAL, SAVE :: coef_clos_ls_omp
144	REAL, SAVE :: supcrit1_omp, supcrit2_omp
145	INTEGER, SAVE :: iflag_mix_omp
146	real, save :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp
147	REAL, SAVE :: tmax_fonte_cv_omp
148
149	REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
150	REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
151	REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
152	REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp
153	REAL,SAVE :: CFC11_ppt_omp,RCFC11_omp,CFC11_ppt_per_omp,RCFC11_per_omp
154	REAL,SAVE :: CFC12_ppt_omp,RCFC12_omp,CFC12_ppt_per_omp,RCFC12_per_omp
155	LOGICAL,SAVE :: ok_4xCO2atm_omp
156	REAL,SAVE :: epmax_omp
157	LOGICAL,SAVE :: ok_adj_ema_omp
158	INTEGER,SAVE :: iflag_clw_omp
159	REAL,SAVE :: cld_lc_lsc_omp,cld_lc_con_omp,cld_tau_lsc_omp,cld_tau_con_omp
160	REAL,SAVE :: ffallv_lsc_omp, ffallv_con_omp,coef_eva_omp
161	LOGICAL,SAVE :: reevap_ice_omp
162	INTEGER,SAVE :: iflag_pdf_omp
163	INTEGER,SAVE :: iflag_ice_thermo_omp
164	INTEGER,SAVE :: iflag_t_glace_omp
165	REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
166	REAL,SAVE :: t_glace_min_omp, t_glace_max_omp
167	REAL,SAVE :: exposant_glace_omp
168	REAL,SAVE :: rei_min_omp, rei_max_omp
169	REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
170	REAL,SAVE :: qsol0_omp
171	REAL,SAVE :: evap0_omp
172	REAL,SAVE :: albsno0_omp
173	REAL :: solarlong0
174	REAL,SAVE :: solarlong0_omp
175	INTEGER,SAVE :: top_height_omp,overlap_omp
176	REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp,f_ri_cd_min_omp
177	LOGICAL,SAVE :: ok_kzmin_omp
178	REAL, SAVE :: fmagic_omp, pmagic_omp
179	INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
180	INTEGER,SAVE :: iflag_pbl_split_omp
181	Integer, save :: lev_histins_omp, lev_histLES_omp
182	INTEGER, SAVE :: lev_histdayNMC_omp
183	INTEGER, SAVE :: levout_histNMC_omp(3)
184	LOGICAL, SAVE :: ok_histNMC_omp(3)
185	REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
186	CHARACTER*4, SAVE :: type_run_omp
187	LOGICAL,SAVE :: ok_cosp_omp
188	LOGICAL,SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
189	REAL,SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
190	REAL,SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
191	REAL,SAVE :: ecrit_ins_omp
192	REAL,SAVE :: ecrit_LES_omp
193	REAL,SAVE :: ecrit_tra_omp
194	REAL,SAVE :: cvl_comp_threshold_omp
195	REAL,SAVE :: cvl_sig2feed_omp
196	REAL,SAVE :: cvl_corr_omp
197	LOGICAL,SAVE :: ok_lic_melt_omp
198	!
199	LOGICAL,SAVE :: cycle_diurne_omp,soil_model_omp,new_oliq_omp
200	LOGICAL,SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp
201	INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp
202	INTEGER, SAVE :: iflag_ener_conserv_omp
203	LOGICAL, SAVE :: ok_conserv_q_omp
204	INTEGER, SAVE :: iflag_fisrtilp_qsat_omp
205	LOGICAL,SAVE :: ok_strato_omp
206	LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
207	real, SAVE:: GWD_RANDO_RUWMAX_omp, gwd_rando_sat_omp
208	LOGICAL,SAVE :: ok_qch4_omp
209	LOGICAL,SAVE :: carbon_cycle_tr_omp
210	LOGICAL,SAVE :: carbon_cycle_cpl_omp
211
212	integer, intent(out):: read_climoz ! read ozone climatology, OpenMP shared
213	! Allowed values are 0, 1 and 2
214	! 0: do not read an ozone climatology
215	! 1: read a single ozone climatology that will be used day and night
216	! 2: read two ozone climatologies, the average day and night
217	! climatology and the daylight climatology
218
219	!-----------------------------------------------------------------
220
221	print*,'CONFPHYS ENTREE'
222	!$OMP MASTER
223	!Config Key = type_ocean
224	!Config Desc = Type d'ocean
225	!Config Def = force
226	!Config Help = Type d'ocean utilise: force, slab,couple
227	!
228	type_ocean_omp = 'force '
229	call getin('type_ocean', type_ocean_omp)
230	!
231	!Config Key = version_ocean
232	!Config Desc = Version d'ocean
233	!Config Def = xxxxxx
234	!Config Help = Version d'ocean utilise: opa8/nemo/sicOBS/xxxxxx
235	!
236	version_ocean_omp = 'xxxxxx'
237	call getin('version_ocean', version_ocean_omp)
238
239	!Config Key = OCEAN
240	!Config Desc = Old parameter name for type_ocean
241	!Config Def = yyyyyy
242	!Config Help = This is only for testing purpose
243	!
244	ocean_omp = 'yyyyyy'
245	call getin('OCEAN', ocean_omp)
246	IF (ocean_omp /= 'yyyyyy') THEN
247	WRITE(lunout,*)'ERROR! Old variable name OCEAN used in parmeter file.'
248	WRITE(lunout,*)'Variable OCEAN has been replaced by the variable type_ocean.'
249	WRITE(lunout,*)'You have to update your parameter file physiq.def to succed running'
250	CALL abort_physic('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1)
251	END IF
252
253	!Config Key = t_coupl
254	!Config Desc = Pas de temps du couplage atm/oce en sec.
255	!Config Def = 86400
256	!Config Help = This is only for testing purpose
257	!
258	t_coupl_omp = 86400.
259	call getin('t_coupl', t_coupl_omp)
260	IF (t_coupl_omp == 0) THEN
261	WRITE(lunout,*)'ERROR! Timestep of coupling between atmosphere and ocean'
262	WRITE(lunout,*)'cannot be zero.'
263	CALL abort_physic('conf_phys','t_coupl = 0.',1)
264	END IF
265
266	!
267	!Config Key = ok_all_xml
268	!Config Desc = utiliser les xml pourles définitions des champs pour xios
269	!Config Def = .false.
270	!Config Help =
271	!
272	ok_all_xml_omp = .false.
273	call getin('ok_all_xml', ok_all_xml_omp)
274	!
275
276	!
277	!Config Key = VEGET
278	!Config Desc = Type de modele de vegetation
279	!Config Def = .false.
280	!Config Help = Type de modele de vegetation utilise
281	!
282	type_veget_omp ='orchidee'
283	call getin('VEGET', type_veget_omp)
284	!
285
286	! Martin
287	!Config Key = ok_snow
288	!Config Desc = Flag to activate snow model SISVAT
289	!Config Def = .false.
290	ok_snow_omp = .false.
291	call getin('ok_snow', ok_snow_omp)
292	! Martin
293
294	!Config Key = OK_journe
295	!Config Desc = Pour des sorties journalieres
296	!Config Def = .false.
297	!Config Help = Pour creer le fichier histday contenant les sorties
298	! journalieres
299	!
300	ok_journe_omp = .false.
301	call getin('OK_journe', ok_journe_omp)
302	!
303	!Config Key = ok_hf
304	!Config Desc = Pour des sorties haute frequence
305	!Config Def = .false.
306	!Config Help = Pour creer le fichier histhf contenant les sorties
307	! haute frequence ( 3h ou 6h)
308	!
309	ok_hf_omp = .false.
310	call getin('ok_hf', ok_hf_omp)
311	!
312	!Config Key = OK_mensuel
313	!Config Desc = Pour des sorties mensuelles
314	!Config Def = .true.
315	!Config Help = Pour creer le fichier histmth contenant les sorties
316	! mensuelles
317	!
318	ok_mensuel_omp = .true.
319	call getin('OK_mensuel', ok_mensuel_omp)
320	!
321	!Config Key = OK_instan
322	!Config Desc = Pour des sorties instantanees
323	!Config Def = .false.
324	!Config Help = Pour creer le fichier histins contenant les sorties
325	! instantanees
326	!
327	ok_instan_omp = .false.
328	call getin('OK_instan', ok_instan_omp)
329	!
330	!Config Key = ok_ade
331	!Config Desc = Aerosol direct effect or not?
332	!Config Def = .false.
333	!Config Help = Used in radlwsw.F
334	!
335	ok_ade_omp = .false.
336	call getin('ok_ade', ok_ade_omp)
337
338	!
339	!Config Key = ok_aie
340	!Config Desc = Aerosol indirect effect or not?
341	!Config Def = .false.
342	!Config Help = Used in nuage.F and radlwsw.F
343	!
344	ok_aie_omp = .false.
345	call getin('ok_aie', ok_aie_omp)
346
347	!
348	!Config Key = ok_cdnc
349	!Config Desc = ok cloud droplet number concentration
350	!Config Def = .false.
351	!Config Help = Used in newmicro.F
352	!
353	ok_cdnc_omp = .false.
354	call getin('ok_cdnc', ok_cdnc_omp)
355	!
356	!Config Key = aerosol_couple
357	!Config Desc = read aerosol in file or calcul by inca
358	!Config Def = .false.
359	!Config Help = Used in physiq.F
360	!
361	aerosol_couple_omp = .false.
362	CALL getin('aerosol_couple',aerosol_couple_omp)
363	!
364	!Config Key = flag_aerosol
365	!Config Desc = which aerosol is use for coupled model
366	!Config Def = 1
367	!Config Help = Used in physiq.F
368	!
369	! - flag_aerosol=0 => no aerosol
370	! - flag_aerosol=1 => so4 only (defaut)
371	! - flag_aerosol=2 => bc only
372	! - flag_aerosol=3 => pom only
373	! - flag_aerosol=4 => seasalt only
374	! - flag_aerosol=5 => dust only
375	! - flag_aerosol=6 => all aerosol
376
377	flag_aerosol_omp = 0
378	CALL getin('flag_aerosol',flag_aerosol_omp)
379	!
380	!Config Key = flag_aerosol_strat
381	!Config Desc = use stratospheric aerosols T/F
382	!Config Def = false
383	!Config Help = Used in physiq.F
384	!
385	!
386	flag_aerosol_strat_omp = .false.
387	CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
388
389	! Temporary variable for testing purpose!
390	!Config Key = new_aod
391	!Config Desc = which calcul of aeropt
392	!Config Def = false
393	!Config Help = Used in physiq.F
394	!
395	new_aod_omp = .true.
396	CALL getin('new_aod',new_aod_omp)
397
398	!
399	!Config Key = aer_type
400	!Config Desc = Use a constant field for the aerosols
401	!Config Def = scenario
402	!Config Help = Used in readaerosol.F90
403	!
404	aer_type_omp = 'scenario'
405	call getin('aer_type', aer_type_omp)
406
407	!
408	!Config Key = bl95_b0
409	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
410	!Config Def = .false.
411	!Config Help = Used in nuage.F
412	!
413	bl95_b0_omp = 2.
414	call getin('bl95_b0', bl95_b0_omp)
415
416	!Config Key = bl95_b1
417	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
418	!Config Def = .false.
419	!Config Help = Used in nuage.F
420	!
421	bl95_b1_omp = 0.2
422	call getin('bl95_b1', bl95_b1_omp)
423
424	!Config Key = freq_ISCCP
425	!Config Desc = Frequence d'appel du simulateur ISCCP en secondes;
426	! par defaut 10800, i.e. 3 heures
427	!Config Def = 10800.
428	!Config Help = Used in ini_histISCCP.h
429	!
430	freq_ISCCP_omp = 10800.
431	call getin('freq_ISCCP', freq_ISCCP_omp)
432	!
433	!Config Key = ecrit_ISCCP
434	!Config Desc = Frequence d'ecriture des resultats du simulateur ISCCP en nombre de jours;
435	! par defaut 1., i.e. 1 jour
436	!Config Def = 1.
437	!Config Help = Used in ini_histISCCP.h
438	!
439	!
440	ecrit_ISCCP_omp = 1.
441	call getin('ecrit_ISCCP', ecrit_ISCCP_omp)
442
443	!Config Key = freq_COSP
444	!Config Desc = Frequence d'appel du simulateur COSP en secondes;
445	! par defaut 10800, i.e. 3 heures
446	!Config Def = 10800.
447	!Config Help = Used in ini_histdayCOSP.h
448	!
449	freq_COSP_omp = 10800.
450	call getin('freq_COSP', freq_COSP_omp)
451
452	!
453	!Config Key = ip_ebil_phy
454	!Config Desc = Niveau de sortie pour les diags bilan d'energie
455	!Config Def = 0
456	!Config Help =
457	!
458	ip_ebil_phy_omp = 0
459	call getin('ip_ebil_phy', ip_ebil_phy_omp)
460	!
461	!Config Key = seuil_inversion
462	!Config Desc = Seuil ur dTh pour le choix entre les schemas de CL
463	!Config Def = -0.1
464	!Config Help =
465	!
466	seuil_inversion_omp = -0.1
467	call getin('seuil_inversion', seuil_inversion_omp)
468
469	!
470	! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
471	!
472	!Config Key = R_ecc
473	!Config Desc = Excentricite
474	!Config Def = 0.016715
475	!Config Help =
476	!
477	!valeur AMIP II
478	R_ecc_omp = 0.016715
479	call getin('R_ecc', R_ecc_omp)
480	!
481	!Config Key = R_peri
482	!Config Desc = Equinoxe
483	!Config Def =
484	!Config Help =
485	!
486	!
487	!valeur AMIP II
488	R_peri_omp = 102.7
489	call getin('R_peri', R_peri_omp)
490	!
491	!Config Key = R_incl
492	!Config Desc = Inclinaison
493	!Config Def =
494	!Config Help =
495	!
496	!
497	!valeur AMIP II
498	R_incl_omp = 23.441
499	call getin('R_incl', R_incl_omp)
500	!
501	!Config Key = solaire
502	!Config Desc = Constante solaire en W/m2
503	!Config Def = 1365.
504	!Config Help =
505	!
506	!
507	!valeur AMIP II
508	solaire_omp = 1365.
509	call getin('solaire', solaire_omp)
510	!
511	!Config Key = co2_ppm
512	!Config Desc = concentration du gaz carbonique en ppmv
513	!Config Def = 348.
514	!Config Help =
515	!
516	!
517	!valeur AMIP II
518	co2_ppm_omp = 348.
519	call getin('co2_ppm', co2_ppm_omp)
520	!
521	!Config Key = RCO2
522	!Config Desc = Concentration du CO2
523	!Config Def = co2_ppm * 1.0e-06 * 44.011/28.97
524	!Config Def = 348. * 1.0e-06 * 44.011/28.97
525	!Config Help =
526	!
527	! RCO2 = 5.286789092164308E-04
528	!ancienne valeur
529	RCO2_omp = co2_ppm_omp * 1.0e-06 * 44.011/28.97 ! pour co2_ppm=348.
530
531	! call getin('RCO2', RCO2)
532	!
533	!Config Key = RCH4
534	!Config Desc = Concentration du CH4
535	!Config Def = 1.65E-06* 16.043/28.97
536	!Config Help =
537	!
538	!
539	!valeur AMIP II
540	!OK RCH4 = 1.65E-06* 16.043/28.97
541	! RCH4 = 9.137366240938903E-07
542	!
543	!ancienne valeur
544	! RCH4 = 1.72E-06* 16.043/28.97
545	!OK call getin('RCH4', RCH4)
546	zzz = 1650.
547	call getin('CH4_ppb', zzz)
548	CH4_ppb_omp = zzz
549	RCH4_omp = CH4_ppb_omp * 1.0E-09 * 16.043/28.97
550	!
551	!Config Key = RN2O
552	!Config Desc = Concentration du N2O
553	!Config Def = 306.E-09* 44.013/28.97
554	!Config Help =
555	!
556	!
557	!valeur AMIP II
558	!OK RN2O = 306.E-09* 44.013/28.97
559	! RN2O = 4.648939592682085E-07
560	!
561	!ancienne valeur
562	! RN2O = 310.E-09* 44.013/28.97
563	!OK call getin('RN2O', RN2O)
564	zzz=306.
565	call getin('N2O_ppb', zzz)
566	N2O_ppb_omp = zzz
567	RN2O_omp = N2O_ppb_omp * 1.0E-09 * 44.013/28.97
568	!
569	!Config Key = RCFC11
570	!Config Desc = Concentration du CFC11
571	!Config Def = 280.E-12* 137.3686/28.97
572	!Config Help =
573	!
574	!
575	!OK RCFC11 = 280.E-12* 137.3686/28.97
576	zzz = 280.
577	call getin('CFC11_ppt',zzz)
578	CFC11_ppt_omp = zzz
579	RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * 137.3686/28.97
580	! RCFC11 = 1.327690990680013E-09
581	!OK call getin('RCFC11', RCFC11)
582	!
583	!Config Key = RCFC12
584	!Config Desc = Concentration du CFC12
585	!Config Def = 484.E-12* 120.9140/28.97
586	!Config Help =
587	!
588	!
589	!OK RCFC12 = 484.E-12* 120.9140/28.97
590	zzz = 484.
591	call getin('CFC12_ppt',zzz)
592	CFC12_ppt_omp = zzz
593	RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * 120.9140/28.97
594	! RCFC12 = 2.020102726958923E-09
595	!OK call getin('RCFC12', RCFC12)
596
597	!ajout CFMIP begin
598	!
599	!Config Key = co2_ppm_per
600	!Config Desc = concentration du co2_ppm_per
601	!Config Def = 348.
602	!Config Help =
603	!
604	co2_ppm_per_omp = co2_ppm_omp
605	call getin('co2_ppm_per', co2_ppm_per_omp)
606	!
607	!Config Key = RCO2_per
608	!Config Desc = Concentration du CO2_per
609	!Config Def = co2_ppm_per * 1.0e-06 * 44.011/28.97
610	!Config Def = 348. * 1.0e-06 * 44.011/28.97
611	!Config Help =
612	!
613	RCO2_per_omp = co2_ppm_per_omp * 1.0e-06 * 44.011/28.97
614
615	!Config Key = ok_4xCO2atm
616	!Config Desc = Calcul ou non effet radiatif 4xco2
617	!Config Def = .false.
618	!Config Help =
619
620	ok_4xCO2atm_omp = .false.
621	call getin('ok_4xCO2atm',ok_4xCO2atm_omp)
622
623	!Config Key = RCH4_per
624	!Config Desc = Concentration du CH4_per
625	!Config Def = 1.65E-06* 16.043/28.97
626	!Config Help =
627	!
628	zzz = CH4_ppb_omp
629	call getin('CH4_ppb_per', zzz)
630	CH4_ppb_per_omp = zzz
631	RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * 16.043/28.97
632	!
633	!Config Key = RN2O_per
634	!Config Desc = Concentration du N2O_per
635	!Config Def = 306.E-09* 44.013/28.97
636	!Config Help =
637	!
638	zzz = N2O_ppb_omp
639	call getin('N2O_ppb_per', zzz)
640	N2O_ppb_per_omp = zzz
641	RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * 44.013/28.97
642	!
643	!Config Key = RCFC11_per
644	!Config Desc = Concentration du CFC11_per
645	!Config Def = 280.E-12* 137.3686/28.97
646	!Config Help =
647	!
648	zzz = CFC11_ppt_omp
649	call getin('CFC11_ppt_per',zzz)
650	CFC11_ppt_per_omp = zzz
651	RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * 137.3686/28.97
652	!
653	!Config Key = RCFC12_per
654	!Config Desc = Concentration du CFC12_per
655	!Config Def = 484.E-12* 120.9140/28.97
656	!Config Help =
657	!
658	zzz = CFC12_ppt_omp
659	call getin('CFC12_ppt_per',zzz)
660	CFC12_ppt_per_omp = zzz
661	RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * 120.9140/28.97
662	!ajout CFMIP end
663
664	!
665	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
666	! Constantes precedemment dans dyn3d/conf_gcm
667
668	!Config Key = cycle_diurne
669	!Config Desc = Cycle ddiurne
670	!Config Def = y
671	!Config Help = Cette option permet d'eteidre le cycle diurne.
672	!Config Peut etre util pour accelerer le code !
673	cycle_diurne_omp = .TRUE.
674	CALL getin('cycle_diurne',cycle_diurne_omp)
675
676	!Config Key = soil_model
677	!Config Desc = Modele de sol
678	!Config Def = y
679	!Config Help = Choix du modele de sol (Thermique ?)
680	!Config Option qui pourait un string afin de pouvoir
681	!Config plus de choix ! Ou meme une liste d'options !
682	soil_model_omp = .TRUE.
683	CALL getin('soil_model',soil_model_omp)
684
685	!Config Key = new_oliq
686	!Config Desc = Nouvelle eau liquide
687	!Config Def = y
688	!Config Help = Permet de mettre en route la
689	!Config nouvelle parametrisation de l'eau liquide !
690	new_oliq_omp = .TRUE.
691	CALL getin('new_oliq',new_oliq_omp)
692
693	!Config Key = ok_orodr
694	!Config Desc = Orodr ???
695	!Config Def = y
696	!Config Help = Y en a pas comprendre !
697	!Config
698	ok_orodr_omp = .TRUE.
699	CALL getin('ok_orodr',ok_orodr_omp)
700
701	!Config Key = ok_orolf
702	!Config Desc = Orolf ??
703	!Config Def = y
704	!Config Help = Connais pas !
705	ok_orolf_omp = .TRUE.
706	CALL getin('ok_orolf', ok_orolf_omp)
707
708	!Config Key = ok_limitvrai
709	!Config Desc = Force la lecture de la bonne annee
710	!Config Def = n
711	!Config Help = On peut forcer le modele a lire le
712	!Config fichier SST de la bonne annee. C'est une tres bonne
713	!Config idee, pourquoi ne pas mettre toujours a y ???
714	ok_limitvrai_omp = .FALSE.
715	CALL getin('ok_limitvrai',ok_limitvrai_omp)
716
717	!Config Key = nbapp_rad
718	!Config Desc = Frequence d'appel au rayonnement
719	!Config Def = 12
720	!Config Help = Nombre d'appels des routines de rayonnements
721	!Config par jour.
722	nbapp_rad_omp = 12
723	CALL getin('nbapp_rad',nbapp_rad_omp)
724
725	!Config Key = iflag_con
726	!Config Desc = Flag de convection
727	!Config Def = 2
728	!Config Help = Flag pour la convection les options suivantes existent :
729	!Config 1 pour LMD,
730	!Config 2 pour Tiedtke,
731	!Config 3 pour CCM(NCAR)
732	iflag_con_omp = 2
733	CALL getin('iflag_con',iflag_con_omp)
734
735	!Config Key = iflag_ener_conserv
736	!Config Desc = Flag de convection
737	!Config Def = 1
738	!Config Help = Flag pour la convection les options suivantes existent :
739	!Config -1 pour Kinetic energy correction
740	!Config 1 conservation kinetic and enthalpy
741	iflag_ener_conserv_omp = -1
742	CALL getin('iflag_ener_conserv',iflag_ener_conserv_omp)
743
744	!Config Key = ok_conserv_q
745	!Config Desc = Switch des corrections de conservation de l'eau
746	!Config Def = y
747	!Config Help = Switch des corrections de conservation de l'eau
748	!Config y -> corrections activees
749	!Config n -> conformite avec versions anterieures au 1/4/2014
750	ok_conserv_q_omp = .FALSE.
751	CALL getin('ok_conserv_q',ok_conserv_q_omp)
752
753	!Config Key = iflag_fisrtilp_qsat
754	!Config Desc = Flag de convection
755	!Config Def = 1
756	!Config Help = Flag pour la convection les options suivantes existent :
757	!Config -1 pour Kinetic energy correction
758	!Config 1 conservation kinetic and enthalpy
759	iflag_fisrtilp_qsat_omp = 0
760	CALL getin('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat_omp)
761
762	!
763	!
764	!
765	! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
766	!
767	! KE
768	!
769
770	!Config key = cvl_comp_threshold
771	!Config Desc = maximum fraction of convective points enabling compression
772	!Config Def = 1.00
773	!Config Help = fields are compressed when less than a fraction cvl_comp_threshold
774	!Config Help = of the points is convective.
775	cvl_comp_threshold_omp = 1.00
776	CALL getin('cvl_comp_threshold', cvl_comp_threshold_omp)
777
778	!Config key = cvl_sig2feed
779	!Config Desc = sigma coordinate at top of feeding layer
780	!Config Def = 0.97
781	!Config Help = deep convection is fed by the layer extending from the surface (pressure ps)
782	!Config Help = and cvl_sig2feed*ps.
783	cvl_sig2feed_omp = 0.97
784	CALL getin('cvl_sig2feed', cvl_sig2feed_omp)
785
786	!Config key = cvl_corr
787	!Config Desc = Facteur multiplication des precip convectives dans KE
788	!Config Def = 1.00
789	!Config Help = 1.02 pour un moderne ou un pre-ind. A ajuster pour un glaciaire
790	cvl_corr_omp = 1.00
791	CALL getin('cvl_corr', cvl_corr_omp)
792
793
794	!Config Key = epmax
795	!Config Desc = Efficacite precip
796	!Config Def = 0.993
797	!Config Help =
798	!
799	epmax_omp = .993
800	call getin('epmax', epmax_omp)
801	!
802	!Config Key = ok_adj_ema
803	!Config Desc =
804	!Config Def = false
805	!Config Help =
806	!
807	ok_adj_ema_omp = .false.
808	call getin('ok_adj_ema',ok_adj_ema_omp)
809	!
810	!Config Key = iflag_clw
811	!Config Desc =
812	!Config Def = 0
813	!Config Help =
814	!
815	iflag_clw_omp = 0
816	call getin('iflag_clw',iflag_clw_omp)
817	!
818	!Config Key = cld_lc_lsc
819	!Config Desc =
820	!Config Def = 2.6e-4
821	!Config Help =
822	!
823	cld_lc_lsc_omp = 2.6e-4
824	call getin('cld_lc_lsc',cld_lc_lsc_omp)
825	!
826	!Config Key = cld_lc_con
827	!Config Desc =
828	!Config Def = 2.6e-4
829	!Config Help =
830	!
831	cld_lc_con_omp = 2.6e-4
832	call getin('cld_lc_con',cld_lc_con_omp)
833	!
834	!Config Key = cld_tau_lsc
835	!Config Desc =
836	!Config Def = 3600.
837	!Config Help =
838	!
839	cld_tau_lsc_omp = 3600.
840	call getin('cld_tau_lsc',cld_tau_lsc_omp)
841	!
842	!Config Key = cld_tau_con
843	!Config Desc =
844	!Config Def = 3600.
845	!Config Help =
846	!
847	cld_tau_con_omp = 3600.
848	call getin('cld_tau_con',cld_tau_con_omp)
849	!
850	!Config Key = ffallv_lsc
851	!Config Desc =
852	!Config Def = 1.
853	!Config Help =
854	!
855	ffallv_lsc_omp = 1.
856	call getin('ffallv_lsc',ffallv_lsc_omp)
857	!
858	!Config Key = ffallv_con
859	!Config Desc =
860	!Config Def = 1.
861	!Config Help =
862	!
863	ffallv_con_omp = 1.
864	call getin('ffallv_con',ffallv_con_omp)
865	!
866	!Config Key = coef_eva
867	!Config Desc =
868	!Config Def = 2.e-5
869	!Config Help =
870	!
871	coef_eva_omp = 2.e-5
872	call getin('coef_eva',coef_eva_omp)
873	!
874	!Config Key = reevap_ice
875	!Config Desc =
876	!Config Def = .false.
877	!Config Help =
878	!
879	reevap_ice_omp = .false.
880	call getin('reevap_ice',reevap_ice_omp)
881
882	!Config Key = iflag_ratqs
883	!Config Desc =
884	!Config Def = 1
885	!Config Help =
886	!
887	iflag_ratqs_omp = 1
888	call getin('iflag_ratqs',iflag_ratqs_omp)
889
890	!
891	!Config Key = iflag_radia
892	!Config Desc =
893	!Config Def = 1
894	!Config Help =
895	!
896	iflag_radia_omp = 1
897	call getin('iflag_radia',iflag_radia_omp)
898
899	!
900	!Config Key = iflag_rrtm
901	!Config Desc =
902	!Config Def = 0
903	!Config Help =
904	!
905	iflag_rrtm_omp = 0
906	call getin('iflag_rrtm',iflag_rrtm_omp)
907
908	!
909	!Config Key = NSW
910	!Config Desc =
911	!Config Def = 0
912	!Config Help =
913	!
914	NSW_omp = 6
915	call getin('NSW',NSW_omp)
916	!albedo SB >>>
917	iflag_albedo_omp = 0
918	call getin('iflag_albedo',iflag_albedo_omp)
919
920	ok_chlorophyll_omp=.false.
921	call getin('ok_chlorophyll',ok_chlorophyll_omp)
922	!albedo SB <<<
923
924	!
925	!Config Key = iflag_cld_th
926	!Config Desc =
927	!Config Def = 1
928	!Config Help =
929	!
930	iflag_cld_th_omp = 1
931	! On lit deux fois avec l'ancien et le nouveau nom
932	! pour assurer une retrocompatiblite.
933	! A abandonner un jour
934	call getin('iflag_cldcon',iflag_cld_th_omp)
935	call getin('iflag_cld_th',iflag_cld_th_omp)
936	iflag_cld_cv_omp = 0
937	call getin('iflag_cld_cv',iflag_cld_cv_omp)
938
939	!
940	!Config Key = tau_cld_cv
941	!Config Desc =
942	!Config Def = 10.
943	!Config Help =
944	!
945	tau_cld_cv_omp = 10.
946	call getin('tau_cld_cv',tau_cld_cv_omp)
947
948	!
949	!Config Key = coefw_cld_cv
950	!Config Desc =
951	!Config Def = 0.1
952	!Config Help =
953	!
954	coefw_cld_cv_omp = 0.1
955	call getin('coefw_cld_cv',coefw_cld_cv_omp)
956
957
958
959
960	!
961	!Config Key = iflag_pdf
962	!Config Desc =
963	!Config Def = 0
964	!Config Help =
965	!
966	iflag_pdf_omp = 0
967	call getin('iflag_pdf',iflag_pdf_omp)
968	!
969	!Config Key = fact_cldcon
970	!Config Desc =
971	!Config Def = 0.375
972	!Config Help =
973	!
974	fact_cldcon_omp = 0.375
975	call getin('fact_cldcon',fact_cldcon_omp)
976
977	!
978	!Config Key = facttemps
979	!Config Desc =
980	!Config Def = 1.e-4
981	!Config Help =
982	!
983	facttemps_omp = 1.e-4
984	call getin('facttemps',facttemps_omp)
985
986	!
987	!Config Key = ok_newmicro
988	!Config Desc =
989	!Config Def = .true.
990	!Config Help =
991	!
992	ok_newmicro_omp = .true.
993	call getin('ok_newmicro',ok_newmicro_omp)
994	!
995	!Config Key = ratqsbas
996	!Config Desc =
997	!Config Def = 0.01
998	!Config Help =
999	!
1000	ratqsbas_omp = 0.01
1001	call getin('ratqsbas',ratqsbas_omp)
1002	!
1003	!Config Key = ratqshaut
1004	!Config Desc =
1005	!Config Def = 0.3
1006	!Config Help =
1007	!
1008	ratqshaut_omp = 0.3
1009	call getin('ratqshaut',ratqshaut_omp)
1010
1011	!Config Key = tau_ratqs
1012	!Config Desc =
1013	!Config Def = 1800.
1014	!Config Help =
1015	!
1016	tau_ratqs_omp = 1800.
1017	call getin('tau_ratqs',tau_ratqs_omp)
1018
1019	!
1020	!-----------------------------------------------------------------------
1021	! Longitude solaire pour le calcul de l'ensoleillement en degre
1022	! si on veut imposer la saison. Sinon, solarlong0=-999.999
1023	!Config Key = solarlong0
1024	!Config Desc =
1025	!Config Def = -999.999
1026	!Config Help =
1027	!
1028	solarlong0_omp = -999.999
1029	call getin('solarlong0',solarlong0_omp)
1030	!
1031	!-----------------------------------------------------------------------
1032	! Valeur imposee pour configuration idealisees
1033	!Config Key = qsol0 pour le bucket, evap0 pour aquaplanetes, albsno0
1034	! Default value -1 to activate the full computation
1035	qsol0_omp = -1.
1036	call getin('qsol0',qsol0_omp)
1037	evap0_omp = -1.
1038	call getin('evap0',evap0_omp)
1039	albsno0_omp = -1.
1040	call getin('albsno0',albsno0_omp)
1041	!
1042	!-----------------------------------------------------------------------
1043	!
1044	!Config Key = inertie_ice
1045	!Config Desc =
1046	!Config Def = 2000.
1047	!Config Help =
1048	!
1049	inertie_ice_omp = 2000.
1050	call getin('inertie_ice',inertie_ice_omp)
1051	!
1052	!Config Key = inertie_sno
1053	!Config Desc =
1054	!Config Def = 2000.
1055	!Config Help =
1056	!
1057	inertie_sno_omp = 2000.
1058	call getin('inertie_sno',inertie_sno_omp)
1059	!
1060	!Config Key = inertie_sol
1061	!Config Desc =
1062	!Config Def = 2000.
1063	!Config Help =
1064	!
1065	inertie_sol_omp = 2000.
1066	call getin('inertie_sol',inertie_sol_omp)
1067
1068	!
1069	!Config Key = rad_froid
1070	!Config Desc =
1071	!Config Def = 35.0
1072	!Config Help =
1073	!
1074	rad_froid_omp = 35.0
1075	call getin('rad_froid',rad_froid_omp)
1076
1077	!
1078	!Config Key = rad_chau1
1079	!Config Desc =
1080	!Config Def = 13.0
1081	!Config Help =
1082	!
1083	rad_chau1_omp = 13.0
1084	call getin('rad_chau1',rad_chau1_omp)
1085
1086	!
1087	!Config Key = rad_chau2
1088	!Config Desc =
1089	!Config Def = 9.0
1090	!Config Help =
1091	!
1092	rad_chau2_omp = 9.0
1093	call getin('rad_chau2',rad_chau2_omp)
1094
1095	!
1096	!Config Key = t_glace_min
1097	!Config Desc =
1098	!Config Def = 258.
1099	!Config Help =
1100	!
1101	t_glace_min_omp = 258.
1102	call getin('t_glace_min',t_glace_min_omp)
1103
1104	!
1105	!Config Key = t_glace_max
1106	!Config Desc =
1107	!Config Def = 273.13
1108	!Config Help =
1109	!
1110	t_glace_max_omp = 273.13
1111	call getin('t_glace_max',t_glace_max_omp)
1112
1113	!
1114	!Config Key = exposant_glace
1115	!Config Desc =
1116	!Config Def = 2.
1117	!Config Help =
1118	!
1119	exposant_glace_omp = 1.
1120	call getin('exposant_glace',exposant_glace_omp)
1121
1122	!
1123	!Config Key = iflag_t_glace
1124	!Config Desc =
1125	!Config Def = 0
1126	!Config Help =
1127	!
1128	iflag_t_glace_omp = 0
1129	call getin('iflag_t_glace',iflag_t_glace_omp)
1130
1131	!
1132	!Config Key = iflag_ice_thermo
1133	!Config Desc =
1134	!Config Def = 0
1135	!Config Help =
1136	!
1137	iflag_ice_thermo_omp = 0
1138	call getin('iflag_ice_thermo',iflag_ice_thermo_omp)
1139
1140	!Config Key = rei_min
1141	!Config Desc =
1142	!Config Def = 3.5
1143	!Config Help =
1144	!
1145	rei_min_omp = 3.5
1146	call getin('rei_min',rei_min_omp)
1147
1148	!
1149	!Config Key = rei_max
1150	!Config Desc =
1151	!Config Def = 61.29
1152	!Config Help =
1153	!
1154	rei_max_omp = 61.29
1155	call getin('rei_max',rei_max_omp)
1156
1157	!
1158	!Config Key = top_height
1159	!Config Desc =
1160	!Config Def = 3
1161	!Config Help =
1162	!
1163	top_height_omp = 3
1164	call getin('top_height',top_height_omp)
1165
1166	!
1167	!Config Key = overlap
1168	!Config Desc =
1169	!Config Def = 3
1170	!Config Help =
1171	!
1172	overlap_omp = 3
1173	call getin('overlap',overlap_omp)
1174
1175
1176	!
1177	!
1178	!Config Key = cdmmax
1179	!Config Desc =
1180	!Config Def = 1.3E-3
1181	!Config Help =
1182	!
1183	cdmmax_omp = 1.3E-3
1184	call getin('cdmmax',cdmmax_omp)
1185
1186	!
1187	!Config Key = cdhmax
1188	!Config Desc =
1189	!Config Def = 1.1E-3
1190	!Config Help =
1191	!
1192	cdhmax_omp = 1.1E-3
1193	call getin('cdhmax',cdhmax_omp)
1194
1195	!261103
1196	!
1197	!Config Key = ksta
1198	!Config Desc =
1199	!Config Def = 1.0e-10
1200	!Config Help =
1201	!
1202	ksta_omp = 1.0e-10
1203	call getin('ksta',ksta_omp)
1204
1205	!
1206	!Config Key = ksta_ter
1207	!Config Desc =
1208	!Config Def = 1.0e-10
1209	!Config Help =
1210	!
1211	ksta_ter_omp = 1.0e-10
1212	call getin('ksta_ter',ksta_ter_omp)
1213
1214	!Config Key = f_ri_cd_min
1215	!Config Desc =
1216	!Config Def = 0.1
1217	!Config Help =
1218	!
1219	f_ri_cd_min_omp = 0.1
1220	call getin('f_ri_cd_min',f_ri_cd_min_omp)
1221
1222	!
1223	!Config Key = ok_kzmin
1224	!Config Desc =
1225	!Config Def = .true.
1226	!Config Help =
1227	!
1228	ok_kzmin_omp = .true.
1229	call getin('ok_kzmin',ok_kzmin_omp)
1230
1231	!
1232	!Config Key = fmagic
1233	!Config Desc = additionnal multiplicator factor used for albedo
1234	!Config Def = 1.
1235	!Config Help = additionnal multiplicator factor used in albedo.F
1236	!
1237	fmagic_omp = 1.
1238	call getin('fmagic',fmagic_omp)
1239
1240	!
1241	!Config Key = pmagic
1242	!Config Desc = additional factor used for albedo
1243	!Config Def = 0.
1244	!Config Help = additional factor used in albedo.F
1245	!
1246	pmagic_omp = 0.
1247	call getin('pmagic',pmagic_omp)
1248
1249
1250	!Config Key = ok_lic_melt
1251	!Config Desc = Prise en compte de la fonte de la calotte dans le bilan d'eau
1252	!Config Def = .false.
1253	!Config Help = mettre a .false. pour assurer la conservation en eau
1254	ok_lic_melt_omp = .false.
1255	call getin('ok_lic_melt', ok_lic_melt_omp)
1256
1257	!
1258	! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
1259	!
1260
1261	!Config Key = iflag_pbl
1262	!Config Desc =
1263	!Config Def = 1
1264	!Config Help =
1265	!
1266	iflag_pbl_omp = 1
1267	call getin('iflag_pbl',iflag_pbl_omp)
1268	!
1269	!Config Key = iflag_pbl_split
1270	!Config Desc = binary flag: least signif bit = split vdf; next bit = split thermals
1271	!Config Def = 0
1272	!Config Help = 0-> no splitting; 1-> vdf splitting; 2-> thermals splitting; 3-> full splitting
1273	!
1274	iflag_pbl_split_omp = 0
1275	call getin('iflag_pbl_split',iflag_pbl_split_omp)
1276	!
1277	!Config Key = iflag_thermals
1278	!Config Desc =
1279	!Config Def = 0
1280	!Config Help =
1281	!
1282	iflag_thermals_omp = 0
1283	call getin('iflag_thermals',iflag_thermals_omp)
1284	!
1285	!Config Key = iflag_thermals_ed
1286	!Config Desc =
1287	!Config Def = 0
1288	!Config Help =
1289	!
1290	fact_thermals_ed_dz_omp = 0.1
1291
1292	call getin('fact_thermals_ed_dz',fact_thermals_ed_dz_omp)
1293	!
1294	!
1295	!Config Key = iflag_thermals_ed
1296	!Config Desc =
1297	!Config Def = 0
1298	!Config Help =
1299	!
1300	iflag_thermals_ed_omp = 0
1301	call getin('iflag_thermals_ed',iflag_thermals_ed_omp)
1302	!
1303	!
1304	!Config Key = iflag_thermals_optflux
1305	!Config Desc =
1306	!Config Def = 0
1307	!Config Help =
1308	!
1309	iflag_thermals_optflux_omp = 0
1310	call getin('iflag_thermals_optflux',iflag_thermals_optflux_omp)
1311	!
1312	!Config Key = iflag_thermals_closure
1313	!Config Desc =
1314	!Config Def = 0
1315	!Config Help =
1316	!
1317	iflag_thermals_closure_omp = 1
1318	call getin('iflag_thermals_closure',iflag_thermals_closure_omp)
1319	!
1320	!
1321	!
1322	!
1323	!Config Key = nsplit_thermals
1324	!Config Desc =
1325	!Config Def = 1
1326	!Config Help =
1327	!
1328	nsplit_thermals_omp = 1
1329	call getin('nsplit_thermals',nsplit_thermals_omp)
1330
1331	!Config Key = alp_bl_k
1332	!Config Desc =
1333	!Config Def = 0.
1334	!Config Help =
1335	!
1336	alp_bl_k_omp = 1.
1337	call getin('alp_bl_k',alp_bl_k_omp)
1338
1339	! nrlmd le 10/04/2012
1340
1341	!Config Key = iflag_trig_bl
1342	!Config Desc =
1343	!Config Def = 0
1344	!Config Help =
1345	!
1346	iflag_trig_bl_omp = 0
1347	call getin('iflag_trig_bl',iflag_trig_bl_omp)
1348
1349	!Config Key = s_trig_bl
1350	!Config Desc =
1351	!Config Def = 0
1352	!Config Help =
1353	!
1354	s_trig_omp = 2e7
1355	call getin('s_trig',s_trig_omp)
1356
1357	!Config Key = tau_trig_shallow
1358	!Config Desc =
1359	!Config Def = 0
1360	!Config Help =
1361	!
1362	tau_trig_shallow_omp = 600
1363	call getin('tau_trig_shallow',tau_trig_shallow_omp)
1364
1365	!Config Key = tau_trig_deep
1366	!Config Desc =
1367	!Config Def = 0
1368	!Config Help =
1369	!
1370	tau_trig_deep_omp = 1800
1371	call getin('tau_trig_deep',tau_trig_deep_omp)
1372
1373	!Config Key = iflag_clos_bl
1374	!Config Desc =
1375	!Config Def = 0
1376	!Config Help =
1377	!
1378	iflag_clos_bl_omp = 0
1379	call getin('iflag_clos_bl',iflag_clos_bl_omp)
1380
1381	! fin nrlmd le 10/04/2012
1382
1383	!
1384	!Config Key = tau_thermals
1385	!Config Desc =
1386	!Config Def = 0.
1387	!Config Help =
1388	!
1389	tau_thermals_omp = 0.
1390	call getin('tau_thermals',tau_thermals_omp)
1391
1392	!
1393	!Config Key = iflag_coupl
1394	!Config Desc =
1395	!Config Def = 0
1396	!Config Help =
1397	!
1398	iflag_coupl_omp = 0
1399	call getin('iflag_coupl',iflag_coupl_omp)
1400
1401	!
1402	!Config Key = iflag_clos
1403	!Config Desc =
1404	!Config Def = 0
1405	!Config Help =
1406	!
1407	iflag_clos_omp = 1
1408	call getin('iflag_clos',iflag_clos_omp)
1409	!
1410	!Config Key = coef_clos_ls
1411	!Config Desc =
1412	!Config Def = 0
1413	!Config Help =
1414	!
1415	coef_clos_ls_omp = 0.
1416	call getin('coef_clos_ls',coef_clos_ls_omp)
1417
1418	!
1419	!Config Key = iflag_cvl_sigd
1420	!Config Desc =
1421	!Config Def = 0
1422	!Config Help =
1423	!
1424	iflag_cvl_sigd_omp = 0
1425	call getin('iflag_cvl_sigd',iflag_cvl_sigd_omp)
1426
1427	!Config Key = iflag_wake
1428	!Config Desc =
1429	!Config Def = 0
1430	!Config Help =
1431	!
1432	iflag_wake_omp = 0
1433	call getin('iflag_wake',iflag_wake_omp)
1434
1435	!Config Key = alp_offset
1436	!Config Desc =
1437	!Config Def = 0
1438	!Config Help =
1439	!
1440	alp_offset_omp = 0.
1441	call getin('alp_offset',alp_offset_omp)
1442
1443	!
1444	!Config Key = lev_histhf
1445	!Config Desc =
1446	!Config Def = 1
1447	!Config Help =
1448	!
1449	lev_histhf_omp = 1
1450	call getin('lev_histhf',lev_histhf_omp)
1451
1452	!
1453	!Config Key = lev_histday
1454	!Config Desc =
1455	!Config Def = 1
1456	!Config Help =
1457	!
1458	lev_histday_omp = 1
1459	call getin('lev_histday',lev_histday_omp)
1460
1461	!
1462	!Config Key = lev_histmth
1463	!Config Desc =
1464	!Config Def = 2
1465	!Config Help =
1466	!
1467	lev_histmth_omp = 2
1468	call getin('lev_histmth',lev_histmth_omp)
1469	!
1470	!Config Key = lev_histins
1471	!Config Desc =
1472	!Config Def = 1
1473	!Config Help =
1474	!
1475	lev_histins_omp = 1
1476	call getin('lev_histins',lev_histins_omp)
1477	!
1478	!Config Key = lev_histLES
1479	!Config Desc =
1480	!Config Def = 1
1481	!Config Help =
1482	!
1483	lev_histLES_omp = 1
1484	call getin('lev_histLES',lev_histLES_omp)
1485	!
1486	!Config Key = lev_histdayNMC
1487	!Config Desc =
1488	!Config Def = 8
1489	!Config Help =
1490	!
1491	lev_histdayNMC_omp = 8
1492	call getin('lev_histdayNMC',lev_histdayNMC_omp)
1493	!
1494	!Config Key = levout_histNMC
1495	!Config Desc =
1496	!Config Def = 5
1497	!Config Help =
1498	!
1499	levout_histNMC_omp(1) = 5
1500	levout_histNMC_omp(2) = 5
1501	levout_histNMC_omp(3) = 5
1502	call getin('levout_histNMC',levout_histNMC_omp)
1503	!
1504	!histNMC BEG
1505	!Config Key = ok_histNMC
1506	!Config Desc = ok_histNMC(1) = frequence de sortie fichiers histmthNMC
1507	!Config Desc = ok_histNMC(2) = frequence de sortie fichiers histdayNMC
1508	!Config Desc = ok_histNMC(3) = frequence de sortie fichiers histhfNMC
1509	!Config Def = n, n, n
1510	!Config Help =
1511	!
1512	ok_histNMC_omp(1) = .false.
1513	ok_histNMC_omp(2) = .false.
1514	ok_histNMC_omp(3) = .false.
1515	call getin('ok_histNMC',ok_histNMC_omp)
1516	!
1517	!Config Key = freq_outNMC
1518	!Config Desc = freq_outNMC(1) = frequence de sortie fichiers histmthNMC
1519	!Config Desc = freq_outNMC(2) = frequence de sortie fichiers histdayNMC
1520	!Config Desc = freq_outNMC(3) = frequence de sortie fichiers histhfNMC
1521	!Config Def = 2592000., 86400., 21600. (1mois, 1jour, 6h)
1522	!Config Help =
1523	!
1524	freq_outNMC_omp(1) = mth_len
1525	freq_outNMC_omp(2) = 1.
1526	freq_outNMC_omp(3) = 1./4.
1527	call getin('freq_outNMC',freq_outNMC_omp)
1528	!
1529	!Config Key = freq_calNMC
1530	!Config Desc = freq_calNMC(1) = frequence de calcul fichiers histmthNMC
1531	!Config Desc = freq_calNMC(2) = frequence de calcul fichiers histdayNMC
1532	!Config Desc = freq_calNMC(3) = frequence de calcul fichiers histhfNMC
1533	!Config Def = pasphys
1534	!Config Help =
1535	!
1536	freq_calNMC_omp(1) = pasphys
1537	freq_calNMC_omp(2) = pasphys
1538	freq_calNMC_omp(3) = pasphys
1539	call getin('freq_calNMC',freq_calNMC_omp)
1540	!
1541	!Config Key = type_run
1542	!Config Desc =
1543	!Config Def = 'AMIP'/'CFMIP' ou 'CLIM'/'ENSP'
1544	!Config Help =
1545	!
1546	type_run_omp = 'AMIP'
1547	call getin('type_run',type_run_omp)
1548
1549	!
1550	!Config Key = ok_cosp
1551	!Config Desc =
1552	!Config Def = .false.
1553	!Config Help =
1554	!
1555	ok_cosp_omp = .false.
1556	call getin('ok_cosp',ok_cosp_omp)
1557
1558	!
1559	!Config Key = ok_mensuelCOSP
1560	!Config Desc =
1561	!Config Def = .true.
1562	!Config Help =
1563	!
1564	ok_mensuelCOSP_omp = .true.
1565	call getin('ok_mensuelCOSP',ok_mensuelCOSP_omp)
1566
1567	!
1568	!Config Key = ok_journeCOSP
1569	!Config Desc =
1570	!Config Def = .true.
1571	!Config Help =
1572	!
1573	ok_journeCOSP_omp = .true.
1574	call getin('ok_journeCOSP',ok_journeCOSP_omp)
1575
1576	!
1577	!Config Key = ok_hfCOSP
1578	!Config Desc =
1579	!Config Def = .false.
1580	!Config Help =
1581	!
1582	ok_hfCOSP_omp = .false.
1583	call getin('ok_hfCOSP',ok_hfCOSP_omp)
1584
1585	!
1586	! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone
1587	! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
1588	!
1589	!Config Key = lonmin_ins
1590	!Config Desc = 100.
1591	!Config Def = longitude minimale sorties "bilKP_ins"
1592	!Config Help =
1593	!
1594	lonmin_ins_omp = 100.
1595	call getin('lonmin_ins',lonmin_ins_omp)
1596	!
1597	!Config Key = lonmax_ins
1598	!Config Desc = 130.
1599	!Config Def = longitude maximale sorties "bilKP_ins"
1600	!Config Help =
1601	!
1602	lonmax_ins_omp = 130.
1603	call getin('lonmax_ins',lonmax_ins_omp)
1604	!
1605	!Config Key = latmin_ins
1606	!Config Desc = -20.
1607	!Config Def = latitude minimale sorties "bilKP_ins"
1608	!Config Help =
1609	!
1610	latmin_ins_omp = -20.
1611	call getin('latmin_ins',latmin_ins_omp)
1612	!
1613	!Config Key = latmax_ins
1614	!Config Desc = 20.
1615	!Config Def = latitude maximale sorties "bilKP_ins"
1616	!Config Help =
1617	!
1618	latmax_ins_omp = 20.
1619	call getin('latmax_ins',latmax_ins_omp)
1620	!
1621	!Config Key = ecrit_hf
1622	!Config Desc =
1623	!Config Def = 1./8. !toutes les 3h
1624	!Config Help =
1625	!
1626	ecrit_hf_omp = 1./8.
1627	call getin('ecrit_hf',ecrit_hf_omp)
1628	!
1629	!Config Key = ecrit_ins
1630	!Config Desc =
1631	!Config Def = 1./48. ! toutes les 1/2 h
1632	!Config Help =
1633	!
1634	ecrit_ins_omp = 1./48.
1635	call getin('ecrit_ins',ecrit_ins_omp)
1636	!
1637	!Config Key = ecrit_day
1638	!Config Desc =
1639	!Config Def = 1.0 !tous les jours
1640	!Config Help = nombre de jours pour ecriture fichier histday.nc
1641	!
1642	ecrit_day_omp = 1.0
1643	call getin('ecrit_day',ecrit_day_omp)
1644	!
1645	!Config Key = ecrit_mth
1646	!Config Desc =
1647	!Config Def = 30. !tous les 30jours (1 fois par mois)
1648	!Config Help =
1649	!
1650	ecrit_mth_omp = 30.
1651	call getin('ecrit_mth',ecrit_mth_omp)
1652	!
1653	!Config Key = ecrit_tra
1654	!Config Desc =
1655	!Config Def = 30. !tous les 30jours (1 fois par mois)
1656	!Config Help =
1657	!
1658	ecrit_tra_omp = 0.
1659	call getin('ecrit_tra',ecrit_tra_omp)
1660	!
1661	!Config Key = ecrit_reg
1662	!Config Desc =
1663	!Config Def = 0.25 !4 fois par jour
1664	!Config Help =
1665	!
1666	ecrit_reg_omp = 0.25 !4 fois par jour
1667	call getin('ecrit_reg',ecrit_reg_omp)
1668	!
1669	!
1670	print*,'CONFPHYS OOK avant drag_ter'
1671	!
1672	! PARAMETRES CDRAG
1673	!
1674	f_cdrag_ter_omp = 0.8
1675	call getin('f_cdrag_ter',f_cdrag_ter_omp)
1676	!
1677	f_cdrag_oce_omp = 0.8
1678	call getin('f_cdrag_oce',f_cdrag_oce_omp)
1679	!
1680
1681	! Gustiness flags
1682	f_z0qh_oce_omp = 1.
1683	call getin('f_z0qh_oce',f_z0qh_oce_omp)
1684	!
1685	f_qsat_oce_omp = 1.
1686	call getin('f_qsat_oce',f_qsat_oce_omp)
1687	!
1688	f_gust_bl_omp = 0.
1689	call getin('f_gust_bl',f_gust_bl_omp)
1690	!
1691	f_gust_wk_omp = 0.
1692	call getin('f_gust_wk',f_gust_wk_omp)
1693	!
1694	iflag_z0_oce_omp=0
1695	call getin('iflag_z0_oce',iflag_z0_oce_omp)
1696	!
1697	iflag_gusts_omp=0
1698	call getin('iflag_gusts',iflag_gusts_omp)
1699	!
1700	min_wind_speed_omp = 1.
1701	call getin('min_wind_speed',min_wind_speed_omp)
1702
1703	z0m_seaice_omp = 0.002 ; call getin('z0m_seaice',z0m_seaice_omp)
1704	z0h_seaice_omp = 0.002 ; call getin('z0h_seaice',z0h_seaice_omp)
1705
1706	f_rugoro_omp = 0.
1707	call getin('f_rugoro',f_rugoro_omp)
1708
1709	z0min_omp = 0.000015
1710	call getin('z0min',z0min_omp)
1711
1712
1713	! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
1714	!
1715	!Config Key = supcrit1
1716	!Config Desc =
1717	!Config Def = .540
1718	!Config Help =
1719	!
1720	supcrit1_omp = .540
1721	call getin('supcrit1',supcrit1_omp)
1722
1723	!
1724	!Config Key = supcrit2
1725	!Config Desc =
1726	!Config Def = .600
1727	!Config Help =
1728	!
1729	supcrit2_omp = .600
1730	call getin('supcrit2',supcrit2_omp)
1731
1732	!
1733	! PARAMETERS FOR THE MIXING DISTRIBUTION
1734	! iflag_mix: 0=OLD,
1735	! 1=NEW (JYG),
1736	! 2=NEW + conv. depth inhib. by tropos. dryness
1737	! '2' is NOT operationnal and should not be used.
1738	!
1739	!Config Key = iflag_mix
1740	!Config Desc =
1741	!Config Def = 1
1742	!Config Help =
1743	!
1744	iflag_mix_omp = 1
1745	call getin('iflag_mix',iflag_mix_omp)
1746
1747	!
1748	!Config Key = scut
1749	!Config Desc =
1750	!Config Def = 0.95
1751	!Config Help =
1752	!
1753	scut_omp = 0.95
1754	call getin('scut',scut_omp)
1755
1756	!
1757	!Config Key = qqa1
1758	!Config Desc =
1759	!Config Def = 1.0
1760	!Config Help =
1761	!
1762	qqa1_omp = 1.0
1763	call getin('qqa1',qqa1_omp)
1764
1765	!
1766	!Config Key = qqa2
1767	!Config Desc =
1768	!Config Def = 0.0
1769	!Config Help =
1770	!
1771	qqa2_omp = 0.0
1772	call getin('qqa2',qqa2_omp)
1773
1774	!
1775	!Config Key = gammas
1776	!Config Desc =
1777	!Config Def = 0.05
1778	!Config Help =
1779	!
1780	gammas_omp = 0.05
1781	call getin('gammas',gammas_omp)
1782
1783	!
1784	!Config Key = Fmax
1785	!Config Desc =
1786	!Config Def = 0.65
1787	!Config Help =
1788	!
1789	Fmax_omp = 0.65
1790	call getin('Fmax',Fmax_omp)
1791
1792	!
1793	!Config Key = tmax_fonte_cv
1794	!Config Desc =
1795	!Config Def = 275.15
1796	!Config Help =
1797	!
1798	tmax_fonte_cv_omp = 275.15
1799	call getin('tmax_fonte_cv',tmax_fonte_cv_omp)
1800
1801	!
1802	!Config Key = alphas
1803	!Config Desc =
1804	!Config Def = -5.
1805	!Config Help =
1806	!
1807	alphas_omp = -5.
1808	call getin('alphas',alphas_omp)
1809
1810	!Config key = ok_strato
1811	!Config Desc = activation de la version strato
1812	!Config Def = .FALSE.
1813	!Config Help = active la version stratosph\'erique de LMDZ de F. Lott
1814
1815	ok_strato_omp=.FALSE.
1816	CALL getin('ok_strato',ok_strato_omp)
1817
1818	!Config key = ok_hines
1819	!Config Desc = activation de la parametrisation de hines
1820	!Config Def = .FALSE.
1821	!Config Help = Clefs controlant la parametrization de Hines
1822	! Et la sponge layer (Runs Stratospheriques)
1823
1824	ok_hines_omp=.FALSE.
1825	CALL getin('ok_hines',ok_hines_omp)
1826
1827	! Random gravity waves:
1828
1829	ok_gwd_rando_omp = .FALSE.
1830	IF ( klon_glo == 1 ) THEN
1831	print*,'La parametrisation des ondes de gravites non orographiques'
1832	print*,'ne fonctionne pas en 1D'
1833	ELSE
1834	CALL getin('ok_gwd_rando', ok_gwd_rando_omp)
1835	ENDIF
1836
1837	GWD_RANDO_RUWMAX_omp = 0.01
1838	CALL getin('gwd_rando_ruwmax', GWD_RANDO_RUWMAX_omp)
1839
1840	gwd_rando_sat_omp = 0.25
1841	CALL getin('gwd_rando_sat', gwd_rando_sat_omp)
1842
1843	!Config key = ok_qch4
1844	!Config Desc = activation de la parametrisation du methane
1845	!Config Def = .FALSE.
1846	!Config Help = Clef controlant l'activation de la parametrisation
1847	! de l'humidite due a oxydation+photolyse du methane strato
1848
1849	ok_qch4_omp=.FALSE.
1850	CALL getin('ok_qch4',ok_qch4_omp)
1851
1852	!Config Key = OK_LES
1853	!Config Desc = Pour des sorties LES
1854	!Config Def = .false.
1855	!Config Help = Pour creer le fichier histLES contenant les sorties
1856	! LES
1857	!
1858	ok_LES_omp = .false.
1859	call getin('OK_LES', ok_LES_omp)
1860
1861	!Config Key = callstats
1862	!Config Desc = Pour des sorties callstats
1863	!Config Def = .false.
1864	!Config Help = Pour creer le fichier stats contenant les sorties
1865	! stats
1866	!
1867	callstats_omp = .false.
1868	call getin('callstats', callstats_omp)
1869	!
1870	!Config Key = ecrit_LES
1871	!Config Desc = Frequence d'ecriture des resultats du LES en nombre de jours;
1872	! par defaut 1., i.e. 1 jour
1873	!Config Def = 1./8.
1874	!Config Help = ...
1875	!
1876	!
1877	ecrit_LES_omp = 1./8.
1878	call getin('ecrit_LES', ecrit_LES_omp)
1879	!
1880	read_climoz = 0 ! default value
1881	call getin('read_climoz', read_climoz)
1882
1883	carbon_cycle_tr_omp=.FALSE.
1884	CALL getin('carbon_cycle_tr',carbon_cycle_tr_omp)
1885
1886	carbon_cycle_cpl_omp=.FALSE.
1887	CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp)
1888
1889	!$OMP END MASTER
1890	!$OMP BARRIER
1891
1892	R_ecc = R_ecc_omp
1893	R_peri = R_peri_omp
1894	R_incl = R_incl_omp
1895	solaire = solaire_omp
1896	co2_ppm = co2_ppm_omp
1897	RCO2 = RCO2_omp
1898	CH4_ppb = CH4_ppb_omp
1899	RCH4 = RCH4_omp
1900	N2O_ppb = N2O_ppb_omp
1901	RN2O = RN2O_omp
1902	CFC11_ppt = CFC11_ppt_omp
1903	RCFC11 = RCFC11_omp
1904	CFC12_ppt = CFC12_ppt_omp
1905	RCFC12 = RCFC12_omp
1906	RCO2_act = RCO2
1907	RCH4_act = RCH4
1908	RN2O_act = RN2O
1909	RCFC11_act = RCFC11
1910	RCFC12_act = RCFC12
1911	ok_4xCO2atm = ok_4xCO2atm_omp
1912	RCO2_per = RCO2_per_omp
1913	RCH4_per = RCH4_per_omp
1914	RN2O_per = RN2O_per_omp
1915	RCFC11_per = RCFC11_per_omp
1916	RCFC12_per = RCFC12_per_omp
1917
1918	cycle_diurne = cycle_diurne_omp
1919	soil_model = soil_model_omp
1920	new_oliq = new_oliq_omp
1921	ok_orodr = ok_orodr_omp
1922	ok_orolf = ok_orolf_omp
1923	ok_limitvrai = ok_limitvrai_omp
1924	nbapp_rad = nbapp_rad_omp
1925	iflag_con = iflag_con_omp
1926	iflag_ener_conserv = iflag_ener_conserv_omp
1927	ok_conserv_q = ok_conserv_q_omp
1928	iflag_fisrtilp_qsat = iflag_fisrtilp_qsat_omp
1929
1930	epmax = epmax_omp
1931	ok_adj_ema = ok_adj_ema_omp
1932	iflag_clw = iflag_clw_omp
1933	cld_lc_lsc = cld_lc_lsc_omp
1934	cld_lc_con = cld_lc_con_omp
1935	cld_tau_lsc = cld_tau_lsc_omp
1936	cld_tau_con = cld_tau_con_omp
1937	ffallv_lsc = ffallv_lsc_omp
1938	ffallv_con = ffallv_con_omp
1939	coef_eva = coef_eva_omp
1940	reevap_ice = reevap_ice_omp
1941	iflag_pdf = iflag_pdf_omp
1942	solarlong0 = solarlong0_omp
1943	qsol0 = qsol0_omp
1944	evap0 = evap0_omp
1945	albsno0 = albsno0_omp
1946	inertie_sol = inertie_sol_omp
1947	inertie_ice = inertie_ice_omp
1948	inertie_sno = inertie_sno_omp
1949	rad_froid = rad_froid_omp
1950	rad_chau1 = rad_chau1_omp
1951	rad_chau2 = rad_chau2_omp
1952	t_glace_min = t_glace_min_omp
1953	t_glace_max = t_glace_max_omp
1954	exposant_glace = exposant_glace_omp
1955	iflag_t_glace = iflag_t_glace_omp
1956	iflag_ice_thermo = iflag_ice_thermo_omp
1957	rei_min = rei_min_omp
1958	rei_max = rei_max_omp
1959	top_height = top_height_omp
1960	overlap = overlap_omp
1961	cdmmax = cdmmax_omp
1962	cdhmax = cdhmax_omp
1963	ksta = ksta_omp
1964	ksta_ter = ksta_ter_omp
1965	f_ri_cd_min = f_ri_cd_min_omp
1966	ok_kzmin = ok_kzmin_omp
1967	fmagic = fmagic_omp
1968	pmagic = pmagic_omp
1969	iflag_pbl = iflag_pbl_omp
1970	iflag_pbl_split = iflag_pbl_split_omp
1971	lev_histhf = lev_histhf_omp
1972	lev_histday = lev_histday_omp
1973	lev_histmth = lev_histmth_omp
1974	lev_histins = lev_histins_omp
1975	lev_histLES = lev_histLES_omp
1976	lev_histdayNMC = lev_histdayNMC_omp
1977	levout_histNMC = levout_histNMC_omp
1978	ok_histNMC(:) = ok_histNMC_omp(:)
1979	freq_outNMC(:) = freq_outNMC_omp(:)
1980	freq_calNMC(:) = freq_calNMC_omp(:)
1981
1982	type_ocean = type_ocean_omp
1983	version_ocean = version_ocean_omp
1984	t_coupl = t_coupl_omp
1985
1986	ok_veget=.true.
1987	type_veget=type_veget_omp
1988	if (type_veget=='n' .or. type_veget=='bucket' .or. type_veget=='betaclim') &
1989	then
1990	ok_veget=.false.
1991	endif
1992	! Martin
1993	ok_snow = ok_snow_omp
1994	! Martin
1995
1996	ok_all_xml = ok_all_xml_omp
1997	ok_newmicro = ok_newmicro_omp
1998	ok_journe = ok_journe_omp
1999	ok_hf = ok_hf_omp
2000	ok_mensuel = ok_mensuel_omp
2001	ok_instan = ok_instan_omp
2002	freq_ISCCP = freq_ISCCP_omp
2003	ecrit_ISCCP = ecrit_ISCCP_omp
2004	freq_COSP = freq_COSP_omp
2005	ok_ade = ok_ade_omp
2006	ok_aie = ok_aie_omp
2007	ok_cdnc = ok_cdnc_omp
2008	aerosol_couple = aerosol_couple_omp
2009	flag_aerosol=flag_aerosol_omp
2010	flag_aerosol_strat=flag_aerosol_strat_omp
2011	new_aod=new_aod_omp
2012	aer_type = aer_type_omp
2013	bl95_b0 = bl95_b0_omp
2014	bl95_b1 = bl95_b1_omp
2015	fact_cldcon = fact_cldcon_omp
2016	facttemps = facttemps_omp
2017	ratqsbas = ratqsbas_omp
2018	ratqshaut = ratqshaut_omp
2019	tau_ratqs = tau_ratqs_omp
2020
2021	iflag_radia = iflag_radia_omp
2022	iflag_rrtm = iflag_rrtm_omp
2023	iflag_albedo = iflag_albedo_omp
2024	ok_chlorophyll = ok_chlorophyll_omp
2025	NSW = NSW_omp
2026	iflag_cld_th = iflag_cld_th_omp
2027	iflag_cld_cv = iflag_cld_cv_omp
2028	tau_cld_cv = tau_cld_cv_omp
2029	coefw_cld_cv = coefw_cld_cv_omp
2030	iflag_ratqs = iflag_ratqs_omp
2031	ip_ebil_phy = ip_ebil_phy_omp
2032	iflag_thermals = iflag_thermals_omp
2033	iflag_thermals_ed = iflag_thermals_ed_omp
2034	fact_thermals_ed_dz = fact_thermals_ed_dz_omp
2035	iflag_thermals_optflux = iflag_thermals_optflux_omp
2036	iflag_thermals_closure = iflag_thermals_closure_omp
2037	nsplit_thermals = nsplit_thermals_omp
2038	tau_thermals = tau_thermals_omp
2039	alp_bl_k = alp_bl_k_omp
2040	! nrlmd le 10/04/2012
2041	iflag_trig_bl = iflag_trig_bl_omp
2042	s_trig = s_trig_omp
2043	tau_trig_shallow = tau_trig_shallow_omp
2044	tau_trig_deep = tau_trig_deep_omp
2045	iflag_clos_bl = iflag_clos_bl_omp
2046	! fin nrlmd le 10/04/2012
2047	iflag_coupl = iflag_coupl_omp
2048	iflag_clos = iflag_clos_omp
2049	iflag_wake = iflag_wake_omp
2050	coef_clos_ls = coef_clos_ls_omp
2051	alp_offset = alp_offset_omp
2052	iflag_cvl_sigd = iflag_cvl_sigd_omp
2053	type_run = type_run_omp
2054	ok_cosp = ok_cosp_omp
2055	ok_mensuelCOSP = ok_mensuelCOSP_omp
2056	ok_journeCOSP = ok_journeCOSP_omp
2057	ok_hfCOSP = ok_hfCOSP_omp
2058	seuil_inversion=seuil_inversion_omp
2059	lonmin_ins = lonmin_ins_omp
2060	lonmax_ins = lonmax_ins_omp
2061	latmin_ins = latmin_ins_omp
2062	latmax_ins = latmax_ins_omp
2063	ecrit_hf = ecrit_hf_omp
2064	ecrit_ins = ecrit_ins_omp
2065	ecrit_day = ecrit_day_omp
2066	ecrit_mth = ecrit_mth_omp
2067	ecrit_tra = ecrit_tra_omp
2068	ecrit_reg = ecrit_reg_omp
2069	cvl_comp_threshold = cvl_comp_threshold_omp
2070	cvl_sig2feed = cvl_sig2feed_omp
2071	cvl_corr = cvl_corr_omp
2072	ok_lic_melt = ok_lic_melt_omp
2073	f_cdrag_ter=f_cdrag_ter_omp
2074	f_cdrag_oce=f_cdrag_oce_omp
2075
2076	f_gust_wk=f_gust_wk_omp
2077	f_gust_bl=f_gust_bl_omp
2078	f_qsat_oce=f_qsat_oce_omp
2079	f_z0qh_oce=f_z0qh_oce_omp
2080	min_wind_speed=min_wind_speed_omp
2081	iflag_gusts=iflag_gusts_omp
2082	iflag_z0_oce=iflag_z0_oce_omp
2083
2084
2085	z0m_seaice=z0m_seaice_omp
2086	z0h_seaice=z0h_seaice_omp
2087
2088	f_rugoro=f_rugoro_omp
2089
2090	z0min=z0min_omp
2091	supcrit1 = supcrit1_omp
2092	supcrit2 = supcrit2_omp
2093	iflag_mix = iflag_mix_omp
2094	scut = scut_omp
2095	qqa1 = qqa1_omp
2096	qqa2 = qqa2_omp
2097	gammas = gammas_omp
2098	Fmax = Fmax_omp
2099	tmax_fonte_cv = tmax_fonte_cv_omp
2100	alphas = alphas_omp
2101	ok_strato = ok_strato_omp
2102	ok_hines = ok_hines_omp
2103	ok_gwd_rando = ok_gwd_rando_omp
2104	GWD_RANDO_RUWMAX = GWD_RANDO_RUWMAX_omp
2105	gwd_rando_sat = gwd_rando_sat_omp
2106	ok_qch4 = ok_qch4_omp
2107	ok_LES = ok_LES_omp
2108	callstats = callstats_omp
2109	ecrit_LES = ecrit_LES_omp
2110	carbon_cycle_tr = carbon_cycle_tr_omp
2111	carbon_cycle_cpl = carbon_cycle_cpl_omp
2112
2113	! Test of coherence between type_ocean and version_ocean
2114	IF (type_ocean=='couple' .AND. (version_ocean/='opa8' .AND. version_ocean/='nemo') ) THEN
2115	WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration'
2116	CALL abort_physic('conf_phys','version_ocean not valid',1)
2117	END IF
2118
2119	IF (type_ocean=='slab' .AND. version_ocean=='xxxxxx') THEN
2120	version_ocean='sicOBS'
2121	ELSE IF (type_ocean=='slab' .AND. version_ocean/='sicOBS' &
2122	.AND. version_ocean/='sicINT' .AND. version_ocean/='sicNO') THEN
2123	WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
2124	CALL abort_physic('conf_phys','version_ocean not valid',1)
2125	END IF
2126
2127	! Test sur new_aod. Ce flag permet de retrouver les resultats de l'AR4
2128	! il n'est utilisable que lors du couplage avec le SO4 seul
2129	IF (ok_ade .OR. ok_aie) THEN
2130	IF ( flag_aerosol .EQ. 0 ) THEN
2131	CALL abort_physic('conf_phys','flag_aerosol=0 not compatible avec ok_ade ou ok_aie=.TRUE.',1)
2132	END IF
2133	IF ( .NOT. new_aod .AND. flag_aerosol .NE. 1) THEN
2134	CALL abort_physic('conf_phys','new_aod=.FALSE. not compatible avec flag_aerosol=1',1)
2135	END IF
2136	END IF
2137
2138	! ok_cdnc must be set to y if ok_aie is activated
2139	IF (ok_aie .AND. .NOT. ok_cdnc) THEN
2140	CALL abort_physic('conf_phys', 'ok_cdnc must be set to y if ok_aie is activated',1)
2141	ENDIF
2142
2143	!$OMP MASTER
2144
2145	write(lunout,*)' ##############################################'
2146	write(lunout,*)' Configuration des parametres de la physique: '
2147	write(lunout,*)' Type ocean = ', type_ocean
2148	write(lunout,*)' Version ocean = ', version_ocean
2149	write(lunout,*)' Config veget = ', ok_veget,type_veget
2150	write(lunout,*)' Snow model SISVAT : ok_snow = ', ok_snow
2151	write(lunout,*)' Config xml pour XIOS : ok_all_xml = ', ok_all_xml
2152	write(lunout,*)' Sortie journaliere = ', ok_journe
2153	write(lunout,*)' Sortie haute frequence = ', ok_hf
2154	write(lunout,*)' Sortie mensuelle = ', ok_mensuel
2155	write(lunout,*)' Sortie instantanee = ', ok_instan
2156	write(lunout,*)' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
2157	write(lunout,*)' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
2158	write(lunout,*)' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
2159	write(lunout,*)' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
2160	write(lunout,*)' Excentricite = ',R_ecc
2161	write(lunout,*)' Equinoxe = ',R_peri
2162	write(lunout,*)' Inclinaison =',R_incl
2163	write(lunout,*)' Constante solaire =',solaire
2164	write(lunout,*)' co2_ppm =',co2_ppm
2165	write(lunout,*)' RCO2_act = ',RCO2_act
2166	write(lunout,*)' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
2167	write(lunout,*)' N2O_ppb =',N2O_ppb,' RN2O_act= ',RN2O_act
2168	write(lunout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11_act= ',RCFC11_act
2169	write(lunout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12_act= ',RCFC12_act
2170	write(lunout,*)' ok_4xCO2atm=',ok_4xCO2atm
2171	write(lunout,*)' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
2172	write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
2173	write(lunout,*)' RCFC12_per = ',RCFC12_per
2174	write(lunout,*)' cvl_comp_threshold=', cvl_comp_threshold
2175	write(lunout,*)' cvl_sig2feed=', cvl_sig2feed
2176	write(lunout,*)' cvl_corr=', cvl_corr
2177	write(lunout,*)'ok_lic_melt=', ok_lic_melt
2178	write(lunout,*)'cycle_diurne=',cycle_diurne
2179	write(lunout,*)'soil_model=',soil_model
2180	write(lunout,*)'new_oliq=',new_oliq
2181	write(lunout,*)'ok_orodr=',ok_orodr
2182	write(lunout,*)'ok_orolf=',ok_orolf
2183	write(lunout,*)'ok_limitvrai=',ok_limitvrai
2184	write(lunout,*)'nbapp_rad=',nbapp_rad
2185	write(lunout,*)'iflag_con=',iflag_con
2186	write(lunout,*)'iflag_ener_conserv=',iflag_ener_conserv
2187	write(lunout,*)'ok_conserv_q=',ok_conserv_q
2188	write(lunout,*)'iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
2189	write(lunout,*)' epmax = ', epmax
2190	write(lunout,*)' ok_adj_ema = ', ok_adj_ema
2191	write(lunout,*)' iflag_clw = ', iflag_clw
2192	write(lunout,*)' cld_lc_lsc = ', cld_lc_lsc
2193	write(lunout,*)' cld_lc_con = ', cld_lc_con
2194	write(lunout,*)' cld_tau_lsc = ', cld_tau_lsc
2195	write(lunout,*)' cld_tau_con = ', cld_tau_con
2196	write(lunout,*)' ffallv_lsc = ', ffallv_lsc
2197	write(lunout,*)' ffallv_con = ', ffallv_con
2198	write(lunout,*)' coef_eva = ', coef_eva
2199	write(lunout,*)' reevap_ice = ', reevap_ice
2200	write(lunout,*)' iflag_pdf = ', iflag_pdf
2201	write(lunout,*)' iflag_cld_th = ', iflag_cld_th
2202	write(lunout,*)' iflag_cld_cv = ', iflag_cld_cv
2203	write(lunout,*)' tau_cld_cv = ', tau_cld_cv
2204	write(lunout,*)' coefw_cld_cv = ', coefw_cld_cv
2205	write(lunout,*)' iflag_radia = ', iflag_radia
2206	write(lunout,*)' iflag_rrtm = ', iflag_rrtm
2207	write(lunout,*)' NSW = ', NSW
2208	write(lunout,*)' iflag_albedo = ', iflag_albedo !albedo SB
2209	write(lunout,*)' ok_chlorophyll =',ok_chlorophyll ! albedo SB
2210	write(lunout,*)' iflag_ratqs = ', iflag_ratqs
2211	write(lunout,*)' seuil_inversion = ', seuil_inversion
2212	write(lunout,*)' fact_cldcon = ', fact_cldcon
2213	write(lunout,*)' facttemps = ', facttemps
2214	write(lunout,*)' ok_newmicro = ',ok_newmicro
2215	write(lunout,*)' ratqsbas = ',ratqsbas
2216	write(lunout,*)' ratqshaut = ',ratqshaut
2217	write(lunout,*)' tau_ratqs = ',tau_ratqs
2218	write(lunout,*)' top_height = ',top_height
2219	write(lunout,*)' rad_froid = ',rad_froid
2220	write(lunout,*)' rad_chau1 = ',rad_chau1
2221	write(lunout,*)' rad_chau2 = ',rad_chau2
2222	write(lunout,*)' t_glace_min = ',t_glace_min
2223	write(lunout,*)' t_glace_max = ',t_glace_max
2224	write(lunout,*)' exposant_glace = ',exposant_glace
2225	write(lunout,*)' iflag_t_glace = ',iflag_t_glace
2226	write(lunout,*)' iflag_ice_thermo = ',iflag_ice_thermo
2227	write(lunout,*)' rei_min = ',rei_min
2228	write(lunout,*)' rei_max = ',rei_max
2229	write(lunout,*)' overlap = ',overlap
2230	write(lunout,*)' cdmmax = ',cdmmax
2231	write(lunout,*)' cdhmax = ',cdhmax
2232	write(lunout,*)' ksta = ',ksta
2233	write(lunout,*)' ksta_ter = ',ksta_ter
2234	write(lunout,*)' f_ri_cd_min = ',f_ri_cd_min
2235	write(lunout,*)' ok_kzmin = ',ok_kzmin
2236	write(lunout,*)' fmagic = ',fmagic
2237	write(lunout,*)' pmagic = ',pmagic
2238	write(lunout,*)' ok_ade = ',ok_ade
2239	write(lunout,*)' ok_aie = ',ok_aie
2240	write(lunout,*)' aerosol_couple = ', aerosol_couple
2241	write(lunout,*)' flag_aerosol = ', flag_aerosol
2242	write(lunout,*)' flag_aerosol_strat = ', flag_aerosol_strat
2243	write(lunout,*)' new_aod = ', new_aod
2244	write(lunout,*)' aer_type = ',aer_type
2245	write(lunout,*)' bl95_b0 = ',bl95_b0
2246	write(lunout,*)' bl95_b1 = ',bl95_b1
2247	write(lunout,*)' lev_histhf = ',lev_histhf
2248	write(lunout,*)' lev_histday = ',lev_histday
2249	write(lunout,*)' lev_histmth = ',lev_histmth
2250	write(lunout,*)' lev_histins = ',lev_histins
2251	write(lunout,*)' lev_histLES = ',lev_histLES
2252	write(lunout,*)' lev_histdayNMC = ',lev_histdayNMC
2253	write(lunout,*)' levout_histNMC = ',levout_histNMC
2254	write(lunout,*)' ok_histNMC = ',ok_histNMC
2255	write(lunout,*)' freq_outNMC = ',freq_outNMC
2256	write(lunout,*)' freq_calNMC = ',freq_calNMC
2257	write(lunout,*)' iflag_pbl = ', iflag_pbl
2258	write(lunout,*)' iflag_pbl_split = ', iflag_pbl_split
2259	write(lunout,*)' iflag_thermals = ', iflag_thermals
2260	write(lunout,*)' iflag_thermals_ed = ', iflag_thermals_ed
2261	write(lunout,*)' fact_thermals_ed_dz = ', fact_thermals_ed_dz
2262	write(lunout,*)' iflag_thermals_optflux = ', iflag_thermals_optflux
2263	write(lunout,*)' iflag_thermals_closure = ', iflag_thermals_closure
2264	write(lunout,*)' iflag_clos = ', iflag_clos
2265	write(lunout,*)' coef_clos_ls = ', coef_clos_ls
2266	write(lunout,*)' type_run = ',type_run
2267	write(lunout,*)' ok_cosp = ',ok_cosp
2268	write(lunout,*)' ok_mensuelCOSP = ',ok_mensuelCOSP
2269	write(lunout,*)' ok_journeCOSP = ',ok_journeCOSP
2270	write(lunout,*)' ok_hfCOSP =',ok_hfCOSP
2271	write(lunout,*)' solarlong0 = ', solarlong0
2272	write(lunout,*)' qsol0 = ', qsol0
2273	write(lunout,*)' evap0 = ', evap0
2274	write(lunout,*)' albsno0 = ', albsno0
2275	write(lunout,*)' inertie_sol = ', inertie_sol
2276	write(lunout,*)' inertie_ice = ', inertie_ice
2277	write(lunout,*)' inertie_sno = ', inertie_sno
2278	write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter
2279	write(lunout,*)' f_cdrag_oce = ',f_cdrag_oce
2280	write(lunout,*)' f_rugoro = ',f_rugoro
2281	write(lunout,*)' z0min = ',z0min
2282	write(lunout,*)' supcrit1 = ', supcrit1
2283	write(lunout,*)' supcrit2 = ', supcrit2
2284	write(lunout,*)' iflag_mix = ', iflag_mix
2285	write(lunout,*)' scut = ', scut
2286	write(lunout,*)' qqa1 = ', qqa1
2287	write(lunout,*)' qqa2 = ', qqa2
2288	write(lunout,*)' gammas = ', gammas
2289	write(lunout,*)' Fmax = ', Fmax
2290	write(lunout,*)' tmax_fonte_cv = ', tmax_fonte_cv
2291	write(lunout,*)' alphas = ', alphas
2292	write(lunout,*)' iflag_wake = ', iflag_wake
2293	write(lunout,*)' alp_offset = ', alp_offset
2294	! nrlmd le 10/04/2012
2295	write(lunout,*)' iflag_trig_bl = ', iflag_trig_bl
2296	write(lunout,*)' s_trig = ', s_trig
2297	write(lunout,*)' tau_trig_shallow = ', tau_trig_shallow
2298	write(lunout,*)' tau_trig_deep = ', tau_trig_deep
2299	write(lunout,*)' iflag_clos_bl = ', iflag_clos_bl
2300	! fin nrlmd le 10/04/2012
2301
2302	write(lunout,*)' lonmin lonmax latmin latmax bilKP_ins =',&
2303	lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
2304	write(lunout,*)' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
2305	ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
2306
2307	write(lunout,*) 'ok_strato = ', ok_strato
2308	write(lunout,*) 'ok_hines = ', ok_hines
2309	write(lunout,*) 'ok_gwd_rando = ', ok_gwd_rando
2310	write(lunout,*) 'ok_qch4 = ', ok_qch4
2311	write(lunout,*) 'GWD_RANDO_RUWMAX = ', GWD_RANDO_RUWMAX
2312	write(lunout,*) 'gwd_rando_sat = ', gwd_rando_sat
2313	write(lunout,*) 'read_climoz = ', read_climoz
2314	write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
2315	write(lunout,*) 'carbon_cycle_cpl = ', carbon_cycle_cpl
2316
2317	!$OMP END MASTER
2318
2319	return
2320
2321	end subroutine conf_phys
2322
2323	end module conf_phys_m
2324	!
2325	!#################################################################
2326	!
2327
2328	subroutine conf_interface(tau_calv)
2329
2330	use IOIPSL
2331	USE print_control_mod, ONLY: lunout
2332	implicit none
2333	! Configuration de l'interace atm/surf
2334	!
2335	! tau_calv: temps de relaxation pour la fonte des glaciers
2336
2337	REAL :: tau_calv
2338	REAL,SAVE :: tau_calv_omp
2339
2340	!
2341	!Config Key = tau_calv
2342	!Config Desc = temps de relaxation pour fonte des glaciers en jours
2343	!Config Def = 1 an
2344	!Config Help =
2345	!
2346	tau_calv_omp = 360.*10.
2347	!$OMP MASTER
2348	call getin('tau_calv',tau_calv_omp)
2349	!$OMP END MASTER
2350	!$OMP BARRIER
2351
2352	tau_calv=tau_calv_omp
2353
2354	!$OMP MASTER
2355	write(lunout,*)' ##############################################'
2356	WRITE(lunout,*)' Configuration de l''interface atm/surfaces : '
2357	WRITE(lunout,*)' tau_calv = ',tau_calv
2358	!$OMP END MASTER
2359
2360	return
2361
2362	end subroutine conf_interface

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format