Context Navigation

conf_phys_m.F90 @ 5360

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

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

Download in other formats:

Original Format