source: LMDZ6/branches/LMDZ_ECRad/DefLists/physiq.def_LUDOBest2 @ 4482

Last change on this file since 4482 was 4482, checked in by lguez, 21 months ago

Sync latest trunk changes to branch LMDZ_ECRad

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