Index: /trunk/LMDZ.MARS/libf/phymars/vdifc.F
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/vdifc.F	(revision 741)
+++ /trunk/LMDZ.MARS/libf/phymars/vdifc.F	(revision 742)
@@ -812,15 +812,15 @@
               zq(ig,1,igcm_h2o_vap)=zq1temp(ig)
               
-c             Take into account H2O latent heat in surface energy budget
-c             We solve dT/dt = (2834.3-0.28*(T-To)-0.004*(T-To)^2)*1e3*iceflux/cpp
-              tsrf_lw(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
-             
-              alpha = exp(-4*abs(T1-T2)*pdqsdif(ig,igcm_h2o_ice)
-     &            *ptimestep/pcapcal(ig))
-
-              tsrf_lw(ig) = (tsrf_lw(ig)*(T2-alpha*T1)+T1*T2*(alpha-1))
-     &                      /(tsrf_lw(ig)*(1-alpha)+alpha*T2-T1)  ! surface temperature at t+1
-
-              pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
+!c             Take into account H2O latent heat in surface energy budget
+!c             We solve dT/dt = (2834.3-0.28*(T-To)-0.004*(T-To)^2)*1e3*iceflux/cpp
+!              tsrf_lw(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
+!             
+!              alpha = exp(-4*abs(T1-T2)*pdqsdif(ig,igcm_h2o_ice)
+!     &            *ptimestep/pcapcal(ig))
+!
+!              tsrf_lw(ig) = (tsrf_lw(ig)*(T2-alpha*T1)+T1*T2*(alpha-1))
+!     &                      /(tsrf_lw(ig)*(1-alpha)+alpha*T2-T1)  ! surface temperature at t+1
+!
+!              pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
 
                if(pqsurf(ig,igcm_h2o_ice)