| 1 | ############################################################################## |
|---|
| 2 | # PARAMETRES ET CLES DE LA PHYSIQUE LMDZ |
|---|
| 3 | ############################################################################## |
|---|
| 4 | # D: Valeur par default |
|---|
| 5 | # |
|---|
| 6 | #Parametres Orographiques et cdrags |
|---|
| 7 | #********************************** |
|---|
| 8 | ## Si=.T. , lecture du fichier limit avec la bonne annee |
|---|
| 9 | ok_limitvrai=n |
|---|
| 10 | iflag_albedo=1 |
|---|
| 11 | #Cdrags |
|---|
| 12 | f_cdrag_stable=1. |
|---|
| 13 | f_cdrag_ter=1. |
|---|
| 14 | f_cdrag_oce=0.6 |
|---|
| 15 | f_ri_cd_min=0.1 |
|---|
| 16 | #cdmmax |
|---|
| 17 | cdmmax=2.5E-3 |
|---|
| 18 | #cdhmax |
|---|
| 19 | cdhmax=2.0E-3 |
|---|
| 20 | ## Orodr ou non pour l orographie |
|---|
| 21 | ok_orodr=y |
|---|
| 22 | ## Orolf ou non pour l orographie |
|---|
| 23 | ok_orolf=y |
|---|
| 24 | #Rugoro |
|---|
| 25 | f_rugoro=0. |
|---|
| 26 | # |
|---|
| 27 | ### Random gravity waves |
|---|
| 28 | #Activitate the stochastic parameterization of convective GWs (should be equal to ok_strato) |
|---|
| 29 | ok_gwd_rando=y |
|---|
| 30 | # Scale the momentum fluxes due to convective GWs |
|---|
| 31 | gwd_rando_ruwmax=2.5 |
|---|
| 32 | #Scale the saturated momentum flux of convective GWs |
|---|
| 33 | gwd_rando_sat=0.2 |
|---|
| 34 | |
|---|
| 35 | sso_gkdrag=0.2 |
|---|
| 36 | sso_gfrcri=0.70 |
|---|
| 37 | sso_gkwake=0.40 |
|---|
| 38 | # |
|---|
| 39 | # Rayonnement |
|---|
| 40 | #********************* |
|---|
| 41 | # activation du nouveau code de rayonnement RRTM |
|---|
| 42 | # 0 : Ancien code et 1 : RRTM (D=0) |
|---|
| 43 | iflag_rrtm=1 |
|---|
| 44 | NSW=6 |
|---|
| 45 | # |
|---|
| 46 | # Parametres nuages |
|---|
| 47 | #****************** |
|---|
| 48 | #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) |
|---|
| 49 | cld_lc_lsc=0.0004 |
|---|
| 50 | cld_lc_con=0.0004 |
|---|
| 51 | #constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.) |
|---|
| 52 | cld_tau_lsc=900. |
|---|
| 53 | cld_tau_con=900. |
|---|
| 54 | #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) |
|---|
| 55 | ffallv_lsc=.5 |
|---|
| 56 | ffallv_con=.5 |
|---|
| 57 | #coefficient sur la reevaporation de la pluie (D: 2.e-5 n) |
|---|
| 58 | coef_eva=2e-4 |
|---|
| 59 | reevap_ice=y |
|---|
| 60 | #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) |
|---|
| 61 | iflag_cldcon=6 |
|---|
| 62 | fact_cldcon=1. |
|---|
| 63 | facttemps=0. |
|---|
| 64 | #calcul eau condensee et fraction nuageuse a partir des PDF (D:0, 0:version avec ratqs sinon nouvelles PDFs) |
|---|
| 65 | iflag_pdf=1 |
|---|
| 66 | #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) |
|---|
| 67 | ok_newmicro=y |
|---|
| 68 | iflag_ratqs=4 |
|---|
| 69 | ratqsbas=0.002 |
|---|
| 70 | ratqshaut=0.40 |
|---|
| 71 | ratqsp0=50000. |
|---|
| 72 | ratqsdp=20000. |
|---|
| 73 | #rayon effectif pour les nuages de glace (D:35) |
|---|
| 74 | rad_froid=35 |
|---|
| 75 | #rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9) |
|---|
| 76 | rad_chau1=12 |
|---|
| 77 | rad_chau2=11 |
|---|
| 78 | ## Choix ou non de New oliq |
|---|
| 79 | new_oliq=y |
|---|
| 80 | #Rayons min et max effectifs des particules des nuages |
|---|
| 81 | rei_min=16. |
|---|
| 82 | rei_max=61.29 |
|---|
| 83 | # |
|---|
| 84 | # Flags Convection |
|---|
| 85 | #***************** |
|---|
| 86 | #Convection scheme switch (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4) |
|---|
| 87 | iflag_con=3 |
|---|
| 88 | #output level of energy conservation diagnostics |
|---|
| 89 | if_ebil=0 |
|---|
| 90 | #maximum efficiency of cld water->precipitation conversion (D: 0.993) |
|---|
| 91 | epmax=0.999 |
|---|
| 92 | #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0) |
|---|
| 93 | ok_adj_ema=n |
|---|
| 94 | iflag_clw=0 |
|---|
| 95 | #Convective closure (D:1, 0:AR4, 2:ALE and ALP) |
|---|
| 96 | iflag_clos=2 |
|---|
| 97 | #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF) |
|---|
| 98 | iflag_mix=1 |
|---|
| 99 | #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0) |
|---|
| 100 | qqa1=1. |
|---|
| 101 | qqa2=0. |
|---|
| 102 | #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05) |
|---|
| 103 | cvl_corr=1.0 |
|---|
| 104 | #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.) |
|---|
| 105 | Fmax=0.65 |
|---|
| 106 | # |
|---|
| 107 | # Flags Couche limite |
|---|
| 108 | #******************** |
|---|
| 109 | #shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) |
|---|
| 110 | iflag_thermals=18 |
|---|
| 111 | #pas de temps du spliting pour les thermiques |
|---|
| 112 | nsplit_thermals=1 |
|---|
| 113 | tau_thermals=0. |
|---|
| 114 | iflag_thermals_ed=8 |
|---|
| 115 | fact_thermals_ed_dz=0.1 |
|---|
| 116 | iflag_thermals_optflux=0 |
|---|
| 117 | #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) |
|---|
| 118 | iflag_pbl=11 |
|---|
| 119 | #diffusion turbulente |
|---|
| 120 | ksta_ter=1.e-7 |
|---|
| 121 | ksta=1.e-10 |
|---|
| 122 | #ok_kzmin : calcul Kzmin dans la CL de surface |
|---|
| 123 | ok_kzmin=y |
|---|
| 124 | #couplage ou non avec la convection |
|---|
| 125 | iflag_coupl=5 (0:AR4 et 1:nouvelle physique) |
|---|
| 126 | seuil_inversion=-0.08 |
|---|
| 127 | #rendement sur la puissance fournie par les thermiques a la convection |
|---|
| 128 | alp_bl_k=0.5 |
|---|
| 129 | # |
|---|
| 130 | ## temperature dependency of qsat for iterative procedure |
|---|
| 131 | ## temporary flag to switch between the old (CMIP5, D=0) and new versions (=1) versions |
|---|
| 132 | iflag_fisrtilp_qsat=4 |
|---|
| 133 | ## Reactivation du calcul d'un zmax continu pour les thermiques |
|---|
| 134 | ## reactivation d'une fermeture humide des thermiques iflag_thermals_closure=2 (D=1) |
|---|
| 135 | iflag_thermals_closure=1 |
|---|
| 136 | # |
|---|
| 137 | # Flags wakes |
|---|
| 138 | #************ |
|---|
| 139 | #Wake scheme switch (D:0, 0:AR4, 1:New Physics) |
|---|
| 140 | iflag_wake=1 |
|---|
| 141 | #Constant value to be added to Alp (D: 0.) |
|---|
| 142 | alp_offset=0. |
|---|
| 143 | |
|---|
| 144 | ### nrlmd le 10/04/2012: Declenchement ### |
|---|
| 145 | #Déclelenchement par couche-limite (0: deterministe, 1: stochastique) |
|---|
| 146 | iflag_trig_bl=1 |
|---|
| 147 | #Surface critique au LCL |
|---|
| 148 | s_trig=1.2e7 |
|---|
| 149 | #Temps de decorrelation entre scenes de cumulus |
|---|
| 150 | tau_trig_shallow=1200 |
|---|
| 151 | #Temps de decorrelation entre scenes de cumulonimbus |
|---|
| 152 | tau_trig_deep=1200 |
|---|
| 153 | # Pour la stabilisation |
|---|
| 154 | random_notrig_max=0.99 |
|---|
| 155 | ok_adjwk=y |
|---|
| 156 | tmax_fonte_cv=278.15 |
|---|
| 157 | |
|---|
| 158 | |
|---|
| 159 | ### nrlmd le 10/04/2012: Fermeture ### |
|---|
| 160 | #Fermeture par couche-limite (0: deterministe, 1: statistique, 2: statistique+convergence grande échelle) |
|---|
| 161 | iflag_clos_bl=1 |
|---|
| 162 | ### fin nrlmd le 10/04/2012 |
|---|
| 163 | |
|---|
| 164 | ## Modification de la plage de temperature pour la phase mixte liquide/glace des nuages |
|---|
| 165 | ## Controle par les t_glace_min/max, exposant_glace, iflag_t_glace (D=0) |
|---|
| 166 | t_glace_min=233.13 |
|---|
| 167 | t_glace_max=273.13 |
|---|
| 168 | exposant_glace=1. |
|---|
| 169 | iflag_t_glace=2 |
|---|
| 170 | iflag_ice_thermo=1 |
|---|
| 171 | |
|---|
| 172 | # |
|---|
| 173 | ## The source terms of the TKE prognostic equation are diagnosed from tendencies (du, dv, dT) |
|---|
| 174 | ## associated with subrgrid scaleamotions and treated as an additional heat source. |
|---|
| 175 | ## Controled by a key, iflag_ener_conserv (D=-1): |
|---|
| 176 | ## 0 : no conservation, -1 : old adhoc correction for kinetic E only (used for CMIP5) |
|---|
| 177 | ## 1 : conservation, 101 : conversion from kinetic to heat only, 110 conversion from potential to heat only |
|---|
| 178 | iflag_ener_conserv=1 |
|---|
| 179 | |
|---|
| 180 | #Activation of corrections for water conservation (D: n). |
|---|
| 181 | ##The default is n for compatibility with earlier versions. |
|---|
| 182 | ##generally it should be y. |
|---|
| 183 | ok_conserv_q=y |
|---|
| 184 | |
|---|
| 185 | iflag_cloudth_vert=0 |
|---|
| 186 | |
|---|
| 187 | # |
|---|
| 188 | lmixmin=10. |
|---|
| 189 | |
|---|
| 190 | #Bergeron |
|---|
| 191 | iflag_bergeron=1 |
|---|
| 192 | ### Convection |
|---|
| 193 | #altitude, in hPa, above LCL at which buoybase is computed, def = original: -40 |
|---|
| 194 | dpbase=-40. |
|---|
| 195 | #altitude, in hPa, above LCL below which Ep=0, def = original: 150, used only if flag_epKEorig=1 |
|---|
| 196 | pbcrit=150.0 |
|---|
| 197 | #altitude, in hPa, above LCL above which Ep=Epmax, def = original: 500, used only if flag_epKEorig=1 |
|---|
| 198 | ptcrit=500.0 |
|---|
| 199 | #reference fractional area of precipitating downdraughts, def = original: 0.01 |
|---|
| 200 | sigdz=0.01 |
|---|
| 201 | #unsaturated fraction of precipitating downdraughts, def = original: 0.15 |
|---|
| 202 | spfac=0.15 |
|---|
| 203 | #relaxation time of cloud base mass flux, def = original: 8000 |
|---|
| 204 | tau=8000. |
|---|
| 205 | #flag for wb (= vert velocity at LFC); 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, def = 1 |
|---|
| 206 | flag_wb=1 |
|---|
| 207 | #parameter used in the computation of wb, def = 6. |
|---|
| 208 | wbmax=6. |
|---|
| 209 | # |
|---|
| 210 | ok_convstop=.False. |
|---|
| 211 | # |
|---|
| 212 | tau_stop=15000. |
|---|
| 213 | # |
|---|
| 214 | ok_intermittent=.False. |
|---|
| 215 | #shedding coefficient (used when iflag_mix_adiab=1) |
|---|
| 216 | coef_peel=0.25 |
|---|
| 217 | # |
|---|
| 218 | ### Clouds |
|---|
| 219 | #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) if 1, def=1 |
|---|
| 220 | flag_epKEorig=1 |
|---|
| 221 | #max condensed water concentration in kg/kg, def=0.0003, original: 0.0011, used only if flag_epKEorig=1 |
|---|
| 222 | elcrit=0.00025 |
|---|
| 223 | #temperature in degre Celsius, at which Ep=Epmax, def=-55., original: -55, used only if flag_epKEorig=1 |
|---|
| 224 | tlcrit=-40. |
|---|
| 225 | |
|---|
