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

Timestamp:

Dec 2, 2022, 7:59:07 PM (2 years ago)

Author:

jnaar

Message:

Changed units of sw_htrt and lw_htrt from W.m-2 to K.s-1 in writediagfi.
Created output diagfi var "surf_h2o_lh" which contains instantaneous latent heat flux from ground ice water condensation / sublimation in W.m-2
JN

File:

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

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

@@ -136 +136 @@
       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 +244 @@
       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 +963 @@
      &                       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 +983 @@
 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 +1057 @@
 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