Index: LMDZ6/trunk/DefLists/physiq.def_NPv6.5 =================================================================== --- LMDZ6/trunk/DefLists/physiq.def_NPv6.5 (revision 5780) +++ LMDZ6/trunk/DefLists/physiq.def_NPv6.5 (revision 5780) @@ -0,0 +1,556 @@ +# +# $Id: physiq.def_NPv6.1 3964 2018-07-06 11:44:14Z acosce $ +# + +############################################################## +# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ +############################################################## +# D: Valeur par default +# + +########################################################### +# divers. Devrait aller dans config.def ? +########################################################### +## Si=.T. , lecture du fichier limit avec la bonne annee +ok_limitvrai=n + +########################################################### +# Parametres couche limite de surface +########################################################### + + # ... +f_cdrag_stable=1 + + # Facteur sur les coefficients de trainee terrestres +f_cdrag_ter=1 + + # Vitesse minimum dans les calculs de trainee Cd ( || V || + Vmin ) +min_wind_speed=0.1 + + # Modelisation des bosquets +ifl_pbltree=1 +Cd_frein=0.016 + +########################################################### +# parameter des flux air-mer +########################################################### + + # Facteur sur les flux sensibles et latent + # Introduit pour reduire l evaporation quand on avait z0h=z0m + # CONSEILLE : 1 +f_cdrag_oce=1 + + # Flag controlant le calcul des z0 +iflag_z0_oce=1 + + # Rapport entre les Z0 latent et sensible +f_z0qh_oce=1.55 + + # Activation des rafales : Oui 1 ou non 0 +iflag_gusts=2 + + # facteur pour les rafales des thermiques (bounary layer) +f_gust_bl=0.845 + + # facteur pour les poches (wake) +f_gust_wk=0.2 + + # facteur sur qsat pour tenir compte de la salinite +f_qsat_oce=0.98 + + #cdmmax : valeur max de Cd moment. Sur ocean seulement +cdmmax=0.2 + + #cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement +cdhmax=0.2 + +#################################################################### +# Couche limite turbulente +#################################################################### + + # Nouvelle version de Mellor et Yamada +new_yamada4=y + + # Choix du schema numerique pour new_yamada4=y + # 1 schema de MAR. Bien pour CL stables mais detruit les stratocu + # 5 schema de MAR modifie. Precaunise. +yamada4_num=5 + + # Flag de controle de la couche limite stable +iflag_corr_sta=4 + + # min sur les fonction de stabilite a la surface +f_ri_cd_min=0.01 + + # max du Ric pour Kz. Decouplage plus grand pour Ric plus grand. +yamada4_ric=0.18 + + # Longueur de melange minimum pour Kz +lmixmin=0 + + #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) +iflag_pbl=12 + + # Seuils pour la diffusion turbulente +ksta_ter=1e-07 +ksta=1e-10 + + #ok_kzmin : calcul Kzmin dans la CL de surface (D: y) +ok_kzmin=n + + # Creation de TKE par orographie sous-maille +addtkeoro=2 +alphatkeoro=1 +smallscales_tkeoro=y + + +#################################################################### +# Ondes de gravite / Orographie sous maille +#################################################################### + + # Orodr ou non pour l orographie +ok_orodr=y + + # Orolf ou non pour l orographie +ok_orolf=y + + #Rugoro : inclusion de l orographie sous maille dans le Cd de surface + # Obsolete avec l activation des ondes de gravite +f_rugoro=0 + + # Parametres controlant l effet de l orographie sous maille +sso_gfrcri=0.7 +sso_gkdrag=0.6 +sso_gkwake=0.4 +sso_gklift=0.1 + + # Random gravity waves + #Activitate the stochastic parameterization of convective GWs + # (should be equal to ok_strato) +ok_gwd_rando=y + + # Scale the momentum fluxes due to convective GWs +gwd_rando_ruwmax=0.5 + + #Scale the saturated momentum flux of convective GWs +gwd_rando_sat=0.6 + + # Important pour controler la periode de la QBO +gwd_front_ruwmax=3 + + +############################################################### +# Rayonnement +############################################################### + # activation du nouveau code de rayonnement RRTM + # 0 : Ancien code et 1 : RRTM (D=0) +iflag_rrtm=1 + + # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0 +NSW=6 + + # Albedo et inertie des calottes +alb_vis_sno_lic=0.96 +alb_nir_sno_lic=0.68 + # Controle de l'inertie des glaciers et banquises + # inertie_sno : pour la neige (y compris sur les points "ter" + # quand orchidee n'est pas actif utilisé) + # inertie_lic/sic : pour les "land ice" et "sea ice" +inertie_sno=350 + # iflag_sic=1 pour forcer à garder l'inertie des "sea ice" + # a inertie_sic meme si la surface est couverte de neige +iflag_sic=1 +inertie_lic=2000 +inertie_sic=2000 + # tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours) + # # Utilise celon le forcage : tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2 +tau_gl= _AUTO_ : DEFAULT = 5. + + # ??? +iflag_albedo=1 + + # Frequence appel rayonnement. Nombre appels par jour +nbapp_cv=48 + + # Pour pouvoir revenir a un bug purement informatique de duplication + # des constantes thermodynamiques entre LMDZ et RRTM + # n : FORTEMENT CONSEILLE +ok_bad_ecmwf_thermo=n + +############################################################### +# Parametres nuages +############################################################### + + #seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4) +cld_lc_lsc=0.001 +cld_lc_con=0.001 + + #constante de temps pour eleminer eau lsc et convective + # (D: 3600. 3600.) +cld_tau_lsc=900 +cld_tau_con=900 + + #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1) +ffallv_lsc=0.8 +ffallv_con=0.8 + + # coefficient sur la reevaporation de la pluie (D: 2.e-5 n) + # regle a 3.e-5 sur des cas de cumulus en 1D +coef_eva=0.0001 + + # flag pour la reevaporation de la pluie + # 0: rien + # 1: a l ancienne + # 2: prend en compte la surface nuageuse max au dessus pour calculer + # le max de reevaporation, comme ecart a la saturation. CONSEILLE +iflag_evap_prec=2 + + # Modification de la plage de temperature pour la phase mixte + # liquide/glace des nuages + # Controle par les t_glace_min/max, exposant_glace, + # iflag_t_glace (D=0) +t_glace_min=243.15 +t_glace_max=273.15 +exposant_glace=0.5 +iflag_t_glace=3 + + # Activation de la thermodynamique de la glace +iflag_ice_thermo=1 + + # normalement elle ne sert plus a rien dans fisrtilp + # sinon elle est un peu toxique: + # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste) +reevap_ice=n + + #Bergeron +iflag_bergeron=2 + + # temperature dependency of qsat for iterative procedure + # temporary flag to switch between the old (CMIP5, D=0) + # and new versions (=1) versions +iflag_fisrtilp_qsat=4 + + #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) +iflag_cld_th=6 +fact_cldcon=1 + + # ????? +facttemps=0 + + #calcul eau condensee et fraction nuageuse a partir des PDF + # (D:0, 0:version avec ratqs sinon nouvelles PDFs) +iflag_pdf=1 + + #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3) +ok_newmicro=y + + # Controle des ratqs grande echelle +iflag_ratqs=4 +ratqsp0=45000 +ratqsdp=10000 +ratqsbas=0.002 +ratqshaut=0.4 + + #rayon effectif pour les nuages de glace (D:35) +rad_froid=35 + + # Rayons effectif des gouttes nuageuses + # ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns) +rad_chau1=12 +rad_chau2=11 + + # Rayon effectif des gouttes si ok_cdnc=0 + # Parameters in CDNC-maer link (Boucher&Lohmann 1995) +bl95_b0=1.3 +bl95_b1=0.2 + + # Choix ou non de New oliq +new_oliq=y + + # Seuils oliqmax, oicemax (en g/kg). + # Au dela on precipite directement au sol +oliqmax=0.0015 +oicemax=0.0015 + + #Rayons min et max effectifs des particules des nuages +rei_min=16 +rei_max=61.29 + + # Distribution vetricale sous maille pour le schema de de nuages. + # 0: rien + # 1: Premiere version A. Jam, + # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5 + # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32 + # cloudth_vert_alpha_th C=2:91 +iflag_cloudth_vert=3 +cloudth_vert_alpha=0.5 +cloudth_vert_alpha_th=0.5 + + +###################################################################### +# Convection +###################################################################### + + #Convection scheme switch + # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4) +iflag_con=3 + + #output level of energy conservation diagnostics +if_ebil=0 + + #maximum efficiency of cld water->precipitation conversion (D: 0.993) +epmax=0.999 + + #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0) +ok_adj_ema=n +iflag_clw=0 + + #Convective closure (D:1, 0:AR4, 2:ALE and ALP) +iflag_clos=2 + + #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF) +iflag_mix=1 + + #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0) +qqa1=1 +qqa2=0 + + #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05) +cvl_corr=1 + + #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.) +Fmax=0.65 + + #altitude, in hPa, above LCL at which buoybase is computed, + # def = original: -40 +dpbase=-40 + + #reference fractional area of precipitating downdraughts, + # def = original: 0.01 +sigdz=0.003 + + #unsaturated fraction of precipitating downdraughts, def = original: 0.15 +spfac=0.15 + + #relaxation time of cloud base mass flux, def = original: 8000 +tau=8000 + + #flag for wb (= vert velocity at LFC); + # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1 + # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1 + # wbmax : assymptotic value +flag_wb=50 +wbmax=2.8 + + # Divers tests +ok_convstop=n +tau_stop=15000 +ok_intermittent=n + + # nouvel epluchage de C.Rio + #shedding coefficient (used when iflag_mix_adiab=1) +iflag_mix_adiab=0 +coef_peel=0.25 + + #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) + #if 1, def=1 +flag_epKEorig=1 + + # Case : flag_epKEorig=1 + # elcrit : max condensed water concentration in kg/kg, D=0.0003, + # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original) + # original: 0.0011, used only if flag_epKEorig=1 +elcrit=0.0003 +tlcrit=-55 + + # Case : flag_epKEorig=0 + # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original) + # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original) +pbcrit=150 +ptcrit=500 + + # Prevent some crashes + # Filter out bad wakes +flag_wk_check_trgl=n +iflag_wk_check_trgl=2 +iflag_alp_wk_cond=1 + + + # Check whether parcel level temperature and specific humidity + # are reasonable and demand that LCL be high enough +cv_flag_feed=2 + + # save 20% CPU +ok_optim_yield=y + + # Convective clouds +iflag_cld_cv=0 + + # To get back to bug where convective output variables + # were set to 0 at the start of the physics +ok_bug_cv_trac=_AUTO_: DEFAULT = n + +########################################################## +# Fermeture convective et couplages +########################################################## + + # Appel separe des thermiques +iflag_pbl_split=10 + + #rendement sur la puissance fournie par les thermiques a la convection +alp_bl_k=0.5 + + #Constant value to be added to Alp (D: 0.) +alp_offset=0 + + #Convection must reach a level at which TC_star formula: C_star=stark*sqrt(2*Wape), + # def: 0.33 +stark=0.33 + + #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided + # by wakes, def: 0.25 +alpk=0.25 + +# +########################################################## +# Couche limite convective / modele du thermique +########################################################## + + # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique) +iflag_thermals=18 + + # pas de temps du spliting pour les thermiques + # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD +nsplit_thermals=1 + + # tau_thermals pour avoir une constante de temps sur les thermiques. + # Non valide +tau_thermals=0 + + # Flag controlant le detrainement et l entrainement +iflag_thermals_ed=8 + + # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer + # le detrainement (A. Jam) +fact_thermals_ed_dz=0.07 + + # Une option de plus sur les thermiques. +iflag_thermals_optflux=0 + + #couplage ou non avec la convection +iflag_coupl=5 + + # Seuil pour desactiver les thermiques en cas de forte inversion + # en sommet de couche limite. Actif si iflag_thermals impair. +seuil_inversion=-0.08 + + # Reactivation du calcul d un zmax continu pour les thermiques + # reactivation d une fermeture humide des thermiques + # iflag_thermals_closure=2 (D=1) +iflag_thermals_closure=2 + + +########################################################## +# Flags wakes +########################################################## + + #Wake scheme switch (D:0, 0:AR4, 1:New Physics) +iflag_wake=1 + + #multiplicative factor of the damping by gravity waves, def: 4. +coefgw=4 + + #wake density = number of wake centers per m2, def: 8.E-12 + #wdens_ref=8.E-12 +wdens_ref_o=1e-09 +wdens_ref_l=8e-12 + + # Ajustement convectif prealable au calcul des poches +ok_adjwk=y + + +#################################################################### +# Conservation +#################################################################### + # The source terms of the TKE prognostic equation are diagnosed +# from tendencies (du, dv, dT) + # associated with subrgrid scaleamotions and treated as an + # additional heat source. + # Controled by a key, iflag_ener_conserv (D=-1): + # 0 : no conservation, -1 : old adhoc correction for kinetic E only + # (used for CMIP5) + # 1 : conservation, 101 : conversion from kinetic to heat only, + # 110 conversion from potential to heat only +iflag_ener_conserv=10 + + #Activation of corrections for water conservation (D: n). + ##The default is n for compatibility with earlier versions. + ##generally it should be y. +ok_conserv_q=y + +fl_cor_ebil=1 + +#################################################################### +# Ozone +#################################################################### + # new treatment for ozone into the tropopause +adjust_tropopause=y + +#################################################################### +# Conservation de l'eau +##################################################################### + # conservation de l eau sur les calottes +ok_lic_cond=y +#################################################################### + +cdnc_min=10. + #min cloud droplet number concentration +ok_bug_split_th=n + # correct a bug in thermal splitting +keepbug_ice_frac=n + # correct convection bug +################################################# +# A ranger +################################################# +iflag_wk_profile=1 +ok_bug_ajs_cv=n +wk_pupper = 3.5 +wk_frac_int_delta_t=0.75 +iflag_wk_new_ptop=2 +ok_conserv_q=y +ok_bug_gfl=n +ok_bug_zg_wk_pbl=n +stark=0.66 +alpk=0.25 +nm_oro_t=0. +zstd_orodr_t=0. +zpmm_orodr_t=0. +zpmm_orolf_t=0. +