Index: LMDZ4/trunk/libf/phylmd/cv3_routines.F
===================================================================
--- LMDZ4/trunk/libf/phylmd/cv3_routines.F	(revision 1043)
+++ LMDZ4/trunk/libf/phylmd/cv3_routines.F	(revision 1044)
@@ -55,5 +55,5 @@
       pbcrit = 150.0
       ptcrit = 500.0
-      epmax  = 0.993
+cIM lu dans physiq.def via conf_phys.F90     epmax  = 0.993
 
       omtrain = 45.0 ! used also for snow (no disctinction rain/snow)
@@ -2093,7 +2093,15 @@
        water(il,i)=revap*revap
       else
-       evap(il,i)=-evap(il,i+1)
-     :            +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
-     :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*50.)
+cJYG/IM : ci-dessous formulation originale de KE
+c      evap(il,i)=-evap(il,i+1)
+c    :            +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
+c    :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*50.)
+c
+cJYG/IM : ci-dessous modification formulation originale de KE
+c        pour eliminer oscillations verticales de pluie se produisant
+c        lorsqu'il y a evaporation totale de la pluie
+c
+       evap(il,i)= +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
+     :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
       end if
 ccc
@@ -2249,5 +2257,5 @@
      :                    ,icb,inb,delt
      :                    ,t,rr,t_wake,rr_wake,u,v,tra
-     :                    ,gz,p,ph,h,hp,lv,cpn,th
+     :                    ,gz,p,ph,h,hp,lv,cpn,th,th_wake
      :                    ,ep,clw,m,tp,mp,rp,up,vp,trap
      :                    ,wt,water,evap,b,sigd
@@ -3016,9 +3024,9 @@
       enddo
 c
-      do i=1,nl
-       do il=1,ncum
-          th_wake(il,i)=t_wake(il,i)*(1000.0/p(il,i))**rdcp
-       enddo
-      enddo
+c     do i=1,nl
+c      do il=1,ncum
+c         th_wake(il,i)=t_wake(il,i)*(1000.0/p(il,i))**rdcp
+c      enddo
+c     enddo
 c
       do i=1,nl
Index: LMDZ4/trunk/libf/phylmd/cva_driver.F
===================================================================
--- LMDZ4/trunk/libf/phylmd/cva_driver.F	(revision 1043)
+++ LMDZ4/trunk/libf/phylmd/cva_driver.F	(revision 1044)
@@ -857,5 +857,5 @@
      :                     ,icb,inb,delt
      :                     ,t,q,t_wake,q_wake,u,v,tra
-     :                     ,gz,p,ph,h,hp,lv,cpn,th
+     :                     ,gz,p,ph,h,hp,lv,cpn,th,th_wake
      :                     ,ep,clw,m,tp,mp,qp,up,vp,trap
      :                     ,wt,water,evap,b,sigd