# # $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 ########################################################### # 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 ########################################################### # 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 nm_oro_t=0 zstd_orodr_t=0 zpmm_orodr_t=0 zpmm_orolf_t=0 #################################################################### # 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=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.00085 #cld_lc_con=0.00085 # pour refroidir un peu la simu valeur de 1850 à 1975 #cld_lc_lsc=0.00095 #cld_lc_con=0.00095 # pour refroidir un peu la simu à partir de 1975 environ (2025/10/24) #cld_lc_lsc=0.00105 #cld_lc_con=0.00105 # pour refroidir un peu la simu à partir de 2106 environ (2025/11/06) cld_lc_lsc=0.00120 cld_lc_con=0.00120 #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.00027144 # 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 #min cloud droplet number concentration cdnc_min=10 ###################################################################### # 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=n ok_bug_ajs_cv=n # correct convection bug keepbug_ice_frac=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 # Controls detrainement / entrainement / nuages thermals_afact=1.11999558474861 thermals_fact_epsilon=0.00246546968311137 thermals_betalpha=0.395730326872031 thermals_detr_q_coef=0.0113315217193354 cloudth_sigma1s_factor=1.34224269825002 cloudth_sigma2s_factor=0.140971937884714 # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer # le detrainement (A. Jam) fact_thermals_ed_dz=0.115742 # 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 ########################################################## # Poches froides ########################################################## #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 # Bug concerning the gust front length computation ok_bug_gfl=n # New specification of wake vertical profile iflag_wk_profile=1 wk_pupper=3.5 wk_frac_int_delta_t=0.75 wk_delta_t_min=0 iflag_wk_new_ptop=2 ok_bug_zg_wk_pbl=n # correct a bug in thermal splitting ok_bug_split_th=n #################################################################### # 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 # conservation de l eau sur les calottes ok_lic_cond=y #################################################################### # Ozone #################################################################### # new treatment for ozone into the tropopause adjust_tropopause=y