Index: /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_plume_6A.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_plume_6A.f90	(revision 5619)
+++ /LMDZ6/trunk/libf/phylmd/lmdz_thermcell_plume_6A.f90	(revision 5620)
@@ -224,6 +224,5 @@
         zta_est(ig,l)=ztva_est(ig,l)
         ztva_est(ig,l) = ztva_est(ig,l)/zpspsk(ig,l)
-        ztva_est(ig,l) = ztva_est(ig,l)*(1.+RETV*(zqta(ig,l-1)  &
-     &      -zqla_est(ig,l))-zqla_est(ig,l))
+        ztva_est(ig,l) = ztva_est(ig,l)*(1.+RETV*(zqta(ig,l-1)-zqla_est(ig,l)))
  
 
@@ -566,6 +565,5 @@
 !on rajoute le calcul de zha pour diagnostiques (temp potentielle)
            zha(ig,l) = ztva(ig,l)
-           ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)  &
-     &              -zqla(ig,l))-zqla(ig,l))
+           ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)-zqla(ig,l)))
            zbuoy(ig,l)=RG*(ztva(ig,l)-ztv(ig,l))/ztv(ig,l)
            zdz=zlev(ig,l+1)-zlev(ig,l)