Changeset 742

Timestamp:

Jul 25, 2012, 6:01:21 PM (12 years ago)

Author:

tnavarro

Message:

latent heat for sublimating ice turned off: at high obliquity it gives horrible, frightful, nightmarish ground temperatures

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/vdifc.F

@@ -812 +812 @@
               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)