Index: LMDZ6/branches/LMDZ-ECRAD/DefLists/physiq.def_NPv6.1
===================================================================
--- LMDZ6/branches/LMDZ-ECRAD/DefLists/physiq.def_NPv6.1	(revision 3879)
+++ LMDZ6/branches/LMDZ-ECRAD/DefLists/physiq.def_NPv6.1	(revision 4171)
@@ -1,4 +1,4 @@
 #
-# $Id: physiq.def_NPv6.0.15_S 3593 2018-02-20 10:07:02Z jgipsl $
+# $Id: physiq.def_NPv6.1 4020 2018-09-05 17:28:38Z jgipsl $
 #
 
@@ -163,11 +163,10 @@
 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= 2.
+
+
           # ???
 iflag_albedo=1
 
-          # Frequence appel convection. Nombre appels par jour
+          # Frequence appel rayonnement. Nombre appels par jour
 nbapp_cv=48
 
@@ -369,149 +368,4 @@
 ptcrit=500
 
-          # Precipitation ejection from adiabatic ascents
-          # cvflag_prec_eject: when TRUE ejection is activated. Default: False
-          # 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.
-          #   Default : cvflag_prec_eject=n, ejectliq=0., ejectice=0.
-cvflag_prec_eject=n
-ejectliq=0.
-ejectice=0.
-
-
-          # 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
-
-           # 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
-
-           # taking into account the dependence of qv on qt at saturation
-           # Default: n
-qsat_depends_on_qt=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 T<t_top_max to be active
-t_top_max=1000
-
-          ### nrlmd le 10/04/2012: Declenchement ###
-          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
-iflag_trig_bl=1
-
-          #Surface critique au LCL
-s_trig=12000000
-
-          #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
-
-          # Temperature a la quelle on fond les descentes glacees
-tmax_fonte_cv=278.15
-
-          #Fermeture par couche-limite (Nicolas Rochetin)
-          #(0: deterministe, 1: statistique, 2: statistique+convergence
-          # grande echelle)
-iflag_clos_bl=1
-
-          #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
-
-#
-##########################################################
-# 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
-
           # Prevent some crashes
           # Filter out bad wakes
@@ -520,4 +374,124 @@
 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
+#libIGCM version : ok_bug_cv_trac=_AUTO_: DEFAULT = n
+ok_bug_cv_trac=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 T<t_top_max to be active
+t_top_max=1000
+
+          ### nrlmd le 10/04/2012: Declenchement ###
+          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
+iflag_trig_bl=1
+
+          #Surface critique au LCL
+s_trig=12000000
+
+          #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
+
+          # Temperature a la quelle on fond les descentes glacees
+tmax_fonte_cv=278.15
+
+          #Fermeture par couche-limite (Nicolas Rochetin)
+          #(0: deterministe, 1: statistique, 2: statistique+convergence
+          # grande echelle)
+iflag_clos_bl=1
+
+          #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
+
+#
+##########################################################
+# 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
 
 
@@ -557,4 +531,2 @@
 
 cdnc_min=10.
-
-
