Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2833)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2840)
@@ -3025,7 +3025,7 @@
 c    &                  zstress)
         call WRITEDIAGFI(ngrid,'sw_htrt','sw heat. rate',
-     &                   'w.m-2',3,zdtsw)
+     &                   'K/s',3,zdtsw)
         call WRITEDIAGFI(ngrid,'lw_htrt','lw heat. rate',
-     &                   'w.m-2',3,zdtlw)
+     &                   'K/s',3,zdtlw)
         call writediagfi(ngrid,"local_time","Local time",
      &                   'sol',2,local_time)
Index: trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 2833)
+++ trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F	(revision 2840)
@@ -136,4 +136,6 @@
       REAL ztsrf(ngrid) ! temporary surface temperature in tsub
       REAL zdtsrf(ngrid) ! surface temperature tendancy in tsub
+      REAL surf_h2o_lh(ngrid) ! Surface h2o latent heat flux
+      REAL zsurf_h2o_lh(ngrid) ! Tsub surface h2o latent heat flux
 
 c     For latent heat release from ground water ice sublimation    
@@ -242,4 +244,6 @@
       pdtsrf(1:ngrid)=0
       zdtsrf(1:ngrid)=0
+      surf_h2o_lh(1:ngrid)=0
+      zsurf_h2o_lh(1:ngrid)=0
       pdqdif(1:ngrid,1:nlay,1:nq)=0
       pdqsdif(1:ngrid,1:nq)=0
@@ -959,4 +963,8 @@
      &                       zdqsdif(ig))*subtimestep
 
+c             Monitoring instantaneous latent heat flux in W.m-2 :
+              zsurf_h2o_lh(ig) = zsurf_h2o_lh(ig)+
+     &                               (zdtsrf(ig)*pcapcal(ig))
+     &                               *subtimestep
 
 c             We ensure that surface temperature can't rise above the solid domain if there
@@ -975,4 +983,5 @@
 c             (btw could also compute the post timestep temp and ice
 c             by simply adding the subtimestep trend instead of this)
+            surf_h2o_lh(ig)= zsurf_h2o_lh(ig)/ptimestep
             pdtsrf(ig)= (ztsrf(ig) -
      &                     ptsrf(ig))/ptimestep
@@ -1048,4 +1057,7 @@
 c --------- end of tracers  ----------------------------
 
+         call WRITEDIAGFI(ngrid,"surf_h2o_lh",
+     &                          "Ground ice latent heat flux",
+     &                               "W.m-2",2,surf_h2o_lh(:))
 
 C       Diagnostic output for HDO