source: LMDZ6/trunk/libf/phylmd/conf_phys_m.F90 @ 4380

Last change on this file since 4380 was 4380, checked in by evignon, 18 months ago

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