Index: trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F	(revision 2285)
+++ trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F	(revision 2293)
@@ -11,5 +11,5 @@
      &                     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,5 +472,5 @@
       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,4 +736,5 @@
 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,5 +769,5 @@
      &              totcloudfrac,wstar,
      &              mem_Mccn_co2,mem_Nccn_co2,
-     &              mem_Mh2o_co2)
+     &              mem_Mh2o_co2,watercap)
 
 c=======================================================================