Index: LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90	(revision 4412)
+++ LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90	(revision 4413)
@@ -8,4 +8,5 @@
    integer        :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure
    integer        :: iflag_thermals_down
+   real           :: entrdn_coef, detrdn_coef
 
 !$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout)
@@ -14,4 +15,5 @@
 !$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure)
 !$OMP THREADPRIVATE(iflag_thermals_down)
+!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef)
 
 
@@ -91,4 +93,11 @@
 
 
+   entrdn_coef = 0.5
+   CALL getin_p('entrdn_coef',entrdn_coef)
+
+   detrdn_coef = 0.5
+   CALL getin_p('detrdn_coef',detrdn_coef)
+
+
      CALL getin_p('thermals_fact_epsilon',fact_epsilon)
      CALL getin_p('thermals_betalpha',betalpha)
Index: LMDZ6/trunk/libf/phylmd/thermcell_main.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_main.F90	(revision 4412)
+++ LMDZ6/trunk/libf/phylmd/thermcell_main.F90	(revision 4413)
@@ -19,4 +19,5 @@
       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: RD,RG
 
@@ -134,4 +135,5 @@
       real, dimension(ngrid,nlay) :: zqsat,csc
       real, dimension(ngrid) :: zcon,zcon2,alim_star_tot,f
+      real, dimension(ngrid,nlay) :: entrdn,detrdn
 
       character (len=20) :: modname='thermcell_main'
@@ -502,9 +504,11 @@
         print*,'AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA'
         print*,'WARNING !!! routine thermcell_down en cours de developpement'
-        ! on veut transporter la temperature potentielle, l'eau totale, qté de mouvement
-        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zthl,zdthladj)
-        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,po,pdoadj)
-        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zu,pduadj)
-        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zv,pdvadj)
+        entrdn=entrdn_coef*detr0
+        detrdn=detrdn_coef*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)
+        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,po,pdoadj)
+        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zu,pduadj)
+        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zv,pdvadj)
       ELSE
       !--------------------------------------------------------------