Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 3094)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 3118)
@@ -34,5 +34,5 @@
       use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
      & nsoilmx,tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var,
-     & flux_geo
+     & flux_geo,qsoil,nqsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
       use geometry_mod, only: longitude,latitude,cell_area
@@ -465,6 +465,6 @@
         fichnom = 'startfi.nc'
         CALL phyetat0 (fichnom,tab0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     &        day_ini,time,tsurf,tsoil,albedo,emis,
-     &        q2,qsurf,tauscaling,totcloudfrac,
+     &        nqsoil,day_ini,time,tsurf,tsoil,albedo,emis,
+     &        q2,qsurf,qsoil,tauscaling,totcloudfrac,
      &        wstar,watercap,perenial_co2ice,
      &        def_slope,def_slope_mean,subslope_dist)
@@ -2007,7 +2007,8 @@
      &              albfi,ithfi,def_slope,
      &              subslope_dist)
-      call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
+      call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,nqsoil,
      &              dtphys,hour_ini,
-     &              tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,
+     &              tsurf,tsoil,inertiesoil,albedo,
+     &              emis,q2,qsurf,qsoil,
      &              tauscaling,totcloudfrac,wstar,watercap,
      &              perenial_co2ice)
Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F	(revision 3094)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F	(revision 3118)
@@ -20,5 +20,6 @@
 
       use infotrac, only: infotrac_init, nqtot, tname
-      use comsoil_h, only: nsoilmx, inertiedat, inertiesoil
+      use comsoil_h, only: nsoilmx, inertiedat, inertiesoil,
+     &                     nqsoil, qsoil
       use surfdat_h, only: ini_surfdat_h, qsurf,watercaptag
       use comsoil_h, only: ini_comsoil_h
@@ -186,6 +187,6 @@
       allocate(subslope_distS(ip1jmp1,nslope))
 
-      CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
+      CALL phyetat0(fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,nqsoil,
+     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,qsoil,
      &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
      &      def_slope, def_slope_mean,subslope_dist)