Changeset 2739

Timestamp:

Jun 30, 2022, 10:31:16 AM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Minor bug fix in co2condens, erroneous use of surface geopotential phisfi(),
wheras geoptential pphi() is already the geopotential relative to the surface.
Tests showed that this specific erroneous computation was only done in rare
cases and had overall little impact on the simulations.
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/co2condens_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3719 +3719 @@
 and top of atmosphere from physiq. Required adding some output arguments
 to callradite.
+
+== 30/06/2022 == EM
+Minor bug fix in co2condens, erroneous use of surface geopotential phisfi(),
+wheras geoptential pphi() is already the geopotential relative to the surface.
+Tests showed that this specific erroneous computation was only done in rare
+cases and had overall little impact on the simulations.

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

@@ -24 +24 @@
      &                      qperemin,masseqmin,
      &                      igcm_co2
-       use surfdat_h, only: emissiv, phisfi
+       use surfdat_h, only: emissiv
        use geometry_mod, only: latitude, ! grid point latitudes (rad)
      &                         longitude_deg, latitude_deg
@@ -400 +400 @@
 
             IF (zfallice(ig,1).gt.0 .and. .not. co2clouds) then
-              zfallheat = zfallice(ig,1) * (pphi(ig,1)- phisfi(ig) +
+              ! potential eneregy release due to the impact of the snowfall
+              ! appendix of forget et al. 1999
+              zfallheat = zfallice(ig,1) * (pphi(ig,1) +
      &                    cpice*(ztcond(ig,1)-ztcondsol(ig)))/latcond
             ELSE