source: CONFIG/UNIFORM/v6/IPSLCM6.5/GENERAL/PARAM/physiq.def_NPv6.0.10 @ 5066

Last change on this file since 5066 was 5066, checked in by cetlod, 4 years ago

First step towards IPSLCM6.5 with the use of NEMOv4 model

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