Changeset 742 for trunk/LMDZ.MARS/libf


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/vdifc.F

    r734 r742  
    812812              zq(ig,1,igcm_h2o_vap)=zq1temp(ig)
    813813             
    814 c             Take into account H2O latent heat in surface energy budget
    815 c             We solve dT/dt = (2834.3-0.28*(T-To)-0.004*(T-To)^2)*1e3*iceflux/cpp
    816               tsrf_lw(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
    817              
    818               alpha = exp(-4*abs(T1-T2)*pdqsdif(ig,igcm_h2o_ice)
    819      &            *ptimestep/pcapcal(ig))
    820 
    821               tsrf_lw(ig) = (tsrf_lw(ig)*(T2-alpha*T1)+T1*T2*(alpha-1))
    822      &                      /(tsrf_lw(ig)*(1-alpha)+alpha*T2-T1)  ! surface temperature at t+1
    823 
    824               pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
     814!c             Take into account H2O latent heat in surface energy budget
     815!c             We solve dT/dt = (2834.3-0.28*(T-To)-0.004*(T-To)^2)*1e3*iceflux/cpp
     816!              tsrf_lw(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
     817!             
     818!              alpha = exp(-4*abs(T1-T2)*pdqsdif(ig,igcm_h2o_ice)
     819!     &            *ptimestep/pcapcal(ig))
     820!
     821!              tsrf_lw(ig) = (tsrf_lw(ig)*(T2-alpha*T1)+T1*T2*(alpha-1))
     822!     &                      /(tsrf_lw(ig)*(1-alpha)+alpha*T2-T1)  ! surface temperature at t+1
     823!
     824!              pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
    825825
    826826               if(pqsurf(ig,igcm_h2o_ice)
Note: See TracChangeset for help on using the changeset viewer.