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

Last change on this file since 1990 was 1989, checked in by Laurent Fairhead, 11 years ago

Inclusion du code RRTM


Adding RRTM code

MPL

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