source: CMIP6/ScenarioMIP/CM61-LR-scen-ssp245-r6/PARAM/physiq.def_NPv6.0.9 @ 4784

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

Added experiment ssp245-r6.

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