Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 2951)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F	(revision 2952)
@@ -32,5 +32,6 @@
      &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var
       use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
-     & nsoilmx, tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var
+     & nsoilmx,tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var,
+     & flux_geo
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
       use geometry_mod, only: longitude,latitude,cell_area
@@ -196,5 +197,7 @@
       REAL,ALLOCATABLE :: watercap_old_slope(:,:) ! Surface water ice (kg.m-2)
       REAL,ALLOCATABLE :: tsoil_old_slope(:,:,:) 
+      REAL,ALLOCATABLE :: inertiesoil_old_slope(:,:,:) 
       REAL,ALLOCATABLE :: albedo_old_slope(:,:,:) ! Surface albedo in each solar band
+      REAL,ALLOCATABLE :: flux_geo_old_slope(:,:)
       integer :: iflat
       integer :: nslope_old, nslope_new
@@ -1798,4 +1801,5 @@
         allocate(watercap_old_slope(ngridmx,nslope_old)) 
 
+
         tsurf_old_slope(:,:)=tsurf(:,:)
         qsurf_old_slope(:,:,:)=qsurf(:,:,:)
@@ -1809,6 +1813,10 @@
 
         allocate(tsoil_old_slope(ngridmx,nsoilmx,nslope_old))
-
+        allocate(inertiesoil_old_slope(ngridmx,nsoilmx,nslope_old))
+        allocate(flux_geo_old_slope(ngridmx,nslope_old))
+
+        inertiesoil_old_slope(:,:,:)=inertiesoil(:,:,:)
         tsoil_old_slope(:,:,:)=tsoil(:,:,:)
+        flux_geo_old_slope(:,:)=flux_geo(:,:)
 
         call end_comsoil_h_slope_var
@@ -1831,4 +1839,6 @@
              tsoil(:,:,islope)=tsoil_old_slope(:,:,1)
              albedo(:,:,islope)=albedo_old_slope(:,:,1)
+             inertiesoil(:,:,islope)=inertiesoil_old_slope(:,:,1)
+             flux_geo(:,islope)=flux_geo_old_slope(:,1)
           enddo
         elseif(nslope_new.eq.1) then
@@ -1839,4 +1849,6 @@
              tsoil(:,:,1)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,1)=albedo_old_slope(:,:,iflat)
+             inertiesoil(:,:,1)=inertiesoil_old_slope(:,:,iflat)
+             flux_geo(:,1)=flux_geo_old_slope(:,iflat)
         elseif(nslope_old.eq.5 .and. nslope_new.eq.7) then
           do islope=1,nslope_new
@@ -1847,4 +1859,6 @@
              tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
+             inertiesoil(:,:,islope)=inertiesoil_old_slope(:,:,iflat)
+             flux_geo(:,islope)=flux_geo_old_slope(:,iflat)
           enddo
         elseif(nslope_old.eq.7 .and. nslope_new.eq.5) then
@@ -1856,4 +1870,6 @@
              tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
+             inertiesoil(:,:,islope)=inertiesoil_old_slope(:,:,iflat)
+             flux_geo(:,islope)=flux_geo_old_slope(:,iflat)
           enddo
         else
Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F	(revision 2951)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F	(revision 2952)
@@ -126,5 +126,5 @@
 	ierr=NF_INQ_VARID(nid,nom,varid)
 
-        if(nom.eq."tsoil") then
+        if(nom.eq."tsoil" .or. nom.eq."inertiesoil") then
 
 	 if (ierr.ne.NF_NOERR) then ! variable not defined yet
Index: trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90	(revision 2951)
+++ trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90	(revision 2952)
@@ -84,4 +84,5 @@
     allocate(alph(ngrid,nsoilmx-1,nslope))
     allocate(beta(ngrid,nsoilmx-1,nslope))
+    allocate(flux_geo(ngrid,nslope))
  
   end subroutine ini_comsoil_h_slope_var
@@ -99,4 +100,5 @@
     if (allocated(alph)) deallocate(alph)
     if (allocated(beta)) deallocate(beta)
+    if (allocated(flux_geo)) deallocate(flux_geo)
 
   end subroutine end_comsoil_h_slope_var
Index: trunk/LMDZ.MARS/libf/phymars/iostart.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/iostart.F90	(revision 2951)
+++ trunk/LMDZ.MARS/libf/phymars/iostart.F90	(revision 2952)
@@ -836,5 +836,5 @@
             ierr=NF90_REDEF(nid_restart)
 #ifdef NC_DOUBLE
-            if(field_name.eq. "tsoil") then
+            if(field_name.eq. "tsoil".or. field_name.eq. "inertiesoil") then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
                               (/idim2,idim3,idim8,idim7/),nvarid)
@@ -844,5 +844,5 @@
             endif
 #else
-            if(field_name.eq. "tsoil") then
+            if(field_name.eq. "tsoil".or. field_name.eq. "inertiesoil") then
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
                               (/idim2,idim3,idim8,idim7/),nvarid)
@@ -1015,5 +1015,4 @@
           IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
           ierr=NF90_ENDDEF(nid_restart)
-          
           firsttime=.false.
           if(first_time1d) first_time1d=.false.