Index: LMDZ5/branches/testing/DefLists/README
===================================================================
--- LMDZ5/branches/testing/DefLists/README (revision 2669)
+++ LMDZ5/branches/testing/DefLists/README (revision 2669)
@@ -0,0 +1,12 @@
+physiq.def Couples model version
+
+physiq.def_L39_AR40.0 IPSL-CM5A
+physiq.def_L39_NPv3.1 IPSL-CM5B
+physiq.def_L39_NPv3.2 Stabilized version
+physiq.def_L59_NPv4.12 First L59 version of the new new physics
+physiq.def_L79_NPv5.17h (5.17h=5.3) IPSL-CM6.0.1
+physiq.def_NPv5.4 IPSL-CM6.0.2
+physiq.def_NPv5.5 IPSL-CM6.0.3
+physiq.def_NPv5.65 IPSL-CM6.0.4
+physiq.def_NPv5.70 IPSL-CM6.0.5
+physiq.def_NPv5.67 IPSL-CM6.0.5 / THC1 (test de sensibilité)
Index: LMDZ5/branches/testing/DefLists/context_lmdz.xml
===================================================================
--- LMDZ5/branches/testing/DefLists/context_lmdz.xml (revision 2641)
+++ LMDZ5/branches/testing/DefLists/context_lmdz.xml (revision 2669)
@@ -31,5 +31,10 @@
-->
-
+
+
+
+
+
+
Index: LMDZ5/branches/testing/DefLists/field_def_lmdz.xml
===================================================================
--- LMDZ5/branches/testing/DefLists/field_def_lmdz.xml (revision 2641)
+++ LMDZ5/branches/testing/DefLists/field_def_lmdz.xml (revision 2669)
@@ -867,4 +867,24 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
Index: LMDZ5/branches/testing/DefLists/file_def_histmth_lmdz.xml
===================================================================
--- LMDZ5/branches/testing/DefLists/file_def_histmth_lmdz.xml (revision 2641)
+++ LMDZ5/branches/testing/DefLists/file_def_histmth_lmdz.xml (revision 2669)
@@ -10,4 +10,26 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: LMDZ5/branches/testing/DefLists/gcm.def_144x142
===================================================================
--- LMDZ5/branches/testing/DefLists/gcm.def_144x142 (revision 2641)
+++ LMDZ5/branches/testing/DefLists/gcm.def_144x142 (revision 2669)
@@ -1,5 +1,5 @@
-## $Id: gcm.def 1403 2010-07-01 09:02:53Z fairhead $
+## $Id$
## nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )
-day_step=960
+day_step=720
## periode pour le pas Matsuno (en pas)
iperiod=5
@@ -13,10 +13,9 @@
niterh=2
## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)
-tetagdiv=1800.
+tetagdiv=3600.
## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
-#tetagrot=5400.
tetagrot=7200.
## temps de dissipation des plus petites long.d ondes pour h ( divgrad)
-tetatemp=5400.
+tetatemp=7200.
## coefficient pour gamdissip
coefdis=0.
@@ -33,8 +32,6 @@
read_start=y
## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1)
+## iphysiq est modifie pour la nouvelle physique par le lmdz.driver pour les autres cas le valeur par defaut est utilise
iphysiq=_AUTO_ : DEFAULT = 5
-#
-#
-### Parametres du zoom
## longitude en degres du centre du zoom
clon=0.
@@ -57,2 +54,4 @@
## Fonction f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.
ysinus=y
+##
+ngroup=4
Index: LMDZ5/branches/testing/DefLists/iodef_dev.xml
===================================================================
--- LMDZ5/branches/testing/DefLists/iodef_dev.xml (revision 2641)
+++ LMDZ5/branches/testing/DefLists/iodef_dev.xml (revision 2669)
@@ -1,5 +1,8 @@
-
+
+
+
+
+ domain_ref="dom_glo" axis_ref="presnivs">
@@ -46,6 +50,4 @@
-
@@ -58,6 +60,4 @@
-
@@ -65,4 +65,5 @@
name="Xdiurnalave"
output_freq="1d"
+ type="one_file"
enabled=".true.">
@@ -93,10 +94,4 @@
wxios_add_vaxis -->
-
- true
+ false
10
Index: LMDZ5/branches/testing/DefLists/physiq.def_L79_NPv5.4
===================================================================
--- LMDZ5/branches/testing/DefLists/physiq.def_L79_NPv5.4 (revision 2641)
+++ (revision )
@@ -1,229 +1,0 @@
-##############################################################################RAMETRES 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
-iflag_albedo=1
-#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
-
-ok_strato=y
-ok_hines=n
-sso_gkdrag=0.2
-sso_gfrcri=0.70
-sso_gkwake=0.40
-#
-# 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.0001
-cld_lc_con=0.0001
-#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=1.
-ffallv_con=1.
-#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.9995
-#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.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.1
-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
-#
-# 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)
-iflag_trig_bl=1
-#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
-# Pour la stabilisation
-random_notrig_max=0.99
-ok_adjwk=y
-tmax_fonte_cv=288.15
-
-
-### 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=233.13
-t_glace_max=273.13
-exposant_glace=1.
-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
-
-### distribution sous-maille sur la verticale
-iflag_cloudth_vert=1
-
-### Convection
-#altitude, in hPa, above LCL at which buoybase is computed, def = original: -40
-dpbase=-40.
-#altitude, in hPa, above LCL below which Ep=0, def = original: 150, used only if flag_epKEorig=0
-pbcrit=150.0
-#altitude, in hPa, above LCL above which Ep=Epmax, def = original: 500, used only if flag_epKEorig=0
-ptcrit=500.0
-#reference fractional area of precipitating downdraughts, def = original: 0.01
-sigdz=0.01
-#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, def = 1
-flag_wb=1
-#parameter used in the computation of wb, def = 6.
-wbmax=6.
-#
-ok_convstop=.False.
-#
-tau_stop=15000.
-#
-ok_intermittent=.False.
-#shedding coefficient (used when iflag_mix_adiab=1)
-coef_peel=0.25
-#
-### Clouds
-#Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) if 1, def=1
-flag_epKEorig=1
-#max condensed water concentration in kg/kg, def=0.0003, original: 0.0011, used only if flag_epKEorig=1
-elcrit=0.0003
-#temperature in degre Celsius, at which Ep=Epmax, def=-55., original: -55, used only if flag_epKEorig=1
-tlcrit=-55.0
-
-### Wakes
-#coefficient in Wape->C_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
-#wake density = number of wake centers per m2, def: 8.E-12
-wdens_ref=8.E-12
-#multiplicative factor of the damping by gravity waves, def: 4.
-coefgw=4.
Index: LMDZ5/branches/testing/DefLists/physiq.def_L79_NPv5.5
===================================================================
--- LMDZ5/branches/testing/DefLists/physiq.def_L79_NPv5.5 (revision 2641)
+++ (revision )
@@ -1,236 +1,0 @@
-##############################################################################
-# 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
-iflag_albedo=1
-#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
-
-ok_strato=y
-ok_hines=n
-sso_gkdrag=0.2
-sso_gfrcri=0.70
-sso_gkwake=0.40
-#
-# Rayonnement
-#*********************
-# activation du nouveau code de rayonnement RRTM
-# 0 : Ancien code et 1 : RRTM (D=0)
-iflag_rrtm=1
-NSW=6
-#
-# Parametres nuages
-#******************
-#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
-cld_lc_lsc=0.00022
-cld_lc_con=0.00022
-#constante de temps pour eleminer eau lsc et convective (D: 3600. 3600.)
-cld_tau_lsc=1200.
-cld_tau_con=1200.
-#facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
-ffallv_lsc=0.67
-ffallv_con=0.67
-#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.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.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.1
-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
-#
-# 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)
-iflag_trig_bl=1
-#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
-# Pour la stabilisation
-random_notrig_max=0.99
-ok_adjwk=y
-tmax_fonte_cv=278.15
-
-
-### 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=233.13
-t_glace_max=273.13
-exposant_glace=1.
-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
-
-iflag_cloudth_vert=1
-
-#
-lmixmin=10.
-
-#Bergeron
-iflag_bergeron=1
-
-
-### Convection
-#altitude, in hPa, above LCL at which buoybase is computed, def = original: -40
-dpbase=-40.
-#altitude, in hPa, above LCL below which Ep=0, def = original: 150, used only if flag_epKEorig=0
-pbcrit=150.0
-#altitude, in hPa, above LCL above which Ep=Epmax, def = original: 500, used only if flag_epKEorig=0
-ptcrit=500.0
-#reference fractional area of precipitating downdraughts, def = original: 0.01
-sigdz=0.01
-#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, def = 1
-flag_wb=1
-#parameter used in the computation of wb, def = 6.
-wbmax=6.
-#
-ok_convstop=.False.
-#
-tau_stop=15000.
-#
-ok_intermittent=.False.
-#shedding coefficient (used when iflag_mix_adiab=1)
-coef_peel=0.25
-#
-### Clouds
-#Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) if 1, def=1
-flag_epKEorig=1
-#max condensed water concentration in kg/kg, def=0.0003, original: 0.0011, used only if flag_epKEorig=1
-elcrit=0.0003
-#temperature in degre Celsius, at which Ep=Epmax, def=-55., original: -55, used only if flag_epKEorig=1
-tlcrit=-55.0
-
-### Wakes
-#coefficient in Wape->C_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
-#wake density = number of wake centers per m2, def: 8.E-12
-wdens_ref=8.E-12
-#multiplicative factor of the damping by gravity waves, def: 4.
-coefgw=4.
Index: LMDZ5/branches/testing/DefLists/physiq.def_NPv5.65
===================================================================
--- LMDZ5/branches/testing/DefLists/physiq.def_NPv5.65 (revision 2641)
+++ LMDZ5/branches/testing/DefLists/physiq.def_NPv5.65 (revision 2669)
@@ -36,6 +36,4 @@
gwd_rando_sat=0.2
-ok_strato=y
-ok_hines=n
sso_gkdrag=0.2
sso_gfrcri=0.70
Index: LMDZ5/branches/testing/DefLists/physiq.def_NPv5.67
===================================================================
--- LMDZ5/branches/testing/DefLists/physiq.def_NPv5.67 (revision 2641)
+++ LMDZ5/branches/testing/DefLists/physiq.def_NPv5.67 (revision 2669)
@@ -36,6 +36,4 @@
gwd_rando_sat=0.2
-ok_strato=y
-ok_hines=n
sso_gkdrag=0.2
sso_gfrcri=0.70
Index: LMDZ5/branches/testing/DefLists/physiq.def_NPv5.70
===================================================================
--- LMDZ5/branches/testing/DefLists/physiq.def_NPv5.70 (revision 2669)
+++ LMDZ5/branches/testing/DefLists/physiq.def_NPv5.70 (revision 2669)
@@ -0,0 +1,240 @@
+#
+# $Id: physiq.def_NPv5.70 2929 2016-07-26 06:49:39Z jgipsl $
+#
+##############################################################################
+# 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
+iflag_albedo=1
+#Cdrags
+f_cdrag_stable=1.
+f_cdrag_ter=1.
+f_cdrag_oce=0.6
+f_ri_cd_min=0.1
+#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
+
+sso_gkdrag=0.2
+sso_gfrcri=0.70
+sso_gkwake=0.40
+#
+# Rayonnement
+#*********************
+# activation du nouveau code de rayonnement RRTM
+# 0 : Ancien code et 1 : RRTM (D=0)
+iflag_rrtm=1
+NSW=6
+#
+# Parametres nuages
+#******************
+#seuils de la precipitation des nuages strateformes (D: 2.6e-4 2.6e-4)
+cld_lc_lsc=0.0002
+cld_lc_con=0.0002
+#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.5
+ffallv_con=0.5
+#coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
+coef_eva=2e-4
+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
+ratqsp0=37500.
+ratqsdp=20000.
+ratqsbas=0.002
+ratqshaut=0.4
+#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.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.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.1
+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
+#
+# 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 ###
+#Declenchement par couche-limite (0: deterministe, 1: stochastique)
+iflag_trig_bl=1
+#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
+# Pour la stabilisation
+random_notrig_max=0.99
+ok_adjwk=y
+tmax_fonte_cv=278.15
+
+
+### nrlmd le 10/04/2012: Fermeture ###
+#Fermeture par couche-limite (0: deterministe, 1: statistique, 2: statistique+convergence grande echelle)
+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=233.13
+t_glace_max=273.13
+exposant_glace=1.
+iflag_t_glace=2
+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
+
+iflag_cloudth_vert=2
+
+#
+lmixmin=10.
+#pbl_lmixmin_alpha=0.05
+
+#Bergeron
+iflag_bergeron=1
+
+
+### Convection
+#altitude, in hPa, above LCL at which buoybase is computed, def = original: -40
+dpbase=-40.
+#altitude, in hPa, above LCL below which Ep=0, def = original: 150, used only if flag_epKEorig=0
+pbcrit=150.0
+#altitude, in hPa, above LCL above which Ep=Epmax, def = original: 500, used only if flag_epKEorig=0
+ptcrit=500.0
+#reference fractional area of precipitating downdraughts, def = original: 0.01
+sigdz=0.01
+#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, def = 1
+flag_wb=1
+#parameter used in the computation of wb, def = 6.
+wbmax=6.
+#
+ok_convstop=.False.
+#
+tau_stop=15000.
+#
+ok_intermittent=.False.
+#shedding coefficient (used when iflag_mix_adiab=1)
+coef_peel=0.25
+#
+### Clouds
+#Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature) if 1, def=1
+flag_epKEorig=1
+#max condensed water concentration in kg/kg, def=0.0003, original: 0.0011, used only if flag_epKEorig=1
+elcrit=0.00025
+#temperature in degre Celsius, at which Ep=Epmax, def=-55., original: -55, used only if flag_epKEorig=1
+tlcrit=-40.
+
+### Wakes
+#coefficient in Wape->C_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
+#wake density = number of wake centers per m2, def: 8.E-12
+wdens_ref=8.E-12
+#multiplicative factor of the damping by gravity waves, def: 4.
+coefgw=4.