source: CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/physiq.def_NPv6.1_strataer @ 4672

Last change on this file since 4672 was 4672, checked in by tlurton, 2 years ago

Added new STRATAER historical experiment.

File size: 17.7 KB
Line 
1#
2# $Id: physiq.def_NPv6.1 3964 2018-07-06 11:44:14Z acosce $
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=1.75
134#--added OB according to FL
135mode_top_bound=3
136
137          #Scale the saturated momentum flux of  convective GWs
138gwd_rando_sat=0.6
139
140          # Important pour controler la periode de la QBO
141gwd_front_ruwmax=3
142
143
144###############################################################
145# Rayonnement
146###############################################################
147          # activation du nouveau code de rayonnement RRTM
148          # 0 : Ancien code et 1 : RRTM (D=0)
149iflag_rrtm=1
150
151          # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0
152NSW=6
153
154    # Albedo et inertie des calottes
155alb_vis_sno_lic=0.96
156alb_nir_sno_lic=0.68
157    # Controle de l'inertie des glaciers et banquises
158    # inertie_sno : pour la neige (y compris sur les points "ter"
159    #      quand orchidee n'est pas actif utilisé)
160    # inertie_lic/sic : pour les "land ice" et "sea ice"
161inertie_sno=350
162    # iflag_sic=1 pour forcer à garder l'inertie des "sea ice"
163    # a inertie_sic meme si la surface est couverte de neige
164iflag_sic=1
165inertie_lic=2000
166inertie_sic=2000
167
168
169          # ???
170iflag_albedo=1
171
172          # Frequence appel rayonnement. Nombre appels par jour
173nbapp_cv=48
174
175          # Pour pouvoir revenir a un bug purement informatique de duplication
176          # des constantes thermodynamiques entre LMDZ et RRTM
177          # n : FORTEMENT CONSEILLE
178ok_bad_ecmwf_thermo=n
179
180###############################################################
181# Parametres nuages
182###############################################################
183
184          #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
185cld_lc_lsc=0.00065
186cld_lc_con=0.00065
187
188          #constante de temps pour eleminer eau lsc et convective
189          # (D: 3600. 3600.)
190cld_tau_lsc=900
191cld_tau_con=900
192
193          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
194ffallv_lsc=0.8
195ffallv_con=0.8
196
197          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
198          # regle a 3.e-5 sur des cas de cumulus en 1D
199coef_eva=0.0001
200
201          # flag pour la reevaporation de la pluie
202          # 0: rien
203          # 1: a l ancienne
204          # 2: prend en compte la surface nuageuse max au dessus pour calculer
205          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE
206iflag_evap_prec=2
207
208          # Modification de la plage de temperature pour la phase mixte
209          # liquide/glace des nuages
210          # Controle par les t_glace_min/max, exposant_glace,
211          # iflag_t_glace (D=0)
212t_glace_min=243.15
213t_glace_max=273.15
214exposant_glace=0.5
215iflag_t_glace=3
216
217          # Activation de la thermodynamique de la glace
218iflag_ice_thermo=1
219
220          # normalement elle ne sert plus a rien dans fisrtilp
221          # sinon elle est un peu toxique:
222          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste)
223reevap_ice=n
224
225          #Bergeron
226iflag_bergeron=2
227
228          # temperature dependency of qsat for iterative procedure
229          # temporary flag to switch between the old (CMIP5, D=0)
230          # and new versions (=1) versions
231iflag_fisrtilp_qsat=4
232
233          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
234iflag_cldcon=6
235fact_cldcon=1
236
237          # ?????
238facttemps=0
239
240          #calcul eau condensee et fraction nuageuse a partir des PDF
241          # (D:0, 0:version avec ratqs sinon nouvelles PDFs) 
242iflag_pdf=1
243
244          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
245ok_newmicro=y
246
247          # Controle des ratqs grande echelle
248iflag_ratqs=4
249ratqsp0=45000
250ratqsdp=10000
251ratqsbas=0.002
252ratqshaut=0.4
253
254          #rayon effectif pour les nuages de glace (D:35)
255rad_froid=35
256
257          # Rayons effectif des gouttes nuageuses
258          # ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns)
259rad_chau1=12
260rad_chau2=11
261
262          # Rayon effectif des gouttes si ok_cdnc=0
263          # Parameters in CDNC-maer link (Boucher&Lohmann 1995)
264bl95_b0=1.3
265bl95_b1=0.2
266
267          #  Choix ou non  de  New oliq               
268new_oliq=y
269
270          # Seuils oliqmax, oicemax (en g/kg).
271          # Au dela on precipite directement au sol
272oliqmax=0.0015
273oicemax=0.0015
274
275          #Rayons min et max effectifs des particules des nuages
276rei_min=16
277rei_max=61.29
278
279          # Distribution vetricale sous maille pour le schema de de nuages.
280          # 0: rien
281          # 1: Premiere version A. Jam,
282          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5
283          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32
284          #                           cloudth_vert_alpha_th C=2:91
285iflag_cloudth_vert=3
286cloudth_vert_alpha=0.5
287cloudth_vert_alpha_th=0.5
288
289
290######################################################################
291# Convection
292######################################################################
293
294          #Convection scheme switch
295          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
296iflag_con=3
297
298          #output level of energy conservation diagnostics
299if_ebil=0
300
301          #maximum efficiency of cld water->precipitation conversion (D: 0.993)
302epmax=0.999
303
304          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0)
305ok_adj_ema=n
306iflag_clw=0
307
308          #Convective closure (D:1, 0:AR4, 2:ALE and ALP)
309iflag_clos=2
310
311          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
312iflag_mix=1
313
314          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0)
315qqa1=1
316qqa2=0
317
318          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
319cvl_corr=1
320
321          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
322Fmax=0.65
323
324          #altitude, in hPa, above LCL at which buoybase is computed,
325          # def = original: -40
326dpbase=-40
327
328          #reference fractional area of precipitating downdraughts,
329          # def = original: 0.01
330sigdz=0.003
331
332          #unsaturated fraction of precipitating downdraughts, def = original: 0.15
333spfac=0.15
334
335          #relaxation time of cloud base mass flux, def = original: 8000
336tau=8000
337
338          #flag for wb (= vert velocity at LFC);
339          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1
340          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1
341          # wbmax : assymptotic value
342flag_wb=50
343wbmax=2.8
344
345          # Divers tests
346ok_convstop=n
347tau_stop=15000
348ok_intermittent=n
349
350          # nouvel epluchage de C.Rio
351          #shedding coefficient (used when iflag_mix_adiab=1)
352iflag_mix_adiab=0
353coef_peel=0.25
354
355          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
356          #if 1, def=1
357flag_epKEorig=1
358
359          # Case : flag_epKEorig=1
360          # elcrit : max condensed water concentration in kg/kg, D=0.0003,
361          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
362          # original: 0.0011, used only if flag_epKEorig=1
363elcrit=0.0003
364tlcrit=-55
365
366          # Case : flag_epKEorig=0
367          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original)
368          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original)
369pbcrit=150
370ptcrit=500
371
372          # Prevent some crashes
373          # Filter out bad wakes
374flag_wk_check_trgl=n
375iflag_wk_check_trgl=2
376iflag_alp_wk_cond=1
377
378
379          # Check whether parcel level temperature and specific humidity
380          # are reasonable and demand that LCL be high enough
381cv_flag_feed=2
382
383          # save 20% CPU
384ok_optim_yield=y
385
386          # Convective clouds
387iflag_cld_cv=0
388
389          # To get back to bug where convective output variables
390          # were set to 0 at the start of the physics
391ok_bug_cv_trac=_AUTO_: DEFAULT = n
392
393##########################################################
394# Fermeture convective et couplages
395##########################################################
396
397          # Appel separe des thermiques
398iflag_pbl_split=10
399
400          #rendement sur la puissance fournie par les thermiques a la convection
401alp_bl_k=0.5
402
403          #Constant value to be added to Alp (D: 0.)
404alp_offset=0
405
406          #Convection must reach a level at which T<t_top_max to be active
407t_top_max=1000
408
409          ### nrlmd le 10/04/2012: Declenchement ###
410          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
411iflag_trig_bl=1
412
413          #Surface critique au LCL
414s_trig=12000000
415
416          #Temps de decorrelation entre scenes de cumulus
417tau_trig_shallow=1200
418
419          #Temps de decorrelation entre scenes de cumulonimbus
420tau_trig_deep=1200
421
422          # Pour la stabilisation
423random_notrig_max=0.99
424
425          # Temperature a la quelle on fond les descentes glacees
426tmax_fonte_cv=278.15
427
428          #Fermeture par couche-limite (Nicolas Rochetin)
429          #(0: deterministe, 1: statistique, 2: statistique+convergence
430          # grande echelle)
431iflag_clos_bl=1
432
433          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
434          # def: 0.33
435stark=0.33
436
437          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
438          # by wakes, def: 0.25
439alpk=0.25
440
441#
442##########################################################
443# Couche limite convective / modele du thermique
444##########################################################
445
446          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
447iflag_thermals=18
448
449          # pas de temps du spliting pour les thermiques
450          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
451nsplit_thermals=1
452
453          # tau_thermals pour avoir une constante de temps sur les thermiques.
454          # Non valide
455tau_thermals=0
456
457          # Flag controlant le detrainement et l entrainement
458iflag_thermals_ed=8
459
460          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
461          # le detrainement (A. Jam)
462fact_thermals_ed_dz=0.07
463
464          # Une option de plus sur les thermiques.
465iflag_thermals_optflux=0
466
467          #couplage ou non avec la convection
468iflag_coupl=5
469
470          # Seuil pour desactiver les thermiques en cas de forte inversion
471          # en sommet de couche limite. Actif si iflag_thermals impair.
472seuil_inversion=-0.08
473
474          # Reactivation du calcul d un zmax continu pour les thermiques
475          # reactivation d une fermeture humide des thermiques
476          # iflag_thermals_closure=2 (D=1)
477iflag_thermals_closure=2
478
479
480##########################################################
481# Flags wakes
482##########################################################
483
484          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
485iflag_wake=1
486
487          #multiplicative factor of the damping by gravity waves, def: 4.
488coefgw=4
489
490          #wake density = number of wake centers per m2, def: 8.E-12
491          #wdens_ref=8.E-12
492wdens_ref_o=1e-09
493wdens_ref_l=8e-12
494
495          # Ajustement convectif prealable au calcul des poches
496ok_adjwk=y
497
498
499####################################################################
500#  Conservation
501####################################################################
502
503          # The source terms of the TKE prognostic equation are diagnosed
504          # from tendencies (du, dv, dT)
505          # associated with subrgrid scaleamotions and treated as an
506          # additional heat source.
507          # Controled by a key, iflag_ener_conserv (D=-1):
508          # 0 : no conservation, -1 : old adhoc correction for kinetic E only
509          # (used for CMIP5)
510          # 1 : conservation, 101 : conversion from kinetic to heat only,
511          # 110 conversion from potential to heat only
512iflag_ener_conserv=10
513
514          #Activation of corrections for water conservation (D: n).
515          ##The default is n for compatibility with earlier versions.
516          ##generally it should be y.
517ok_conserv_q=y
518
519fl_cor_ebil=1
520
521####################################################################
522# Ozone
523####################################################################
524          # new treatment for ozone into the tropopause
525adjust_tropopause=y
526
527####################################################################
528#  Conservation de l'eau
529#####################################################################
530          # conservation de l eau sur les calottes
531ok_lic_cond=y
532
533cdnc_min=10.
534
535####################################################################
536###  StratAer
537######################################################################
538flag_aerosol_strat=2
539convscav=y
540lessivage=y
541iflag_lscav=4
542ccntrAA=1
543ccntrENV=1.
544coefcoli=0.001
545
546####################################################################
547##  Injection
548####################################################################
549# Nucleation rate in lat/pres box
550flag_nuc_rate_box=n         # Nucleation rate limit or not to a lat/pres
551nuclat_min=-15.0            # min lat to activate nuc rate (degN)
552nuclat_max=15.0             # max lat to activate nuc rate (degN)
553nucpres_min=8000.0          # min pres to activate nuc rate (Pa)
554nucpres_max=11000.0         # max pres to activate nuc rate (Pa)
555# Aerosol radiative feedback (0: off, 1: on)
556flag_aer_feedback = y
557# Coag traceur
558flag_sulf_emit = 1
559# Duree de l'injection (en jours)
560injdur = 1
561# Flag either Gaussian distribution (0) or uniform distribution (1)
562flag_sulf_emit_distrib = 0
563#--flag_sulf_emit=1 --example Pinatubo
564nErupt=1                     # eruption nb during simulation
565year_emit_vol=1991          # year of emission date
566mth_emit_vol=6               # month of emission date
567day_emit_vol=15             # day of emission date
568m_aer_emiss_vol=5.0e9        # emitted sulfur mass in kgS, e.g. 7Tg(S)=14Tg(SO2)
569altemiss_vol=19.e3           # emission altitude in m
570sigma_alt_vol=1.0e3         # standard deviation of emission altitude in m
571xlat_min_vol=15.0         # min latitude of volcano in degree
572xlat_max_vol=15.0      # max latitude of volcano in degree
573xlon_min_vol=120.0   # min longitude of volcano in degree
574xlon_max_vol=120.0   # max longitude of volcano in degree
575#--flag_sulf_emit=2 --SAI
576m_aer_emiss_sai=1.e10  # emitted sulfur mass in kgS, eg 1e9=1TgS, 1e10=10TgS (kg/yr in case flag_sulf_emit=3)
577altemiss_sai=19.e3     # emission altitude in m
578sigma_alt_sai=1.e3     # standard deviation of emission altitude in m
579xlat_sai=0.0          # latitude of SAI in degree
580xlon_sai=0.01         # longitude of SAI in degree
581#--flag_sulf_emit=3 -- SAI between two latitudes
582xlat_min_sai=-10.0    # maximum latitude of SAI in degrees
583xlat_max_sai=10.0     # maximum latitude of SAI in degrees
Note: See TracBrowser for help on using the repository browser.