source: CMIP6/GeoMIP/CM61-LR-G6sulfur/PARAM/physiq.def_NPv6.0.7 @ 5128

Last change on this file since 5128 was 5128, checked in by tlurton, 4 years ago

Added CMIP6 experiment G6sulfur.

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