source: LMDZ5/branches/AI-cosp/libf/phylmd/conf_phys_m.F90 @ 4049

Last change on this file since 4049 was 2420, checked in by crio, 9 years ago

Nouvelle option d'epluchage de l'ascendance adiabatique dans le schema d'Emanuel: epluchage fonction de B/w2 au lieu de w. S'active avec iflag_mix_adiab=1 (valeur par defaut iflag_mix_adiab=0). Fonctionne avec iflag_mix=0 et iflag_mix=1.
Correction de bugs dans le schema de convection pour le calcul de inb, cape et buoy (sous le meme flag pour l'instant).
New option for the erosion of the adiabatic ascent in the Emanuel scheme: erosion function of B/w2 instead of w. Activated by iflag_mix_adiab=1 (standard value iflag_mix_adiab=0). Should work with iflag_mix=0 and iflag_mix=1.
Various bug corrections in the convection scheme for the computation of inb, cape, buoy (protected by the same flag for now).

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