Changeset 2218

Timestamp:

Jan 16, 2020, 10:44:19 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Change "latentheat" flag to a more descriptive "latentheat_surfwater" and
set its default value to .true.
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/callkeys.h (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/vdifc_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -2825 +2825 @@
 Bug fix in newstart "composition" option; added option "q=factor" in
 newstart to multiply a tracer mixing ratio by a constant factor.
+
+== 16/01/2020 == EM
+Change "latentheat" flag to a more descriptive "latentheat_surfwater" and
+set its default value to .true.

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -15 +15 @@
      &   ,calltherm,callrichsl,callslope,tituscap,callyamada4,co2clouds &
      &   ,co2useh2o,meteo_flux,CLFvaryingCO2,spantCO2,CLFvarying        &
-     &   ,satindexco2,rdstorm,slpwind,calllott_nonoro,latentheat
+     &   ,satindexco2,rdstorm,slpwind,calllott_nonoro                   &
+     &   ,latentheat_surfwater
      
       COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
@@ -67 +68 @@
       logical rdstorm ! rocket dust storm parametrization
       logical slpwind ! entrainment by slope wind parametrization
-      logical latentheat ! latent heat release from ground water ice sublimation/condensation
+      logical latentheat_surfwater ! latent heat release from ground water ice sublimation/condensation
       logical sedimentation
       logical activice,tifeedback,supersat,caps

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

@@ -307 +307 @@
          write(*,*)"latent heat release during sublimation", 
      &              " /condensation of ground water ice"
-         latentheat=.false. ! default value
-         call getin("latentheat",latentheat)
-         write(*,*)" latentheat = ",latentheat
+         latentheat_surfwater=.true. ! default value
+         call getin("latentheat_surfwater",latentheat_surfwater)
+         write(*,*)" latentheat_surfwater = ",latentheat_surfwater
 
          write(*,*)"rad.transfer is computed every iradia",

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

@@ -867 +867 @@
 !              pdtsrf(ig) = (tsrf_lw(ig)-ptsrf(ig))/ptimestep
 c             Take into account the H2O latent heat impact on the surface temperature
-              if (latentheat) then
+              if (latentheat_surfwater) then
                 tsrf_lh(ig) = ptsrf(ig) + pdtsrf(ig) *ptimestep
                 lh=(2834.3-0.28*(tsrf_lh(ig)-To)-