Ignore:
Timestamp:
Mar 20, 2014, 10:57:19 AM (11 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes r1920:1997 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/phylmd/conf_phys_m.F90

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