source: LMDZ6/trunk/DefLists/physiq.def_NPv5.80b @ 4432

Last change on this file since 4432 was 3474, checked in by Laurent Fairhead, 6 years ago

Adding proper hybrid coordinates coefficients as per the CMIP6 data request
r3284 of IPSLCM6.0.15 branch
Also added Cecill license information to a bunch of DefLists? files

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