source: CMIP6/ScenarioMIP/CM61-LR-scen-ssp126/PARAM/physiq.def_NPv6.0.11 @ 4220

Last change on this file since 4220 was 4220, checked in by oboucher, 5 years ago

ScenarioMIP ssp126

File size: 14.2 KB
Line 
1#
2# $Id: physiq.def_NPv6.0.11 3312 2017-07-12 09:43:27Z aclsce $
3#
4
5##############################################################
6# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ
7##############################################################
8# D: Valeur par default
9#
10
11###########################################################
12# divers. Devrait aller dans config.def ?
13###########################################################
14##   Si=.T. ,  lecture du fichier limit avec la bonne annee
15ok_limitvrai=n
16
17###########################################################
18# Parametres couche limite de surface
19###########################################################
20
21          #  ...
22f_cdrag_stable=1.
23
24          # Facteur sur les coefficients de trainee terrestres
25f_cdrag_ter=1.
26
27          # Vitesse minimum dans les calculs de trainee Cd ( || V || + Vmin )
28min_wind_speed=0.1
29
30###########################################################
31# parameter des flux air-mer
32###########################################################
33
34          # Facteur sur les flux sensibles et latent
35          # Introduit pour reduire l evaporation quand on avait z0h=z0m
36          # CONSEILLE : 1
37f_cdrag_oce=1.
38
39          # Flag controlant le calcul des z0
40iflag_z0_oce=1
41
42          # Rapport entre les Z0 latent et sensible
43f_z0qh_oce=1.55
44
45          # Activation des rafales : Oui 1 ou non 0
46iflag_gusts=0
47
48          # facteur pour les rafales des thermiques (bounary layer)
49f_gust_bl=0.845
50
51          # facteur pour les poches (wake)
52f_gust_wk=0.2
53
54          # facteur sur qsat pour tenir compte de la salinite
55f_qsat_oce=1.
56
57          #cdmmax : valeur max de Cd moment. Sur ocean seulement
58cdmmax=2.5E-3
59
60          #cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement
61cdhmax=2.0E-3
62
63####################################################################
64# Couche limite turbulente
65####################################################################
66
67          # Nouvelle version de Mellor et Yamada
68new_yamada4=y
69
70          # Choix du schema numerique pour new_yamada4=y
71          #  1 schema de MAR. Bien pour CL stables mais detruit les stratocu
72          #  5 schema de MAR modifie. Precaunise.
73yamada4_num=5
74
75          # Flag de controle de la couche limite stable
76iflag_corr_sta=4
77
78          # min sur les fonction de stabilite a la surface
79f_ri_cd_min=0.01
80
81          # max du Ric pour Kz. Decouplage plus grand pour Ric plus grand.
82yamada4_ric=0.18
83
84          # Longueur de melange minimum pour Kz
85lmixmin=0.
86
87          #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)
88iflag_pbl=12
89
90          # Seuils pour la diffusion turbulente
91ksta_ter=1.e-7
92ksta=1.e-10
93
94          #ok_kzmin : calcul Kzmin dans la CL de surface (D: y)
95ok_kzmin=n
96
97          # Creation de TKE par orographie sous-maille
98addtkeoro=0
99alphatkeoro=1.
100smallscales_tkeoro=n
101
102
103####################################################################
104# Ondes de gravite / Orographie sous maille
105####################################################################
106
107          #  Orodr  ou  non   pour l orographie
108ok_orodr=y
109
110          #  Orolf  ou  non   pour l orographie             
111ok_orolf=y
112
113          #Rugoro : inclusion de l orographie sous maille dans le Cd de surface
114          # Obsolete avec l activation des ondes de gravite
115f_rugoro=0.
116
117          # Parametres controlant l effet de l orographie sous maille
118sso_gfrcri=0.70
119sso_gkdrag=0.6
120sso_gkwake=0.4
121sso_gklift=0.1
122
123          # Random gravity waves
124          #Activitate the stochastic parameterization of convective GWs
125          # (should be equal to ok_strato)
126ok_gwd_rando=y
127
128          # Scale the momentum fluxes due to convective GWs
129gwd_rando_ruwmax=0.23
130
131          #Scale the saturated momentum flux of  convective GWs
132gwd_rando_sat=0.60
133
134          # Important pour controler la periode de la QBO
135gwd_front_ruwmax = 3.00
136
137
138###############################################################
139# Rayonnement
140###############################################################
141          # activation du nouveau code de rayonnement RRTM
142          # 0 : Ancien code et 1 : RRTM (D=0)
143iflag_rrtm=1
144
145          # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0
146NSW=6
147
148    # Albedo et inertie des calottes
149alb_vis_sno_lic=0.96
150alb_nir_sno_lic=0.68
151inertie_sno=350.
152
153          # ???
154iflag_albedo=1
155
156          # Frequence appel rayonnement. Nombre appels par jour
157nbapp_cv=48
158
159          # Pour pouvoir revenir a un bug purement informatique de duplication
160          # des constantes thermodynamiques entre LMDZ et RRTM
161          # n : FORTEMENT CONSEILLE
162ok_bad_ecmwf_thermo=n
163
164###############################################################
165# Parametres nuages
166###############################################################
167
168          #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
169cld_lc_lsc=0.000075
170cld_lc_con=0.000075
171
172          #constante de temps pour eleminer eau lsc et convective
173          # (D: 3600. 3600.)
174cld_tau_lsc=900.
175cld_tau_con=900.
176
177          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
178ffallv_lsc=0.6
179ffallv_con=0.6
180
181          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
182          # regle a 3.e-5 sur des cas de cumulus en 1D
183coef_eva=1.e-4
184
185          # flag pour la reevaporation de la pluie
186          # 0: rien
187          # 1: a l ancienne
188          # 2: prend en compte la surface nuageuse max au dessus pour calculer
189          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE
190iflag_evap_prec=1
191
192          # Modification de la plage de temperature pour la phase mixte
193          # liquide/glace des nuages
194          # Controle par les t_glace_min/max, exposant_glace,
195          # iflag_t_glace (D=0)
196t_glace_min=243.15
197t_glace_max=273.15
198exposant_glace=2.
199iflag_t_glace=2
200
201          # Activation de la thermodynamique de la glace
202iflag_ice_thermo=1
203
204          # normalement elle ne sert plus a rien dans fisrtilp
205          # sinon elle est un peu toxique:
206          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste)
207reevap_ice=y
208
209          #Bergeron
210iflag_bergeron=2
211
212          # temperature dependency of qsat for iterative procedure
213          # temporary flag to switch between the old (CMIP5, D=0)
214          # and new versions (=1) versions
215iflag_fisrtilp_qsat=4
216
217          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
218iflag_cldcon=6
219fact_cldcon=1.
220
221          # ?????
222facttemps=0.
223
224          #calcul eau condensee et fraction nuageuse a partir des PDF
225          # (D:0, 0:version avec ratqs sinon nouvelles PDFs) 
226iflag_pdf=1
227
228          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
229ok_newmicro=y
230
231          # Controle des ratqs grande echelle
232iflag_ratqs=4
233ratqsp0=45000.
234ratqsdp=20000.
235ratqsbas=0.002
236ratqshaut=0.4
237
238          #rayon effectif pour les nuages de glace (D:35)
239rad_froid=35
240
241          #rayons effectifs pour les nuages eau liq (tailles des gouttes
242          # d eau nuageuse, D: 13 9)
243rad_chau1=12
244rad_chau2=11
245
246          #  Choix ou non  de  New oliq               
247new_oliq=y
248
249          # Seuils oliqmax, oicemax (en g/kg).
250          # Au dela on precipite directement au sol
251oliqmax=1.5E-3
252oicemax=1.5E-3
253
254          #Rayons min et max effectifs des particules des nuages
255rei_min=16.
256rei_max=61.29
257
258          # Distribution vetricale sous maille pour le schema de de nuages.
259          # 0: rien
260          # 1: Premiere version A. Jam,
261          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5
262          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32
263          #                           cloudth_vert_alpha_th C=2:91
264iflag_cloudth_vert=3
265cloudth_vert_alpha=0.5
266cloudth_vert_alpha_th=0.5
267
268
269######################################################################
270# Convection
271######################################################################
272
273          #Convection scheme switch
274          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
275iflag_con=3
276
277          #output level of energy conservation diagnostics
278if_ebil=0
279
280          #maximum efficiency of cld water->precipitation conversion (D: 0.993)
281epmax=0.9985
282
283          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0)
284ok_adj_ema=n
285iflag_clw=0
286
287          #Convective closure (D:1, 0:AR4, 2:ALE and ALP)
288iflag_clos=2
289
290          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
291iflag_mix=1
292
293          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0)
294qqa1=1.
295qqa2=0.
296
297          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
298cvl_corr=1.0
299
300          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
301Fmax=0.65
302          #
303t_top_max=255.
304#
305### Convection
306          #altitude, in hPa, above LCL at which buoybase is computed,
307          # def = original: -40
308dpbase=-40.
309
310          #reference fractional area of precipitating downdraughts,
311          # def = original: 0.01
312sigdz=0.01
313
314          #unsaturated fraction of precipitating downdraughts, def = original: 0.15
315spfac=0.15
316
317          #relaxation time of cloud base mass flux, def = original: 8000
318tau=8000.
319
320          #flag for wb (= vert velocity at LFC);
321          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1
322          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1
323          # wbmax : assymptotic value
324flag_wb=1
325wbmax=6.
326
327          # Divers tests
328ok_convstop=.False.
329tau_stop=15000.
330ok_intermittent=.False.
331
332          # nouvel epluchage de C.Rio
333          #shedding coefficient (used when iflag_mix_adiab=1)
334iflag_mix_adiab=0
335coef_peel=0.25
336
337          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
338          #if 1, def=1
339flag_epKEorig=1
340
341          # Case : flag_epKEorig=1
342          # elcrit : max condensed water concentration in kg/kg, D=0.0003,
343          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
344          # original: 0.0011, used only if flag_epKEorig=1
345elcrit=0.00025
346tlcrit=-40.
347
348          # Case : flag_epKEorig=0
349          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original)
350          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original)
351pbcrit=150.0
352ptcrit=500.0
353
354          # Prevent some crashes
355          # Filter out bad wakes
356flag_wk_check_trgl=y
357
358          # Check whether parcel level temperature and specific humidity
359          # are reasonable and demand that LCL be high enough
360cv_flag_feed=2
361
362          # save 20% CPU
363ok_optim_yield=y
364
365
366##########################################################
367# Fermeture convective et couplages
368##########################################################
369
370          # Appel separe des thermiques
371iflag_pbl_split=0
372
373          #rendement sur la puissance fournie par les thermiques a la convection
374alp_bl_k=0.5
375
376          #Constant value to be added to Alp (D: 0.)
377alp_offset=0.
378
379          ### nrlmd le 10/04/2012: Declenchement ###
380          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
381iflag_trig_bl=1
382
383          #Surface critique au LCL
384s_trig=1.2e7
385
386          #Temps de decorrelation entre scenes de cumulus
387tau_trig_shallow=1200
388
389          #Temps de decorrelation entre scenes de cumulonimbus
390tau_trig_deep=1200
391
392          # Pour la stabilisation
393random_notrig_max=0.99
394
395          # Temperature a la quelle on fond les descentes glacees
396tmax_fonte_cv=278.15
397
398          #Fermeture par couche-limite (Nicolas Rochetin)
399          #(0: deterministe, 1: statistique, 2: statistique+convergence
400          # grande echelle)
401iflag_clos_bl=1
402
403          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
404          # def: 0.33
405stark=0.33
406
407          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
408          # by wakes, def: 0.25
409alpk=0.25
410
411#
412##########################################################
413# Couche limite convective / modele du thermique
414##########################################################
415
416          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
417iflag_thermals=18
418
419          # pas de temps du spliting pour les thermiques
420          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
421nsplit_thermals=1
422
423          # tau_thermals pour avoir une constante de temps sur les thermiques.
424          # Non valide
425tau_thermals=0.
426
427          # Flag controlant le detrainement et l entrainement
428iflag_thermals_ed=8
429
430          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
431          # le detrainement (A. Jam)
432fact_thermals_ed_dz=0.07
433
434          # Une option de plus sur les thermiques.
435iflag_thermals_optflux=0
436
437          #couplage ou non avec la convection
438iflag_coupl=5 (0:AR4 et 1:nouvelle physique)
439
440          # Seuil pour desactiver les thermiques en cas de forte inversion
441          # en sommet de couche limite. Actif si iflag_thermals impair.
442seuil_inversion=-0.08
443
444          # Reactivation du calcul d un zmax continu pour les thermiques
445          # reactivation d une fermeture humide des thermiques
446          # iflag_thermals_closure=2 (D=1)
447iflag_thermals_closure=2
448
449
450##########################################################
451# Flags wakes
452##########################################################
453
454          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
455iflag_wake=1
456
457          #multiplicative factor of the damping by gravity waves, def: 4.
458coefgw=4.
459
460          #wake density = number of wake centers per m2, def: 8.E-12
461          #wdens_ref=8.E-12
462wdens_ref_o=1.e-9
463wdens_ref_l=8.e-12
464
465          # Ajustement convectif prealable au calcul des poches
466ok_adjwk=y
467
468
469####################################################################
470#  Conservation
471####################################################################
472
473          # The source terms of the TKE prognostic equation are diagnosed
474          # from tendencies (du, dv, dT)
475          # associated with subrgrid scaleamotions and treated as an
476          # additional heat source.
477          # Controled by a key, iflag_ener_conserv (D=-1):
478          # 0 : no conservation, -1 : old adhoc correction for kinetic E only
479          # (used for CMIP5)
480          # 1 : conservation, 101 : conversion from kinetic to heat only,
481          # 110 conversion from potential to heat only
482iflag_ener_conserv=1
483
484          #Activation of corrections for water conservation (D: n).
485          ##The default is n for compatibility with earlier versions.
486          ##generally it should be y.
487ok_conserv_q=y
488
489fl_cor_ebil=1
Note: See TracBrowser for help on using the repository browser.