source: LMDZ6/trunk/DefLists/physiq.def_NPv7.0.1b @ 5374

Last change on this file since 5374 was 4574, checked in by musat, 18 months ago

NPv7.0.1b : 1er reglage CMIP7 b(efore tuning) [ou b(eta)]
de la physique LUDO corrigée avec

  • activation lscp
  • modifs sur les poches
  • activation des drags/lifts partout (seuils à 0).

Ionela Musat

  • Property svn:executable set to *
File size: 20.3 KB
Line 
1#
2# $Id: physiq.def_NPv6.2_MR1 3867 2021-03-23 16:49:48Z fairhead $
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          #  Modelisation des bosquets
31ifl_pbltree=1
32Cd_frein=0.016
33
34###########################################################
35# parameter des flux air-mer
36###########################################################
37
38          # Facteur sur les flux sensibles et latent
39          # Introduit pour reduire l evaporation quand on avait z0h=z0m
40          # CONSEILLE : 1
41f_cdrag_oce=1
42
43          # Flag controlant le calcul des z0
44iflag_z0_oce=1
45
46          # Rapport entre les Z0 latent et sensible
47f_z0qh_oce=1.55
48
49          # Activation des rafales : Oui 1 ou non 0
50iflag_gusts=2
51
52          # facteur pour les rafales des thermiques (bounary layer)
53f_gust_bl=0.845
54
55          # facteur pour les poches (wake)
56f_gust_wk=0.2
57
58          # facteur sur qsat pour tenir compte de la salinite
59f_qsat_oce=0.98
60
61          #cdmmax : valeur max de Cd moment. Sur ocean seulement
62cdmmax=0.2
63
64          #cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement
65cdhmax=0.2
66
67####################################################################
68# Couche limite turbulente
69####################################################################
70
71          # Nouvelle version de Mellor et Yamada
72new_yamada4=y
73
74          # Choix du schema numerique pour new_yamada4=y
75          #  1 schema de MAR. Bien pour CL stables mais detruit les stratocu
76          #  5 schema de MAR modifie. Precaunise.
77yamada4_num=5
78
79          # Flag de controle de la couche limite stable
80iflag_corr_sta=4
81
82          # min sur les fonction de stabilite a la surface
83f_ri_cd_min=0.01
84
85          # max du Ric pour Kz. Decouplage plus grand pour Ric plus grand.
86yamada4_ric=0.18
87
88          # Longueur de melange minimum pour Kz
89lmixmin=0
90
91          #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)
92iflag_pbl=12
93
94          # Seuils pour la diffusion turbulente
95ksta_ter=1e-07
96ksta=1e-10
97
98          #ok_kzmin : calcul Kzmin dans la CL de surface (D: y)
99ok_kzmin=n
100
101          # Creation de TKE par orographie sous-maille
102addtkeoro=2
103alphatkeoro=1
104smallscales_tkeoro=y
105
106
107####################################################################
108# Ondes de gravite / Orographie sous maille
109####################################################################
110
111          #  Orodr  ou  non   pour l orographie
112ok_orodr=y
113
114          #  Orolf  ou  non   pour l orographie             
115ok_orolf=y
116
117          #Rugoro : inclusion de l orographie sous maille dans le Cd de surface
118          # Obsolete avec l activation des ondes de gravite
119f_rugoro=0
120
121          # Parametres controlant l effet de l orographie sous maille
122sso_gfrcri=0.7
123sso_gkdrag=0.6
124sso_gkwake=0.4
125sso_gklift=0.1
126
127          # Random gravity waves
128          #Activitate the stochastic parameterization of convective GWs
129          # (should be equal to ok_strato)
130ok_gwd_rando=y
131
132          # Scale the momentum fluxes due to convective GWs
133gwd_rando_ruwmax=0.5
134
135          #Scale the saturated momentum flux of  convective GWs
136gwd_rando_sat=0.6
137
138          # Important pour controler la periode de la QBO
139gwd_front_ruwmax=3
140
141          # zstd_orodr_t : minimum subgrid orography standard deviation for activation of drag parameterisation (D:10, 0.:CMIP7)
142zstd_orodr_t = 0.
143
144          # zpmm_orodr_t : minimum zpic-zval value for activation of drag parameterisation (D:100., 0.:CMIP7)
145zpmm_orodr_t = 0.
146
147          # zpmm_orolf_t : minimum zpic-zval value for activation of lift parameterisation (D:100., 0.:CMIP7)
148zpmm_orolf_t =0.
149
150###############################################################
151# Rayonnement
152###############################################################
153          # activation du nouveau code de rayonnement RRTM
154          # 0 : Ancien code, 1 : RRTM et 2 : ECRAD (D=0)
155iflag_rrtm=1
156
157          # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0
158NSW=6
159
160    # Albedo et inertie des calottes
161alb_vis_sno_lic=0.96
162alb_nir_sno_lic=0.68
163    # Controle de l'inertie des glaciers et banquises
164    # inertie_sno : pour la neige (y compris sur les points "ter"
165    #      quand orchidee n'est pas actif)
166    # inertie_lic/sic : pour les "land ice" et "sea ice"
167inertie_sno=350
168    # iflag_sic=1 pour forcer a garder l'inertie des "sea ice"
169    # a inertie_sic meme si la surface est couverte de neige
170iflag_sic=1
171inertie_lic=2000
172inertie_sic=2000
173    # tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours)
174    # Utilise selon le forcage : tau_gl=5 pour Clim, tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2
175tau_gl=2
176
177          # ???
178iflag_albedo=1
179
180          # Frequence appel convection. Nombre appels par jour
181          # (nbapp_cv=0 ==> appel tous les pas de temps)
182nbapp_cv=48
183
184          # Pour pouvoir revenir a un bug purement informatique de duplication
185          # des constantes thermodynamiques entre LMDZ et RRTM
186          # n : FORTEMENT CONSEILLE
187ok_bad_ecmwf_thermo=n
188
189###############################################################
190# Parametres nuages
191###############################################################
192
193          # Activation condensation de grande echelle epuree par E.Vignon : ok_new_lscp=y + ok_icefra_lscp=y
194ok_new_lscp=y
195ok_icefra_lscp=y
196
197
198          #seuils de la precipitation des nuages stratiformes (D: 2.6e-4 2.6e-4)
199cld_lc_lsc=0.00096841923606032
200cld_lc_con=0.00096841923606032
201
202          #constante de temps pour eliminer eau lsc et convective
203          # (D: 3600. 3600.)
204cld_tau_lsc=8151.30533480977
205cld_tau_con=8151.30533480977
206
207          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
208ffallv_lsc=0.451743342889991
209ffallv_con=0.451743342889991
210
211          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
212          # regle a 3.e-5 sur des cas de cumulus en 1D
213coef_eva=0.00490135133370158
214          # coefficient sur la sublimation de la glace (D: 2.e-5 n)
215coef_eva_i=0.00490135133370158
216
217### iflag_autoconversion= 0 avec reste de bogue sur l autoconversion (default)
218### iflag_autoconversion= 1, sans bogue sur l autoconversion
219### iflag_autoconversion= 2 : 1 + new numerical integration of autoconversion
220###                         : have to add niter_lscp=1
221iflag_autoconversion=1
222
223          # flag pour la reevaporation de la pluie
224          # 0: rien
225          # 1: a l ancienne
226          # 2: prend en compte la surface nuageuse max au dessus pour calculer
227          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE
228          # 4: Ludo original
229          # 5: Ludo + correction Etienne
230iflag_evap_prec=5
231
232          # min local intensity for rainfall before decrease of the associated precipitation fraction (mm/s)
233          # rain_int_min is used if iflag_evap_prec>=4
234rain_int_min=1.44716149682657e-5
235
236          # Modification de la plage de temperature pour la phase mixte
237          # liquide/glace des nuages
238          # Controle par les t_glace_min/max, exposant_glace,
239          # iflag_t_glace (D=0)
240t_glace_min=243.15
241t_glace_max=273.15
242exposant_glace=0.5
243iflag_t_glace=3
244
245          # Activation de la thermodynamique de la glace
246iflag_ice_thermo=1
247
248          # normalement elle ne sert plus a rien dans fisrtilp
249          # sinon elle est un peu toxique:
250          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste)
251reevap_ice=n
252
253          #Bergeron
254iflag_bergeron=2
255
256          # temperature dependency of qsat for iterative procedure
257          # temporary flag to switch between the old (CMIP5, D=0)
258          # and new versions (=1) versions
259iflag_fisrtilp_qsat=4
260
261
262          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
263iflag_cld_th=6
264fact_cldcon=1
265
266          # ?????
267facttemps=0
268
269          #calcul eau condensee et fraction nuageuse a partir des PDF
270          # (D:0, 0:version avec ratqs sinon nouvelles PDFs) 
271iflag_pdf=1
272
273          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
274ok_newmicro=y
275
276          # Controle des ratqs grande echelle
277iflag_ratqs=4
278ratqsp0=48254.6771101272
279ratqsdp=16383.2736877474
280ratqsbas=0.002
281ratqshaut=0.342148444633161
282
283          #rayon effectif pour les nuages de glace (D:35)
284rad_froid=35
285
286          # Rayons effectif des gouttes nuageuses
287          # ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns)
288rad_chau1=12
289rad_chau2=11
290
291          # Rayon effectif des gouttes si ok_cdnc=0
292          # Parameters in CDNC-maer link (Boucher&Lohmann 1995)
293bl95_b0=1.3
294bl95_b1=0.2
295
296          #  Choix ou non  de  New oliq               
297new_oliq=y
298
299          # Seuils oliqmax, oicemax (en g/kg).
300          # Au dela on precipite directement au sol
301oliqmax=0.0015
302oicemax=0.0015
303
304          #Rayons min et max effectifs des particules des nuages
305rei_min=18.1205
306rei_max=69.4129
307
308          # Distribution vetricale sous maille pour le schema de de nuages.
309          # 0: rien
310          # 1: Premiere version A. Jam,
311          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5
312          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32
313          #                           cloudth_vert_alpha_th C=2:91
314iflag_cloudth_vert=3
315cloudth_vert_alpha=0.5
316cloudth_vert_alpha_th=0.5
317
318
319######################################################################
320# Convection
321######################################################################
322
323          #Convection scheme switch
324          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
325iflag_con=3
326
327          # correct convection bug
328          #(default = y)
329keepbug_ice_frac=n
330
331          #output level of energy conservation diagnostics
332if_ebil=0
333
334          #maximum efficiency of cld water->precipitation conversion (D: 0.993)
335epmax=0.999071
336
337          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0)
338ok_adj_ema=n
339iflag_clw=0
340
341          #Convective closure (D:1, 0:AR4, 2:ALE and ALP)
342iflag_clos=2
343
344          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
345iflag_mix=1
346
347          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0)
348qqa1=1
349qqa2=0
350
351          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
352cvl_corr=1
353
354          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
355Fmax=0.65
356
357          #altitude, in hPa, above LCL at which buoybase is computed,
358          # def = original: -40
359dpbase=-40
360
361          #reference fractional area of precipitating downdraughts,
362          # def = original: 0.01
363sigdz=0.003
364
365          #unsaturated fraction of precipitating downdraughts, def = original: 0.15
366spfac=0.15
367
368          #relaxation time of cloud base mass flux, def = original: 8000
369tau=8000
370
371          #flag for wb (= vert velocity at LFC);
372          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1
373          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1
374          # wbmax : assymptotic value
375flag_wb=50
376wbmax=2.8
377
378          # Divers tests
379ok_convstop=n
380tau_stop=15000
381ok_intermittent=n
382
383          # nouvel epluchage de C.Rio
384          #shedding coefficient (used when iflag_mix_adiab=1)
385iflag_mix_adiab=0
386coef_peel=0.25
387
388          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
389          #if 1, def=1
390flag_epKEorig=1
391
392          # Case : flag_epKEorig=1
393          # elcrit : max condensed water concentration in kg/kg, D=0.0003,
394          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
395          # original: 0.0011, used only if flag_epKEorig=1
396elcrit=0.0003
397tlcrit=-55
398
399          # Case : flag_epKEorig=0
400          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original)
401          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original)
402pbcrit=150
403ptcrit=500
404
405          # Precipitation ejection from adiabatic ascents
406          # cvflag_prec_eject: when TRUE ejection is activated. Default: False
407          # WARNING ! : cvflag_prec_eject=y is effective only if icvflag_Tpa=1
408          # ejectliq : 0. =>liquid precipitation is carried up in the adiab updraughts
409          #            1. =>liquid precipitation is ejected from adiab updraughts
410          # ejectice   0. =>solid precipitation is carried up in the adiab updraughts
411          #            1. =>solid precipitation is ejected from adiab updraughts
412          # ejectliq may take only 0. and 1. values
413          # ejectice may take any value between 0. and 1.
414          #   Default : cvflag_prec_eject=n, ejectliq=0., ejectice=0.
415cvflag_prec_eject=n
416ejectliq=0.
417ejectice=0.
418
419          # Check whether parcel level temperature and specific humidity
420          # are reasonable and demand that LCL be high enough
421cv_flag_feed=2
422
423          # save 20% CPU
424ok_optim_yield=y
425
426          # Convective clouds
427iflag_cld_cv=0
428
429          # To get back to bug where convective output variables
430          # were set to 0 at the start of the physics
431ok_bug_cv_trac=n
432
433           # solving of the temperature equation for the adiabatic ascent
434           # 0 -> AJ scheme in two loops ; ice fraction given by envt. temp.
435           # 1 -> JYG scheme in one loop ; ice fraction given by envt. temp.
436           # 2 -> JYG scheme in one loop ; ice fraction given by adiab. asc. temp.
437           # Default: 0
438           # 2 is costly and ill tested; experts only.
439icvflag_Tpa=0
440
441           # taking into account the dependence of qv on qt at saturation
442           # Default: n
443qsat_depends_on_qt=n
444
445##########################################################
446# Fermeture convective et couplages
447##########################################################
448
449          # Decimal code:  + first digit drives the splitting of thermals:
450          #                                0 = thermals see the average profiles of T and q
451          #                                1 = thermals are only outside wakes
452          #                + second digit drives the splitting of diffuse PBL:
453          #                                0 = The diffuse PBL sees the average profiles of T and q
454          #                                1 = Two BLs, one in (w) and one in (x), with the same
455          #                                    surface temperature and humidity.
456          #                                2 = Two BLs, one in (w) and one in (x), with separate
457          #                                    surface temperature and humidity.
458          #                                3 = same as 1 over ocean; same as 0 everywhere else.
459          # Example: iflag_pbl_split=12 => thermals outside wakes, split BL with Tsurf and qsurf difference
460          # iflag_pbl_split=10 => thermals outside wakes
461iflag_pbl_split=10
462
463          #rendement sur la puissance fournie par les thermiques a la convection
464alp_bl_k=0.5
465
466          #Constant value to be added to Alp (D: 0.)
467alp_offset=0
468
469          #Convection must reach a level at which T<t_top_max to be active
470t_top_max=1000
471
472          ### nrlmd le 10/04/2012: Declenchement ###
473          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
474iflag_trig_bl=1
475
476          #Surface critique au LCL
477s_trig=12000000
478
479          #Temps de decorrelation entre scenes de cumulus
480tau_trig_shallow=1200
481
482          #Temps de decorrelation entre scenes de cumulonimbus
483tau_trig_deep=1200
484
485          # Pour la stabilisation
486random_notrig_max=0.99
487
488          # Temperature a la quelle on fond les descentes glacees
489tmax_fonte_cv=278.15
490
491          #Fermeture par couche-limite (Nicolas Rochetin)
492          #(0: deterministe, 1: statistique, 2: statistique+convergence
493          # grande echelle)
494iflag_clos_bl=1
495
496          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
497          # def: 0.33
498stark=0.33
499
500          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
501          # by wakes, def: 0.25
502alpk=0.25
503
504#
505##########################################################
506# Couche limite convective / modele du thermique
507##########################################################
508
509          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
510iflag_thermals=18
511
512          # pas de temps du spliting pour les thermiques
513          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
514nsplit_thermals=1
515
516          # tau_thermals pour avoir une constante de temps sur les thermiques.
517          # Non valide
518tau_thermals=0
519
520          # Flag controlant le detrainement et l entrainement
521iflag_thermals_ed=8
522
523          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
524          # le detrainement (A. Jam)
525fact_thermals_ed_dz=0.103105610491144
526
527          # Une option de plus sur les thermiques.
528iflag_thermals_optflux=0
529
530          #couplage ou non avec la convection
531iflag_coupl=5
532
533          # Seuil pour desactiver les thermiques en cas de forte inversion
534          # en sommet de couche limite. Actif si iflag_thermals impair.
535seuil_inversion=-0.08
536
537          # Reactivation du calcul d un zmax continu pour les thermiques
538          # reactivation d une fermeture humide des thermiques
539          # iflag_thermals_closure=2 (D=1)
540iflag_thermals_closure=2
541
542          # Flag permettant de garder le bug induisant des wake_deltaT non nuls
543          # a partir des tendances des thermiques, meme en l'absence de poches.
544          # Default : ok_bug_split_th=y ==> le bug est maintenu
545ok_bug_split_th=n
546
547          # contribution of buoyancy to the plume acceleration
548thermals_afact=1.05641487083682
549
550          # drag term in the plume acceleration
551thermals_fact_epsilon=0.00205282352527457
552
553          # scaling factor for entrainment and detrainment
554thermals_betalpha=0.411960019760327
555
556          # width of the environment subgrid scale water distribution
557cloudth_sigma1s_factor=0.673674753185742
558
559##########################################################
560# Flags wakes
561##########################################################
562
563          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
564iflag_wake=1
565
566          #multiplicative factor of the damping by gravity waves, def: 4.
567coefgw=4
568
569          #wake density = number of wake centers per m2, def: 8.E-12
570          #wdens_ref=8.E-12
571wdens_ref_o=1e-09
572wdens_ref_l=8e-12
573
574          # Ajustement convectif prealable au calcul des poches
575ok_adjwk=y
576
577          # Prevent some crashes
578          # Filter out bad wakes
579flag_wk_check_trgl=n
580iflag_wk_check_trgl=2
581iflag_alp_wk_cond=1
582
583          #Wake population dynamic switch (D:0):
584          # iflag_wk_pop_dyn=0 ==> wake_dens is set to wdens_ref
585          # iflag_wk_pop_dyn=1 ==> wake_dens is a prognostic variable
586iflag_wk_pop_dyn=0
587
588iflag_wk_profile=1
589wk_pupper=0.6
590####################################################################
591#  Conservation
592####################################################################
593          # The source terms of the TKE prognostic equation are diagnosed
594          # from tendencies (du, dv, dT)
595          # associated with subrgrid scaleamotions and treated as an
596          # additional heat source.
597          # Controled by a key, iflag_ener_conserv (D=-1):
598          # 0 : no conservation, -1 : old adhoc correction for kinetic E only
599          # (used for CMIP5)
600          # 1 : conservation, 101 : conversion from kinetic to heat only,
601          # 110 conversion from potential to heat only
602iflag_ener_conserv=10
603
604          #Activation of corrections for water conservation (D: n).
605          ##The default is n for compatibility with earlier versions.
606          ##generally it should be y.
607ok_conserv_q=y
608
609fl_cor_ebil=1
610
611####################################################################
612# Ozone
613####################################################################
614          # new treatment for ozone into the tropopause
615adjust_tropopause=y
616
617####################################################################
618#  Conservation de l'eau
619#####################################################################
620          # conservation de l eau sur les calottes
621ok_lic_cond=y
622####################################################################
623          #min cloud droplet number concentration
624cdnc_min=10.
Note: See TracBrowser for help on using the repository browser.