source: CMIP6/GeoMIP/CM61-LR-G6sulfur/PARAM/physiq.def_NPv6.0.13 @ 5128

Last change on this file since 5128 was 5128, checked in by tlurton, 4 years ago

Added CMIP6 experiment G6sulfur.

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