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