Index: trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90	(revision 3974)
+++ trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90	(revision 3996)
@@ -431,5 +431,5 @@
 call ini_fillgeom(1,latitude,longitude,(/1.0/))
 call conf_phys(1,llm,nq)
-call initracer(ngrid,nq,qsurf)
+call initracer(ngrid,nq)
 
 ! In 1D model physics are called every time step
Index: trunk/LMDZ.MARS/libf/phymars/initracer.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/initracer.F90	(revision 3974)
+++ trunk/LMDZ.MARS/libf/phymars/initracer.F90	(revision 3996)
@@ -5,5 +5,5 @@
 contains
 
-SUBROUTINE initracer(ngrid,nq,qsurf)
+SUBROUTINE initracer(ngrid,nq)
 
 use tracer_mod
@@ -40,5 +40,4 @@
 integer, intent(in) :: ngrid ! number of atmospheric columns
 integer, intent(in) :: nq    ! number of tracers
-real, intent(out) :: qsurf(ngrid,nq) ! tracer on surface (e.g.  kg.m-2)
 
 integer :: iq,ig,count
@@ -548,19 +547,4 @@
 endif
 
-! if water cycle but iceparty=.false., there will nevertheless be
-! water ice at the surface (iceparty is not used anymore, but this
-! part is still relevant, as we want to stay compatible with the
-! older versions).
-if (water.and.(igcm_h2o_ice == 0)) then
-    igcm_h2o_ice=igcm_h2o_vap ! so that qsurf(i_h2o_ice) is identified
-                              ! even though there is no q(i_h2o_ice)
-else
-    ! surface ice qsurf(i_h2o_ice) was loaded twice by phyetat0,
-    ! as qsurf(i_h2o_vap) & as qsurf(i_h2o_ice), so to be clean:
-    if (igcm_h2o_vap /= 0) then
-        qsurf(1:ngrid,igcm_h2o_vap)=0
-    endif
-endif
-
 ! Additional test required for HDO
 ! We need to compute some things for H2O before HDO
Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 3974)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 3996)
@@ -711,5 +711,5 @@
 c        initialize tracers
 c        ~~~~~~~~~~~~~~~~~~
-         call initracer(ngrid,nq,qsurf)
+         call initracer(ngrid,nq)
 
 c        Initialize albedo and orbital calculation
Index: trunk/LMDZ.MARS/libf/phymars/planete_h.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/planete_h.F90	(revision 3974)
+++ trunk/LMDZ.MARS/libf/phymars/planete_h.F90	(revision 3996)
@@ -17,5 +17,5 @@
 REAL,SAVE :: lsperi    ! Solar longitude of the perihelion, angle in rad
 REAL,SAVE :: e_elips   ! Orbit eccentricity
-REAL,SAVE :: p_elips   ! Ellipse semi-latus rectum
+REAL,SAVE :: p_elips   ! Ellipse parameter (semi-latus rectum)
 !$OMP THREADPRIVATE(lsperi,e_elips,p_elips)
 REAL,PARAMETER :: unitastr=149.597927 ! Astronomical unit AU, in Mkm
Index: trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90	(revision 3974)
+++ trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90	(revision 3996)
@@ -62,7 +62,7 @@
 
   implicit none
-  integer,intent(in) :: ngrid ! number of atmospheric columns
-  integer,intent(in) :: nq ! number of tracers
-  integer,intent(in) :: nslope ! number of sub-grid scale slope
+  integer,intent(in) :: ngrid  ! number of atmospheric columns
+  integer,intent(in) :: nq     ! number of tracers
+  integer,intent(in) :: nslope ! number of sub-grid scale slopes
 
     allocate(albedodat(ngrid))