source: LMDZ6/trunk/DefLists/physiq.def_NPv6.0.14splith @ 3399

Last change on this file since 3399 was 3162, checked in by Laurent Fairhead, 7 years ago

Official physiq.def files for IPSLCM6.0.14 configuration

  • Property svn:keywords set to Id
File size: 15.3 KB
Line 
1#
2# $Id: physiq.def_NPv6.0.14splith 3162 2018-01-26 07:22:57Z 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 utilisé)
158    # inertie_lic/sic : pour les "land ice" et "sea ice"
159inertie_sno=350
160    # iflag_sic=1 pour forcer à 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
166
167          # ???
168iflag_albedo=1
169
170          # Frequence appel rayonnement. Nombre appels par jour
171nbapp_cv=48
172
173          # Pour pouvoir revenir a un bug purement informatique de duplication
174          # des constantes thermodynamiques entre LMDZ et RRTM
175          # n : FORTEMENT CONSEILLE
176ok_bad_ecmwf_thermo=n
177
178###############################################################
179# Parametres nuages
180###############################################################
181
182
183          #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
184cld_lc_lsc=0.0004
185cld_lc_con=0.0004
186
187          #constante de temps pour eleminer eau lsc et convective
188          # (D: 3600. 3600.)
189cld_tau_lsc=900
190cld_tau_con=900
191
192          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
193ffallv_lsc=0.8
194ffallv_con=0.8
195
196          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
197          # regle a 3.e-5 sur des cas de cumulus en 1D
198coef_eva=0.0001
199
200          # flag pour la reevaporation de la pluie
201          # 0: rien
202          # 1: a l ancienne
203          # 2: prend en compte la surface nuageuse max au dessus pour calculer
204          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE
205iflag_evap_prec=2
206
207          # Modification de la plage de temperature pour la phase mixte
208          # liquide/glace des nuages
209          # Controle par les t_glace_min/max, exposant_glace,
210          # iflag_t_glace (D=0)
211t_glace_min=243.15
212t_glace_max=273.15
213exposant_glace=0.5
214iflag_t_glace=3
215
216          # Activation de la thermodynamique de la glace
217iflag_ice_thermo=1
218
219          # normalement elle ne sert plus a rien dans fisrtilp
220          # sinon elle est un peu toxique:
221          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste)
222reevap_ice=n
223
224          #Bergeron
225iflag_bergeron=2
226
227          # temperature dependency of qsat for iterative procedure
228          # temporary flag to switch between the old (CMIP5, D=0)
229          # and new versions (=1) versions
230iflag_fisrtilp_qsat=4
231
232          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
233iflag_cldcon=6
234fact_cldcon=1
235
236          # ?????
237facttemps=0
238
239          #calcul eau condensee et fraction nuageuse a partir des PDF
240          # (D:0, 0:version avec ratqs sinon nouvelles PDFs) 
241iflag_pdf=1
242
243          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
244ok_newmicro=y
245
246          # Controle des ratqs grande echelle
247iflag_ratqs=4
248ratqsp0=45000
249ratqsdp=10000
250ratqsbas=0.002
251ratqshaut=0.4
252
253          #rayon effectif pour les nuages de glace (D:35)
254rad_froid=35
255
256          # Rayons effectif des gouttes nuageuses
257          # ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns)
258rad_chau1=12
259rad_chau2=11
260
261          # Rayon effectif des gouttes si ok_cdnc=0
262          # Parameters in CDNC-maer link (Boucher&Lohmann 1995)
263bl95_b0=1.3
264bl95_b1=0.2
265
266          #  Choix ou non  de  New oliq               
267new_oliq=y
268
269          # Seuils oliqmax, oicemax (en g/kg).
270          # Au dela on precipite directement au sol
271oliqmax=0.0015
272oicemax=0.0015
273
274          #Rayons min et max effectifs des particules des nuages
275rei_min=16
276rei_max=61.29
277
278          # Distribution vetricale sous maille pour le schema de de nuages.
279          # 0: rien
280          # 1: Premiere version A. Jam,
281          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5
282          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32
283          #                           cloudth_vert_alpha_th C=2:91
284iflag_cloudth_vert=3
285cloudth_vert_alpha=0.5
286cloudth_vert_alpha_th=0.5
287
288
289######################################################################
290# Convection
291######################################################################
292
293          #Convection scheme switch
294          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
295iflag_con=3
296
297          #output level of energy conservation diagnostics
298if_ebil=0
299
300          #maximum efficiency of cld water->precipitation conversion (D: 0.993)
301epmax=0.999
302
303          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0)
304ok_adj_ema=n
305iflag_clw=0
306
307          #Convective closure (D:1, 0:AR4, 2:ALE and ALP)
308iflag_clos=2
309
310          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
311iflag_mix=1
312
313          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0)
314qqa1=1
315qqa2=0
316
317          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
318cvl_corr=1
319
320          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
321Fmax=0.65
322
323          #altitude, in hPa, above LCL at which buoybase is computed,
324          # def = original: -40
325dpbase=-40
326
327          #reference fractional area of precipitating downdraughts,
328          # def = original: 0.01
329sigdz=0.003
330
331          #unsaturated fraction of precipitating downdraughts, def = original: 0.15
332spfac=0.15
333
334          #relaxation time of cloud base mass flux, def = original: 8000
335tau=8000
336
337          #flag for wb (= vert velocity at LFC);
338          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1
339          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1
340          # wbmax : assymptotic value
341flag_wb=50
342wbmax=2.8
343
344          # Divers tests
345ok_convstop=n
346tau_stop=15000
347ok_intermittent=n
348
349          # nouvel epluchage de C.Rio
350          #shedding coefficient (used when iflag_mix_adiab=1)
351iflag_mix_adiab=0
352coef_peel=0.25
353
354          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
355          #if 1, def=1
356flag_epKEorig=1
357
358          # Case : flag_epKEorig=1
359          # elcrit : max condensed water concentration in kg/kg, D=0.0003,
360          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
361          # original: 0.0011, used only if flag_epKEorig=1
362elcrit=0.0003
363tlcrit=-55
364
365          # Case : flag_epKEorig=0
366          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original)
367          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original)
368pbcrit=150
369ptcrit=500
370
371          # Prevent some crashes
372          # Filter out bad wakes
373flag_wk_check_trgl=n
374iflag_wk_check_trgl=2
375iflag_alp_wk_cond=1
376
377
378          # Check whether parcel level temperature and specific humidity
379          # are reasonable and demand that LCL be high enough
380cv_flag_feed=2
381
382          # save 20% CPU
383ok_optim_yield=y
384
385          # Convective clouds
386iflag_cld_cv=0
387
388##########################################################
389# Fermeture convective et couplages
390##########################################################
391
392          # Appel separe des thermiques
393iflag_pbl_split=10
394
395          #rendement sur la puissance fournie par les thermiques a la convection
396alp_bl_k=0.5
397
398          #Constant value to be added to Alp (D: 0.)
399alp_offset=0
400
401          #Convection must reach a level at which T<t_top_max to be active
402t_top_max=1000
403
404          ### nrlmd le 10/04/2012: Declenchement ###
405          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
406iflag_trig_bl=1
407
408          #Surface critique au LCL
409s_trig=12000000
410
411          #Temps de decorrelation entre scenes de cumulus
412tau_trig_shallow=1200
413
414          #Temps de decorrelation entre scenes de cumulonimbus
415tau_trig_deep=1200
416
417          # Pour la stabilisation
418random_notrig_max=0.99
419
420          # Temperature a la quelle on fond les descentes glacees
421tmax_fonte_cv=278.15
422
423          #Fermeture par couche-limite (Nicolas Rochetin)
424          #(0: deterministe, 1: statistique, 2: statistique+convergence
425          # grande echelle)
426iflag_clos_bl=1
427
428          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
429          # def: 0.33
430stark=0.33
431
432          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
433          # by wakes, def: 0.25
434alpk=0.25
435
436#
437##########################################################
438# Couche limite convective / modele du thermique
439##########################################################
440
441          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique)
442iflag_thermals=18
443
444          # pas de temps du spliting pour les thermiques
445          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
446nsplit_thermals=1
447
448          # tau_thermals pour avoir une constante de temps sur les thermiques.
449          # Non valide
450tau_thermals=0
451
452          # Flag controlant le detrainement et l entrainement
453iflag_thermals_ed=8
454
455          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
456          # le detrainement (A. Jam)
457fact_thermals_ed_dz=0.07
458
459          # Une option de plus sur les thermiques.
460iflag_thermals_optflux=0
461
462          #couplage ou non avec la convection
463iflag_coupl=5
464
465          # Seuil pour desactiver les thermiques en cas de forte inversion
466          # en sommet de couche limite. Actif si iflag_thermals impair.
467seuil_inversion=-0.08
468
469          # Reactivation du calcul d un zmax continu pour les thermiques
470          # reactivation d une fermeture humide des thermiques
471          # iflag_thermals_closure=2 (D=1)
472iflag_thermals_closure=2
473
474
475##########################################################
476# Flags wakes
477##########################################################
478
479          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
480iflag_wake=1
481
482          #multiplicative factor of the damping by gravity waves, def: 4.
483coefgw=4
484
485          #wake density = number of wake centers per m2, def: 8.E-12
486          #wdens_ref=8.E-12
487wdens_ref_o=1e-09
488wdens_ref_l=8e-12
489
490          # Ajustement convectif prealable au calcul des poches
491ok_adjwk=y
492
493
494####################################################################
495#  Conservation
496####################################################################
497
498          # The source terms of the TKE prognostic equation are diagnosed
499          # from tendencies (du, dv, dT)
500          # associated with subrgrid scaleamotions and treated as an
501          # additional heat source.
502          # Controled by a key, iflag_ener_conserv (D=-1):
503          # 0 : no conservation, -1 : old adhoc correction for kinetic E only
504          # (used for CMIP5)
505          # 1 : conservation, 101 : conversion from kinetic to heat only,
506          # 110 conversion from potential to heat only
507iflag_ener_conserv=10
508
509          #Activation of corrections for water conservation (D: n).
510          ##The default is n for compatibility with earlier versions.
511          ##generally it should be y.
512ok_conserv_q=y
513
514fl_cor_ebil=1
515
516####################################################################
517# Ozone
518####################################################################
519          # new treatment for ozone into the tropopause
520adjust_tropopause=y
521
522####################################################################
523#  Conservation de l'eau
524#####################################################################
525          # conservation de l eau sur les calottes
526ok_lic_cond=y
527
Note: See TracBrowser for help on using the repository browser.