Index: trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3400)
+++ trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 3404)
@@ -1121,6 +1121,9 @@
      &                .and.lag_layer)) then
 
-     
-!                     zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
+                    if ((-zdqsdif_tot(ig)*subtimestep)
+     &                .ge.(zqsurf(ig))) then
+                          zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
+                    endif
+
                     if((h2o_ice_depth(ig,islope).gt.0)
      &                  .and. (zqsurf(ig) .lt. tol_frost)