Context Navigation

conf_phys_m.F90 @ 5408

Last change on this file since 5408 was 4669, checked in by Laurent Fairhead, 16 months ago
Merged with trunk revision 4586 corresponding to june 2023 testing
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: 91.6 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_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat, aerosol_couple, &
20	chemistry_couple, flag_aerosol, flag_aerosol_strat, flag_aer_feedback, &
21	flag_bc_internal_mixture, 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, carbon_cycle_rad, level_coupling_esm
29	USE carbon_cycle_mod, ONLY: read_fco2_ocean_cor,var_fco2_ocean_cor
30	USE carbon_cycle_mod, ONLY: read_fco2_land_cor,var_fco2_land_cor
31	USE mod_grid_phy_lmdz, ONLY: klon_glo
32	USE print_control_mod, ONLY: lunout
33	use config_ocean_skin_m, only: config_ocean_skin
34	USE phys_state_var_mod, ONLY: phys_tstep
35
36	INCLUDE "conema3.h"
37	INCLUDE "nuage.h"
38	INCLUDE "YOMCST.h"
39	INCLUDE "YOMCST2.h"
40	INCLUDE "alpale.h"
41
42	!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
43	INCLUDE "clesphys.h"
44	INCLUDE "compbl.h"
45	INCLUDE "comsoil.h"
46	INCLUDE "YOEGWD.h"
47	!
48	! Configuration de la "physique" de LMDZ a l'aide de la fonction
49	! GETIN de IOIPSL
50	!
51	! LF 05/2001
52	!
53	! type_ocean: type d'ocean (force, slab, couple)
54	! version_ocean: version d'ocean (opa8/nemo pour type_ocean=couple ou
55	! sicOBS,sicINT,sicNO pour type_ocean=slab)
56	! ok_veget: type de modele de vegetation
57	! ok_journe: sorties journalieres
58	! ok_hf: sorties haute frequence
59	! ok_mensuel: sorties mensuelles
60	! ok_instan: sorties instantanees
61	! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
62	! ok_alw: activate aerosol LW effect
63	! ok_cdnc, ok cloud droplet number concentration
64	! flag_aerosol_strat : flag pour les aerosols stratos
65	! flag_bc_internal_mixture : use BC internal mixture if true
66	! bl95_b*: parameters in the formula to link CDNC to aerosol mass conc
67	! ok_volcan: activate volcanic diags (SW heat & LW cool rate, SW & LW flux)
68	! flag_volc_surfstrat: VolMIP flag, activate forcing surface cooling rate (=1), strato heating rate (=2) or nothing (=0, default)
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_alw, ok_cdnc, ok_volcan
78	LOGICAL :: aerosol_couple, chemistry_couple
79	INTEGER :: flag_aerosol
80	INTEGER :: flag_aerosol_strat
81	INTEGER :: flag_volc_surfstrat
82	LOGICAL :: flag_aer_feedback
83	LOGICAL :: flag_bc_internal_mixture
84	REAL :: bl95_b0, bl95_b1
85	REAL :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
86	INTEGER :: iflag_cld_th
87	INTEGER :: iflag_ratqs
88
89	CHARACTER (len = 6), SAVE :: type_ocean_omp, version_ocean_omp, ocean_omp
90	CHARACTER (len = 10),SAVE :: type_veget_omp
91	CHARACTER (len = 8), SAVE :: aer_type_omp
92	INTEGER, SAVE :: landice_opt_omp
93	INTEGER, SAVE :: iflag_tsurf_inlandsis_omp,iflag_temp_inlandsis_omp
94	INTEGER, SAVE :: iflag_albcalc_omp,iflag_z0m_snow_omp
95	LOGICAL, SAVE :: SnoMod_omp,BloMod_omp,ok_outfor_omp,ok_zsn_ii_omp
96	LOGICAL, SAVE :: discret_xf_omp,opt_runoff_ac_omp
97	LOGICAL, SAVE :: is_ok_slush_omp,is_ok_z0h_rn_omp,is_ok_density_kotlyakov_omp
98	REAL, SAVE :: prescribed_z0m_snow_omp,correc_alb_omp
99	REAL, SAVE :: buf_sph_pol_omp,buf_siz_pol_omp
100	LOGICAL, SAVE :: ok_newmicro_omp
101	LOGICAL, SAVE :: ok_all_xml_omp
102	LOGICAL, SAVE :: ok_lwoff_omp
103	LOGICAL, SAVE :: ok_journe_omp, ok_mensuel_omp, ok_instan_omp, ok_hf_omp
104	LOGICAL, SAVE :: ok_LES_omp
105	LOGICAL, SAVE :: callstats_omp
106	LOGICAL, SAVE :: ok_ade_omp, ok_aie_omp, ok_alw_omp, ok_cdnc_omp, ok_volcan_omp
107	LOGICAL, SAVE :: aerosol_couple_omp, chemistry_couple_omp
108	INTEGER, SAVE :: flag_aerosol_omp
109	INTEGER, SAVE :: flag_aerosol_strat_omp
110	INTEGER, SAVE :: flag_volc_surfstrat_omp
111	LOGICAL, SAVE :: flag_aer_feedback_omp
112	LOGICAL, SAVE :: flag_bc_internal_mixture_omp
113	REAL,SAVE :: bl95_b0_omp, bl95_b1_omp
114	REAL,SAVE :: freq_ISCCP_omp, ecrit_ISCCP_omp
115	REAL,SAVE :: freq_COSP_omp, freq_AIRS_omp
116	REAL,SAVE :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
117	REAL,SAVE :: tau_cld_cv_omp, coefw_cld_cv_omp
118	INTEGER, SAVE :: iflag_cld_cv_omp
119
120	REAL, SAVE :: ratqshaut_omp
121	REAL, SAVE :: tau_ratqs_omp
122	REAL, SAVE :: t_coupl_omp
123	INTEGER, SAVE :: iflag_radia_omp
124	INTEGER, SAVE :: iflag_rrtm_omp
125	INTEGER, SAVE :: iflag_albedo_omp !albedo SB
126	LOGICAL, SAVE :: ok_chlorophyll_omp ! albedo SB
127	INTEGER, SAVE :: NSW_omp
128	INTEGER, SAVE :: iflag_cld_th_omp, ip_ebil_phy_omp
129	INTEGER, SAVE :: iflag_ratqs_omp
130
131	REAL, SAVE :: f_cdrag_ter_omp,f_cdrag_oce_omp
132	REAL, SAVE :: f_rugoro_omp , z0min_omp
133	REAL, SAVE :: z0m_seaice_omp,z0h_seaice_omp
134	REAL, SAVE :: z0m_landice_omp,z0h_landice_omp
135	REAL, SAVE :: min_wind_speed_omp,f_gust_wk_omp,f_gust_bl_omp,f_qsat_oce_omp, f_z0qh_oce_omp
136	INTEGER, SAVE :: iflag_gusts_omp,iflag_z0_oce_omp
137
138	REAL :: seuil_inversion
139	REAL,SAVE :: seuil_inversion_omp
140
141	INTEGER,SAVE :: iflag_thermals_omp,nsplit_thermals_omp
142	REAL,SAVE :: tau_thermals_omp,alp_bl_k_omp
143	! nrlmd le 10/04/2012
144	INTEGER,SAVE :: iflag_trig_bl_omp,iflag_clos_bl_omp
145	INTEGER,SAVE :: tau_trig_shallow_omp,tau_trig_deep_omp
146	REAL,SAVE :: s_trig_omp
147	! fin nrlmd le 10/04/2012
148	REAL :: alp_offset
149	REAL, SAVE :: alp_offset_omp
150	INTEGER,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
151	INTEGER,SAVE :: iflag_cvl_sigd_omp
152	REAL, SAVE :: coef_clos_ls_omp
153	REAL, SAVE :: supcrit1_omp, supcrit2_omp
154	INTEGER, SAVE :: iflag_mix_omp
155	INTEGER, SAVE :: iflag_mix_adiab_omp
156	REAL, SAVE :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp
157	REAL, SAVE :: tmax_fonte_cv_omp
158
159	REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp
160	REAL,SAVE :: solaire_omp_init
161	LOGICAL,SAVE :: ok_suntime_rrtm_omp
162	REAL,SAVE :: co2_ppm_omp, RCO2_omp, co2_ppm_per_omp, RCO2_per_omp
163	REAL,SAVE :: co2_ppm0_omp
164	REAL,SAVE :: CH4_ppb_omp, RCH4_omp, CH4_ppb_per_omp, RCH4_per_omp
165	REAL,SAVE :: N2O_ppb_omp, RN2O_omp, N2O_ppb_per_omp, RN2O_per_omp
166	REAL,SAVE :: CFC11_ppt_omp,RCFC11_omp,CFC11_ppt_per_omp,RCFC11_per_omp
167	REAL,SAVE :: CFC12_ppt_omp,RCFC12_omp,CFC12_ppt_per_omp,RCFC12_per_omp
168	REAL,SAVE :: epmax_omp
169	REAL,SAVE :: coef_epmax_cape_omp
170	LOGICAL,SAVE :: ok_adj_ema_omp
171	INTEGER,SAVE :: iflag_clw_omp
172	INTEGER,SAVE :: iflag_ice_thermo_omp
173	LOGICAL,SAVE :: ok_ice_sursat_omp
174	LOGICAL,SAVE :: ok_plane_h2o_omp, ok_plane_contrail_omp
175	INTEGER,SAVE :: iflag_rei_omp
176	REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
177	REAL,SAVE :: rei_min_omp, rei_max_omp
178	INTEGER,SAVE :: iflag_sic_omp, iflag_inertie_omp
179	REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
180	REAL,SAVE :: inertie_lic_omp
181	REAL,SAVE :: qsol0_omp
182	REAL,SAVE :: evap0_omp
183	REAL,SAVE :: albsno0_omp
184	REAL :: solarlong0
185	REAL,SAVE :: solarlong0_omp
186	INTEGER,SAVE :: top_height_omp,overlap_omp
187	REAL,SAVE :: cdmmax_omp,cdhmax_omp,ksta_omp,ksta_ter_omp,f_ri_cd_min_omp
188	LOGICAL,SAVE :: ok_kzmin_omp
189	REAL, SAVE :: fmagic_omp, pmagic_omp
190	INTEGER,SAVE :: iflag_pbl_omp,lev_histhf_omp,lev_histday_omp,lev_histmth_omp
191	INTEGER,SAVE :: iflag_pbl_split_omp,iflag_physiq_omp
192	!FC
193	INTEGER,SAVE :: ifl_pbltree_omp
194	REAL,SAVE :: Cd_frein_omp
195	!FC
196	INTEGER,SAVE :: iflag_order2_sollw_omp
197	INTEGER, SAVE :: lev_histins_omp, lev_histLES_omp
198	INTEGER, SAVE :: lev_histdayNMC_omp
199	INTEGER, SAVE :: levout_histNMC_omp(3)
200	LOGICAL, SAVE :: ok_histNMC_omp(3)
201	REAL, SAVE :: freq_outNMC_omp(3), freq_calNMC_omp(3)
202	CHARACTER*4, SAVE :: type_run_omp
203	LOGICAL, SAVE :: ok_cosp_omp, ok_airs_omp
204	LOGICAL, SAVE :: ok_mensuelCOSP_omp,ok_journeCOSP_omp,ok_hfCOSP_omp
205	REAL, SAVE :: lonmin_ins_omp, lonmax_ins_omp, latmin_ins_omp, latmax_ins_omp
206	REAL, SAVE :: ecrit_hf_omp, ecrit_day_omp, ecrit_mth_omp, ecrit_reg_omp
207	REAL, SAVE :: ecrit_ins_omp
208	REAL, SAVE :: ecrit_LES_omp
209	REAL, SAVE :: ecrit_tra_omp
210	REAL, SAVE :: cvl_comp_threshold_omp
211	REAL, SAVE :: cvl_sig2feed_omp
212	REAL, SAVE :: cvl_corr_omp
213	LOGICAL, SAVE :: ok_lic_melt_omp
214	LOGICAL, SAVE :: ok_lic_cond_omp
215	!
216	REAL, SAVE :: zrel_oro_t_omp, zstd_orodr_t_omp
217	REAL, SAVE :: zpmm_orodr_t_omp, zpmm_orolf_t_omp
218	INTEGER, SAVE :: iflag_cycle_diurne_omp
219	LOGICAL, SAVE :: soil_model_omp,new_oliq_omp
220	LOGICAL, SAVE :: ok_orodr_omp, ok_orolf_omp, ok_limitvrai_omp
221	INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp
222	INTEGER, SAVE :: nbapp_cv_omp, nbapp_wk_omp
223	INTEGER, SAVE :: iflag_ener_conserv_omp
224	LOGICAL, SAVE :: ok_conserv_q_omp
225	LOGICAL,SAVE :: ok_strato_omp
226	LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
227	REAL, SAVE :: gwd_rando_ruwmax_omp, gwd_rando_sat_omp
228	REAL, SAVE :: gwd_front_ruwmax_omp, gwd_front_sat_omp
229	REAL, SAVE :: sso_gkdrag_omp,sso_grahil_omp,sso_grcrit_omp
230	REAL, SAVE :: sso_gfrcri_omp,sso_gkwake_omp,sso_gklift_omp
231	LOGICAL, SAVE :: ok_qch4_omp
232	LOGICAL, SAVE :: carbon_cycle_tr_omp
233	LOGICAL, SAVE :: carbon_cycle_cpl_omp
234	LOGICAL, SAVE :: carbon_cycle_rad_omp
235	INTEGER, SAVE :: level_coupling_esm_omp
236	LOGICAL, SAVE :: read_fco2_ocean_cor_omp
237	REAL, SAVE :: var_fco2_ocean_cor_omp
238	LOGICAL, SAVE :: read_fco2_land_cor_omp
239	REAL, SAVE :: var_fco2_land_cor_omp
240	LOGICAL, SAVE :: adjust_tropopause_omp
241	LOGICAL, SAVE :: ok_daily_climoz_omp
242	LOGICAL, SAVE :: ok_new_lscp_omp
243	LOGICAL, SAVE :: ok_icefra_lscp_omp
244	!rajout de choix_bulk et nit_bulk par Olivier Torres
245	INTEGER,SAVE :: choix_bulk_omp
246	INTEGER,SAVE :: nit_bulk_omp
247	INTEGER,SAVE :: kz0_omp
248	LOGICAL, SAVE :: ok_bs_omp, ok_rad_bs_omp
249
250
251	INTEGER, INTENT(OUT):: read_climoz ! read ozone climatology, OpenMP shared
252	! Allowed values are 0, 1 and 2
253	! 0: do not read an ozone climatology
254	! 1: read a single ozone climatology that will be used day and night
255	! 2: read two ozone climatologies, the average day and night
256	! climatology and the daylight climatology
257
258	!-----------------------------------------------------------------
259
260	print*,'CONFPHYS ENTREE'
261	!$OMP MASTER
262	!Config Key = type_ocean
263	!Config Desc = Type d'ocean
264	!Config Def = force
265	!Config Help = Type d'ocean utilise: force, slab,couple
266	!
267	type_ocean_omp = 'force '
268	CALL getin('type_ocean', type_ocean_omp)
269	!
270	!Config Key = version_ocean
271	!Config Desc = Version d'ocean
272	!Config Def = xxxxxx
273	!Config Help = Version d'ocean utilise: opa8/nemo/sicOBS/xxxxxx
274	!
275	version_ocean_omp = 'xxxxxx'
276	CALL getin('version_ocean', version_ocean_omp)
277
278	!Config Key = OCEAN
279	!Config Desc = Old parameter name for type_ocean
280	!Config Def = yyyyyy
281	!Config Help = This is only for testing purpose
282	!
283	ocean_omp = 'yyyyyy'
284	CALL getin('OCEAN', ocean_omp)
285	IF (ocean_omp /= 'yyyyyy') THEN
286	WRITE(lunout,*)'ERROR! Old variable name OCEAN used in parmeter file.'
287	WRITE(lunout,*)'Variable OCEAN has been replaced by the variable type_ocean.'
288	WRITE(lunout,*)'You have to update your parameter file physiq.def to succed running'
289	CALL abort_physic('conf_phys','Variable OCEAN no longer existing, use variable name type_ocean',1)
290	ENDIF
291
292	!Config Key = t_coupl
293	!Config Desc = Pas de temps du couplage atm/oce en sec.
294	!Config Def = 86400
295	!Config Help = This is only for testing purpose
296	!
297	t_coupl_omp = 86400.
298	CALL getin('t_coupl', t_coupl_omp)
299	IF (t_coupl_omp == 0) THEN
300	WRITE(lunout,*)'ERROR! Timestep of coupling between atmosphere and ocean'
301	WRITE(lunout,*)'cannot be zero.'
302	CALL abort_physic('conf_phys','t_coupl = 0.',1)
303	ENDIF
304
305	!
306	!Config Key = ok_all_xml
307	!Config Desc = utiliser les xml pourles définitions des champs pour xios
308	!Config Def = .FALSE.
309	!Config Help =
310	!
311	ok_all_xml_omp = .FALSE.
312	CALL getin('ok_all_xml', ok_all_xml_omp)
313
314	!
315	!Config Key = ok_lwoff
316	!Config Desc = inhiber l effet radiatif LW des nuages
317	!Config Def = .FALSE.
318	!Config Help =
319	!
320	ok_lwoff_omp = .FALSE.
321	CALL getin('ok_lwoff', ok_lwoff_omp)
322	!
323
324	!
325	!Config Key = VEGET
326	!Config Desc = Type de modele de vegetation
327	!Config Def = .FALSE.
328	!Config Help = Type de modele de vegetation utilise
329	!
330	type_veget_omp ='orchidee'
331	CALL getin('VEGET', type_veget_omp)
332	!
333
334	! INLANDSIS
335	!==================================================================
336	! Martin et Etienne
337	!Config Key = landice_opt
338	!Config Desc = which landice snow model (BULK, or INLANDSIS)
339	!Config Def = 0
340	landice_opt_omp = 0
341	CALL getin('landice_opt', landice_opt_omp)
342	! Martin et Etienne
343
344	!Etienne
345	!Config Key = iflag_temp_inlandsis
346	!Config Desc = which method to calculate temp within the soil in INLANDSIS
347	!Config Def = 0
348	iflag_temp_inlandsis_omp = 0
349	CALL getin('iflag_temp_inlandsis', iflag_temp_inlandsis_omp)
350
351	!Etienne
352	!Config Key = iflag_tsurf_inlandsis
353	!Config Desc = which method to calculate tsurf in INLANDSIS
354	!Config Def = 0
355	iflag_tsurf_inlandsis_omp = 1
356	CALL getin('iflag_tsurf_inlandsis', iflag_tsurf_inlandsis_omp)
357
358
359	!Etienne
360	!Config Key = iflag_albcalc
361	!Config Desc = method to calculate snow albedo in INLANDSIS
362	!Config Def = 0
363	iflag_albcalc_omp = 0
364	CALL getin('iflag_albcalc', iflag_albcalc_omp)
365
366
367	!Etienne
368	!Config Key = SnoMod
369	!Config Desc = activation of snow modules in inlandsis
370	!Config Def = .TRUE.
371	SnoMod_omp = .TRUE.
372	CALL getin('SnoMod', SnoMod_omp)
373
374	!Etienne
375	!Config Key = BloMod
376	!Config Desc = activation of blowing snow in inlandsis
377	!Config Def = .FALSE.
378	BloMod_omp = .FALSE.
379	CALL getin('BloMod', BloMod_omp)
380
381	!Etienne
382	!Config Key = ok_outfor
383	!Config Desc = activation of output ascii file in inlandsis
384	!Config Def = .FALSE.
385	ok_outfor_omp = .FALSE.
386	CALL getin('ok_outfor', ok_outfor_omp)
387
388
389	!Etienne
390	!Config Key = ok_sn_ii
391	!Config Desc = activation of ice/snow detection
392	!Config Def = .TRUE.
393	ok_zsn_ii_omp = .TRUE.
394	CALL getin('ok_zsn_ii', ok_zsn_ii_omp)
395
396
397	!Etienne
398	!Config Key = discret_xf
399	!Config Desc = snow discretization following XF
400	!Config Def = .TRUE.
401	discret_xf_omp = .TRUE.
402	CALL getin('discret_xf', discret_xf_omp)
403
404
405	!Etienne
406	!Config Key = is_ok_slush
407	!Config Desc = activation of the slush option
408	!Config Def = .TRUE.
409	is_ok_slush_omp = .TRUE.
410	CALL getin('is_ok_slush', is_ok_slush_omp)
411
412	!Etienne
413	!Config Key = opt_runoff_ac
414	!Config Desc = option runoff AC
415	!Config Def = .TRUE.
416	opt_runoff_ac_omp = .TRUE.
417	CALL getin('opt_runoff_ac', opt_runoff_ac_omp)
418
419	!Etienne
420	!Config Key = is_ok_z0h_rn
421	!Config Desc = z0h calculation following RN method
422	!Config Def = .TRUE.
423	is_ok_z0h_rn_omp = .TRUE.
424	CALL getin('is_ok_z0h_rn', is_ok_z0h_rn_omp)
425
426
427	!Etienne
428	!Config Key = is_ok_density_kotlyakov
429	!Config Desc = snow density calculation following kotlyakov
430	!Config Def = .FALSE.
431	is_ok_density_kotlyakov_omp = .FALSE.
432	CALL getin('is_ok_density_kotlyakov', is_ok_density_kotlyakov_omp)
433
434
435	!Etienne
436	!Config Key = prescribed_z0m_snow
437	!Config Desc = prescribed snow z0m
438	!Config Def = 0.005
439	prescribed_z0m_snow_omp = 0.005
440	CALL getin('prescribed_z0m_snow', prescribed_z0m_snow_omp)
441
442
443	!Etienne
444	!Config Key = iflag_z0m_snow
445	!Config Desc = method to calculate snow z0m
446	!Config Def = 0
447	iflag_z0m_snow_omp = 0
448	CALL getin('iflag_z0m_snow', iflag_z0m_snow_omp)
449
450
451	!Etienne
452	!Config Key = correc_alb
453	!Config Desc = correction term for albedo
454	!Config Def = 1.01
455	correc_alb_omp=1.01
456	CALL getin('correc_alb', correc_alb_omp)
457
458
459	!Etienne
460	!Config Key = buf_sph_pol
461	!Config Desc = sphericity of buffer layer in polar regions
462	!Config Def = 99.
463	buf_sph_pol_omp=99.
464	CALL getin('buf_sph_pol', buf_sph_pol_omp)
465
466	!Etienne
467	!Config Key = buf_siz_pol
468	!Config Desc = grain size of buffer layer in polar regions in e-4m
469	!Config Def = 4.
470	buf_siz_pol_omp=4.
471	CALL getin('buf_siz_pol', buf_siz_pol_omp)
472
473	!==================================================================
474
475	!Config Key = OK_journe
476	!Config Desc = Pour des sorties journalieres
477	!Config Def = .FALSE.
478	!Config Help = Pour creer le fichier histday contenant les sorties
479	! journalieres
480	!
481	ok_journe_omp = .FALSE.
482	CALL getin('OK_journe', ok_journe_omp)
483	!
484	!Config Key = ok_hf
485	!Config Desc = Pour des sorties haute frequence
486	!Config Def = .FALSE.
487	!Config Help = Pour creer le fichier histhf contenant les sorties
488	! haute frequence ( 3h ou 6h)
489	!
490	ok_hf_omp = .FALSE.
491	CALL getin('ok_hf', ok_hf_omp)
492	!
493	!Config Key = OK_mensuel
494	!Config Desc = Pour des sorties mensuelles
495	!Config Def = .TRUE.
496	!Config Help = Pour creer le fichier histmth contenant les sorties
497	! mensuelles
498	!
499	ok_mensuel_omp = .TRUE.
500	CALL getin('OK_mensuel', ok_mensuel_omp)
501	!
502	!Config Key = OK_instan
503	!Config Desc = Pour des sorties instantanees
504	!Config Def = .FALSE.
505	!Config Help = Pour creer le fichier histins contenant les sorties
506	! instantanees
507	!
508	ok_instan_omp = .FALSE.
509	CALL getin('OK_instan', ok_instan_omp)
510	!
511	!Config Key = ok_ade
512	!Config Desc = Aerosol direct effect or not?
513	!Config Def = .FALSE.
514	!Config Help = Used in radlwsw.F
515	!
516	ok_ade_omp = .FALSE.
517	CALL getin('ok_ade', ok_ade_omp)
518
519	!Config Key = ok_alw
520	!Config Desc = Aerosol longwave effect or not?
521	!Config Def = .FALSE.
522	!Config Help = Used in radlwsw.F
523	!
524	ok_alw_omp = .FALSE.
525	CALL getin('ok_alw', ok_alw_omp)
526
527	!
528	!Config Key = ok_aie
529	!Config Desc = Aerosol indirect effect or not?
530	!Config Def = .FALSE.
531	!Config Help = Used in nuage.F and radlwsw.F
532	!
533	ok_aie_omp = .FALSE.
534	CALL getin('ok_aie', ok_aie_omp)
535
536	!
537	!Config Key = ok_cdnc
538	!Config Desc = ok cloud droplet number concentration
539	!Config Def = .FALSE.
540	!Config Help = Used in newmicro.F
541	!
542	ok_cdnc_omp = .FALSE.
543	CALL getin('ok_cdnc', ok_cdnc_omp)
544
545	!
546	!Config Key = ok_volcan
547	!Config Desc = ok to generate volcanic diags
548	!Config Def = .FALSE.
549	!Config Help = Used in radlwsw_m.F
550	!
551	ok_volcan_omp = .FALSE.
552	CALL getin('ok_volcan', ok_volcan_omp)
553
554	!
555	!Config Key = flag_volc_surfstrat
556	!Config Desc = impose cooling rate at the surface (=1),
557	! heating rate in the strato (=2), or nothing (=0)
558	!Config Def = 0
559	!Config Help = Used in radlwsw_m.F
560	!
561	flag_volc_surfstrat_omp = 0 ! NL: SURFSTRAT
562	CALL getin('flag_volc_surfstrat', flag_volc_surfstrat_omp)
563
564	!
565	!Config Key = aerosol_couple
566	!Config Desc = read aerosol in file or calcul by inca
567	!Config Def = .FALSE.
568	!Config Help = Used in physiq.F
569	!
570	aerosol_couple_omp = .FALSE.
571	CALL getin('aerosol_couple',aerosol_couple_omp)
572	!
573	!Config Key = chemistry_couple
574	!Config Desc = read O3 chemistry in file or calcul by inca
575	!Config Def = .FALSE.
576	!Config Help = Used in physiq.F
577	!
578	chemistry_couple_omp = .FALSE.
579	CALL getin('chemistry_couple',chemistry_couple_omp)
580	!
581	!Config Key = flag_aerosol
582	!Config Desc = which aerosol is use for coupled model
583	!Config Def = 1
584	!Config Help = Used in physiq.F
585	!
586	! - flag_aerosol=0 => no aerosol
587	! - flag_aerosol=1 => so4 only (defaut)
588	! - flag_aerosol=2 => bc only
589	! - flag_aerosol=3 => pom only
590	! - flag_aerosol=4 => seasalt only
591	! - flag_aerosol=5 => dust only
592	! - flag_aerosol=6 => all aerosol
593	! - flag_aerosol=7 => natural aerosol + MACv2SP
594	! - (in this case aerosols.1980.nc should point to aerosols.nat.nc)
595
596	flag_aerosol_omp = 0
597	CALL getin('flag_aerosol',flag_aerosol_omp)
598
599	!
600	!Config Key = flag_bc_internal_mixture
601	!Config Desc = state of mixture for BC aerosols
602	! - n = external mixture
603	! - y = internal mixture
604	!Config Def = n
605	!Config Help = Used in physiq.F / aeropt
606	!
607	flag_bc_internal_mixture_omp = .FALSE.
608	CALL getin('flag_bc_internal_mixture',flag_bc_internal_mixture_omp)
609
610	!
611	!Config Key = aer_type
612	!Config Desc = Use a constant field for the aerosols
613	!Config Def = scenario
614	!Config Help = Used in readaerosol.F90
615	!
616	aer_type_omp = 'scenario'
617	CALL getin('aer_type', aer_type_omp)
618
619	!
620	!Config Key = bl95_b0
621	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
622	!Config Def = .FALSE.
623	!Config Help = Used in nuage.F
624	!
625	bl95_b0_omp = 2.
626	CALL getin('bl95_b0', bl95_b0_omp)
627
628	!Config Key = bl95_b1
629	!Config Desc = Parameter in CDNC-maer link (Boucher&Lohmann 1995)
630	!Config Def = .FALSE.
631	!Config Help = Used in nuage.F
632	!
633	bl95_b1_omp = 0.2
634	CALL getin('bl95_b1', bl95_b1_omp)
635
636	!Config Key = freq_ISCCP
637	!Config Desc = Frequence d'appel du simulateur ISCCP en secondes;
638	! par defaut 10800, i.e. 3 heures
639	!Config Def = 10800.
640	!Config Help = Used in ini_histISCCP.h
641	!
642	freq_ISCCP_omp = 10800.
643	CALL getin('freq_ISCCP', freq_ISCCP_omp)
644	!
645	!Config Key = ecrit_ISCCP
646	!Config Desc = Frequence d'ecriture des resultats du simulateur ISCCP en nombre de jours;
647	! par defaut 1., i.e. 1 jour
648	!Config Def = 1.
649	!Config Help = Used in ini_histISCCP.h
650	!
651	!
652	ecrit_ISCCP_omp = 1.
653	CALL getin('ecrit_ISCCP', ecrit_ISCCP_omp)
654
655	!Config Key = freq_COSP
656	!Config Desc = Frequence d'appel du simulateur COSP en secondes;
657	! par defaut 10800, i.e. 3 heures
658	!Config Def = 10800.
659	!Config Help = Used in ini_histdayCOSP.h
660	!
661	freq_COSP_omp = 10800.
662	CALL getin('freq_COSP', freq_COSP_omp)
663
664	!Config Key = freq_AIRS
665	!Config Desc = Frequence d'appel du simulateur AIRS en secondes;
666	! par defaut 10800, i.e. 3 heures
667	!Config Def = 10800.
668	!Config Help = Used in ini_histdayAIRS.h
669	!
670	freq_AIRS_omp = 10800.
671	CALL getin('freq_AIRS', freq_AIRS_omp)
672
673	!
674	!Config Key = ip_ebil_phy
675	!Config Desc = Niveau de sortie pour les diags bilan d'energie
676	!Config Def = 0
677	!Config Help =
678	!
679	ip_ebil_phy_omp = 0
680	CALL getin('ip_ebil_phy', ip_ebil_phy_omp)
681	IF (ip_ebil_phy_omp/=0) THEN
682	CALL abort_physic('conf_phys','ip_ebil_phy_omp doit etre 0 sur cette version',1)
683	ENDIF
684
685	!
686	!Config Key = seuil_inversion
687	!Config Desc = Seuil ur dTh pour le choix entre les schemas de CL
688	!Config Def = -0.1
689	!Config Help =
690	!
691	seuil_inversion_omp = -0.1
692	CALL getin('seuil_inversion', seuil_inversion_omp)
693
694	!
695	! Constante solaire & Parametres orbitaux & taux gaz effet de serre BEG
696	!
697	!Config Key = R_ecc
698	!Config Desc = Excentricite
699	!Config Def = 0.016715
700	!Config Help =
701	!
702	!valeur AMIP II
703	R_ecc_omp = 0.016715
704	CALL getin('R_ecc', R_ecc_omp)
705	!
706	!Config Key = R_peri
707	!Config Desc = Equinoxe
708	!Config Def =
709	!Config Help =
710	!
711	!
712	!valeur AMIP II
713	R_peri_omp = 102.7
714	CALL getin('R_peri', R_peri_omp)
715	!
716	!Config Key = R_incl
717	!Config Desc = Inclinaison
718	!Config Def =
719	!Config Help =
720	!
721	!
722	!valeur AMIP II
723	R_incl_omp = 23.441
724	CALL getin('R_incl', R_incl_omp)
725	!
726	!Config Key = solaire
727	!Config Desc = Constante solaire en W/m2
728	!Config Def = 1365.
729	!Config Help =
730	!
731	!
732	!valeur AMIP II
733	solaire_omp = 1365.
734	solaire_omp_init = solaire_omp !--we keep track of the default value
735	CALL getin('solaire', solaire_omp)
736	!
737	!Config Key = co2_ppm
738	!Config Desc = concentration du CO2 en ppmv
739	!Config Def = 348.
740	!Config Help =
741	!
742	!valeur AMIP II
743	co2_ppm_omp = 348.
744	CALL getin('co2_ppm', co2_ppm_omp)
745	!
746	!conversion en rapport de mélange massique
747	RCO2_omp = co2_ppm_omp * 1.0e-06 * RMCO2 / RMD
748
749	!
750	!Config Key = co2_ppm0
751	!Config Desc = concentration initiale du CO2 en ppmv pour la version ESM avec CO2 interactif dans le cas
752	! où cette concentration de figure pas dans l'état de redémarrage de la physique
753	!Config Def = 284.32
754	!Config Help =
755	!
756	co2_ppm0_omp = 284.32
757	CALL getin('co2_ppm0', co2_ppm0_omp)
758	!
759	!Config Key = RCH4
760	!Config Desc = Concentration du CH4
761	!Config Def = 1.65E-06* 16.043/28.97
762	!Config Help =
763	!
764	CH4_ppb_omp = 1650.
765	CALL getin('CH4_ppb', CH4_ppb_omp)
766	!conversion en rapport de mélange massique
767	RCH4_omp = CH4_ppb_omp * 1.0E-09 * RMCH4 / RMD
768	!
769	!Config Key = RN2O
770	!Config Desc = Concentration du N2O
771	!Config Def = 306.E-09* 44.013/28.97
772	!Config Help =
773	!
774	N2O_ppb_omp = 306.
775	CALL getin('N2O_ppb', N2O_ppb_omp)
776	!conversion en rapport de mélange massique
777	RN2O_omp = N2O_ppb_omp * 1.0E-09 * RMN2O / RMD
778	!
779	!Config Key = RCFC11
780	!Config Desc = Concentration du CFC11
781	!Config Def = 280.E-12* 137.3686/28.97
782	!Config Help =
783	!
784	CFC11_ppt_omp = 280.
785	CALL getin('CFC11_ppt',CFC11_ppt_omp)
786	!conversion en rapport de mélange massique
787	RCFC11_omp=CFC11_ppt_omp* 1.0E-12 * RMCFC11 / RMD
788	!
789	!Config Key = RCFC12
790	!Config Desc = Concentration du CFC12
791	!Config Def = 484.E-12* 120.9140/28.97
792	!Config Help =
793	!
794	CFC12_ppt_omp = 484.
795	CALL getin('CFC12_ppt',CFC12_ppt_omp)
796	!conversion en rapport de mélange massique
797	RCFC12_omp = CFC12_ppt_omp * 1.0E-12 * RMCFC12 / RMD
798
799	!
800	!Config Key = co2_ppm_per
801	!Config Desc = concentration du CO2 perturbé en ppmv (CFMIP)
802	!Config Def = 348.
803	!Config Help =
804	!
805	co2_ppm_per_omp = co2_ppm_omp
806	CALL getin('co2_ppm_per', co2_ppm_per_omp)
807	!conversion en rapport de mélange massique
808	RCO2_per_omp = co2_ppm_per_omp * 1.0e-06 * RMCO2 / RMD
809
810	!Config Key = RCH4_per
811	!Config Desc = Concentration du CH4_per
812	!Config Def = 1.65E-06* 16.043/28.97
813	!Config Help =
814	!
815	CH4_ppb_per_omp = CH4_ppb_omp
816	CALL getin('CH4_ppb_per', CH4_ppb_per_omp)
817	!conversion en rapport de mélange massique
818	RCH4_per_omp = CH4_ppb_per_omp * 1.0E-09 * RMCH4 / RMD
819	!
820	!Config Key = RN2O_per
821	!Config Desc = Concentration du N2O_per
822	!Config Def = 306.E-09* 44.013/28.97
823	!Config Help =
824	!
825	N2O_ppb_per_omp = N2O_ppb_omp
826	CALL getin('N2O_ppb_per', N2O_ppb_per_omp)
827	!conversion en rapport de mélange massique
828	RN2O_per_omp = N2O_ppb_per_omp * 1.0E-09 * RMN2O / RMD
829	!
830	!Config Key = RCFC11_per
831	!Config Desc = Concentration du CFC11_per
832	!Config Def = 280.E-12* 137.3686/28.97
833	!Config Help =
834	!
835	CFC11_ppt_per_omp = CFC11_ppt_omp
836	CALL getin('CFC11_ppt_per',CFC11_ppt_per_omp)
837	!conversion en rapport de mélange massique
838	RCFC11_per_omp=CFC11_ppt_per_omp* 1.0E-12 * RMCFC11 / RMD
839	!
840	!Config Key = RCFC12_per
841	!Config Desc = Concentration du CFC12_per
842	!Config Def = 484.E-12* 120.9140/28.97
843	!Config Help =
844	!
845	CFC12_ppt_per_omp = CFC12_ppt_omp
846	CALL getin('CFC12_ppt_per',CFC12_ppt_per_omp)
847	!conversion en rapport de mélange massique
848	RCFC12_per_omp = CFC12_ppt_per_omp * 1.0E-12 * RMCFC12 / RMD
849
850	!
851	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
852	! Constantes precedemment dans dyn3d/conf_gcm
853
854	!Config Key = iflag_cycle_diurne
855	!Config Desc = Cycle diurne
856	!Config Def = 1
857	!Config Help = Cette option permet d'eteidre le cycle diurne.
858	!Config Peut etre util pour accelerer le code !
859	iflag_cycle_diurne_omp = 1
860	CALL getin('iflag_cycle_diurne',iflag_cycle_diurne_omp)
861
862	!Config Key = soil_model
863	!Config Desc = Modele de sol
864	!Config Def = y
865	!Config Help = Choix du modele de sol (Thermique ?)
866	!Config Option qui pourait un string afin de pouvoir
867	!Config plus de choix ! Ou meme une liste d'options !
868	soil_model_omp = .TRUE.
869	CALL getin('soil_model',soil_model_omp)
870
871	!Config Key = new_oliq
872	!Config Desc = Nouvelle eau liquide
873	!Config Def = y
874	!Config Help = Permet de mettre en route la
875	!Config nouvelle parametrisation de l'eau liquide !
876	new_oliq_omp = .TRUE.
877	CALL getin('new_oliq',new_oliq_omp)
878
879	!Config Key = ok_orodr
880	!Config Desc = Orodr ???
881	!Config Def = y
882	!Config Help = Y en a pas comprendre !
883	!Config
884	ok_orodr_omp = .TRUE.
885	CALL getin('ok_orodr',ok_orodr_omp)
886
887	!Config Key = ok_orolf
888	!Config Desc = Orolf ??
889	!Config Def = y
890	!Config Help = Connais pas !
891	ok_orolf_omp = .TRUE.
892	CALL getin('ok_orolf', ok_orolf_omp)
893
894
895	!Config Key = zrel_oro_t
896	!Config Desc = zrel_oro_t
897	!Config Def = 9999.
898	!Config Help = Connais pas !
899	zrel_oro_t_omp = 9999.
900	CALL getin('zrel_oro_t', zrel_oro_t_omp)
901
902	!Config Key = zstd_orodr_t
903	!Config Desc = zstd_orodr_t
904	!Config Def = 0.
905	!Config Help = Connais pas !
906	zstd_orodr_t_omp = 10.
907	CALL getin('zstd_orodr_t', zstd_orodr_t_omp)
908
909
910	!Config Key = zpmm_orodr_t
911	!Config Desc = zpmm_orodr_t
912	!Config Def = 0.
913	!Config Help = Connais pas !
914	zpmm_orodr_t_omp = 100.
915	CALL getin('zpmm_orodr_t', zpmm_orodr_t_omp)
916
917
918	!Config Key = zpmm_orolf_t
919	!Config Desc = zpmm_orolf_t
920	!Config Def = 0.
921	!Config Help = Connais pas !
922	zpmm_orolf_t_omp = 100.
923	CALL getin('zpmm_orolf_t', zpmm_orolf_t_omp)
924
925	!Config Key = ok_limitvrai
926	!Config Desc = Force la lecture de la bonne annee
927	!Config Def = n
928	!Config Help = On peut forcer le modele a lire le
929	!Config fichier SST de la bonne annee. C'est une tres bonne
930	!Config idee, pourquoi ne pas mettre toujours a y ???
931	ok_limitvrai_omp = .FALSE.
932	CALL getin('ok_limitvrai',ok_limitvrai_omp)
933
934	!Config Key = nbapp_rad
935	!Config Desc = Frequence d'appel au rayonnement
936	!Config Def = 12
937	!Config Help = Nombre d'appels des routines de rayonnements
938	!Config par jour.
939	nbapp_rad_omp = 12
940	CALL getin('nbapp_rad',nbapp_rad_omp)
941
942	!rajout Olivier Torres
943	!Config Key = choix_bulk
944	!Config Desc = choix de la formulation bulk a prendre dans clcdrag au-dessus de l'ocean
945	!Config Def = 0
946	!Config 0 -> originale (lmdz/Louis 79)
947	!Config 1 -> COARE
948	!Config 2 -> CORE-"pure" (cf. Large)
949	!Config 3 -> CORE-"mixte" (avec z_0 et C_T^N donnees par Smith 88)
950	choix_bulk_omp = 0
951	CALL getin('choix_bulk',choix_bulk_omp)
952
953	!Config Key = nit_bulk
954	!Config Desc = choix du nombre d'it de pt fixe dans la bulk
955	!Config Def = 5
956	nit_bulk_omp = 1
957	CALL getin('nit_bulk',nit_bulk_omp)
958
959	!Config Key = kz0
960	!Config Desc = choix de la formulation z0 pour la bulk ECUME
961	!Config Def = 1
962	!Config 0 -> ARPEGE formulation
963	!Config 1 -> Smith Formulation
964	!Config 2 -> Direct computation using the stability functions
965	kz0_omp = 0
966	CALL getin('kz0',kz0_omp)
967
968
969	!Config Key = iflag_con
970	!Config Desc = Flag de convection
971	!Config Def = 2
972	!Config Help = Flag pour la convection les options suivantes existent :
973	!Config 1 pour LMD,
974	!Config 2 pour Tiedtke,
975	!Config 3 pour CCM(NCAR)
976	iflag_con_omp = 2
977	CALL getin('iflag_con',iflag_con_omp)
978
979	!Config Key = nbapp_cv
980	!Config Desc = Frequence d'appel a la convection
981	!Config Def = 0
982	!Config Help = Nombre d'appels des routines de convection
983	!Config par jour. Si =0, appel a chaque pas de temps physique.
984	nbapp_cv_omp = 0
985	CALL getin('nbapp_cv',nbapp_cv_omp)
986
987	!Config Key = nbapp_wk
988	!Config Desc = Frequence d'appel aux wakes
989	!Config Def = 0
990	!Config Help = Nombre d'appels des routines de wakes
991	!Config par jour. Si =0, appel a chaque pas de temps physique.
992	nbapp_wk_omp = 0
993	CALL getin('nbapp_wk',nbapp_wk_omp)
994
995	!Config Key = iflag_ener_conserv
996	!Config Desc = Flag de convection
997	!Config Def = 1
998	!Config Help = Flag pour la convection les options suivantes existent :
999	!Config -1 pour Kinetic energy correction
1000	!Config 1 conservation kinetic and enthalpy
1001	iflag_ener_conserv_omp = -1
1002	CALL getin('iflag_ener_conserv',iflag_ener_conserv_omp)
1003
1004	!Config Key = ok_conserv_q
1005	!Config Desc = Switch des corrections de conservation de l'eau
1006	!Config Def = y
1007	!Config Help = Switch des corrections de conservation de l'eau
1008	!Config y -> corrections activees
1009	!Config n -> conformite avec versions anterieures au 1/4/2014
1010	ok_conserv_q_omp = .FALSE.
1011	CALL getin('ok_conserv_q',ok_conserv_q_omp)
1012
1013	!
1014	!
1015	!
1016	! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
1017	!
1018	! KE
1019	!
1020
1021	!Config key = cvl_comp_threshold
1022	!Config Desc = maximum fraction of convective points enabling compression
1023	!Config Def = 1.00
1024	!Config Help = fields are compressed when less than a fraction cvl_comp_threshold
1025	!Config Help = of the points is convective.
1026	cvl_comp_threshold_omp = 1.00
1027	CALL getin('cvl_comp_threshold', cvl_comp_threshold_omp)
1028
1029	!Config key = cvl_sig2feed
1030	!Config Desc = sigma coordinate at top of feeding layer
1031	!Config Def = 0.97
1032	!Config Help = deep convection is fed by the layer extending from the surface (pressure ps)
1033	!Config Help = and cvl_sig2feed*ps.
1034	cvl_sig2feed_omp = 0.97
1035	CALL getin('cvl_sig2feed', cvl_sig2feed_omp)
1036
1037	!Config key = cvl_corr
1038	!Config Desc = Facteur multiplication des precip convectives dans KE
1039	!Config Def = 1.00
1040	!Config Help = 1.02 pour un moderne ou un pre-ind. A ajuster pour un glaciaire
1041	cvl_corr_omp = 1.00
1042	CALL getin('cvl_corr', cvl_corr_omp)
1043
1044
1045	!Config Key = epmax
1046	!Config Desc = Efficacite precip
1047	!Config Def = 0.993
1048	!Config Help =
1049	!
1050	epmax_omp = .993
1051	CALL getin('epmax', epmax_omp)
1052
1053	coef_epmax_cape_omp = 0.0
1054	CALL getin('coef_epmax_cape', coef_epmax_cape_omp)
1055	!
1056	!Config Key = ok_adj_ema
1057	!Config Desc =
1058	!Config Def = FALSE
1059	!Config Help =
1060	!
1061	ok_adj_ema_omp = .FALSE.
1062	CALL getin('ok_adj_ema',ok_adj_ema_omp)
1063	!
1064	!Config Key = iflag_clw
1065	!Config Desc =
1066	!Config Def = 0
1067	!Config Help =
1068	!
1069	iflag_clw_omp = 0
1070	CALL getin('iflag_clw',iflag_clw_omp)
1071
1072
1073	!Config Key = iflag_ratqs
1074	!Config Desc =
1075	!Config Def = 1
1076	!Config Help =
1077	!
1078	iflag_ratqs_omp = 1
1079	CALL getin('iflag_ratqs',iflag_ratqs_omp)
1080
1081	!
1082	!Config Key = iflag_radia
1083	!Config Desc =
1084	!Config Def = 1
1085	!Config Help =
1086	!
1087	iflag_radia_omp = 1
1088	CALL getin('iflag_radia',iflag_radia_omp)
1089
1090	!
1091	!Config Key = iflag_rrtm
1092	!Config Desc =
1093	!Config Def = 0
1094	!Config Help =
1095	!
1096	iflag_rrtm_omp = 0
1097	CALL getin('iflag_rrtm',iflag_rrtm_omp)
1098
1099	!
1100	!Config Key = NSW
1101	!Config Desc =
1102	!Config Def = 0
1103	!Config Help =
1104	!
1105	NSW_omp = 2
1106	CALL getin('NSW',NSW_omp)
1107	!albedo SB >>>
1108	iflag_albedo_omp = 0
1109	CALL getin('iflag_albedo',iflag_albedo_omp)
1110
1111	ok_chlorophyll_omp=.FALSE.
1112	CALL getin('ok_chlorophyll',ok_chlorophyll_omp)
1113	!albedo SB <<<
1114	!
1115	!Config Key = ok_sun_time
1116	!Config Desc = oui ou non variabilite solaire
1117	!Config Def = .FALSE.
1118	!Config Help =
1119	!
1120	!
1121	!valeur AMIP II
1122	ok_suntime_rrtm_omp = .FALSE.
1123	IF (iflag_rrtm_omp==1) THEN
1124	CALL getin('ok_suntime_rrtm',ok_suntime_rrtm_omp)
1125	ENDIF
1126
1127	!Config Key = flag_aerosol_strat
1128	!Config Desc = use stratospheric aerosols 0, 1, 2
1129	! - 0 = no stratospheric aerosols
1130	! - 1 = stratospheric aerosols scaled from 550 nm AOD
1131	! - 2 = stratospheric aerosol properties from CMIP6
1132	!Option 2 is only available with RRTM, this is tested later on
1133	!Config Def = 0
1134	!Config Help = Used in physiq.F
1135	!
1136	flag_aerosol_strat_omp = 0
1137	CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
1138
1139	!Config Key = flag_aer_feedback
1140	!Config Desc = (des)activate aerosol radiative feedback
1141	! - F = no aerosol radiative feedback
1142	! - T = aerosol radiative feedback
1143	!Config Def = T
1144	!Config Help = Used in physiq.F
1145	!
1146	flag_aer_feedback_omp = .TRUE.
1147	IF (iflag_rrtm_omp==1) THEN
1148	CALL getin('flag_aer_feedback',flag_aer_feedback_omp)
1149	ENDIF
1150
1151	!Config Key = iflag_cld_th
1152	!Config Desc =
1153	!Config Def = 1
1154	!Config Help =
1155	!
1156	iflag_cld_th_omp = 1
1157	! On lit deux fois avec l'ancien et le nouveau nom
1158	! pour assurer une retrocompatiblite.
1159	! A abandonner un jour
1160	CALL getin('iflag_cldcon',iflag_cld_th_omp)
1161	CALL getin('iflag_cld_th',iflag_cld_th_omp)
1162	iflag_cld_cv_omp = 0
1163	CALL getin('iflag_cld_cv',iflag_cld_cv_omp)
1164
1165	!
1166	!Config Key = tau_cld_cv
1167	!Config Desc =
1168	!Config Def = 10.
1169	!Config Help =
1170	!
1171	tau_cld_cv_omp = 10.
1172	CALL getin('tau_cld_cv',tau_cld_cv_omp)
1173
1174	!
1175	!Config Key = coefw_cld_cv
1176	!Config Desc =
1177	!Config Def = 0.1
1178	!Config Help =
1179	!
1180	coefw_cld_cv_omp = 0.1
1181	CALL getin('coefw_cld_cv',coefw_cld_cv_omp)
1182
1183
1184
1185
1186	!
1187	!Config Key = fact_cldcon
1188	!Config Desc =
1189	!Config Def = 0.375
1190	!Config Help =
1191	!
1192	fact_cldcon_omp = 0.375
1193	CALL getin('fact_cldcon',fact_cldcon_omp)
1194
1195	!
1196	!Config Key = facttemps
1197	!Config Desc =
1198	!Config Def = 1.e-4
1199	!Config Help =
1200	!
1201	facttemps_omp = 1.e-4
1202	CALL getin('facttemps',facttemps_omp)
1203
1204	!
1205	!Config Key = ok_newmicro
1206	!Config Desc =
1207	!Config Def = .TRUE.
1208	!Config Help =
1209	!
1210	ok_newmicro_omp = .TRUE.
1211	CALL getin('ok_newmicro',ok_newmicro_omp)
1212	!
1213	!Config Key = ratqsbas
1214	!Config Desc =
1215	!Config Def = 0.01
1216	!Config Help =
1217	!
1218	ratqsbas_omp = 0.01
1219	CALL getin('ratqsbas',ratqsbas_omp)
1220	!
1221	!Config Key = ratqshaut
1222	!Config Desc =
1223	!Config Def = 0.3
1224	!Config Help =
1225	!
1226	ratqshaut_omp = 0.3
1227	CALL getin('ratqshaut',ratqshaut_omp)
1228
1229	!Config Key = tau_ratqs
1230	!Config Desc =
1231	!Config Def = 1800.
1232	!Config Help =
1233	!
1234	tau_ratqs_omp = 1800.
1235	CALL getin('tau_ratqs',tau_ratqs_omp)
1236
1237	!
1238	!-----------------------------------------------------------------------
1239	! Longitude solaire pour le calcul de l'ensoleillement en degre
1240	! si on veut imposer la saison. Sinon, solarlong0=-999.999
1241	!Config Key = solarlong0
1242	!Config Desc =
1243	!Config Def = -999.999
1244	!Config Help =
1245	!
1246	solarlong0_omp = -999.999
1247	CALL getin('solarlong0',solarlong0_omp)
1248	!
1249	!-----------------------------------------------------------------------
1250	! Valeur imposee pour configuration idealisees
1251	!Config Key = qsol0 pour le bucket, evap0 pour aquaplanetes, albsno0
1252	! Default value -1 to activate the full computation
1253	qsol0_omp = -1.
1254	CALL getin('qsol0',qsol0_omp)
1255	evap0_omp = -1.
1256	CALL getin('evap0',evap0_omp)
1257	albsno0_omp = -1.
1258	CALL getin('albsno0',albsno0_omp)
1259	!
1260	!-----------------------------------------------------------------------
1261	!
1262	!Config Key = iflag_sic
1263	!Config Desc =
1264	!Config Def = 0
1265	!Config Help =
1266	!
1267	iflag_sic_omp = 0
1268	CALL getin('iflag_sic',iflag_sic_omp)
1269	!
1270	!Config Key = iflag_inertie
1271	!Config Desc =
1272	!Config Def = 0
1273	!Config Help =
1274	!
1275	iflag_inertie_omp = 0
1276	CALL getin('iflag_inertie',iflag_inertie_omp)
1277	!
1278	!Config Key = inertie_sic
1279	!Config Desc =
1280	!Config Def = 2000.
1281	!Config Help =
1282	!
1283	inertie_sic_omp = 2000.
1284	CALL getin('inertie_sic',inertie_sic_omp)
1285	!
1286	!Config Key = inertie_lic
1287	!Config Desc =
1288	!Config Def = 2000.
1289	!Config Help =
1290	!
1291	inertie_lic_omp = 2000.
1292	CALL getin('inertie_lic',inertie_lic_omp)
1293	!
1294	!Config Key = inertie_sno
1295	!Config Desc =
1296	!Config Def = 2000.
1297	!Config Help =
1298	!
1299	inertie_sno_omp = 2000.
1300	CALL getin('inertie_sno',inertie_sno_omp)
1301	!
1302	!Config Key = inertie_sol
1303	!Config Desc =
1304	!Config Def = 2000.
1305	!Config Help =
1306	!
1307	inertie_sol_omp = 2000.
1308	CALL getin('inertie_sol',inertie_sol_omp)
1309
1310	!
1311	!Config Key = rad_froid
1312	!Config Desc =
1313	!Config Def = 35.0
1314	!Config Help =
1315	!
1316	rad_froid_omp = 35.0
1317	CALL getin('rad_froid',rad_froid_omp)
1318
1319	!
1320	!Config Key = rad_chau1
1321	!Config Desc =
1322	!Config Def = 13.0
1323	!Config Help =
1324	!
1325	rad_chau1_omp = 13.0
1326	CALL getin('rad_chau1',rad_chau1_omp)
1327
1328	!
1329	!Config Key = rad_chau2
1330	!Config Desc =
1331	!Config Def = 9.0
1332	!Config Help =
1333	!
1334	rad_chau2_omp = 9.0
1335	CALL getin('rad_chau2',rad_chau2_omp)
1336
1337
1338	!Config Key = iflag_rei
1339	!Config Desc =
1340	!Config Def = 0
1341	!Config Help =
1342	!
1343	iflag_rei_omp = 0
1344	CALL getin('iflag_rei',iflag_rei_omp)
1345
1346
1347	!
1348	!Config Key = iflag_ice_thermo
1349	!Config Desc =
1350	!Config Def = 0
1351	!Config Help =
1352	!
1353	iflag_ice_thermo_omp = 0
1354	CALL getin('iflag_ice_thermo',iflag_ice_thermo_omp)
1355
1356	!
1357	!Config Key = ok_ice_sursat
1358	!Config Desc =
1359	!Config Def = 0
1360	!Config Help =
1361	!
1362	ok_ice_sursat_omp = 0
1363	CALL getin('ok_ice_sursat',ok_ice_sursat_omp)
1364
1365	!Config Key = ok_plane_h2o
1366	!Config Desc =
1367	!Config Def = 0
1368	!Config Help =
1369	!
1370	ok_plane_h2o_omp = .FALSE.
1371	CALL getin('ok_plane_h2o',ok_plane_h2o_omp)
1372
1373	!Config Key = ok_plane_contrail
1374	!Config Desc =
1375	!Config Def = 0
1376	!Config Help =
1377	!
1378	ok_plane_contrail_omp = .FALSE.
1379	CALL getin('ok_plane_contrail',ok_plane_contrail_omp)
1380
1381	!
1382	!Config Key = rei_min
1383	!Config Desc =
1384	!Config Def = 3.5
1385	!Config Help =
1386	!
1387	rei_min_omp = 3.5
1388	CALL getin('rei_min',rei_min_omp)
1389
1390	!
1391	!Config Key = rei_max
1392	!Config Desc =
1393	!Config Def = 61.29
1394	!Config Help =
1395	!
1396	rei_max_omp = 61.29
1397	CALL getin('rei_max',rei_max_omp)
1398
1399	!
1400	!Config Key = top_height
1401	!Config Desc =
1402	!Config Def = 3
1403	!Config Help =
1404	!
1405	top_height_omp = 3
1406	CALL getin('top_height',top_height_omp)
1407
1408	!
1409	!Config Key = overlap
1410	!Config Desc =
1411	!Config Def = 3
1412	!Config Help =
1413	!
1414	overlap_omp = 3
1415	CALL getin('overlap',overlap_omp)
1416
1417	!
1418	!Config Key = cdmmax
1419	!Config Desc =
1420	!Config Def = 1.3E-3
1421	!Config Help =
1422	!
1423	cdmmax_omp = 1.3E-3
1424	CALL getin('cdmmax',cdmmax_omp)
1425
1426	!
1427	!Config Key = cdhmax
1428	!Config Desc =
1429	!Config Def = 1.1E-3
1430	!Config Help =
1431	!
1432	cdhmax_omp = 1.1E-3
1433	CALL getin('cdhmax',cdhmax_omp)
1434
1435	!261103
1436	!
1437	!Config Key = ksta
1438	!Config Desc =
1439	!Config Def = 1.0e-10
1440	!Config Help =
1441	!
1442	ksta_omp = 1.0e-10
1443	CALL getin('ksta',ksta_omp)
1444
1445	!
1446	!Config Key = ksta_ter
1447	!Config Desc =
1448	!Config Def = 1.0e-10
1449	!Config Help =
1450	!
1451	ksta_ter_omp = 1.0e-10
1452	CALL getin('ksta_ter',ksta_ter_omp)
1453
1454	!Config Key = f_ri_cd_min
1455	!Config Desc =
1456	!Config Def = 0.1
1457	!Config Help =
1458	!
1459	f_ri_cd_min_omp = 0.1
1460	CALL getin('f_ri_cd_min',f_ri_cd_min_omp)
1461
1462	!
1463	!Config Key = ok_kzmin
1464	!Config Desc =
1465	!Config Def = .TRUE.
1466	!Config Help =
1467	!
1468	ok_kzmin_omp = .TRUE.
1469	CALL getin('ok_kzmin',ok_kzmin_omp)
1470
1471	!
1472	!Config Key = fmagic
1473	!Config Desc = additionnal multiplicator factor used for albedo
1474	!Config Def = 1.
1475	!Config Help = additionnal multiplicator factor used in albedo.F
1476	!
1477	fmagic_omp = 1.
1478	CALL getin('fmagic',fmagic_omp)
1479
1480	!
1481	!Config Key = pmagic
1482	!Config Desc = additional factor used for albedo
1483	!Config Def = 0.
1484	!Config Help = additional factor used in albedo.F
1485	!
1486	pmagic_omp = 0.
1487	CALL getin('pmagic',pmagic_omp)
1488
1489
1490	!Config Key = ok_lic_melt
1491	!Config Desc = Prise en compte de la fonte de la calotte dans le bilan d'eau
1492	!Config Def = .FALSE.
1493	!Config Help = mettre a .FALSE. pour assurer la conservation en eau
1494	ok_lic_melt_omp = .FALSE.
1495	CALL getin('ok_lic_melt', ok_lic_melt_omp)
1496
1497
1498	!Config Key = ok_lic_cond
1499	!Config Desc = Prise en compte depot de vapeur d'eau sur la calotte dans le bilan d'eau
1500	!Config Def = .FALSE.
1501	!Config Help = mettre a .TRUE. pour assurer la conservation en eau
1502	ok_lic_cond_omp = .FALSE.
1503	CALL getin('ok_lic_cond', ok_lic_cond_omp)
1504
1505	!
1506	! PARAMETER FOR THE PLANETARY BOUNDARY LAYER
1507	!
1508
1509	!Config Key = iflag_pbl
1510	!Config Desc =
1511	!Config Def = 1
1512	!Config Help =
1513	!
1514	iflag_pbl_omp = 1
1515	CALL getin('iflag_pbl',iflag_pbl_omp)
1516
1517	iflag_physiq_omp = 0 ! 0: std, 1: with physiqex output, 2: for physiqex
1518	CALL getin('iflag_physiq',iflag_physiq_omp)
1519
1520	!FC
1521	!Config Key = ifl_pbltree
1522	!Config Desc = drag from trees 0 no activated
1523	!Config Def = 0
1524	!Config Help =
1525	!
1526	ifl_pbltree_omp = 0
1527	CALL getin('ifl_pbltree',ifl_pbltree_omp)
1528	!FC
1529	!Config Key = Cd_frein
1530	!Config Desc = drag from trees
1531	!Config Def = 7.5E-02 (valeur Masson mais fait planter avec des LAI eleves)
1532	!Config Help =
1533	!
1534	Cd_frein_omp = 7.5E-02
1535	CALL getin('Cd_frein',Cd_frein_omp)
1536
1537	!
1538	!Config Key = iflag_pbl_split
1539	!Config Desc = decimal flag: least signif digit = split vdf; next digit = split thermals
1540	!Config Def = 0
1541	!Config Help = 0-> no splitting; 1-> vdf splitting; 10-> thermals splitting; 11-> full splitting
1542	!
1543	iflag_pbl_split_omp = 0
1544	call getin('iflag_pbl_split',iflag_pbl_split_omp)
1545	!
1546	!Config Key = iflag_order2_sollw
1547	!Config Desc =
1548	!Config Def = 0
1549	!Config Help =
1550	!
1551	iflag_order2_sollw_omp = 0
1552	CALL getin('iflag_order2_sollw',iflag_order2_sollw_omp)
1553	!
1554	!Config Key = iflag_thermals
1555	!Config Desc =
1556	!Config Def = 0
1557	!Config Help =
1558	!
1559	iflag_thermals_omp = 0
1560	CALL getin('iflag_thermals',iflag_thermals_omp)
1561	!
1562	!Config Key = nsplit_thermals
1563	!Config Desc =
1564	!Config Def = 0
1565	!Config Help =
1566	!
1567	nsplit_thermals_omp = 1
1568	CALL getin('nsplit_thermals',nsplit_thermals_omp)
1569	!
1570	!Config Key = nsplit_thermals
1571	!Config Desc =
1572	!Config Def = 0
1573	!Config Help =
1574	alp_bl_k_omp = 1.
1575	CALL getin('alp_bl_k',alp_bl_k_omp)
1576
1577	! nrlmd le 10/04/2012
1578
1579	!Config Key = iflag_trig_bl
1580	!Config Desc =
1581	!Config Def = 0
1582	!Config Help =
1583	!
1584	iflag_trig_bl_omp = 0
1585	CALL getin('iflag_trig_bl',iflag_trig_bl_omp)
1586
1587	!Config Key = s_trig_bl
1588	!Config Desc =
1589	!Config Def = 0
1590	!Config Help =
1591	!
1592	s_trig_omp = 2e7
1593	CALL getin('s_trig',s_trig_omp)
1594
1595	!Config Key = tau_trig_shallow
1596	!Config Desc =
1597	!Config Def = 0
1598	!Config Help =
1599	!
1600	tau_trig_shallow_omp = 600
1601	CALL getin('tau_trig_shallow',tau_trig_shallow_omp)
1602
1603	!Config Key = tau_trig_deep
1604	!Config Desc =
1605	!Config Def = 0
1606	!Config Help =
1607	!
1608	tau_trig_deep_omp = 1800
1609	CALL getin('tau_trig_deep',tau_trig_deep_omp)
1610
1611	!Config Key = iflag_clos_bl
1612	!Config Desc =
1613	!Config Def = 0
1614	!Config Help =
1615	!
1616	iflag_clos_bl_omp = 0
1617	CALL getin('iflag_clos_bl',iflag_clos_bl_omp)
1618
1619	! fin nrlmd le 10/04/2012
1620
1621	!
1622	!Config Key = tau_thermals
1623	!Config Desc =
1624	!Config Def = 0.
1625	!Config Help =
1626	!
1627	tau_thermals_omp = 0.
1628	CALL getin('tau_thermals',tau_thermals_omp)
1629
1630	!
1631	!Config Key = iflag_coupl
1632	!Config Desc =
1633	!Config Def = 0
1634	!Config Help =
1635	!
1636	iflag_coupl_omp = 0
1637	CALL getin('iflag_coupl',iflag_coupl_omp)
1638
1639	!
1640	!Config Key = iflag_clos
1641	!Config Desc =
1642	!Config Def = 0
1643	!Config Help =
1644	!
1645	iflag_clos_omp = 1
1646	CALL getin('iflag_clos',iflag_clos_omp)
1647	!
1648	!Config Key = coef_clos_ls
1649	!Config Desc =
1650	!Config Def = 0
1651	!Config Help =
1652	!
1653	coef_clos_ls_omp = 0.
1654	CALL getin('coef_clos_ls',coef_clos_ls_omp)
1655
1656	!
1657	!Config Key = iflag_cvl_sigd
1658	!Config Desc =
1659	!Config Def = 0
1660	!Config Help =
1661	!
1662	iflag_cvl_sigd_omp = 0
1663	CALL getin('iflag_cvl_sigd',iflag_cvl_sigd_omp)
1664
1665	!Config Key = iflag_wake
1666	!Config Desc =
1667	!Config Def = 0
1668	!Config Help =
1669	!
1670	iflag_wake_omp = 0
1671	CALL getin('iflag_wake',iflag_wake_omp)
1672
1673	!Config Key = alp_offset
1674	!Config Desc =
1675	!Config Def = 0
1676	!Config Help =
1677	!
1678	alp_offset_omp = 0.
1679	CALL getin('alp_offset',alp_offset_omp)
1680
1681	!
1682	!Config Key = lev_histhf
1683	!Config Desc =
1684	!Config Def = 1
1685	!Config Help =
1686	!
1687	lev_histhf_omp = 1
1688	CALL getin('lev_histhf',lev_histhf_omp)
1689
1690	!
1691	!Config Key = lev_histday
1692	!Config Desc =
1693	!Config Def = 1
1694	!Config Help =
1695	!
1696	lev_histday_omp = 1
1697	CALL getin('lev_histday',lev_histday_omp)
1698
1699	!
1700	!Config Key = lev_histmth
1701	!Config Desc =
1702	!Config Def = 2
1703	!Config Help =
1704	!
1705	lev_histmth_omp = 2
1706	CALL getin('lev_histmth',lev_histmth_omp)
1707	!
1708	!Config Key = lev_histins
1709	!Config Desc =
1710	!Config Def = 1
1711	!Config Help =
1712	!
1713	lev_histins_omp = 1
1714	CALL getin('lev_histins',lev_histins_omp)
1715	!
1716	!Config Key = lev_histLES
1717	!Config Desc =
1718	!Config Def = 1
1719	!Config Help =
1720	!
1721	lev_histLES_omp = 1
1722	CALL getin('lev_histLES',lev_histLES_omp)
1723	!
1724	!Config Key = lev_histdayNMC
1725	!Config Desc =
1726	!Config Def = 8
1727	!Config Help =
1728	!
1729	lev_histdayNMC_omp = 8
1730	CALL getin('lev_histdayNMC',lev_histdayNMC_omp)
1731	!
1732	!Config Key = levout_histNMC
1733	!Config Desc =
1734	!Config Def = 5
1735	!Config Help =
1736	!
1737	levout_histNMC_omp(1) = 5
1738	levout_histNMC_omp(2) = 5
1739	levout_histNMC_omp(3) = 5
1740	CALL getin('levout_histNMC',levout_histNMC_omp)
1741	!
1742	!histNMC BEG
1743	!Config Key = ok_histNMC
1744	!Config Desc = ok_histNMC(1) = frequence de sortie fichiers histmthNMC
1745	!Config Desc = ok_histNMC(2) = frequence de sortie fichiers histdayNMC
1746	!Config Desc = ok_histNMC(3) = frequence de sortie fichiers histhfNMC
1747	!Config Def = n, n, n
1748	!Config Help =
1749	!
1750	ok_histNMC_omp(1) = .FALSE.
1751	ok_histNMC_omp(2) = .FALSE.
1752	ok_histNMC_omp(3) = .FALSE.
1753	CALL getin('ok_histNMC',ok_histNMC_omp)
1754	!
1755	!Config Key = freq_outNMC
1756	!Config Desc = freq_outNMC(1) = frequence de sortie fichiers histmthNMC
1757	!Config Desc = freq_outNMC(2) = frequence de sortie fichiers histdayNMC
1758	!Config Desc = freq_outNMC(3) = frequence de sortie fichiers histhfNMC
1759	!Config Def = 2592000., 86400., 21600. (1mois, 1jour, 6h)
1760	!Config Help =
1761	!
1762	freq_outNMC_omp(1) = mth_len
1763	freq_outNMC_omp(2) = 1.
1764	freq_outNMC_omp(3) = 1./4.
1765	CALL getin('freq_outNMC',freq_outNMC_omp)
1766	!
1767	!Config Key = freq_calNMC
1768	!Config Desc = freq_calNMC(1) = frequence de calcul fichiers histmthNMC
1769	!Config Desc = freq_calNMC(2) = frequence de calcul fichiers histdayNMC
1770	!Config Desc = freq_calNMC(3) = frequence de calcul fichiers histhfNMC
1771	!Config Def = phys_tstep
1772	!Config Help =
1773	!
1774	freq_calNMC_omp(1) = phys_tstep
1775	freq_calNMC_omp(2) = phys_tstep
1776	freq_calNMC_omp(3) = phys_tstep
1777	CALL getin('freq_calNMC',freq_calNMC_omp)
1778	!
1779	!Config Key = type_run
1780	!Config Desc =
1781	!Config Def = 'AMIP'/'CFMIP' ou 'CLIM'/'ENSP'
1782	!Config Help =
1783	!
1784	type_run_omp = 'AMIP'
1785	CALL getin('type_run',type_run_omp)
1786
1787	!
1788	!Config Key = ok_cosp
1789	!Config Desc =
1790	!Config Def = .FALSE.
1791	!Config Help =
1792	!
1793	ok_cosp_omp = .FALSE.
1794	CALL getin('ok_cosp',ok_cosp_omp)
1795
1796	!
1797	!Config Key = ok_airs
1798	!Config Desc =
1799	!Config Def = .FALSE.
1800	!Config Help =
1801	!
1802	ok_airs_omp = .FALSE.
1803	CALL getin('ok_airs',ok_airs_omp)
1804
1805	!
1806	!Config Key = ok_mensuelCOSP
1807	!Config Desc =
1808	!Config Def = .TRUE.
1809	!Config Help =
1810	!
1811	ok_mensuelCOSP_omp = .TRUE.
1812	CALL getin('ok_mensuelCOSP',ok_mensuelCOSP_omp)
1813
1814	!
1815	!Config Key = ok_journeCOSP
1816	!Config Desc =
1817	!Config Def = .TRUE.
1818	!Config Help =
1819	!
1820	ok_journeCOSP_omp = .TRUE.
1821	CALL getin('ok_journeCOSP',ok_journeCOSP_omp)
1822
1823	!
1824	!Config Key = ok_hfCOSP
1825	!Config Desc =
1826	!Config Def = .FALSE.
1827	!Config Help =
1828	!
1829	ok_hfCOSP_omp = .FALSE.
1830	CALL getin('ok_hfCOSP',ok_hfCOSP_omp)
1831
1832	!
1833	! coordonnees (lonmin_ins, lonmax_ins, latmin_ins, latmax_ins) pour la zone
1834	! avec sorties instantannees tous les pas de temps de la physique => "histbilKP_ins.nc"
1835	!
1836	!Config Key = lonmin_ins
1837	!Config Desc = 100.
1838	!Config Def = longitude minimale sorties "bilKP_ins"
1839	!Config Help =
1840	!
1841	lonmin_ins_omp = 100.
1842	CALL getin('lonmin_ins',lonmin_ins_omp)
1843	!
1844	!Config Key = lonmax_ins
1845	!Config Desc = 130.
1846	!Config Def = longitude maximale sorties "bilKP_ins"
1847	!Config Help =
1848	!
1849	lonmax_ins_omp = 130.
1850	CALL getin('lonmax_ins',lonmax_ins_omp)
1851	!
1852	!Config Key = latmin_ins
1853	!Config Desc = -20.
1854	!Config Def = latitude minimale sorties "bilKP_ins"
1855	!Config Help =
1856	!
1857	latmin_ins_omp = -20.
1858	CALL getin('latmin_ins',latmin_ins_omp)
1859	!
1860	!Config Key = latmax_ins
1861	!Config Desc = 20.
1862	!Config Def = latitude maximale sorties "bilKP_ins"
1863	!Config Help =
1864	!
1865	latmax_ins_omp = 20.
1866	CALL getin('latmax_ins',latmax_ins_omp)
1867	!
1868	!Config Key = ecrit_hf
1869	!Config Desc =
1870	!Config Def = 1./8. !toutes les 3h
1871	!Config Help =
1872	!
1873	ecrit_hf_omp = 1./8.
1874	CALL getin('ecrit_hf',ecrit_hf_omp)
1875	!
1876	!Config Key = ecrit_ins
1877	!Config Desc =
1878	!Config Def = 1./48. ! toutes les 1/2 h
1879	!Config Help =
1880	!
1881	ecrit_ins_omp = 1./48.
1882	CALL getin('ecrit_ins',ecrit_ins_omp)
1883	!
1884	!Config Key = ecrit_day
1885	!Config Desc =
1886	!Config Def = 1.0 !tous les jours
1887	!Config Help = nombre de jours pour ecriture fichier histday.nc
1888	!
1889	ecrit_day_omp = 1.0
1890	CALL getin('ecrit_day',ecrit_day_omp)
1891	!
1892	!Config Key = ecrit_mth
1893	!Config Desc =
1894	!Config Def = 30. !tous les 30jours (1 fois par mois)
1895	!Config Help =
1896	!
1897	ecrit_mth_omp = 30.
1898	CALL getin('ecrit_mth',ecrit_mth_omp)
1899	!
1900	!Config Key = ecrit_tra
1901	!Config Desc =
1902	!Config Def = 30. !tous les 30jours (1 fois par mois)
1903	!Config Help =
1904	!
1905	ecrit_tra_omp = 0.
1906	CALL getin('ecrit_tra',ecrit_tra_omp)
1907	!
1908	!Config Key = ecrit_reg
1909	!Config Desc =
1910	!Config Def = 0.25 !4 fois par jour
1911	!Config Help =
1912	!
1913	ecrit_reg_omp = 0.25 !4 fois par jour
1914	CALL getin('ecrit_reg',ecrit_reg_omp)
1915	!
1916	!
1917	print*,'CONFPHYS OOK avant drag_ter'
1918	!
1919	! PARAMETRES CDRAG
1920	!
1921	f_cdrag_ter_omp = 0.8
1922	CALL getin('f_cdrag_ter',f_cdrag_ter_omp)
1923	!
1924	f_cdrag_oce_omp = 0.8
1925	CALL getin('f_cdrag_oce',f_cdrag_oce_omp)
1926	!
1927
1928	! Gustiness flags
1929	f_z0qh_oce_omp = 1.
1930	CALL getin('f_z0qh_oce',f_z0qh_oce_omp)
1931	!
1932	f_qsat_oce_omp = 1.
1933	CALL getin('f_qsat_oce',f_qsat_oce_omp)
1934	!
1935	f_gust_bl_omp = 0.
1936	CALL getin('f_gust_bl',f_gust_bl_omp)
1937	!
1938	f_gust_wk_omp = 0.
1939	CALL getin('f_gust_wk',f_gust_wk_omp)
1940	!
1941	!Config Key = iflag_z0_oce
1942	!Config Desc = 0 (z0h=z0m), 1 (diff. equ. for z0h and z0m), -1 (z0m=z0h=z0min)
1943	!Config Def = 0 ! z0h = z0m
1944	!Config Help =
1945	!
1946	iflag_z0_oce_omp=0
1947	CALL getin('iflag_z0_oce',iflag_z0_oce_omp)
1948	!
1949	iflag_gusts_omp=0
1950	CALL getin('iflag_gusts',iflag_gusts_omp)
1951	!
1952	min_wind_speed_omp = 1.
1953	CALL getin('min_wind_speed',min_wind_speed_omp)
1954
1955	z0m_seaice_omp = 0.002 ; CALL getin('z0m_seaice',z0m_seaice_omp)
1956	z0h_seaice_omp = 0.002 ; CALL getin('z0h_seaice',z0h_seaice_omp)
1957
1958
1959	z0m_landice_omp = 0.001 ; CALL getin('z0m_landice',z0m_landice_omp)
1960	z0h_landice_omp = 0.001 ; CALL getin('z0h_landice',z0h_landice_omp)
1961
1962	f_rugoro_omp = 0.
1963	CALL getin('f_rugoro',f_rugoro_omp)
1964
1965	z0min_omp = 0.000015
1966	CALL getin('z0min',z0min_omp)
1967
1968
1969	! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
1970	!
1971	!Config Key = supcrit1
1972	!Config Desc =
1973	!Config Def = .540
1974	!Config Help =
1975	!
1976	supcrit1_omp = .540
1977	CALL getin('supcrit1',supcrit1_omp)
1978
1979	!
1980	!Config Key = supcrit2
1981	!Config Desc =
1982	!Config Def = .600
1983	!Config Help =
1984	!
1985	supcrit2_omp = .600
1986	CALL getin('supcrit2',supcrit2_omp)
1987
1988	!
1989	! PARAMETERS FOR THE MIXING DISTRIBUTION
1990	! iflag_mix: 0=OLD,
1991	! 1=NEW (JYG),
1992	! 2=NEW + conv. depth inhib. by tropos. dryness
1993	! '2' is NOT operationnal and should not be used.
1994	!
1995	!Config Key = iflag_mix
1996	!Config Desc =
1997	!Config Def = 1
1998	!Config Help =
1999	!
2000	iflag_mix_omp = 1
2001	CALL getin('iflag_mix',iflag_mix_omp)
2002
2003	!
2004	! PARAMETERS FOR THE EROSION OF THE ADIABATIC ASCENTS
2005	! iflag_mix_adiab: 0=OLD,
2006	! 1=NEW (CR),
2007	!
2008	!
2009	!Config Key = iflag_mix_adiab
2010	!Config Desc =
2011	!Config Def = 1
2012	!Config Help =
2013	!
2014	iflag_mix_adiab_omp = 0
2015	CALL getin('iflag_mix_adiab',iflag_mix_adiab_omp)
2016
2017	!
2018	!Config Key = scut
2019	!Config Desc =
2020	!Config Def = 0.95
2021	!Config Help =
2022	!
2023	scut_omp = 0.95
2024	CALL getin('scut',scut_omp)
2025
2026	!
2027	!Config Key = qqa1
2028	!Config Desc =
2029	!Config Def = 1.0
2030	!Config Help =
2031	!
2032	qqa1_omp = 1.0
2033	CALL getin('qqa1',qqa1_omp)
2034
2035	!
2036	!Config Key = qqa2
2037	!Config Desc =
2038	!Config Def = 0.0
2039	!Config Help =
2040	!
2041	qqa2_omp = 0.0
2042	CALL getin('qqa2',qqa2_omp)
2043
2044	!
2045	!Config Key = gammas
2046	!Config Desc =
2047	!Config Def = 0.05
2048	!Config Help =
2049	!
2050	gammas_omp = 0.05
2051	CALL getin('gammas',gammas_omp)
2052
2053	!
2054	!Config Key = Fmax
2055	!Config Desc =
2056	!Config Def = 0.65
2057	!Config Help =
2058	!
2059	Fmax_omp = 0.65
2060	CALL getin('Fmax',Fmax_omp)
2061
2062	!
2063	!Config Key = tmax_fonte_cv
2064	!Config Desc =
2065	!Config Def = 275.15
2066	!Config Help =
2067	!
2068	tmax_fonte_cv_omp = 275.15
2069	CALL getin('tmax_fonte_cv',tmax_fonte_cv_omp)
2070
2071	!
2072	!Config Key = alphas
2073	!Config Desc =
2074	!Config Def = -5.
2075	!Config Help =
2076	!
2077	alphas_omp = -5.
2078	CALL getin('alphas',alphas_omp)
2079
2080	!Config key = ok_strato
2081	!Config Desc = activation de la version strato
2082	!Config Def = .FALSE.
2083	!Config Help = active la version stratosph\'erique de LMDZ de F. Lott
2084	! Et la sponge layer (Runs Stratospheriques)
2085
2086	ok_strato_omp=.FALSE.
2087	CALL getin('ok_strato',ok_strato_omp)
2088
2089	!Config key = ok_hines
2090	!Config Desc = activation de la parametrisation de hines
2091	!Config Def = .FALSE.
2092	!Config Help = Clefs controlant la parametrization de Hines
2093
2094	ok_hines_omp=.FALSE.
2095	CALL getin('ok_hines',ok_hines_omp)
2096
2097	! Parametres pour les ondes de gravite
2098	!
2099	! Subgrid Scale Orography (Lott Miller (1997), Lott (1999))
2100
2101	sso_gkdrag_omp = merge(0.1875, 0.2, ok_strato_omp)
2102	CALL getin('sso_gkdrag', sso_gkdrag_omp)
2103
2104	sso_grahil_omp=merge(0.1,1.,ok_strato_omp)
2105	CALL getin('sso_grahil', sso_grahil_omp)
2106
2107	sso_grcrit_omp =merge(1.,0.01,ok_strato_omp)
2108	CALL getin('sso_grcrit', sso_grcrit_omp)
2109
2110	sso_gfrcri_omp = 1.
2111	CALL getin('sso_gfrcri', sso_gfrcri_omp)
2112
2113	sso_gkwake_omp = 0.50
2114	CALL getin('sso_gkwake', sso_gkwake_omp)
2115
2116	sso_gklift_omp = merge(0.25,0.50,ok_strato_omp)
2117	CALL getin('sso_gklift', sso_gklift_omp)
2118
2119	! Random gravity waves:
2120
2121	ok_gwd_rando_omp = .FALSE.
2122	IF ( klon_glo == 1 ) THEN
2123	print*,'La parametrisation des ondes de gravites non orographiques'
2124	print*,'ne fonctionne pas en 1D'
2125	ELSE
2126	CALL getin('ok_gwd_rando', ok_gwd_rando_omp)
2127	ENDIF
2128
2129	gwd_rando_ruwmax_omp = 2.00
2130	CALL getin('gwd_rando_ruwmax', gwd_rando_ruwmax_omp)
2131
2132	gwd_rando_sat_omp = 0.25
2133	CALL getin('gwd_rando_sat', gwd_rando_sat_omp)
2134
2135	gwd_front_ruwmax_omp = 2.50
2136	CALL getin('gwd_front_ruwmax', gwd_front_ruwmax_omp)
2137
2138	gwd_front_sat_omp = 0.60
2139	CALL getin('gwd_front_sat', gwd_front_sat_omp)
2140
2141
2142	!Config key = ok_qch4
2143	!Config Desc = activation de la parametrisation du methane
2144	!Config Def = .FALSE.
2145	!Config Help = Clef controlant l'activation de la parametrisation
2146	! de l'humidite due a oxydation+photolyse du methane strato
2147
2148	ok_qch4_omp=.FALSE.
2149	CALL getin('ok_qch4',ok_qch4_omp)
2150
2151	!Config Key = OK_LES
2152	!Config Desc = Pour des sorties LES
2153	!Config Def = .FALSE.
2154	!Config Help = Pour creer le fichier histLES contenant les sorties
2155	! LES
2156	!
2157	ok_LES_omp = .FALSE.
2158	CALL getin('OK_LES', ok_LES_omp)
2159
2160	!Config Key = callstats
2161	!Config Desc = Pour des sorties callstats
2162	!Config Def = .FALSE.
2163	!Config Help = Pour creer le fichier stats contenant les sorties
2164	! stats
2165	!
2166	callstats_omp = .FALSE.
2167	CALL getin('callstats', callstats_omp)
2168	!
2169	!Config Key = ecrit_LES
2170	!Config Desc = Frequence d'ecriture des resultats du LES en nombre de jours;
2171	! par defaut 1., i.e. 1 jour
2172	!Config Def = 1./8.
2173	!Config Help = ...
2174	!
2175	!
2176	adjust_tropopause = .FALSE.
2177	CALL getin('adjust_tropopause', adjust_tropopause_omp)
2178	!
2179	!Config Key = adjust_tropopause
2180	!Config Desc = Adjust the ozone field from the climoz file by stretching its
2181	! tropopause so that it matches the one of LMDZ.
2182	!Config Def = .FALSE.
2183	!Config Help = Ensure tropospheric ozone column conservation.
2184	!
2185	!
2186	ok_daily_climoz = .FALSE.
2187	CALL getin('ok_daily_climoz', ok_daily_climoz_omp)
2188	!
2189	!Config Key = ok_daily_climoz
2190	!Config Desc = Interpolate in time the ozone forcings within ce0l.
2191	! .TRUE. if backward compatibility is needed.
2192	!Config Def = .TRUE.
2193	!Config Help = .FALSE. ensure much fewer (no calendar dependency)
2194	! and lighter monthly climoz files, inetrpolated in time at gcm run time.
2195
2196	ok_new_lscp_omp = .FALSE.
2197	CALL getin('ok_new_lscp', ok_new_lscp_omp)
2198	!
2199	!Config Key = ok_new_lscp_omp
2200	!Config Desc = new cloud scheme ith ice and mixed phase (Etienne and JB)
2201	!Config Def = .FALSE.
2202	!Config Help = ...
2203
2204
2205
2206	ok_icefra_lscp_omp = .FALSE.
2207	CALL getin('ok_icefra_lscp', ok_icefra_lscp_omp)
2208	!
2209	!Config Key = ok_icefra_lscp_omp
2210	!Config Desc = ice fraction in radiation from lscp
2211	!Config Def = .FALSE.
2212	!Config Help = ...
2213
2214
2215
2216	ok_bs_omp = .FALSE.
2217	CALL getin('ok_bs', ok_bs_omp)
2218	!
2219	!Config Key = ok_bs_omp
2220	!Config Desc = include blowing snow
2221	!Config Def = .FALSE.
2222	!Config Help = ...
2223
2224
2225	ok_rad_bs_omp = .FALSE.
2226	CALL getin('ok_rad_bs', ok_rad_bs_omp)
2227	!
2228	!Config Key = ok_rad_bs_omp
2229	!Config Desc = include blowing snow radiative effect
2230	!Config Def = .FALSE.
2231	!Config Help = ...
2232
2233	ecrit_LES_omp = 1./8.
2234	CALL getin('ecrit_LES', ecrit_LES_omp)
2235	!
2236	read_climoz = 0 ! default value
2237	CALL getin('read_climoz', read_climoz)
2238
2239	carbon_cycle_tr_omp=.FALSE.
2240	CALL getin('carbon_cycle_tr',carbon_cycle_tr_omp)
2241
2242	carbon_cycle_cpl_omp=.FALSE.
2243	CALL getin('carbon_cycle_cpl',carbon_cycle_cpl_omp)
2244
2245	carbon_cycle_rad_omp=.FALSE.
2246	CALL getin('carbon_cycle_rad',carbon_cycle_rad_omp)
2247
2248	read_fco2_ocean_cor_omp=.FALSE.
2249	CALL getin('read_fco2_ocean_cor',read_fco2_ocean_cor_omp)
2250
2251	var_fco2_ocean_cor_omp=0. ! default value
2252	CALL getin('var_fco2_ocean_cor',var_fco2_ocean_cor_omp)
2253
2254	read_fco2_land_cor_omp=.FALSE.
2255	CALL getin('read_fco2_land_cor',read_fco2_land_cor_omp)
2256
2257	var_fco2_land_cor_omp=0. ! default value
2258	CALL getin('var_fco2_land_cor',var_fco2_land_cor_omp)
2259
2260	! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO
2261	! Definitions of level_coupling_esm in physiq.def
2262	! level_coupling_esm = 0 ! No field exchange between LMDZ and ORCHIDEE models
2263	! ! No field exchange between LMDZ and NEMO
2264	! level_coupling_esm = 1 ! Field exchange between LMDZ and ORCHIDEE models
2265	! ! No field exchange between LMDZ and NEMO models
2266	! level_coupling_esm = 2 ! No field exchange between LMDZ and ORCHIDEE models
2267	! ! Field exchange between LMDZ and NEMO models
2268	! level_coupling_esm = 3 ! Field exchange between LMDZ and ORCHIDEE models
2269	! ! Field exchange between LMDZ and NEMO models
2270	level_coupling_esm_omp=0 ! default value
2271	CALL getin('level_coupling_esm',level_coupling_esm_omp)
2272
2273	!$OMP END MASTER
2274	!$OMP BARRIER
2275
2276	R_ecc = R_ecc_omp
2277	R_peri = R_peri_omp
2278	R_incl = R_incl_omp
2279	solaire = solaire_omp
2280	ok_suntime_rrtm = ok_suntime_rrtm_omp
2281	co2_ppm = co2_ppm_omp
2282	co2_ppm0 = co2_ppm0_omp
2283	RCO2 = RCO2_omp
2284	CH4_ppb = CH4_ppb_omp
2285	RCH4 = RCH4_omp
2286	N2O_ppb = N2O_ppb_omp
2287	RN2O = RN2O_omp
2288	CFC11_ppt = CFC11_ppt_omp
2289	RCFC11 = RCFC11_omp
2290	CFC12_ppt = CFC12_ppt_omp
2291	RCFC12 = RCFC12_omp
2292	RCO2_act = RCO2
2293	RCH4_act = RCH4
2294	RN2O_act = RN2O
2295	RCFC11_act = RCFC11
2296	RCFC12_act = RCFC12
2297	RCO2_per = RCO2_per_omp
2298	RCH4_per = RCH4_per_omp
2299	RN2O_per = RN2O_per_omp
2300	RCFC11_per = RCFC11_per_omp
2301	RCFC12_per = RCFC12_per_omp
2302
2303	iflag_cycle_diurne = iflag_cycle_diurne_omp
2304	soil_model = soil_model_omp
2305	new_oliq = new_oliq_omp
2306	ok_orodr = ok_orodr_omp
2307	ok_orolf = ok_orolf_omp
2308	zrel_oro_t=zrel_oro_t_omp
2309	zstd_orodr_t=zstd_orodr_t_omp
2310	zpmm_orodr_t=zpmm_orodr_t_omp
2311	zpmm_orolf_t=zpmm_orolf_t_omp
2312	ok_limitvrai = ok_limitvrai_omp
2313	nbapp_rad = nbapp_rad_omp
2314	iflag_con = iflag_con_omp
2315	nbapp_cv = nbapp_cv_omp
2316	nbapp_wk = nbapp_wk_omp
2317	iflag_ener_conserv = iflag_ener_conserv_omp
2318	ok_conserv_q = ok_conserv_q_omp
2319	epmax = epmax_omp
2320	coef_epmax_cape = coef_epmax_cape_omp
2321	ok_adj_ema = ok_adj_ema_omp
2322	iflag_clw = iflag_clw_omp
2323	solarlong0 = solarlong0_omp
2324	qsol0 = qsol0_omp
2325	evap0 = evap0_omp
2326	albsno0 = albsno0_omp
2327	iflag_sic = iflag_sic_omp
2328	iflag_inertie = iflag_inertie_omp
2329	inertie_sol = inertie_sol_omp
2330	inertie_sic = inertie_sic_omp
2331	inertie_lic = inertie_lic_omp
2332	inertie_sno = inertie_sno_omp
2333	ok_bs = ok_bs_omp
2334	ok_rad_bs=ok_rad_bs_omp
2335	rad_froid = rad_froid_omp
2336	rad_chau1 = rad_chau1_omp
2337	rad_chau2 = rad_chau2_omp
2338	iflag_rei=iflag_rei_omp
2339	iflag_ice_thermo = iflag_ice_thermo_omp
2340	ok_ice_sursat = ok_ice_sursat_omp
2341	ok_plane_h2o = ok_plane_h2o_omp
2342	ok_plane_contrail = ok_plane_contrail_omp
2343	rei_min = rei_min_omp
2344	rei_max = rei_max_omp
2345	top_height = top_height_omp
2346	overlap = overlap_omp
2347	cdmmax = cdmmax_omp
2348	cdhmax = cdhmax_omp
2349	ksta = ksta_omp
2350	ksta_ter = ksta_ter_omp
2351	f_ri_cd_min = f_ri_cd_min_omp
2352	ok_kzmin = ok_kzmin_omp
2353	fmagic = fmagic_omp
2354	pmagic = pmagic_omp
2355	iflag_pbl = iflag_pbl_omp
2356	iflag_physiq = iflag_physiq_omp
2357	iflag_pbl_split = iflag_pbl_split_omp
2358	!FC
2359	ifl_pbltree = ifl_pbltree_omp
2360	Cd_frein =Cd_frein_omp
2361	iflag_order2_sollw = iflag_order2_sollw_omp
2362	lev_histhf = lev_histhf_omp
2363	lev_histday = lev_histday_omp
2364	lev_histmth = lev_histmth_omp
2365	lev_histins = lev_histins_omp
2366	lev_histLES = lev_histLES_omp
2367	lev_histdayNMC = lev_histdayNMC_omp
2368	levout_histNMC = levout_histNMC_omp
2369	ok_histNMC(:) = ok_histNMC_omp(:)
2370	freq_outNMC(:) = freq_outNMC_omp(:)
2371	freq_calNMC(:) = freq_calNMC_omp(:)
2372
2373	type_ocean = type_ocean_omp
2374	version_ocean = version_ocean_omp
2375	t_coupl = t_coupl_omp
2376
2377	ok_veget=.TRUE.
2378	type_veget=type_veget_omp
2379	IF (type_veget=='n' .or. type_veget=='bucket' .or. type_veget=='betaclim') THEN
2380	ok_veget=.FALSE.
2381	ENDIF
2382	! INLANDSIS
2383	!=================================================
2384	landice_opt = landice_opt_omp
2385	iflag_tsurf_inlandsis = iflag_tsurf_inlandsis_omp
2386	iflag_temp_inlandsis = iflag_temp_inlandsis_omp
2387	iflag_albcalc = iflag_albcalc_omp
2388	SnoMod=SnoMod_omp
2389	BloMod=BloMod_omp
2390	ok_outfor=ok_outfor_omp
2391	is_ok_slush=is_ok_slush_omp
2392	opt_runoff_ac=opt_runoff_ac_omp
2393	is_ok_z0h_rn=is_ok_z0h_rn_omp
2394	is_ok_density_kotlyakov=is_ok_density_kotlyakov_omp
2395	prescribed_z0m_snow=prescribed_z0m_snow_omp
2396	correc_alb=correc_alb_omp
2397	iflag_z0m_snow=iflag_z0m_snow_omp
2398	ok_zsn_ii=ok_zsn_ii_omp
2399	discret_xf=discret_xf_omp
2400	buf_sph_pol=buf_sph_pol_omp
2401	buf_siz_pol=buf_siz_pol_omp
2402	!=================================================
2403	ok_all_xml = ok_all_xml_omp
2404	ok_lwoff = ok_lwoff_omp
2405	ok_newmicro = ok_newmicro_omp
2406	ok_journe = ok_journe_omp
2407	ok_hf = ok_hf_omp
2408	ok_mensuel = ok_mensuel_omp
2409	ok_instan = ok_instan_omp
2410	freq_ISCCP = freq_ISCCP_omp
2411	ecrit_ISCCP = ecrit_ISCCP_omp
2412	freq_COSP = freq_COSP_omp
2413	freq_AIRS = freq_AIRS_omp
2414	ok_ade = ok_ade_omp
2415	ok_aie = ok_aie_omp
2416	ok_alw = ok_alw_omp
2417	ok_cdnc = ok_cdnc_omp
2418	ok_volcan = ok_volcan_omp
2419	flag_volc_surfstrat = flag_volc_surfstrat_omp
2420	aerosol_couple = aerosol_couple_omp
2421	chemistry_couple = chemistry_couple_omp
2422	flag_aerosol = flag_aerosol_omp
2423	flag_aerosol_strat = flag_aerosol_strat_omp
2424	flag_aer_feedback = flag_aer_feedback_omp
2425	flag_bc_internal_mixture=flag_bc_internal_mixture_omp
2426	aer_type = aer_type_omp
2427	bl95_b0 = bl95_b0_omp
2428	bl95_b1 = bl95_b1_omp
2429	fact_cldcon = fact_cldcon_omp
2430	facttemps = facttemps_omp
2431	ratqsbas = ratqsbas_omp
2432	ratqshaut = ratqshaut_omp
2433	tau_ratqs = tau_ratqs_omp
2434
2435	iflag_radia = iflag_radia_omp
2436	iflag_rrtm = iflag_rrtm_omp
2437	iflag_albedo = iflag_albedo_omp
2438	ok_chlorophyll = ok_chlorophyll_omp
2439	NSW = NSW_omp
2440	iflag_cld_th = iflag_cld_th_omp
2441	iflag_cld_cv = iflag_cld_cv_omp
2442	tau_cld_cv = tau_cld_cv_omp
2443	coefw_cld_cv = coefw_cld_cv_omp
2444	iflag_ratqs = iflag_ratqs_omp
2445	ip_ebil_phy = ip_ebil_phy_omp
2446	iflag_thermals = iflag_thermals_omp
2447	nsplit_thermals = nsplit_thermals_omp
2448	tau_thermals = tau_thermals_omp
2449	alp_bl_k = alp_bl_k_omp
2450	! nrlmd le 10/04/2012
2451	iflag_trig_bl = iflag_trig_bl_omp
2452	s_trig = s_trig_omp
2453	tau_trig_shallow = tau_trig_shallow_omp
2454	tau_trig_deep = tau_trig_deep_omp
2455	iflag_clos_bl = iflag_clos_bl_omp
2456	! fin nrlmd le 10/04/2012
2457	iflag_coupl = iflag_coupl_omp
2458	iflag_clos = iflag_clos_omp
2459	iflag_wake = iflag_wake_omp
2460	coef_clos_ls = coef_clos_ls_omp
2461	alp_offset = alp_offset_omp
2462	iflag_cvl_sigd = iflag_cvl_sigd_omp
2463	type_run = type_run_omp
2464	ok_cosp = ok_cosp_omp
2465	ok_airs = ok_airs_omp
2466
2467	ok_mensuelCOSP = ok_mensuelCOSP_omp
2468	ok_journeCOSP = ok_journeCOSP_omp
2469	ok_hfCOSP = ok_hfCOSP_omp
2470	seuil_inversion=seuil_inversion_omp
2471	lonmin_ins = lonmin_ins_omp
2472	lonmax_ins = lonmax_ins_omp
2473	latmin_ins = latmin_ins_omp
2474	latmax_ins = latmax_ins_omp
2475	ecrit_hf = ecrit_hf_omp
2476	ecrit_ins = ecrit_ins_omp
2477	ecrit_day = ecrit_day_omp
2478	ecrit_mth = ecrit_mth_omp
2479	ecrit_tra = ecrit_tra_omp
2480	ecrit_reg = ecrit_reg_omp
2481	cvl_comp_threshold = cvl_comp_threshold_omp
2482	cvl_sig2feed = cvl_sig2feed_omp
2483	cvl_corr = cvl_corr_omp
2484	ok_lic_melt = ok_lic_melt_omp
2485	ok_lic_cond = ok_lic_cond_omp
2486	f_cdrag_ter=f_cdrag_ter_omp
2487	f_cdrag_oce=f_cdrag_oce_omp
2488
2489	f_gust_wk=f_gust_wk_omp
2490	f_gust_bl=f_gust_bl_omp
2491	f_qsat_oce=f_qsat_oce_omp
2492	f_z0qh_oce=f_z0qh_oce_omp
2493	min_wind_speed=min_wind_speed_omp
2494	iflag_gusts=iflag_gusts_omp
2495	iflag_z0_oce=iflag_z0_oce_omp
2496
2497	z0m_seaice=z0m_seaice_omp
2498	z0h_seaice=z0h_seaice_omp
2499	z0m_landice=z0m_landice_omp
2500	z0h_landice=z0h_landice_omp
2501
2502	f_rugoro=f_rugoro_omp
2503
2504	z0min=z0min_omp
2505	supcrit1 = supcrit1_omp
2506	supcrit2 = supcrit2_omp
2507	iflag_mix = iflag_mix_omp
2508	iflag_mix_adiab = iflag_mix_adiab_omp
2509	scut = scut_omp
2510	qqa1 = qqa1_omp
2511	qqa2 = qqa2_omp
2512	gammas = gammas_omp
2513	Fmax = Fmax_omp
2514	tmax_fonte_cv = tmax_fonte_cv_omp
2515	alphas = alphas_omp
2516
2517	gkdrag=sso_gkdrag_omp
2518	grahilo=sso_grahil_omp
2519	grcrit=sso_grcrit_omp
2520	gfrcrit=sso_gfrcri_omp
2521	gkwake=sso_gkwake_omp
2522	gklift=sso_gklift_omp
2523
2524	ok_strato = ok_strato_omp
2525	ok_hines = ok_hines_omp
2526	ok_gwd_rando = ok_gwd_rando_omp
2527	gwd_rando_ruwmax = gwd_rando_ruwmax_omp
2528	gwd_rando_sat = gwd_rando_sat_omp
2529	gwd_front_ruwmax = gwd_front_ruwmax_omp
2530	gwd_front_sat = gwd_front_sat_omp
2531	ok_qch4 = ok_qch4_omp
2532	ok_LES = ok_LES_omp
2533	callstats = callstats_omp
2534	ecrit_LES = ecrit_LES_omp
2535	adjust_tropopause = adjust_tropopause_omp
2536	ok_daily_climoz = ok_daily_climoz_omp
2537	carbon_cycle_tr = carbon_cycle_tr_omp
2538	carbon_cycle_cpl = carbon_cycle_cpl_omp
2539	carbon_cycle_rad = carbon_cycle_rad_omp
2540	level_coupling_esm = level_coupling_esm_omp
2541	ok_new_lscp = ok_new_lscp_omp
2542	ok_icefra_lscp=ok_icefra_lscp_omp
2543	read_fco2_ocean_cor = read_fco2_ocean_cor_omp
2544	var_fco2_ocean_cor = var_fco2_ocean_cor_omp
2545	read_fco2_land_cor = read_fco2_land_cor_omp
2546	var_fco2_land_cor = var_fco2_land_cor_omp
2547	!rajout Olivier Torres
2548	kz0=kz0_omp
2549	choix_bulk = choix_bulk_omp
2550	nit_bulk = nit_bulk_omp
2551
2552	! Test of coherence between type_ocean and version_ocean
2553	IF (type_ocean=='couple' .AND. (version_ocean/='opa8' .AND. version_ocean/='nemo') ) THEN
2554	WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid in coupled configuration'
2555	CALL abort_physic('conf_phys','version_ocean not valid',1)
2556	ENDIF
2557
2558	IF (type_ocean=='slab' .AND. version_ocean=='xxxxxx') THEN
2559	version_ocean='sicOBS'
2560	ELSE IF (type_ocean=='slab' .AND. version_ocean/='sicOBS' &
2561	.AND. version_ocean/='sicINT' .AND. version_ocean/='sicNO') THEN
2562	WRITE(lunout,*)' ERROR version_ocean=',version_ocean,' not valid with slab ocean'
2563	CALL abort_physic('conf_phys','version_ocean not valid',1)
2564	ENDIF
2565
2566	!--test on radiative scheme
2567	IF (iflag_rrtm .EQ. 0) THEN
2568	IF (NSW.NE.2) THEN
2569	WRITE(lunout,*) ' ERROR iflag_rrtm=0 and NSW<>2 not possible'
2570	CALL abort_physic('conf_phys','choice NSW not valid',1)
2571	ENDIF
2572	ELSE IF (iflag_rrtm .EQ. 1) THEN
2573	IF (NSW.NE.2.AND.NSW.NE.4.AND.NSW.NE.6) THEN
2574	WRITE(lunout,*) ' ERROR iflag_rrtm=1 and NSW<>2,4,6 not possible'
2575	CALL abort_physic('conf_phys','choice NSW not valid',1)
2576	ENDIF
2577	ELSE IF (iflag_rrtm .EQ. 2) THEN
2578	IF (NSW.NE.2.AND.NSW.NE.4.AND.NSW.NE.6) THEN
2579	WRITE(lunout,*) ' ERROR iflag_rrtm=1 and NSW<>2,4,6 not possible'
2580	CALL abort_physic('conf_phys','choice NSW not valid',1)
2581	ENDIF
2582	ELSE
2583	WRITE(lunout,*) ' ERROR iflag_rrtm<>0,1'
2584	CALL abort_physic('conf_phys','choice iflag_rrtm not valid',1)
2585	ENDIF
2586	!--here we test that solaire has not been changed if ok_suntime_rrtm is activated
2587	! IF (ok_suntime_rrtm.AND.ABS(solaire-solaire_omp_init).GT.1.E-7) THEN
2588	! WRITE(lunout,*) ' ERROR ok_suntime_rrtm=y and solaire is provided in def file'
2589	! CALL abort_physic('conf_phys','ok_suntime_rrtm=y and solaire is provided',1)
2590	! ENDIF
2591	#ifdef CPP_StratAer
2592	IF (iflag_rrtm .NE. 1) THEN
2593	WRITE(lunout,*) ' ERROR iflag_rrtm<>1 but StratAer activated'
2594	CALL abort_physic('conf_phys','iflag_rrtm not valid for StratAer',1)
2595	ENDIF
2596	IF (NSW .NE. 6) THEN
2597	WRITE(lunout,*) ' ERROR NSW<>6 but StratAer activated'
2598	CALL abort_physic('conf_phys','NSW not valid for StratAer',1)
2599	ENDIF
2600	#endif
2601
2602	!--test on ocean surface albedo
2603	IF (iflag_albedo.LT.0.OR.iflag_albedo.GT.2) THEN
2604	WRITE(lunout,*) ' ERROR iflag_albedo<>0,1'
2605	CALL abort_physic('conf_phys','choice iflag_albedo not valid',1)
2606	ENDIF
2607
2608	! Flag_aerosol cannot be set to zero if aerosol direct effect (ade) or aerosol indirect effect (aie) are activated
2609	IF (ok_ade .OR. ok_aie) THEN
2610	IF ( flag_aerosol .EQ. 0 ) THEN
2611	CALL abort_physic('conf_phys','flag_aerosol=0 not compatible avec ok_ade ou ok_aie=.TRUE.',1)
2612	ENDIF
2613	ENDIF
2614
2615	! Flag_aerosol cannot be set to zero if we are in coupled mode for aerosol
2616	IF (aerosol_couple .AND. flag_aerosol .EQ. 0 ) THEN
2617	CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if aerosol_couple=y ', 1)
2618	ENDIF
2619
2620	! Read_climoz needs to be set zero if we are in couple mode for chemistry
2621	IF (chemistry_couple .AND. read_climoz .ne. 0) THEN
2622	CALL abort_physic('conf_phys', 'read_climoz need to be to zero if chemistry_couple=y ', 1)
2623	ENDIF
2624
2625	! flag_aerosol need to be different to zero if ok_cdnc is activated
2626	IF (ok_cdnc .AND. flag_aerosol .EQ. 0) THEN
2627	CALL abort_physic('conf_phys', 'flag_aerosol cannot be to zero if ok_cdnc is activated ', 1)
2628	ENDIF
2629
2630	! ok_cdnc must be set to y if ok_aie is activated
2631	IF (ok_aie .AND. .NOT. ok_cdnc) THEN
2632	CALL abort_physic('conf_phys', 'ok_cdnc must be set to y if ok_aie is activated',1)
2633	ENDIF
2634
2635	! flag_aerosol=7 => MACv2SP climatology
2636	IF (flag_aerosol.EQ.7.AND. iflag_rrtm.NE.1) THEN
2637	CALL abort_physic('conf_phys', 'flag_aerosol=7 (MACv2SP) can only be activated with RRTM',1)
2638	ENDIF
2639	IF (flag_aerosol.EQ.7.AND. NSW.NE.6) THEN
2640	CALL abort_physic('conf_phys', 'flag_aerosol=7 (MACv2SP) can only be activated with NSW=6',1)
2641	ENDIF
2642
2643	! BC internal mixture is only possible with RRTM & NSW=6 & flag_aerosol=6 or aerosol_couple
2644	IF (flag_bc_internal_mixture .AND. NSW.NE.6) THEN
2645	CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with NSW=6',1)
2646	ENDIF
2647	IF (flag_bc_internal_mixture .AND. iflag_rrtm.NE.1) THEN
2648	CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with RRTM',1)
2649	ENDIF
2650	IF (flag_bc_internal_mixture .AND. flag_aerosol.NE.6) THEN
2651	CALL abort_physic('conf_phys', 'flag_bc_internal_mixture can only be activated with flag_aerosol=6',1)
2652	ENDIF
2653
2654	! test sur flag_volc_surfstrat
2655	IF (flag_volc_surfstrat.LT.0.OR.flag_volc_surfstrat.GT.2) THEN
2656	CALL abort_physic('conf_phys', 'flag_volc_surfstrat can only be 0 1 or 2',1)
2657	ENDIF
2658	IF ((.NOT.ok_volcan.OR..NOT.ok_ade.OR..NOT.ok_aie).AND.flag_volc_surfstrat.GT.0) THEN
2659	CALL abort_physic('conf_phys', 'ok_ade, ok_aie, ok_volcan need to be activated if flag_volc_surfstrat is 1 or 2',1)
2660	ENDIF
2661
2662	! Test on carbon cycle
2663	IF (carbon_cycle_tr .AND. .NOT. carbon_cycle_cpl) THEN
2664	CALL abort_physic('conf_phys', 'carbon_cycle_cpl has to be TRUE if carbon_cycle_tr is on',1)
2665	ENDIF
2666	IF (carbon_cycle_rad .AND. .NOT. carbon_cycle_cpl) THEN
2667	CALL abort_physic('conf_phys', 'carbon_cycle_cpl has to be TRUE if carbon_cycle_rad is on',1)
2668	ENDIF
2669
2670	! ORCHIDEE must be activated for ifl_pbltree=1
2671	IF (.NOT. ok_veget .AND. ifl_pbltree==1) THEN
2672	WRITE(lunout,*)'Warning: ORCHIDEE must be activated for ifl_pbltree=1'
2673	WRITE(lunout,*)'ifl_pbltree is now changed to zero'
2674	ifl_pbltree=0
2675	ENDIF
2676
2677	!$OMP MASTER
2678
2679	WRITE(lunout,*) ' ##############################################'
2680	WRITE(lunout,*) ' Configuration des parametres de la physique: '
2681	WRITE(lunout,*) ' Type ocean = ', type_ocean
2682	WRITE(lunout,*) ' Version ocean = ', version_ocean
2683	WRITE(lunout,*) ' Config veget = ', ok_veget,type_veget
2684	WRITE(lunout,*) ' Snow model landice : landice_opt = ', landice_opt
2685	WRITE(lunout,*) ' Config xml pour XIOS : ok_all_xml = ', ok_all_xml
2686	WRITE(lunout,*) ' Sortie journaliere = ', ok_journe
2687	WRITE(lunout,*) ' Sortie haute frequence = ', ok_hf
2688	WRITE(lunout,*) ' Sortie mensuelle = ', ok_mensuel
2689	WRITE(lunout,*) ' Sortie instantanee = ', ok_instan
2690	WRITE(lunout,*) ' Frequence appel simulateur ISCCP, freq_ISCCP =', freq_ISCCP
2691	WRITE(lunout,*) ' Frequence appel simulateur ISCCP, ecrit_ISCCP =', ecrit_ISCCP
2692	WRITE(lunout,*) ' Frequence appel simulateur COSP, freq_COSP =', freq_COSP
2693	WRITE(lunout,*) ' Frequence appel simulateur AIRS, freq_AIRS =', freq_AIRS
2694	WRITE(lunout,*) ' Sortie bilan d''energie, ip_ebil_phy =', ip_ebil_phy
2695	WRITE(lunout,*) ' Excentricite = ',R_ecc
2696	WRITE(lunout,*) ' Equinoxe = ',R_peri
2697	WRITE(lunout,*) ' Inclinaison =',R_incl
2698	WRITE(lunout,*) ' Constante solaire =',solaire
2699	WRITE(lunout,*) ' ok_suntime_rrtm =',ok_suntime_rrtm
2700	WRITE(lunout,*) ' co2_ppm =',co2_ppm
2701	WRITE(lunout,*) ' co2_ppm0 =',co2_ppm0
2702	WRITE(lunout,*) ' RCO2_act = ',RCO2_act
2703	WRITE(lunout,*) ' CH4_ppb =',CH4_ppb,' RCH4_act = ',RCH4_act
2704	WRITE(lunout,*) ' N2O_ppb =',N2O_ppb,' RN2O_act= ',RN2O_act
2705	WRITE(lunout,*) ' CFC11_ppt=',CFC11_ppt,' RCFC11_act= ',RCFC11_act
2706	WRITE(lunout,*) ' CFC12_ppt=',CFC12_ppt,' RCFC12_act= ',RCFC12_act
2707	WRITE(lunout,*) ' RCO2_per = ',RCO2_per,' RCH4_per = ', RCH4_per
2708	WRITE(lunout,*) ' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
2709	WRITE(lunout,*) ' RCFC12_per = ',RCFC12_per
2710	WRITE(lunout,*) ' cvl_comp_threshold=', cvl_comp_threshold
2711	WRITE(lunout,*) ' cvl_sig2feed=', cvl_sig2feed
2712	WRITE(lunout,*) ' cvl_corr=', cvl_corr
2713	WRITE(lunout,*) ' ok_lic_melt=', ok_lic_melt
2714	WRITE(lunout,*) ' ok_lic_cond=', ok_lic_cond
2715	WRITE(lunout,*) ' iflag_cycle_diurne=',iflag_cycle_diurne
2716	WRITE(lunout,*) ' soil_model=',soil_model
2717	WRITE(lunout,*) ' new_oliq=',new_oliq
2718	WRITE(lunout,*) ' ok_orodr=',ok_orodr
2719	WRITE(lunout,*) ' ok_orolf=',ok_orolf
2720	WRITE(lunout,*) ' zrel_oro_t=',zrel_oro_t
2721	WRITE(lunout,*) ' zstd_orodr_t=',zstd_orodr_t
2722	WRITE(lunout,*) ' zpmm_orodr_t=',zpmm_orodr_t
2723	WRITE(lunout,*) ' zpmm_orolf_t=',zpmm_orolf_t
2724	WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai
2725	WRITE(lunout,*) ' nbapp_rad=',nbapp_rad
2726	WRITE(lunout,*) ' iflag_con=',iflag_con
2727	WRITE(lunout,*) ' nbapp_cv=',nbapp_cv
2728	WRITE(lunout,*) ' nbapp_wk=',nbapp_wk
2729	WRITE(lunout,*) ' iflag_ener_conserv=',iflag_ener_conserv
2730	WRITE(lunout,*) ' ok_conserv_q=',ok_conserv_q
2731	WRITE(lunout,*) ' epmax = ', epmax
2732	WRITE(lunout,*) ' coef_epmax_cape = ', coef_epmax_cape
2733	WRITE(lunout,*) ' ok_adj_ema = ', ok_adj_ema
2734	WRITE(lunout,*) ' iflag_clw = ', iflag_clw
2735	WRITE(lunout,*) ' iflag_cld_th = ', iflag_cld_th
2736	WRITE(lunout,*) ' iflag_cld_cv = ', iflag_cld_cv
2737	WRITE(lunout,*) ' tau_cld_cv = ', tau_cld_cv
2738	WRITE(lunout,*) ' coefw_cld_cv = ', coefw_cld_cv
2739	WRITE(lunout,*) ' iflag_radia = ', iflag_radia
2740	WRITE(lunout,*) ' iflag_rrtm = ', iflag_rrtm
2741	WRITE(lunout,*) ' NSW = ', NSW
2742	WRITE(lunout,*) ' iflag_albedo = ', iflag_albedo !albedo SB
2743	WRITE(lunout,*) ' ok_chlorophyll =',ok_chlorophyll ! albedo SB
2744	WRITE(lunout,*) ' iflag_ratqs = ', iflag_ratqs
2745	WRITE(lunout,*) ' seuil_inversion = ', seuil_inversion
2746	WRITE(lunout,*) ' fact_cldcon = ', fact_cldcon
2747	WRITE(lunout,*) ' facttemps = ', facttemps
2748	WRITE(lunout,*) ' ok_newmicro = ',ok_newmicro
2749	WRITE(lunout,*) ' ratqsbas = ',ratqsbas
2750	WRITE(lunout,*) ' ratqshaut = ',ratqshaut
2751	WRITE(lunout,*) ' tau_ratqs = ',tau_ratqs
2752	WRITE(lunout,*) ' top_height = ',top_height
2753	WRITE(lunout,*) ' rad_froid = ',rad_froid
2754	WRITE(lunout,*) ' rad_chau1 = ',rad_chau1
2755	WRITE(lunout,*) ' rad_chau2 = ',rad_chau2
2756	WRITE(lunout,*) ' iflag_rei = ',iflag_rei
2757	WRITE(lunout,*) ' iflag_ice_thermo = ',iflag_ice_thermo
2758	WRITE(lunout,*) ' ok_ice_sursat = ',ok_ice_sursat
2759	WRITE(lunout,*) ' ok_plane_h2o = ',ok_plane_h2o
2760	WRITE(lunout,*) ' ok_plane_contrail = ',ok_plane_contrail
2761	WRITE(lunout,*) ' rei_min = ',rei_min
2762	WRITE(lunout,*) ' rei_max = ',rei_max
2763	WRITE(lunout,*) ' overlap = ',overlap
2764	WRITE(lunout,*) ' cdmmax = ',cdmmax
2765	WRITE(lunout,*) ' cdhmax = ',cdhmax
2766	WRITE(lunout,*) ' ksta = ',ksta
2767	WRITE(lunout,*) ' ksta_ter = ',ksta_ter
2768	WRITE(lunout,*) ' f_ri_cd_min = ',f_ri_cd_min
2769	WRITE(lunout,*) ' ok_kzmin = ',ok_kzmin
2770	WRITE(lunout,*) ' fmagic = ',fmagic
2771	WRITE(lunout,*) ' pmagic = ',pmagic
2772	WRITE(lunout,*) ' ok_ade = ',ok_ade
2773	WRITE(lunout,*) ' ok_volcan = ',ok_volcan
2774	WRITE(lunout,*) ' flag_volc_surfstrat = ',flag_volc_surfstrat
2775	WRITE(lunout,*) ' ok_aie = ',ok_aie
2776	WRITE(lunout,*) ' ok_alw = ',ok_alw
2777	WRITE(lunout,*) ' aerosol_couple = ', aerosol_couple
2778	WRITE(lunout,*) ' chemistry_couple = ', chemistry_couple
2779	WRITE(lunout,*) ' flag_aerosol = ', flag_aerosol
2780	WRITE(lunout,*) ' flag_aerosol_strat= ', flag_aerosol_strat
2781	WRITE(lunout,*) ' flag_aer_feedback= ', flag_aer_feedback
2782	WRITE(lunout,*) ' aer_type = ',aer_type
2783	WRITE(lunout,*) ' bl95_b0 = ',bl95_b0
2784	WRITE(lunout,*) ' bl95_b1 = ',bl95_b1
2785	WRITE(lunout,*) ' lev_histhf = ',lev_histhf
2786	WRITE(lunout,*) ' lev_histday = ',lev_histday
2787	WRITE(lunout,*) ' lev_histmth = ',lev_histmth
2788	WRITE(lunout,*) ' lev_histins = ',lev_histins
2789	WRITE(lunout,*) ' lev_histLES = ',lev_histLES
2790	WRITE(lunout,*) ' lev_histdayNMC = ',lev_histdayNMC
2791	WRITE(lunout,*) ' levout_histNMC = ',levout_histNMC
2792	WRITE(lunout,*) ' ok_histNMC = ',ok_histNMC
2793	WRITE(lunout,*) ' freq_outNMC = ',freq_outNMC
2794	WRITE(lunout,*) ' freq_calNMC = ',freq_calNMC
2795	WRITE(lunout,*) ' iflag_pbl = ', iflag_pbl
2796	WRITE(lunout,*) ' iflag_physiq = ', iflag_physiq
2797	!FC
2798	WRITE(lunout,*) ' ifl_pbltree = ', ifl_pbltree
2799	WRITE(lunout,*) ' Cd_frein = ', Cd_frein
2800	WRITE(lunout,*) ' iflag_pbl_split = ', iflag_pbl_split
2801	WRITE(lunout,*) ' iflag_order2_sollw = ', iflag_order2_sollw
2802	WRITE(lunout,*) ' iflag_thermals = ', iflag_thermals
2803	WRITE(lunout,*) ' iflag_clos = ', iflag_clos
2804	WRITE(lunout,*) ' coef_clos_ls = ', coef_clos_ls
2805	WRITE(lunout,*) ' type_run = ',type_run
2806	WRITE(lunout,*) ' ok_cosp = ',ok_cosp
2807	WRITE(lunout,*) ' ok_airs = ',ok_airs
2808
2809	WRITE(lunout,*) ' ok_mensuelCOSP = ',ok_mensuelCOSP
2810	WRITE(lunout,*) ' ok_journeCOSP = ',ok_journeCOSP
2811	WRITE(lunout,*) ' ok_hfCOSP =',ok_hfCOSP
2812	WRITE(lunout,*) ' solarlong0 = ', solarlong0
2813	WRITE(lunout,*) ' qsol0 = ', qsol0
2814	WRITE(lunout,*) ' evap0 = ', evap0
2815	WRITE(lunout,*) ' albsno0 = ', albsno0
2816	WRITE(lunout,*) ' iflag_sic = ', iflag_sic
2817	WRITE(lunout,*) ' iflag_inertie = ', iflag_inertie
2818	WRITE(lunout,*) ' inertie_sol = ', inertie_sol
2819	WRITE(lunout,*) ' inertie_sic = ', inertie_sic
2820	WRITE(lunout,*) ' inertie_lic = ', inertie_lic
2821	WRITE(lunout,*) ' inertie_sno = ', inertie_sno
2822	WRITE(lunout,*) ' ok_bs = ', ok_bs
2823	WRITE(lunout,*) ' ok_rad_bs = ', ok_rad_bs
2824	WRITE(lunout,*) ' f_cdrag_ter = ',f_cdrag_ter
2825	WRITE(lunout,*) ' f_cdrag_oce = ',f_cdrag_oce
2826	WRITE(lunout,*) ' f_rugoro = ',f_rugoro
2827	WRITE(lunout,*) ' z0min = ',z0min
2828	WRITE(lunout,*) ' supcrit1 = ', supcrit1
2829	WRITE(lunout,*) ' supcrit2 = ', supcrit2
2830	WRITE(lunout,*) ' iflag_mix = ', iflag_mix
2831	WRITE(lunout,*) ' iflag_mix_adiab = ', iflag_mix_adiab
2832	WRITE(lunout,*) ' scut = ', scut
2833	WRITE(lunout,*) ' qqa1 = ', qqa1
2834	WRITE(lunout,*) ' qqa2 = ', qqa2
2835	WRITE(lunout,*) ' gammas = ', gammas
2836	WRITE(lunout,*) ' Fmax = ', Fmax
2837	WRITE(lunout,*) ' tmax_fonte_cv = ', tmax_fonte_cv
2838	WRITE(lunout,*) ' alphas = ', alphas
2839	WRITE(lunout,*) ' iflag_wake = ', iflag_wake
2840	WRITE(lunout,*) ' alp_offset = ', alp_offset
2841	! nrlmd le 10/04/2012
2842	WRITE(lunout,*) ' iflag_trig_bl = ', iflag_trig_bl
2843	WRITE(lunout,*) ' s_trig = ', s_trig
2844	WRITE(lunout,*) ' tau_trig_shallow = ', tau_trig_shallow
2845	WRITE(lunout,*) ' tau_trig_deep = ', tau_trig_deep
2846	WRITE(lunout,*) ' iflag_clos_bl = ', iflag_clos_bl
2847	! fin nrlmd le 10/04/2012
2848
2849	WRITE(lunout,*) ' lonmin lonmax latmin latmax bilKP_ins =',&
2850	lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
2851	WRITE(lunout,*) ' ecrit_ hf, ins, day, mth, reg, tra, ISCCP, LES',&
2852	ecrit_hf, ecrit_ins, ecrit_day, ecrit_mth, ecrit_reg, ecrit_tra, ecrit_ISCCP, ecrit_LES
2853
2854	WRITE(lunout,*) ' ok_strato = ', ok_strato
2855	WRITE(lunout,*) ' ok_hines = ', ok_hines
2856	WRITE(lunout,*) ' ok_gwd_rando = ', ok_gwd_rando
2857	WRITE(lunout,*) ' ok_qch4 = ', ok_qch4
2858	WRITE(lunout,*) ' gwd_rando_ruwmax = ', gwd_rando_ruwmax
2859	WRITE(lunout,*) ' gwd_rando_sat = ', gwd_rando_sat
2860	WRITE(lunout,*) ' gwd_front_ruwmax = ', gwd_front_ruwmax
2861	WRITE(lunout,*) ' gwd_front_sat = ', gwd_front_sat
2862	WRITE(lunout,*) ' SSO gkdrag =',gkdrag
2863	WRITE(lunout,*) ' SSO grahilo=',grahilo
2864	WRITE(lunout,*) ' SSO grcrit=',grcrit
2865	WRITE(lunout,*) ' SSO gfrcrit=',gfrcrit
2866	WRITE(lunout,*) ' SSO gkwake=',gkwake
2867	WRITE(lunout,*) ' SSO gklift=',gklift
2868	WRITE(lunout,*) ' adjust_tropopause = ', adjust_tropopause
2869	WRITE(lunout,*) ' ok_daily_climoz = ',ok_daily_climoz
2870	WRITE(lunout,*) ' ok_new_lscp = ', ok_new_lscp
2871	WRITE(lunout,*) ' ok_icefra_lscp = ', ok_icefra_lscp
2872	WRITE(lunout,*) ' read_climoz = ', read_climoz
2873	WRITE(lunout,*) ' carbon_cycle_tr = ', carbon_cycle_tr
2874	WRITE(lunout,*) ' carbon_cycle_cpl = ', carbon_cycle_cpl
2875	WRITE(lunout,*) ' carbon_cycle_rad = ', carbon_cycle_rad
2876	WRITE(lunout,*) ' level_coupling_esm = ', level_coupling_esm
2877	WRITE(lunout,*) ' read_fco2_ocean_cor = ', read_fco2_ocean_cor
2878	WRITE(lunout,*) ' var_fco2_ocean_cor = ', var_fco2_ocean_cor
2879	WRITE(lunout,*) ' read_fco2_land_cor = ', read_fco2_land_cor
2880	WRITE(lunout,*) ' var_fco2_land_cor = ', var_fco2_land_cor
2881	WRITE(lunout,*) ' iflag_tsurf_inlandsis = ', iflag_tsurf_inlandsis
2882	WRITE(lunout,*) ' iflag_temp_inlandsis = ', iflag_temp_inlandsis
2883	WRITE(lunout,*) ' iflag_albcalc = ', iflag_albcalc
2884	WRITE(lunout,*) ' SnoMod = ', SnoMod
2885	WRITE(lunout,*) ' BloMod = ', BloMod
2886	WRITE(lunout,*) ' ok_outfor = ', ok_outfor
2887	WRITE(lunout,*) ' is_ok_slush = ', is_ok_slush
2888	WRITE(lunout,*) ' opt_runoff_ac = ', opt_runoff_ac
2889	WRITE(lunout,*) ' is_ok_z0h_rn = ', is_ok_z0h_rn
2890	WRITE(lunout,*) ' is_ok_density_kotlyakov = ', is_ok_density_kotlyakov
2891	WRITE(lunout,*) ' prescribed_z0m_snow = ', prescribed_z0m_snow
2892	WRITE(lunout,*) ' iflag_z0m_snow = ', iflag_z0m_snow
2893	WRITE(lunout,*) ' ok_zsn_ii = ', ok_zsn_ii
2894	WRITE(lunout,*) ' discret_xf = ', discret_xf
2895	WRITE(lunout,*) ' correc_alb= ', correc_alb
2896	WRITE(lunout,*) ' buf_sph_pol = ', buf_sph_pol
2897	WRITE(lunout,*) ' buf_siz_pol= ', buf_siz_pol
2898	!rajout Olivier Torres
2899	write(lunout,*) 'choix_bulk = ', choix_bulk
2900	write(lunout,*) 'nit_bulk = ', nit_bulk
2901	write(lunout,*) 'kz0 = ', kz0
2902
2903	!$OMP END MASTER
2904	call config_ocean_skin
2905
2906	END SUBROUTINE conf_phys
2907
2908	END MODULE conf_phys_m
2909	!
2910	!#################################################################
2911	!
2912
2913	SUBROUTINE conf_interface(tau_calv)
2914
2915	USE IOIPSL
2916	USE print_control_mod, ONLY: lunout
2917	IMPLICIT NONE
2918	! Configuration de l'interace atm/surf
2919	!
2920	! tau_calv: temps de relaxation pour la fonte des glaciers
2921	!
2922	REAL :: tau_calv
2923	REAL, SAVE :: tau_calv_omp
2924	!
2925	!Config Key = tau_calv
2926	!Config Desc = temps de relaxation pour fonte des glaciers en jours
2927	!Config Def = 1 an
2928	!Config Help =
2929	!
2930	tau_calv_omp = 360.*10.
2931	!$OMP MASTER
2932	CALL getin('tau_calv',tau_calv_omp)
2933	!$OMP END MASTER
2934	!$OMP BARRIER
2935	!
2936	tau_calv=tau_calv_omp
2937	!
2938	!$OMP MASTER
2939	WRITE(lunout,*)' ##############################################'
2940	WRITE(lunout,*)' Configuration de l''interface atm/surfaces : '
2941	WRITE(lunout,*)' tau_calv = ',tau_calv
2942	!$OMP END MASTER
2943	!
2944	RETURN
2945
2946	END SUBROUTINE conf_interface

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

Download in other formats:

Original Format