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

Last change on this file since 2154 was 2136, checked in by lguez, 10 years ago

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