Ignore:
Timestamp:
07/24/17 15:28:06 (7 years ago)
Author:
jgipsl
Message:

Update physiq.def_NPv6.0.11 as done in IPSLCM6 revision [3310], [3311] and [3312].

File:
1 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/LMDZOR_v6/GENERAL/PARAM/physiq.def_NPv6.0.11

    r3301 r3323  
    22# $Id$ 
    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 
Note: See TracChangeset for help on using the changeset viewer.