Changeset 2293 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Apr 20, 2020, 6:03:41 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Fix for the 1D, follow up from the recent addition of watercap.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -11 +11 @@
      &                     albedice, iceradius, dtemisice, z0,
      &                     zmea, zstd, zsig, zgam, zthe, phisfi,
-     &                     watercaptag, hmons, summit, base
+     &                     watercaptag, watercap, hmons, summit, base
       use slope_mod, only: theta_sl, psi_sl
       use phyredem, only: physdem0,physdem1
@@ -472 +472 @@
       paleomars=.false. ! Default: no water ice reservoir
       call getin("paleomars",paleomars)
-      if (paleomars==.true.) then
+      if (paleomars.eqv..true.) then
         write(*,*) "paleomars=", paleomars
         write(*,*) "Orbital parameters from callphys.def"
@@ -736 +736 @@
 c Check if the surface is a water ice reservoir 
 c --------------------------------------------------
+      watercap(1)=0 ! Initialize watercap
       watercaptag(1)=.false. ! Default: no water ice reservoir
       print *,'Water ice cap on ground ?'
@@ -768 +769 @@
      &              totcloudfrac,wstar,
      &              mem_Mccn_co2,mem_Nccn_co2,
-     &              mem_Mh2o_co2)
+     &              mem_Mh2o_co2,watercap)
 
 c=======================================================================