source: LMDZ5/trunk/libf/phylmd/conf_phys_m.F90 @ 2101

Last change on this file since 2101 was 2075, checked in by Laurent Fairhead, 10 years ago

Paramétrisation de la fréquence de couplage océan/atmosphère

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