source: dynamico_lmdz/aquaplanet/NEW_PHYSIC/RUN_AQUAPLANET_DYNAMICO_LMDZ5/physiq.def @ 4110

##############################################################################
# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ
##############################################################################
# D: Valeur par default
#
#Parametres Orographiques et cdrags
#**********************************
##   Si=.T. ,  lecture du fichier limit avec la bonne annee 
ok_limitvrai=n
#Cdrags 
f_cdrag_stable=1.
f_cdrag_ter=1.
f_cdrag_oce=0.7
f_ri_cd_min=0.6
#cdmmax
cdmmax=2.5E-3
#cdhmax
cdhmax=2.0E-3
##  Orodr  ou  non   pour l orographie
ok_orodr=y
##  Orolf  ou  non   pour l orographie              
ok_orolf=y
#Rugoro
f_rugoro=0.
#
### 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=2.5
#Scale the saturated momentum flux of  convective GWs 
gwd_rando_sat=0.2
#
# Rayonnement
#*********************
# activation du nouveau code de rayonnement RRTM
# 0 : Ancien code et 1 : RRTM (D=0)
iflag_rrtm=0
#NSW=2
#
# Parametres nuages
#******************
#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
#cld_lc_lsc=0.000192
#cld_lc_con=0.000192
cld_lc_lsc=0.0005
cld_lc_con=0.0005
#constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.)
cld_tau_lsc=3600.
cld_tau_con=3600.
#facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
ffallv_lsc=0.66528
ffallv_con=0.66528
#coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
coef_eva=2e-5
reevap_ice=y
#calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) 
iflag_cldcon=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
iflag_ratqs=4
ratqsbas=0.002
ratqshaut=0.312
#rayon effectif pour les nuages de glace (D:35)
rad_froid=35
#rayons effectifs pour les nuages eau liq (tailles des gouttes d eau nuageuse, D: 13 9)
rad_chau1=12
rad_chau2=11
##  Choix ou non  de  New oliq               
new_oliq=y
#Rayons min et max effectifs des particules des nuages
rei_min=16.
rei_max=61.29
#
# Flags 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.9949 
#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=0
#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.0
#Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
Fmax=0.65
#
# Flags Couche limite
#********************
#shema de CL des thermiques (D:0, 0:ajustement sec,=>1:versions thermiques) 
iflag_thermals=18
#pas de temps du spliting pour les thermiques
nsplit_thermals=1
tau_thermals=0.
iflag_thermals_ed=8
fact_thermals_ed_dz=0.04
iflag_thermals_optflux=0
#shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) 
iflag_pbl=11
#diffusion turbulente
ksta_ter=1.e-7
ksta=1.e-10
#ok_kzmin : calcul Kzmin dans la CL de surface
ok_kzmin=y
#couplage ou non avec la convection
iflag_coupl=5 (0:AR4 et 1:nouvelle physique)
seuil_inversion=-0.08
#rendement sur la puissance fournie par les thermiques a la convection
alp_bl_k=0.5
#
## 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 
## 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=1
#
iflag_cld_cv=2
iflag_cloudth_vert=1 
# Flags wakes
#************
#Wake scheme switch (D:0, 0:AR4, 1:New Physics)
iflag_wake=1
#Constant value to be added to Alp (D: 0.)
alp_offset=0.

### nrlmd le 10/04/2012: Declenchement ###
#Déclelenchement par couche-limite (0: deterministe, 1: stochastique)
#ym iflag_trig_bl=1
iflag_trig_bl=0
#Surface critique au LCL
s_trig=1.2e7
#Temps de decorrelation entre scenes de cumulus
tau_trig_shallow=1200
#Temps de decorrelation entre scenes de cumulonimbus
tau_trig_deep=1200
### fin nrlmd le 10/04/2012

### nrlmd le 10/04/2012: Fermeture ###
#Fermeture par couche-limite (0: deterministe, 1: statistique, 2: statistique+convergence grande échelle)
iflag_clos_bl=1
### fin nrlmd le 10/04/2012

## 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=240.
t_glace_max=273.13
exposant_glace=1.7
iflag_t_glace=1
iflag_ice_thermo=1

#
## 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=1

#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 

tmax_fonte_cv=288.15