source: CONFIG/UNIFORM/v6/LMDZOR_v6/GENERAL/PARAM/physiq.def_NPv5.80bz0 @ 3722

Last change on this file since 3722 was 3034, checked in by jgipsl, 7 years ago

Updated LMDZOR_v6 with IPSLCM6:

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