Index: LMDZ6/trunk/libf/phylmd/thermcell_down.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_down.F90	(revision 4440)
+++ LMDZ6/trunk/libf/phylmd/thermcell_down.F90	(revision 4441)
@@ -223,5 +223,5 @@
 
 
-   USE thermcell_ini_mod, ONLY : prt_level,RLvCp,RKAPPA,RETV
+   USE thermcell_ini_mod, ONLY : prt_level,RLvCp,RKAPPA,RETV,fact_thermals_down
    IMPLICIT NONE
 
@@ -271,6 +271,6 @@
       do ig=1,ngrid
          if (ilay.le.lmax(ig).and.lmax(ig)>1) then
-            edn(ig,ilay)=0.5*dup(ig,ilay)
-            ddn(ig,ilay)=0.5*eup(ig,ilay)
+            edn(ig,ilay)=fact_thermals_down*dup(ig,ilay)
+            ddn(ig,ilay)=fact_thermals_down*eup(ig,ilay)
             fdn(ig,ilay)=fdn(ig,ilay+1)+edn(ig,ilay)-ddn(ig,ilay)
             thetad(ig,ilay)=( fdn(ig,ilay+1)*thetad(ig,ilay+1) + edn(ig,ilay)*theta(ig,ilay) ) / (fdn(ig,ilay)+ddn(ig,ilay))
Index: LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90	(revision 4440)
+++ LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90	(revision 4441)
@@ -3,4 +3,5 @@
 
 save
+
    integer :: dvdq=1,dqimpl=-1,prt_level=0,lunout
    real RG,RD,RCPD,RKAPPA,RLVTT,RLvCp,RETV
@@ -8,5 +9,5 @@
    integer        :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure
    integer        :: iflag_thermals_down
-   real           :: entrdn_coef, detrdn_coef
+   real           :: fact_thermals_down
 
 !$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout)
@@ -15,5 +16,5 @@
 !$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure)
 !$OMP THREADPRIVATE(iflag_thermals_down)
-!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef)
+!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef,fact_thermals_down)
 
 
@@ -92,11 +93,6 @@
    CALL getin_p('iflag_thermals_down',iflag_thermals_down)
 
-
-   entrdn_coef = 0.5
-   CALL getin_p('entrdn_coef',entrdn_coef)
-
-   detrdn_coef = 0.5
-   CALL getin_p('detrdn_coef',detrdn_coef)
-
+   fact_thermals_down = 0.5
+   CALL getin_p('fact_thermals_down',fact_thermals_down)
 
      CALL getin_p('thermals_fact_epsilon',fact_epsilon)
Index: LMDZ6/trunk/libf/phylmd/thermcell_main.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_main.F90	(revision 4440)
+++ LMDZ6/trunk/libf/phylmd/thermcell_main.F90	(revision 4441)
@@ -18,6 +18,5 @@
       USE thermcell_ini_mod, ONLY: thermcell_ini,dqimpl,dvdq,prt_level,lunout,prt_level
       USE thermcell_ini_mod, ONLY: iflag_thermals_closure,iflag_thermals_ed,tau_thermals,r_aspect_thermals
-      USE thermcell_ini_mod, ONLY: iflag_thermals_down
-      USE thermcell_ini_mod, ONLY: entrdn_coef, detrdn_coef
+      USE thermcell_ini_mod, ONLY: iflag_thermals_down, fact_thermals_down
       USE thermcell_ini_mod, ONLY: RD,RG
 
@@ -503,6 +502,6 @@
       IF (iflag_thermals_down .GT. 0) THEN
         if (debut) print*,'WARNING !!! routine thermcell_down en cours de developpement'
-        entrdn=entrdn_coef*detr0
-        detrdn=detrdn_coef*entr0 
+        entrdn=fact_thermals_down*detr0
+        detrdn=fact_thermals_down*entr0 
         ! we want to transport potential temperature, total water and momentum
         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zthl,zdthladj)