Changeset 2179 for trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F

Timestamp:

Nov 14, 2019, 6:03:49 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Add the sublimation/condensation latent heat release from the surface water ice in vdifc_mod.F. It can be activated or deactivated with the new flag "latentheat" (for now latentheat=false by default).
MV

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (2 diffs)

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

@@ -122 +122 @@
       REAL,SAVE :: acond,bcond
       
-c     For latent heat release from ground ice sublimation    
+c     For latent heat release from ground water ice sublimation    
+      REAL tsrf_lh(ngrid) ! temporary surface temperature
+      REAL lh ! latent heat, formulation given in the Technical Document: 
+              ! "Modeling water ice sublimation under Phoenix-like conditions", Montmessin et al. 2004  
 !      REAL tsrf_lw(ngrid)
 !      REAL alpha
-      REAL T1,T2
-      SAVE T1,T2
-      DATA T1,T2/-604.3,1080.7/ ! zeros of latent heat equation for ice
+!      REAL T1,T2
+!      SAVE T1,T2
+!      DATA T1,T2/-604.3,1080.7/ ! zeros of latent heat equation for ice
 
 c    Tracers :
@@ -863 +866 @@
 !
 !              pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
+c             Take into account the H2O latent heat impact on the surface temperature
+              if (latentheat) then
+                tsrf_lh(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
+                lh=(2834.3-0.28*(tsrf_lh(ig)-To)-
+     &              0.004*(tsrf_lh(ig)-To)*(tsrf_lh(ig)-To))*1.e+3
+                pdtsrf(ig)= pdtsrf(ig) + pdqsdif(ig,igcm_h2o_ice)*lh
+     &                    /pcapcal(ig)
+              endif
 
                if(pqsurf(ig,igcm_h2o_ice)