Context Navigation

conf_phys_m.F90 @ 5444

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

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

Download in other formats:

Original Format