Index: LMDZ6/trunk/libf/phylmd/yamada4.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/yamada4.F90	(revision 4288)
+++ LMDZ6/trunk/libf/phylmd/yamada4.F90	(revision 4289)
@@ -444,5 +444,5 @@
          tkeprov= tkeprov*                           &
            &  (tkeprov+dt*(shear(ig,k)+max(0.,buoy(ig,k))))/ &
-           &  (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)))
+           &  (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)+drgpro(ig,k)*tkeprov))
          q2(ig,k)=tkeprov*ydeux
         ENDDO
@@ -497,10 +497,4 @@
          DO ig=1,ngrid
          tkeprov=q2(ig,k)/ydeux
-
-!             if(ifl_pbltree .eq. 0) then
-!         disseff=dissip(ig,k)-min(0.,buoy(ig,k))
-!         tkeexp=exp(-dt*disseff/tkeprov)
-!         tkeprov= shear(ig,k)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp
-!           else
 !FC on ajoute la dissipation due aux arbres
          disseff=dissip(ig,k)-min(0.,buoy(ig,k)) + drgpro(ig,k)*tkeprov
@@ -510,8 +504,5 @@
          tkeprov= (shear(ig,k)+ &
           & drgpro(ig,k)*(winds(ig,k))**3)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp
-!               endif
-
          q2(ig,k)=tkeprov*ydeux
-
          ! En cas stable, on traite la flotabilite comme la
          ! dissipation, en supposant que buoy/q2^3 est constant.
@@ -522,14 +513,7 @@
        DO k = 2, klev - 1
          DO ig=1,ngrid
-         !tkeprov=q2(ig,k)/ydeux
-         !tkeprov=tkeprov+max(buoy(ig,k)+shear(ig,k),0.)*dt
-         !disseff=dissip(ig,k)-min(0.,buoy(ig,k)+shear(ig,k)) 
-         !tkeexp=exp(-dt*disseff/tkeprov)
-         !tkeprov= tkeprov*tkeexp
-         !q2(ig,k)=tkeprov*ydeux
          ! En cas stable, on traite la flotabilite comme la
          ! dissipation, en supposant que dissipeff/TKE est constant.
          ! Puis on prend la solution exacte
-         ! 
          ! With drag and dissipation from high vegetation (EV & FC, 05/10/2020)
          winds(ig,k)=sqrt(u(ig,k)**2+v(ig,k)**2)