Index: LMDZ4/trunk/libf/phylmd/cva_driver.F
===================================================================
--- LMDZ4/trunk/libf/phylmd/cva_driver.F	(revision 1061)
+++ LMDZ4/trunk/libf/phylmd/cva_driver.F	(revision 1062)
@@ -733,5 +733,5 @@
       IF (iflag_con .eq. 3) THEN
        IF (iflag_mix .ge. 1 ) THEN
-    
+         CALL zilch(supmax,nloc*klev)    
          CALL cv3p_mixing(nloc,ncum,nd,nd,ntra,icb,nk,inb    ! na->nd
      :                       ,ph,t,q,qs,u,v,tra,h,lv,qnk
Index: LMDZ4/trunk/libf/phylmd/phys_output_mod.F90
===================================================================
--- LMDZ4/trunk/libf/phylmd/phys_output_mod.F90	(revision 1061)
+++ LMDZ4/trunk/libf/phylmd/phys_output_mod.F90	(revision 1062)
@@ -265,4 +265,12 @@
   integer, dimension(nfiles) , save :: flag_dtwak        = (/ 4, 5, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dqwak        = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_wake_h       = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_wake_s       = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_wake_deltat  = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_wake_deltaq  = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_wake_omg     = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_Vprecip      = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_ftd          = (/ 4, 5, 10, 10, 10 /)
+  integer, dimension(nfiles) , save :: flag_fqd          = (/ 4, 5, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtlsc        = (/ 4, 10, 10, 10, 10 /)
   integer, dimension(nfiles) , save :: flag_dtlschr      = (/ 4, 10, 10, 10, 10 /)
@@ -408,5 +416,6 @@
         jmax_ins=jjmp1
 
-        do i=1,iim-1
+! correction abderr        do i=1,iim-1
+        do i=1,iim
            print*,'io_lon(i)=',io_lon(i)
            if (io_lon(i).le.lonmin_reg(iff)) imin_ins=i
@@ -705,8 +714,17 @@
  IF(iflag_con.GE.3) THEN
  IF (iflag_wake == 1) THEN
+   CALL histdef2d(iff,flag_wake_h,"wake_h","wake_h", "-")
+   CALL histdef2d(iff,flag_wake_s,"wake_s","wake_s", "-")
    CALL histdef3d(iff,flag_dtwak,"dtwak","Wake dT","K/s")
    CALL histdef3d(iff,flag_dqwak,"dqwak","Wake dQ","(kg/kg)/s")
+   CALL histdef3d(iff,flag_wake_deltat,"wake_deltat","wake_deltat", " ")
+   CALL histdef3d(iff,flag_wake_deltaq,"wake_deltaq","wake_deltaq", " ")
+   CALL histdef3d(iff,flag_wake_omg,"wake_omg","wake_omg", " ")
  ENDIF
+   CALL histdef3d(iff,flag_Vprecip,"Vprecip","precipitation vertical profile","-")
+   CALL histdef3d(iff,flag_ftd,"ftd","tend temp due aux descentes precip","-")
+   CALL histdef3d(iff,flag_fqd,"fqd","tend vap eau due aux descentes precip","-")
  ENDIF !(iflag_con.GE.3)
+
  CALL histdef3d(iff,flag_dtlsc,"dtlsc","Condensation dT", "K/s")
  CALL histdef3d(iff,flag_dtlschr,"dtlschr","Large-scale condensational heating rate","K/s")