Ignore:
Timestamp:
01/23/18 14:50:55 (6 years ago)
Author:
acosce
Message:

update LMDZORINCA config to version 6.0.14

Location:
CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/config.def_actuel

    r3251 r3542  
    5555### ok_cdnc=y/n   Cloud droplet number concentration 
    5656ok_cdnc=_AUTO_ 
    57 ### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    58 bl95_b0=1.7 
    59 ### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    60 bl95_b1=0.2 
    6157# ok_alw=y flag pour activer l effet LW des poussieres 
    6258ok_alw=y 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/config.def_annuel

    r3158 r3542  
    2222R_incl = 23.441 
    2323### solaire =    Constante solaire 
    24 solaire = _AUTO_  
     24# solaire est modifie avec les valeurs du fichier SOLARANDVOLCANOES.txt 
     25# solaire parameter est utilise uniquement avec si iflag_rrtm=0 et ok_suntime_rrtm=n 
     26solaire = _AUTO_: DEFAULT = 0.0  
    2527### constante solaire lue dans un fichier par defaut 
    2628ok_suntime_rrtm=y 
    … …  
    5456### ok_cdnc=y/n   Cloud droplet number concentration 
    5557ok_cdnc=_AUTO_ 
    56 ### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    57 bl95_b0=1.7 
    58 ### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    59 bl95_b1=0.2 
    6058# ok_alw=y flag pour activer l effet LW des poussieres 
    6159ok_alw=y 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/config.def_preind

    r3158 r3542  
    5656### ok_cdnc=y/n   Cloud droplet number concentration 
    5757ok_cdnc=_AUTO_ 
    58 ### bl95_b0 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    59 bl95_b0=1.7 
    60 ### bl95_b1 =    Parameter in CDNC-maer link (Boucher&Lohmann 1995) 
    61 bl95_b1=0.2 
    6258# ok_alw=y flag pour activer l effet LW des poussieres 
    6359ok_alw=y 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/iodef.xml

    r2954 r3542  
    1515    <variable_definition>        
    1616      <variable id="optimal_buffer_size" type="string">performance</variable> 
    17       <variable id="buffer_size_factor" type="double">2.0</variable> 
     17      <variable id="buffer_size_factor" type="double">1.0</variable> 
    1818      <variable id="min_buffer_size" type="int">10000000</variable> 
    1919      <variable id="info_level"            type="int">1</variable> 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/orchidee.def_CWRR

    r3302 r3542  
    22#************************************************************************** 
    33#                    Parameter file list for ORCHIDEE 
     4# See all optional parameters in modipsl/modeles/ORCHIDEE/orchidee.default 
    45#************************************************************************** 
    56 
    … …  
    2627WRITE_STEP = 0 
    2728 
     29# Print level for text output in out_orchidee_* 
     30# PRINTLEV = 1 : some output in initialization phase 
     31# PRINTLEV = 2 : more output in initialization phase and printing of new date 
     32# PRINTLEV = 3, 4,.. : more output also at each timestep 
     33PRINTLEV=1 
     34 
    2835# Hydrology parameters 
    2936#************************************************************************** 
    … …  
    3643# default = n 
    3744RIVER_ROUTING = y 
     45 
     46# Activate creation of river_desc.nc file 
     47# RIVER_DESC will be activated only the first execution in the simulation.  
     48RIVER_DESC = _AUTO_ 
    3849 
    3950# SOILTYPE_CLASSIF :  Type of classification used for the map of soil types (default zobler) 
    … …  
    4354 
    4455# DO_RSOIL : Reduce soil evaporation with a soil resistance 
    45 DO_RSOIL = y 
     56DO_RSOIL = n 
     57 
     58# To make AVAN and NVAN constant with depth 
     59CWRR_AKS_A0 = 0. 
     60CWRR_AKS_POWER = 0. 
     61CWRR_NKS_A0 = 0. 
     62CWRR_NKS_POWER = 0. 
    4663 
    4764#************************************************************************** 
    4865# Activate Stomate component 
    49 # default = n 
    5066STOMATE_OK_STOMATE = _AUTOBLOCKER_  
    5167 
    … …  
    6177SECHIBA_QSINT = 0.02 
    6278 
    63 # STRESS_GS, STRESS_GM, STRESS_VCMAX : Stress on GS, GM and VCMAX 
    64 STRESS_GS=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. 
    65 STRESS_GM=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. 
    66 STRESS_VCMAX=0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,0.,1. 
     79# Use 15 PFTs. Note that PFTmap.nc must contains 15PFTs. 
     80# PFT=10, 14 and 15 share the same parametrization except for VCMAX and LAI_MAX specified below.   
     81NVM=15 
     82PFT_TO_MTC=1,2,3,4,5,6,7,8,9,10,11,12,13,10,10 
     83NAME__10='temperate C3           grass      ' 
     84NAME__14='tropical  C3           grass      ' 
     85NAME__15='boreal    C3           grass      ' 
     86 
     87# STRESS_GS, STRESS_GM, STRESS_VCMAX : Stress on GS, GM and VCMAX (15 values, without space to have short line readable) 
     88STRESS_GS=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. 
     89STRESS_GM=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. 
     90STRESS_VCMAX=1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. 
    6791 
    6892# Optimized parameters for albedo 
    69 ALB_LEAF_NIR = 0., 0.23, 0.18, 0.18, 0.20, 0.24, 0.15, 0.26, 0.20, 0.24, 0.27, 0.28, 0.26 
    70 ALB_LEAF_VIS = 0., 0.04, 0.04, 0.04, 0.04, 0.03, 0.03, 0.03, 0.03, 0.06, 0.06, 0.06, 0.06 
    71 SNOWA_AGED_NIR = 0.50, 0, 0, 0.10, 0.37, 0.08, 0.16, 0.17, 0.27, 0.44, 0.44, 0.44, 0.44 
    72 SNOWA_DEC_NIR = 0.13, 0, 0, 0.10, 0.10, 0.16, 0.04, 0.07, 0.08, 0.12, 0.12, 0.12, 0.12 
    73 SNOWA_AGED_VIS = 0.74, 0, 0, 0.08, 0.24, 0.07, 0.18, 0.18, 0.33, 0.57, 0.57, 0.57, 0.57 
    74 SNOWA_DEC_VIS = 0.21, 0, 0, 0.14, 0.08, 0.17, 0.05, 0.06, 0.09, 0.15, 0.15, 0.15, 0.15 
     93ALB_LEAF_NIR = 0.,0.23,0.18,0.18,0.20,0.24,0.15,0.26,0.20,0.24,0.27,0.28,0.26,0.24,0.24 
     94ALB_LEAF_VIS = 0.,0.04,0.04,0.04,0.04,0.03,0.03,0.03,0.03,0.06,0.06,0.06,0.06,0.06,0.06 
     95SNOWA_AGED_NIR = 0.50,0,0,0.10,0.37,0.08,0.16,0.17,0.27,0.44,0.44,0.44,0.44,0.44,0.44 
     96SNOWA_DEC_NIR = 0.13,0,0,0.10,0.10,0.16,0.04,0.07,0.08,0.12,0.12,0.12,0.12,0.12,0.12 
     97SNOWA_AGED_VIS = 0.74,0,0,0.08,0.24,0.07,0.18,0.18,0.33,0.57,0.57,0.57,0.57,0.57,0.57 
     98SNOWA_DEC_VIS = 0.21,0,0,0.14,0.08,0.17,0.05,0.06,0.09,0.15,0.15,0.15,0.15,0.15,0.15 
     99 
     100# Optimized parameters for photosynthetic capacity  
     101VCMAX25__02 = 45.0 
     102VCMAX25__03 = 45.0 
     103VCMAX25__04 = 35.0 
     104VCMAX25__05 = 40.0 
     105VCMAX25__06 = 50.0 
     106VCMAX25__07 = 45.0 
     107VCMAX25__08 = 35.0 
     108VCMAX25__09 = 35.0 
     109VCMAX25__10 = 50.0 
     110VCMAX25__11 = 50.0 
     111VCMAX25__12 = 60.0 
     112VCMAX25__13 = 60.0 
     113VCMAX25__14 = 50.0 
     114VCMAX25__15 = 40.0 
     115 
     116# Optimized parameters for LAI 
     117LAI_MAX__02 = 7.0 
     118LAI_MAX__03 = 5.0 
     119LAI_MAX__04 = 5.0 
     120LAI_MAX__05 = 4.0 
     121LAI_MAX__06 = 5.0 
     122LAI_MAX__07 = 3.5 
     123LAI_MAX__08 = 4.0 
     124LAI_MAX__09 = 3.0 
     125LAI_MAX__10 = 2.5 
     126LAI_MAX__11 = 2.0 
     127LAI_MAX__12 = 5.0 
     128LAI_MAX__13 = 5.0 
     129LAI_MAX__14 = 2.5 
     130LAI_MAX__15 = 2.0 
     131 
     132# Activate downregulation for CO2 
     133DOWNREGULATION_CO2=y 
     134DOWNREGULATION_CO2_BASELEVEL=380. 
     135 
     136# Activate harvest of wood 
     137DO_WOOD_HARVEST=y 
     138 
     139# Growth respiration 
     140FRAC_GROWTHRESP__02 = 0.35 
     141FRAC_GROWTHRESP__03 = 0.35 
     142FRAC_GROWTHRESP__04 = 0.28 
     143FRAC_GROWTHRESP__05 = 0.28 
     144FRAC_GROWTHRESP__06 = 0.28 
     145FRAC_GROWTHRESP__07 = 0.35 
     146FRAC_GROWTHRESP__08 = 0.35 
     147FRAC_GROWTHRESP__09 = 0.35 
     148FRAC_GROWTHRESP__10 = 0.28 
     149FRAC_GROWTHRESP__11 = 0.28 
     150FRAC_GROWTHRESP__12 = 0.28 
     151FRAC_GROWTHRESP__13 = 0.28 
     152FRAC_GROWTHRESP__14 = 0.35 
     153FRAC_GROWTHRESP__15 = 0.35 
     154 
     155# Maintenance respiration slope C 
     156MAINT_RESP_SLOPE_C__02 = 0.12 
     157MAINT_RESP_SLOPE_C__03 = 0.12 
     158MAINT_RESP_SLOPE_C__04 = 0.16 
     159MAINT_RESP_SLOPE_C__05 = 0.16 
     160MAINT_RESP_SLOPE_C__06 = 0.16 
     161MAINT_RESP_SLOPE_C__07 = 0.25 
     162MAINT_RESP_SLOPE_C__08 = 0.25 
     163MAINT_RESP_SLOPE_C__09 = 0.25 
     164MAINT_RESP_SLOPE_C__10 = 0.16 
     165MAINT_RESP_SLOPE_C__11 = 0.12 
     166MAINT_RESP_SLOPE_C__12 = 0.16 
     167MAINT_RESP_SLOPE_C__13 = 0.12 
     168MAINT_RESP_SLOPE_C__14 = 0.12 
     169MAINT_RESP_SLOPE_C__15 = 0.25 
     170 
     171# Soil turnover 
     172CARBON_TAU_ISLOW = 7 
     173CARBON_TAU_IPASSIVE = 300 
    75174 
    76175# Parameters related to vegetation map 
    … …  
    96195FIRE_DISABLE = y 
    97196 
    98 # Create river description file 
    99 RIVER_DESC = _AUTO_ 
    100 RIVER_DESC_FILE =  river_desc.nc 
    101  
    102197# Parameters related to surface and thermal physical properties 
    103198#************************************************************************ 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/orchidee.def_Choi

    r2770 r3542  
    6060SECHIBA_QSINT = 0.02 
    6161 
     62# Use 15 PFTs. Note that PFTmap.nc must contains 15PFTs. 
     63# PFT=10, 14 and 15 share the same parametrization.   
     64NVM=15 
     65PFT_TO_MTC=1,2,3,4,5,6,7,8,9,10,11,12,13,10,10 
     66 
     67# Set NVM and PFT_TO_MTC as follow to run with PFTmap files containing 13 pfts 
     68# Do this for example to run with the CMIP5 vegetation files 
     69#NVM=13 
     70#PFT_TO_MTC=1,2,3,4,5,6,7,8,9,10,11,12,13 
     71 
    6272# Parmeters related to vegetation map 
    6373#************************************************************************** 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/physiq.def_NPv6.0.11

    r3316 r3542  
    11# 
    2 # $Id: physiq.def_NPv6.0.11 3301 2017-07-07 09:46:45Z jgipsl $ 
     2# $Id: physiq.def_NPv6.0.11 3323 2017-07-24 13:28:06Z jgipsl $ 
    33# 
    4 ############################################################################## 
     4 
     5############################################################## 
    56# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ 
    6 ############################################################################## 
     7############################################################## 
    78# D: Valeur par default 
    89# 
    9 #Parametres Orographiques et cdrags 
    10 #********************************** 
     10 
     11########################################################### 
     12# divers. Devrait aller dans config.def ? 
     13########################################################### 
    1114##   Si=.T. ,  lecture du fichier limit avec la bonne annee  
    1215ok_limitvrai=n 
    13 iflag_albedo=1 
    14 #Cdrags  
     16 
     17########################################################### 
     18# Parametres couche limite de surface 
     19########################################################### 
     20 
     21          #  ... 
    1522f_cdrag_stable=1. 
     23 
     24          # Facteur sur les coefficients de trainee terrestres 
    1625f_cdrag_ter=1. 
    1726 
    18 #parameter for Z0 
     27          # Vitesse minimum dans les calculs de trainee Cd ( || V || + Vmin ) 
     28min_wind_speed=0.1 
     29 
     30########################################################### 
     31# parameter des flux air-mer 
     32########################################################### 
     33 
     34          # Facteur sur les flux sensibles et latent 
     35          # Introduit pour reduire l evaporation quand on avait z0h=z0m 
     36          # CONSEILLE : 1 
    1937f_cdrag_oce=1. 
     38 
     39          # Flag controlant le calcul des z0 
    2040iflag_z0_oce=1 
     41 
     42          # Rapport entre les Z0 latent et sensible 
    2143f_z0qh_oce=1.55 
    22 min_wind_speed=0.1 
    23  
     44 
     45          # Activation des rafales : Oui 1 ou non 0 
     46iflag_gusts=0 
     47 
     48          # facteur pour les rafales des thermiques (bounary layer) 
     49f_gust_bl=0.845 
     50 
     51          # facteur pour les poches (wake) 
     52f_gust_wk=0.2 
     53 
     54          # facteur sur qsat pour tenir compte de la salinite 
     55f_qsat_oce=1. 
     56 
     57          #cdmmax : valeur max de Cd moment. Sur ocean seulement 
     58cdmmax=2.5E-3 
     59 
     60          #cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement 
     61cdhmax=2.0E-3 
     62 
     63#################################################################### 
     64# Couche limite turbulente 
     65#################################################################### 
     66 
     67          # Nouvelle version de Mellor et Yamada 
    2468new_yamada4=y 
     69 
     70          # Choix du schema numerique pour new_yamada4=y 
     71          #  1 schema de MAR. Bien pour CL stables mais detruit les stratocu 
     72          #  5 schema de MAR modifie. Precaunise. 
    2573yamada4_num=5 
     74 
     75          # Flag de controle de la couche limite stable 
    2676iflag_corr_sta=4 
     77 
     78          # min sur les fonction de stabilite a la surface 
    2779f_ri_cd_min=0.01 
     80 
     81          # max du Ric pour Kz. Decouplage plus grand pour Ric plus grand. 
    2882yamada4_ric=0.18 
     83 
     84          # Longueur de melange minimum pour Kz 
    2985lmixmin=0. 
    30 #cdmmax 
    31 cdmmax=2.5E-3 
    32 #cdhmax 
    33 cdhmax=2.0E-3 
    34 ##  Orodr  ou  non   pour l orographie 
     86 
     87          #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)  
     88iflag_pbl=12 
     89 
     90          # Seuils pour la diffusion turbulente 
     91ksta_ter=1.e-7 
     92ksta=1.e-10 
     93 
     94          #ok_kzmin : calcul Kzmin dans la CL de surface (D: y) 
     95ok_kzmin=n 
     96 
     97          # Creation de TKE par orographie sous-maille 
     98addtkeoro=0 
     99alphatkeoro=1. 
     100smallscales_tkeoro=n 
     101 
     102 
     103#################################################################### 
     104# Ondes de gravite / Orographie sous maille 
     105#################################################################### 
     106 
     107          #  Orodr  ou  non   pour l orographie 
    35108ok_orodr=y 
    36 ##  Orolf  ou  non   pour l orographie               
     109 
     110          #  Orolf  ou  non   pour l orographie               
    37111ok_orolf=y 
    38 #Rugoro 
     112 
     113          #Rugoro : inclusion de l orographie sous maille dans le Cd de surface 
     114          # Obsolete avec l activation des ondes de gravite 
    39115f_rugoro=0. 
    40 # 
    41 ### Random gravity waves 
    42 #Activitate the stochastic parameterization of convective GWs (should be equal to ok_strato) 
    43 ok_gwd_rando=y 
    44 # Scale the momentum fluxes due to convective GWs 
    45 gwd_rando_ruwmax=0.23 
    46 #Scale the saturated momentum flux of  convective GWs  
    47 gwd_rando_sat=0.60 
    48 # 
    49 gwd_front_ruwmax = 4.00 
    50  
     116 
     117          # Parametres controlant l effet de l orographie sous maille 
    51118sso_gfrcri=0.70 
    52119sso_gkdrag=0.6 
    53120sso_gkwake=0.4 
    54121sso_gklift=0.1 
    55 # 
     122 
     123          # Random gravity waves 
     124          #Activitate the stochastic parameterization of convective GWs 
     125          # (should be equal to ok_strato) 
     126ok_gwd_rando=y 
     127 
     128          # Scale the momentum fluxes due to convective GWs 
     129gwd_rando_ruwmax=0.23 
     130 
     131          #Scale the saturated momentum flux of  convective GWs  
     132gwd_rando_sat=0.60 
     133 
     134          # Important pour controler la periode de la QBO 
     135gwd_front_ruwmax = 3.00 
     136 
     137 
     138############################################################### 
    56139# Rayonnement 
    57 #********************* 
    58 # activation du nouveau code de rayonnement RRTM 
    59 # 0 : Ancien code et 1 : RRTM (D=0) 
     140############################################################### 
     141          # activation du nouveau code de rayonnement RRTM 
     142          # 0 : Ancien code et 1 : RRTM (D=0) 
    60143iflag_rrtm=1 
     144 
     145          # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0 
    61146NSW=6 
    62 # 
     147 
     148    # Albedo et inertie des calottes 
     149alb_vis_sno_lic=0.96 
     150alb_nir_sno_lic=0.68 
     151inertie_sno=350. 
     152 
     153          # ??? 
     154iflag_albedo=1 
     155 
     156          # Frequence appel rayonnement. Nombre appels par jour 
     157nbapp_cv=48 
     158 
     159          # Pour pouvoir revenir a un bug purement informatique de duplication 
     160          # des constantes thermodynamiques entre LMDZ et RRTM 
     161          # n : FORTEMENT CONSEILLE 
     162ok_bad_ecmwf_thermo=n 
     163 
     164############################################################### 
    63165# Parametres nuages 
    64 #****************** 
    65 #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) 
    66 cld_lc_lsc=0.0000783 
    67 cld_lc_con=0.0000783 
    68 #constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.) 
     166############################################################### 
     167 
     168          #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) 
     169cld_lc_lsc=0.000075 
     170cld_lc_con=0.000075 
     171 
     172          #constante de temps pour eleminer eau lsc et convective 
     173          # (D: 3600. 3600.) 
    69174cld_tau_lsc=900. 
    70175cld_tau_con=900. 
    71 #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) 
     176 
     177          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) 
    72178ffallv_lsc=0.6 
    73179ffallv_con=0.6 
    74 #coefficient sur la reevaporation de la pluie (D: 2.e-5 n) 
     180 
     181          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n) 
     182          # regle a 3.e-5 sur des cas de cumulus en 1D 
    75183coef_eva=1.e-4 
     184 
     185          # flag pour la reevaporation de la pluie 
     186          # 0: rien 
     187          # 1: a l ancienne 
     188          # 2: prend en compte la surface nuageuse max au dessus pour calculer 
     189          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE 
     190iflag_evap_prec=1 
     191 
     192          # Modification de la plage de temperature pour la phase mixte 
     193          # liquide/glace des nuages 
     194          # Controle par les t_glace_min/max, exposant_glace, 
     195          # iflag_t_glace (D=0) 
     196t_glace_min=243.15 
     197t_glace_max=273.15 
     198exposant_glace=2. 
     199iflag_t_glace=2 
     200 
     201          # Activation de la thermodynamique de la glace 
     202iflag_ice_thermo=1 
     203 
     204          # normalement elle ne sert plus a rien dans fisrtilp 
     205          # sinon elle est un peu toxique: 
     206          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste) 
    76207reevap_ice=y 
    77 #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)  
     208 
     209          #Bergeron 
     210iflag_bergeron=2 
     211 
     212          # temperature dependency of qsat for iterative procedure 
     213          # temporary flag to switch between the old (CMIP5, D=0) 
     214          # and new versions (=1) versions  
     215iflag_fisrtilp_qsat=4  
     216 
     217          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)  
    78218iflag_cldcon=6 
    79219fact_cldcon=1. 
     220 
     221          # ????? 
    80222facttemps=0. 
    81 #calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs)   
     223 
     224          #calcul eau condensee et fraction nuageuse a partir des PDF 
     225          # (D:0, 0:version avec ratqs sinon nouvelles PDFs)   
    82226iflag_pdf=1 
    83 #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) 
     227 
     228          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) 
    84229ok_newmicro=y 
     230 
     231          # Controle des ratqs grande echelle 
    85232iflag_ratqs=4 
    86233ratqsp0=45000. 
    … …  
    88235ratqsbas=0.002 
    89236ratqshaut=0.4 
    90 #rayon effectif pour les nuages de glace (D:35) 
     237 
     238          #rayon effectif pour les nuages de glace (D:35) 
    91239rad_froid=35 
    92 #rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9) 
     240 
     241          #rayons effectifs pour les nuages eau liq (tailles des gouttes 
     242          # d eau nuageuse, D: 13 9) 
    93243rad_chau1=12 
    94244rad_chau2=11 
    95 ##  Choix ou non  de  New oliq                
     245 
     246          #  Choix ou non  de  New oliq                
    96247new_oliq=y 
    97 # Seuils oliqmax, oicemax 
     248 
     249          # Seuils oliqmax, oicemax (en g/kg). 
     250          # Au dela on precipite directement au sol 
    98251oliqmax=1.5E-3 
    99252oicemax=1.5E-3 
    100 #Rayons min et max effectifs des particules des nuages 
     253 
     254          #Rayons min et max effectifs des particules des nuages 
    101255rei_min=16. 
    102256rei_max=61.29 
    103 # 
    104 # Flags Convection  
    105 #***************** 
    106 #Convection scheme switch (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4) 
     257 
     258          # Distribution vetricale sous maille pour le schema de de nuages. 
     259          # 0: rien 
     260          # 1: Premiere version A. Jam, 
     261          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5 
     262          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32 
     263          #                           cloudth_vert_alpha_th C=2:91 
     264iflag_cloudth_vert=3 
     265cloudth_vert_alpha=0.5 
     266cloudth_vert_alpha_th=0.5 
     267 
     268 
     269###################################################################### 
     270# Convection  
     271###################################################################### 
     272 
     273          #Convection scheme switch 
     274          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4) 
    107275iflag_con=3 
    108 #output level of energy conservation diagnostics 
     276 
     277          #output level of energy conservation diagnostics 
    109278if_ebil=0 
    110 #maximum efficiency of cld water->precipitation conversion (D: 0.993) 
     279 
     280          #maximum efficiency of cld water->precipitation conversion (D: 0.993) 
    111281epmax=0.9985 
    112 #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0) 
     282 
     283          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0) 
    113284ok_adj_ema=n 
    114285iflag_clw=0 
    115 #Convective closure (D:1, 0:AR4, 2:ALE and ALP) 
     286 
     287          #Convective closure (D:1, 0:AR4, 2:ALE and ALP) 
    116288iflag_clos=2 
    117 #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF) 
     289 
     290          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF) 
    118291iflag_mix=1 
    119 #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0) 
     292 
     293          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0) 
    120294qqa1=1. 
    121295qqa2=0. 
    122 #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05) 
     296 
     297          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05) 
    123298cvl_corr=1.0 
    124 #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.) 
     299 
     300          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.) 
    125301Fmax=0.65 
    126 # 
     302          # 
    127303t_top_max=255. 
    128304# 
    129 # Flags Couche limite 
    130 #******************** 
    131 #shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques)  
     305### Convection 
     306          #altitude, in hPa, above LCL at which buoybase is computed, 
     307          # def = original: -40 
     308dpbase=-40. 
     309 
     310          #reference fractional area of precipitating downdraughts, 
     311          # def = original: 0.01 
     312sigdz=0.01 
     313 
     314          #unsaturated fraction of precipitating downdraughts, def = original: 0.15 
     315spfac=0.15 
     316 
     317          #relaxation time of cloud base mass flux, def = original: 8000 
     318tau=8000. 
     319 
     320          #flag for wb (= vert velocity at LFC); 
     321          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1 
     322          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1 
     323          # wbmax : assymptotic value  
     324flag_wb=1 
     325wbmax=6. 
     326 
     327          # Divers tests 
     328ok_convstop=.False. 
     329tau_stop=15000. 
     330ok_intermittent=.False. 
     331 
     332          # nouvel epluchage de C.Rio 
     333          #shedding coefficient (used when iflag_mix_adiab=1) 
     334iflag_mix_adiab=0 
     335coef_peel=0.25 
     336 
     337          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) 
     338          #if 1, def=1 
     339flag_epKEorig=1 
     340 
     341          # Case : flag_epKEorig=1 
     342          # elcrit : max condensed water concentration in kg/kg, D=0.0003, 
     343          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original) 
     344          # original: 0.0011, used only if flag_epKEorig=1 
     345elcrit=0.00025 
     346tlcrit=-40. 
     347 
     348          # Case : flag_epKEorig=0 
     349          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original) 
     350          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original) 
     351pbcrit=150.0 
     352ptcrit=500.0 
     353 
     354          # Prevent some crashes 
     355          # Filter out bad wakes 
     356flag_wk_check_trgl=y 
     357 
     358          # Check whether parcel level temperature and specific humidity 
     359          # are reasonable and demand that LCL be high enough 
     360cv_flag_feed=2 
     361 
     362          # save 20% CPU  
     363ok_optim_yield=y 
     364 
     365 
     366########################################################## 
     367# Fermeture convective et couplages 
     368########################################################## 
     369 
     370          # Appel separe des thermiques  
     371iflag_pbl_split=0 
     372 
     373          #rendement sur la puissance fournie par les thermiques a la convection 
     374alp_bl_k=0.5 
     375 
     376          #Constant value to be added to Alp (D: 0.) 
     377alp_offset=0. 
     378 
     379          ### nrlmd le 10/04/2012: Declenchement ### 
     380          #Declenchement par couche-limite (0: deterministe, 1: stochastique) 
     381iflag_trig_bl=1 
     382 
     383          #Surface critique au LCL 
     384s_trig=1.2e7 
     385 
     386          #Temps de decorrelation entre scenes de cumulus 
     387tau_trig_shallow=1200 
     388 
     389          #Temps de decorrelation entre scenes de cumulonimbus 
     390tau_trig_deep=1200 
     391 
     392          # Pour la stabilisation 
     393random_notrig_max=0.99 
     394 
     395          # Temperature a la quelle on fond les descentes glacees 
     396tmax_fonte_cv=278.15 
     397 
     398          #Fermeture par couche-limite (Nicolas Rochetin) 
     399          #(0: deterministe, 1: statistique, 2: statistique+convergence 
     400          # grande echelle) 
     401iflag_clos_bl=1 
     402 
     403          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape), 
     404          # def: 0.33 
     405stark=0.33 
     406 
     407          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided 
     408          # by wakes, def: 0.25 
     409alpk=0.25 
     410 
     411# 
     412########################################################## 
     413# Couche limite convective / modele du thermique 
     414########################################################## 
     415 
     416          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)  
    132417iflag_thermals=18 
    133 #pas de temps du spliting pour les thermiques 
     418 
     419          # pas de temps du spliting pour les thermiques 
     420          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD 
    134421nsplit_thermals=1 
     422 
     423          # tau_thermals pour avoir une constante de temps sur les thermiques. 
     424          # Non valide 
    135425tau_thermals=0. 
     426 
     427          # Flag controlant le detrainement et l entrainement 
    136428iflag_thermals_ed=8 
     429 
     430          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer 
     431          # le detrainement (A. Jam) 
    137432fact_thermals_ed_dz=0.07 
     433 
     434          # Une option de plus sur les thermiques. 
    138435iflag_thermals_optflux=0 
    139 #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)  
    140 iflag_pbl=12 
    141 #diffusion turbulente 
    142 ksta_ter=1.e-7 
    143 ksta=1.e-10 
    144 #ok_kzmin : calcul Kzmin dans la CL de surface 
    145 ### ok_kzmin=y 
    146 #couplage ou non avec la convection 
     436 
     437          #couplage ou non avec la convection 
    147438iflag_coupl=5 (0:AR4 et 1:nouvelle physique) 
     439 
     440          # Seuil pour desactiver les thermiques en cas de forte inversion 
     441          # en sommet de couche limite. Actif si iflag_thermals impair. 
    148442seuil_inversion=-0.08 
    149 #rendement sur la puissance fournie par les thermiques a la convection 
    150 alp_bl_k=0.5 
    151 # 
    152 ## temperature dependency of qsat for iterative procedure 
    153 ## temporary flag to switch between the old (CMIP5, D=0) and new versions (=1) versions  
    154 iflag_fisrtilp_qsat=4  
    155 ## Reactivation du calcul d un zmax continu pour les thermiques 
    156 ## reactivation d une fermeture humide des thermiques  iflag_thermals_closure=2 (D=1) 
     443 
     444          # Reactivation du calcul d un zmax continu pour les thermiques 
     445          # reactivation d une fermeture humide des thermiques 
     446          # iflag_thermals_closure=2 (D=1) 
    157447iflag_thermals_closure=2 
    158 # 
     448 
     449 
     450########################################################## 
    159451# Flags wakes 
    160 #************ 
    161 #Wake scheme switch (D:0, 0:AR4, 1:New Physics) 
     452########################################################## 
     453 
     454          #Wake scheme switch (D:0, 0:AR4, 1:New Physics) 
    162455iflag_wake=1 
    163 #Constant value to be added to Alp (D: 0.) 
    164 alp_offset=0. 
    165  
    166 ### nrlmd le 10/04/2012: Declenchement ### 
    167 #Declenchement par couche-limite (0: deterministe, 1: stochastique) 
    168 iflag_trig_bl=1 
    169 #Surface critique au LCL 
    170 s_trig=1.2e7 
    171 #Temps de decorrelation entre scenes de cumulus 
    172 tau_trig_shallow=1200 
    173 #Temps de decorrelation entre scenes de cumulonimbus 
    174 tau_trig_deep=1200 
    175 # Pour la stabilisation 
    176 random_notrig_max=0.99 
    177 ok_adjwk=y 
    178 tmax_fonte_cv=278.15 
    179 # 
    180  
    181 ### nrlmd le 10/04/2012: Fermeture ### 
    182 #Fermeture par couche-limite (0: deterministe, 1: statistique, 2: statistique+convergence grande echelle) 
    183 iflag_clos_bl=1 
    184 ### fin nrlmd le 10/04/2012 
    185  
    186 ## Modification de la plage de temperature pour la phase mixte liquide/glace des nuages 
    187 ## Controle par les t_glace_min/max, exposant_glace, iflag_t_glace (D=0) 
    188 t_glace_min=243.15 
    189 t_glace_max=273.15 
    190 exposant_glace=2. 
    191 iflag_t_glace=2 
    192 iflag_ice_thermo=1 
    193  
    194 # 
    195 ## The source terms of the TKE prognostic equation are diagnosed from tendencies (du, dv, dT)  
    196 ## associated with subrgrid scaleamotions and treated as an additional heat source. 
    197 ## Controled by a key, iflag_ener_conserv (D=-1): 
    198 ## 0 : no conservation, -1 : old adhoc correction for kinetic E only (used for CMIP5)  
    199 ## 1 : conservation, 101 : conversion from kinetic to heat only, 110 conversion from potential to heat only 
    200 iflag_ener_conserv=1 
    201  
    202 #Activation of corrections for water conservation (D: n).  
    203 ##The default is n for compatibility with earlier versions. 
    204 ##generally it should be y. 
    205 ok_conserv_q=y  
    206  
    207 iflag_cloudth_vert=3 
    208  
    209 # 
    210  
    211 #Bergeron 
    212 iflag_bergeron=2 
    213  
    214  
    215 ### Convection 
    216 #altitude, in hPa, above LCL at which buoybase is computed, def = original: -40 
    217 dpbase=-40. 
    218 #altitude, in hPa, above LCL below which Ep=0, def = original: 150, used only if flag_epKEorig=0 
    219 pbcrit=150.0 
    220 #altitude, in hPa, above LCL above which Ep=Epmax, def = original: 500, used only if flag_epKEorig=0 
    221 ptcrit=500.0 
    222 #reference fractional area of precipitating downdraughts, def = original: 0.01 
    223 sigdz=0.01 
    224 #unsaturated fraction of precipitating downdraughts, def = original: 0.15 
    225 spfac=0.15 
    226 #relaxation time of cloud base mass flux, def = original: 8000 
    227 tau=8000. 
    228 #flag for wb (= vert velocity at LFC); 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, def = 1 
    229 flag_wb=1 
    230 #parameter used in the computation of wb, def = 6. 
    231 wbmax=6. 
    232 # 
    233 ok_convstop=.False. 
    234 # 
    235 tau_stop=15000. 
    236 # 
    237 ok_intermittent=.False. 
    238 # nouvel epluchage de C.Rio 
    239 iflag_mix_adiab=0 
    240 #shedding coefficient (used when iflag_mix_adiab=1) 
    241 coef_peel=0.25 
    242 # 
    243 ### Clouds 
    244 #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) if 1, def=1 
    245 flag_epKEorig=1 
    246 #max condensed water concentration in kg/kg, def=0.0003, original: 0.0011, used only if flag_epKEorig=1 
    247 elcrit=0.00025 
    248 #temperature in degre Celsius, at which Ep=Epmax, def=-55., original: -55, used only if flag_epKEorig=1 
    249 tlcrit=-40. 
    250  
    251 ### Wakes 
    252 #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape), def: 0.33 
    253 stark=0.33 
    254 #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided by wakes, def: 0.25 
    255 alpk=0.25 
    256 #wake density = number of wake centers per m2, def: 8.E-12 
    257 #wdens_ref=8.E-12 
     456 
     457          #multiplicative factor of the damping by gravity waves, def: 4. 
     458coefgw=4. 
     459 
     460          #wake density = number of wake centers per m2, def: 8.E-12 
     461          #wdens_ref=8.E-12 
    258462wdens_ref_o=1.e-9 
    259463wdens_ref_l=8.e-12 
    260 # 
    261 #multiplicative factor of the damping by gravity waves, def: 4. 
    262 coefgw=4. 
    263 # 
    264 alb_vis_sno_lic=0.96 
    265 alb_nir_sno_lic=0.68 
    266 inertie_sno=350. 
    267  
    268 nbapp_cv=48 
    269 # Prevent some crashes 
    270 # Filter out bad wakes 
    271 flag_wk_check_trgl=y 
    272 # Check whether parcel level temperature and specific humidity 
    273 # are reasonable and demand that LCL be high enough 
    274 cv_flag_feed=2 
    275 # 
    276 # 
    277 # Correction sur ctes passees a rrtm 
    278 ok_bad_ecmwf_thermo=n 
    279 # Conservation de l energie 
     464 
     465          # Ajustement convectif prealable au calcul des poches 
     466ok_adjwk=y 
     467 
     468 
     469#################################################################### 
     470#  Conservation  
     471#################################################################### 
     472 
     473          # The source terms of the TKE prognostic equation are diagnosed 
     474          # from tendencies (du, dv, dT)  
     475          # associated with subrgrid scaleamotions and treated as an  
     476          # additional heat source. 
     477          # Controled by a key, iflag_ener_conserv (D=-1): 
     478          # 0 : no conservation, -1 : old adhoc correction for kinetic E only 
     479          # (used for CMIP5)  
     480          # 1 : conservation, 101 : conversion from kinetic to heat only, 
     481          # 110 conversion from potential to heat only 
     482iflag_ener_conserv=1 
     483 
     484          #Activation of corrections for water conservation (D: n).  
     485          ##The default is n for compatibility with earlier versions. 
     486          ##generally it should be y. 
     487ok_conserv_q=y  
     488 
    280489fl_cor_ebil=1 
    281 # 
    282 # save 20% CPU  
    283 ok_optim_yield=y 

  • CONFIG/UNIFORM/v6/LMDZORINCA_v6/GENERAL/PARAM/run.def

    r3395 r3542  
    5454### Only read by ce0l : Do not set land points on Antartic 
    5555no_ter_antartique=y 
    56 ### Choose to run with tracers  
    57 iflag_phytrac=1 
    5856 
    5957 
Note: See TracChangeset for help on using the changeset viewer.