Index: trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3336)
+++ trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3337)
@@ -1016,6 +1016,4 @@
 c           nsubtimestep(ig)=1 !for debug
            subtimestep = ptimestep/nsubtimestep(ig)
-                          call write_output('subtimestep',
-     &                'vdifc substimestep length','s',subtimestep)
             ztsrf(ig)=ptsrf(ig,islope)  !  +pdtsrf(ig)*subtimestep
             zq_tmp_vap(ig,:,:) =zq(ig,:,:)
@@ -1414,4 +1412,7 @@
      &                     "Ground ice latent heat flux",
      &                     "W.m-2",surf_h2o_lh(:,iflat))
+         call write_output("nsubtimestep",
+     &                     "Number of subtimesteps in vdifc",
+     &                     "int",nsubtimestep(:))
 
          call write_output('zdqsdif_ssi_frost_tot',
@@ -1498,5 +1499,5 @@
 
       do i=1,naersize
-        if ((qsurf(i).eq.0).and.
+        if ((qsurf(i).lt.0.005).and.
      &      (.not.watercaptag(i))) then
           ntsub(i) = 1