Index: trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 4169)
+++ trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 4170)
@@ -37,5 +37,5 @@
 use comsoil_h,           only: flux_geo
 use comslope_mod,        only: nslope, major_slope
-use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif, zdqsdif_ssi_tot
+use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif
 use comcstfi_h,          only: pi
 use geometry_mod,        only: latitude
@@ -914,12 +914,4 @@
      h2oice_depth(:,:) = -999.
    endif
-   
-! Total flux with SSI
-   call get_field("flux_ssice",zdqsdif_ssi_tot,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <flux_ssice> : ", &
-                          "<flux_ssice> is set as 0 (no subsurface ice)"
-     zdqsdif_ssi_tot(:,:) = 0.
-   endif
 
 ! Diffusion coeficent
@@ -950,5 +942,4 @@
     h2oice_depth(:,:) = -999.
     co2ice_depth(:,:) = -999.
-    zdqsdif_ssi_tot(:,:) = 0.
     coef_ssdif(:,:) = 4.e-4
     perennial_co2ice(:,:) = 0.
@@ -958,5 +949,4 @@
    h2oice_depth(:,:) = -999.
    co2ice_depth(:,:) = -999.
-   zdqsdif_ssi_tot(:,:) = 0.
    coef_ssdif(:,:) = 4.e-4
    perennial_co2ice(:,:) = 0.
Index: trunk/LMDZ.MARS/libf/phymars/phyredem.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/phyredem.F90	(revision 4169)
+++ trunk/LMDZ.MARS/libf/phymars/phyredem.F90	(revision 4170)
@@ -210,5 +210,5 @@
   use callkeys_mod,        only: calltherm, dustinjection, calllott_nonoro
   use callkeys_mod, only: CLFvarying
-  use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif, zdqsdif_ssi_tot
+  use paleoclimate_mod,    only: paleoclimate, h2oice_depth, co2ice_depth, coef_ssdif
 
   implicit none
@@ -390,5 +390,4 @@
     call put_field("h2oice_depth","Depth of the shallowest H2O ice layer",h2oice_depth)
     call put_field("co2ice_depth","Depth of the shallowest CO2 ice layer",co2ice_depth)
-    call put_field("flux_ssice","Total flux exchanged with subsurface water ice",zdqsdif_ssi_tot)
     call put_field("coef_ssdif","Diffusion coefficent for subsurface water",coef_ssdif)
   endif