Context Navigation

conf_phys_m.F90 @ 3437

Last change on this file since 3437 was 3435, checked in by Laurent Fairhead, 5 years ago

"Historic" :-) commit merging the physics branch used for DYNAMICO with the LMDZ trunk.
The same physics branch can now be used seamlessly with the traditional lon-lat LMDZ
dynamical core and DYNAMICO.
Testing consisted in running a lon-lat LMDZ bucket simulation with the NPv6.1 physics package
with the original trunk sources and the merged sources. Tests were succesful in the sense that
numeric continuity was preserved in the restart files from both simulation. Further tests
included running both versions of the physics codes for one year in a LMDZOR setting in which
the restart files also came out identical.

Caution:

as the physics package now manages unstructured grids, grid information needs to be transmitted

to the surface scheme ORCHIDEE. This means that the interface defined in surf_land_orchidee_mod.F90
is only compatible with ORCHIDEE version orchidee2.1 and later versions. If previous versions of
ORCHIDEE need to be used, the CPP key ORCHIDEE_NOUNSTRUCT needs to be set at compilation time.
This is done automatically if makelmdz/makelmdz_fcm are called with the veget orchidee2.0 switch

due to a limitation in XIOS, the time at which limit conditions will be read in by DYNAMICO will be

delayed by one physic timestep with respect to the time it is read in by the lon-lat model. This is caused
by the line

IF (MOD(itime-1, lmt_pas) == 0 .OR. (jour_lu /= jour .AND. grid_type /= unstructured)) THEN ! time to read

in limit_read_mod.F90

Work still needed on COSP integration and XML files for DYNAMICO

EM, YM, LF

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

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

Download in other formats:

Original Format