Changeset 2840 for trunk/LMDZ.MARS

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (1 diff)
• libf/phymars/vdifc_mod.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -3830 +3830 @@
 Add an example of setup in LMDZ.MARS/util/compile to compile utilitary programs
 on Irene-Rome supercomputer.
+
+== 02/12/2022 == JN
+Add "surf_h2o_lh" output. This is the latent heat term from ground water ice
+sublimation / condensation in W.m-2 
+Also corrected units of shortwave & longwave heat rates in writediagfi

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

@@ -3025 +3025 @@
 c    &                  zstress)
         call WRITEDIAGFI(ngrid,'sw_htrt','sw heat. rate',
-     &                   'w.m-2',3,zdtsw)
+     &                   'K/s',3,zdtsw)
         call WRITEDIAGFI(ngrid,'lw_htrt','lw heat. rate',
-     &                   'w.m-2',3,zdtlw)
+     &                   'K/s',3,zdtlw)
         call writediagfi(ngrid,"local_time","Local time",
      &                   'sol',2,local_time)

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