#
# $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 (empty documentation)
###########################################################
# Parametres couche limite de surface
###########################################################
# ...
f_cdrag_stable =1 (empty documentation)
# Facteur sur les coefficients de trainee terrestres
f_cdrag_ter =1 (empty documentation)
# Vitesse minimum dans les calculs de trainee Cd ( || V || + Vmin )
min_wind_speed =0.1 (empty documentation)
# Modelisation des bosquets
ifl_pbltree =1 (empty documentation)
Cd_frein =0.016 (empty documentation)
###########################################################
# 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 (empty documentation)
# Flag controlant le calcul des z0
iflag_z0_oce =1 (empty documentation)
# Rapport entre les Z0 latent et sensible
f_z0qh_oce =1.55 (empty documentation)
# Activation des rafales : Oui 1 ou non 0
iflag_gusts =2 (empty documentation)
# facteur pour les rafales des thermiques (bounary layer)
f_gust_bl =0.845 (empty documentation)
# facteur pour les poches (wake)
f_gust_wk =0.2 (empty documentation)
# facteur sur qsat pour tenir compte de la salinite
f_qsat_oce =0.98 (empty documentation)
#cdmmax : valeur max de Cd moment. Sur ocean seulement
cdmmax =0.2 (empty documentation)
#cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement
cdhmax =0.2 (empty documentation)
####################################################################
# Couche limite turbulente
####################################################################
# Nouvelle version de Mellor et Yamada
new_yamada4 =y (empty documentation)
new_yamada4 =y (empty documentation)
# 1 schema de MAR. Bien pour CL stables mais detruit les stratocu
# 5 schema de MAR modifie. Precaunise.
yamada4_num =5 (empty documentation)
# Flag de controle de la couche limite stable
iflag_corr_sta =4 (empty documentation)
# min sur les fonction de stabilite a la surface
f_ri_cd_min =0.01 (empty documentation)
# max du Ric pour Kz. Decouplage plus grand pour Ric plus grand.
yamada4_ric =0.18 (empty documentation)
# Longueur de melange minimum pour Kz
lmixmin =0 (empty documentation)
#shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada)
iflag_pbl =12 (empty documentation)
# Seuils pour la diffusion turbulente
ksta_ter =1e-07 (empty documentation)
ksta =1e-10 (empty documentation)
#ok_kzmin : calcul Kzmin dans la CL de surface (D: y)
ok_kzmin =n (empty documentation)
# Creation de TKE par orographie sous-maille
addtkeoro =2 (empty documentation)
alphatkeoro =1 (empty documentation)
smallscales_tkeoro =y (empty documentation)
####################################################################
# Ondes de gravite / Orographie sous maille
####################################################################
# Orodr ou non pour l orographie
ok_orodr =y (empty documentation)
# Orolf ou non pour l orographie
ok_orolf =y (empty documentation)
#Rugoro : inclusion de l orographie sous maille dans le Cd de surface
# Obsolete avec l activation des ondes de gravite
f_rugoro =0 (empty documentation)
# Parametres controlant l effet de l orographie sous maille
sso_gfrcri =0.7 (empty documentation)
sso_gkdrag =0.6 (empty documentation)
sso_gkwake =0.4 (empty documentation)
sso_gklift =0.1 (empty documentation)
# Random gravity waves
#Activitate the stochastic parameterization of convective GWs
# (should be equal to ok_strato)
ok_gwd_rando =y (empty documentation)
# Scale the momentum fluxes due to convective GWs
gwd_rando_ruwmax =0.5 (empty documentation)
#Scale the saturated momentum flux of convective GWs
gwd_rando_sat =0.6 (empty documentation)
# Important pour controler la periode de la QBO
gwd_front_ruwmax =3 (empty documentation)
###############################################################
# Rayonnement
###############################################################
# activation du nouveau code de rayonnement RRTM
# 0 : Ancien code et 1 : RRTM (D=0)
iflag_rrtm =1 (empty documentation)
# Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0
NSW =6 (empty documentation)
# Albedo et inertie des calottes
alb_vis_sno_lic =0.96 (empty documentation)
alb_nir_sno_lic =0.68 (empty documentation)
# 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)
# inertie_lic/sic : pour les "land ice" et "sea ice"
inertie_sno =350 (empty documentation)
iflag_sic =1 (empty documentation)
# a inertie_sic meme si la surface est couverte de neige
iflag_sic =1 (empty documentation)
inertie_lic =2000 (empty documentation)
inertie_sic =2000 (empty documentation)
# tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours)
tau_gl =2 (empty documentation)
tau_gl =2 (empty documentation)
# ???
iflag_albedo =1 (empty documentation)
# Frequence appel convection. Nombre appels par jour
# (nbapp_cv=0 ==> appel tous les pas de temps)
nbapp_cv =48 (empty documentation)
# 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 (empty documentation)
###############################################################
# Parametres nuages
###############################################################
#seuils de la precipitation des nuages stratiformes (D: 2.6e-4 2.6e-4)
cld_lc_lsc =0.00065
cld_lc_con =0.00065
#constante de temps pour eliminer eau lsc et convective
# (D: 3600. 3600.)
cld_tau_lsc =900 (empty documentation)
cld_tau_con =900 (empty documentation)
#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 (empty documentation)
# 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 (empty documentation)
# 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 (empty documentation)
t_glace_max =273.15 (empty documentation)
exposant_glace =0.5 (empty documentation)
iflag_t_glace =3 (empty documentation)
# Activation de la thermodynamique de la glace
iflag_ice_thermo =1 (empty documentation)
# 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 (empty documentation)
#Bergeron
iflag_bergeron =2 (empty documentation)
# 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 (empty documentation)
#calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4)
iflag_cld_th =6 (empty documentation)
fact_cldcon =1 (empty documentation)
# ?????
facttemps =0 (empty documentation)
#calcul eau condensee et fraction nuageuse a partir des PDF
# (D:0, 0:version avec ratqs sinon nouvelles PDFs)
iflag_pdf =1 (empty documentation)
#calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
ok_newmicro =y (empty documentation)
# Controle des ratqs grande echelle
iflag_ratqs =4 (empty documentation)
ratqsp0 =45000
ratqsdp =10000
ratqsbas =0.002 (empty documentation)
ratqshaut =0.4
#rayon effectif pour les nuages de glace (D:35)
rad_froid =35 (empty documentation)
# Rayons effectif des gouttes nuageuses
# ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns)
rad_chau1 =12 (empty documentation)
rad_chau2 =11 (empty documentation)
# Rayon effectif des gouttes si ok_cdnc=0
# Parameters in CDNC-maer link (Boucher&Lohmann 1995)
bl95_b0 =1.3 (empty documentation)
bl95_b1 =0.2 (empty documentation)
# Choix ou non de New oliq
new_oliq =y (empty documentation)
# Seuils oliqmax, oicemax (en g/kg).
# Au dela on precipite directement au sol
oliqmax =0.0015 (empty documentation)
oicemax =0.0015 (empty documentation)
#Rayons min et max effectifs des particules des nuages
rei_min =16 (empty documentation)
rei_max =61.29 (empty documentation)
# 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 (empty documentation)
cloudth_vert_alpha =0.5 (empty documentation)
cloudth_vert_alpha_th =0.5 (empty documentation)
######################################################################
# Convection
######################################################################
#Convection scheme switch
# (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
iflag_con =3 (empty documentation)
# correct convection bug
#(default = y)
keepbug_ice_frac =n (empty documentation)
#output level of energy conservation diagnostics
if_ebil =0 (empty documentation)
#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 (empty documentation)
iflag_clw =0 (empty documentation)
#Convective closure (D:1, 0:AR4, 2:ALE and ALP)
iflag_clos =2 (empty documentation)
#Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
iflag_mix =1 (empty documentation)
iflag_mix =1 (empty documentation)
qqa1 =1 (empty documentation)
qqa2 =0 (empty documentation)
#Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
cvl_corr =1 (empty documentation)
#Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
Fmax =0.65 (empty documentation)
#altitude, in hPa, above LCL at which buoybase is computed,
# def = original: -40
dpbase =-40 (empty documentation)
#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 (empty documentation)
#relaxation time of cloud base mass flux, def = original: 8000
tau =8000 (empty documentation)
#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 (empty documentation)
tau_stop =15000 (empty documentation)
ok_intermittent =n (empty documentation)
# nouvel epluchage de C.Rio
#shedding coefficient (used when iflag_mix_adiab=1)
iflag_mix_adiab =0 (empty documentation)
coef_peel =0.25 (empty documentation)
#Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
#if 1, def=1
flag_epKEorig =1 (empty documentation)
flag_epKEorig =1 (empty documentation)
# elcrit : max condensed water concentration in kg/kg, D=0.0003,
# tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
flag_epKEorig =1 (empty documentation)
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 (empty documentation)
ptcrit =500 (empty documentation)
# Precipitation ejection from adiabatic ascents
# cvflag_prec_eject: when TRUE ejection is activated. Default: False
# WARNING ! : cvflag_prec_eject=y is effective only if icvflag_Tpa=1
# ejectliq : 0. =>liquid precipitation is carried up in the adiab updraughts
# 1. =>liquid precipitation is ejected from adiab updraughts
# ejectice 0. =>solid precipitation is carried up in the adiab updraughts
# 1. =>solid precipitation is ejected from adiab updraughts
# ejectliq may take only 0. and 1. values
# ejectice may take any value between 0. and 1.
cvflag_prec_eject =n (empty documentation)
cvflag_prec_eject =n (empty documentation)
ejectliq =0. (empty documentation)
ejectice =0. (empty documentation)
# Check whether parcel level temperature and specific humidity
# are reasonable and demand that LCL be high enough
cv_flag_feed =2 (empty documentation)
# save 20% CPU
ok_optim_yield =y (empty documentation)
# Convective clouds
iflag_cld_cv =0 (empty documentation)
# To get back to bug where convective output variables
# were set to 0 at the start of the physics
ok_bug_cv_trac =n (empty documentation)
# solving of the temperature equation for the adiabatic ascent
# 0 -> AJ scheme in two loops ; ice fraction given by envt. temp.
# 1 -> JYG scheme in one loop ; ice fraction given by envt. temp.
# 2 -> JYG scheme in one loop ; ice fraction given by adiab. asc. temp.
# Default: 0
# 2 is costly and ill tested; experts only.
icvflag_Tpa =0 (empty documentation)
# taking into account the dependence of qv on qt at saturation
# Default: n
qsat_depends_on_qt =n (empty documentation)
##########################################################
# Fermeture convective et couplages
##########################################################
# Decimal code: + first digit drives the splitting of thermals:
# 0 = thermals see the average profiles of T and q
# 1 = thermals are only outside wakes
# + second digit drives the splitting of diffuse PBL:
# 0 = The diffuse PBL sees the average profiles of T and q
# 1 = Two BLs, one in (w) and one in (x), with the same
# surface temperature and humidity.
# 2 = Two BLs, one in (w) and one in (x), with separate
# surface temperature and humidity.
# 3 = same as 1 over ocean; same as 0 everywhere else.
# Example: iflag_pbl_split=12 => thermals outside wakes, split BL with Tsurf and qsurf difference
iflag_pbl_split =10 (empty documentation)
iflag_pbl_split =10 (empty documentation)
#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 (empty documentation)
#Convection must reach a level at which T
t_top_max =1000 (empty documentation)
### nrlmd le 10/04/2012: Declenchement ###
#Declenchement par couche-limite (0: deterministe, 1: stochastique)
iflag_trig_bl =1 (empty documentation)
#Surface critique au LCL
s_trig =12000000 (empty documentation)
#Temps de decorrelation entre scenes de cumulus
tau_trig_shallow =1200 (empty documentation)
#Temps de decorrelation entre scenes de cumulonimbus
tau_trig_deep =1200 (empty documentation)
# Pour la stabilisation
random_notrig_max =0.99 (empty documentation)
# Temperature a la quelle on fond les descentes glacees
tmax_fonte_cv =278.15 (empty documentation)
#Fermeture par couche-limite (Nicolas Rochetin)
#(0: deterministe, 1: statistique, 2: statistique+convergence
# grande echelle)
iflag_clos_bl =1 (empty documentation)
#coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
# def: 0.33
stark =0.33 (empty documentation)
#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 (empty documentation)
# pas de temps du spliting pour les thermiques
# TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
nsplit_thermals =1 (empty documentation)
# tau_thermals pour avoir une constante de temps sur les thermiques.
# Non valide
tau_thermals =0 (empty documentation)
# Flag controlant le detrainement et l entrainement
iflag_thermals_ed =8 (empty documentation)
# 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 (empty documentation)
#couplage ou non avec la convection
iflag_coupl =5 (empty documentation)
# Seuil pour desactiver les thermiques en cas de forte inversion
# en sommet de couche limite. Actif si iflag_thermals impair.
seuil_inversion =-0.08 (empty documentation)
# Reactivation du calcul d un zmax continu pour les thermiques
# reactivation d une fermeture humide des thermiques
iflag_thermals_closure =2 (empty documentation)
iflag_thermals_closure =2 (empty documentation)
# Flag permettant de garder le bug induisant des wake_deltaT non nuls
# a partir des tendances des thermiques, meme en l'absence de poches.
# Default : ok_bug_split_th=y ==> le bug est maintenu
ok_bug_split_th =n (empty documentation)
# contribution of buoyancy to the plume acceleration
thermals_afact =.66666666666666666666 (empty documentation)
# drag term in the plume acceleration
thermals_fact_epsilon =0.002 (empty documentation)
# scaling factor for entrainment and detrainment
thermals_betalpha =0.9 (empty documentation)
# width of the environment subgrid scale water distribution
cloudth_sigma1s_factor =1.1 (empty documentation)
##########################################################
# Flags wakes
##########################################################
#Wake scheme switch (D:0, 0:AR4, 1:New Physics)
iflag_wake =1 (empty documentation)
#multiplicative factor of the damping by gravity waves, def: 4.
coefgw =4 (empty documentation)
#wake density = number of wake centers per m2, def: 8.E-12
#wdens_ref=8.E-12
wdens_ref_o =1e-09 (empty documentation)
wdens_ref_l =8e-12 (empty documentation)
# Ajustement convectif prealable au calcul des poches
ok_adjwk =y (empty documentation)
# Prevent some crashes
# Filter out bad wakes
flag_wk_check_trgl =n (empty documentation)
iflag_wk_check_trgl =2 (empty documentation)
iflag_alp_wk_cond =1 (empty documentation)
#Wake population dynamic switch (D:0):
iflag_wk_pop_dyn =0 (empty documentation)
# iflag_wk_pop_dyn=1 ==> wake_dens is a prognostic variable
iflag_wk_pop_dyn =0 (empty documentation)
####################################################################
# 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 (empty documentation)
#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 (empty documentation)
fl_cor_ebil =1 (empty documentation)
####################################################################
# Ozone
####################################################################
# new treatment for ozone into the tropopause
adjust_tropopause =y (empty documentation)
####################################################################
# Conservation de l'eau
#####################################################################
# conservation de l eau sur les calottes
ok_lic_cond =y (empty documentation)
####################################################################
#min cloud droplet number concentration
cdnc_min =10. (empty documentation)
# Activation condensation de grande echelle epuree par E.Vignon : ok_new_lscp:y + ok_icefra_lscp:y
ok_new_lscp =y
ok_icefra_lscp =y