Index: trunk/LMDZ.COMMON/libf/evolution/planet.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/planet.F90	(revision 4157)
+++ trunk/LMDZ.COMMON/libf/evolution/planet.F90	(revision 4170)
@@ -38,6 +38,6 @@
 
 ! Ice-related:
-real(dp),    dimension(:,:), allocatable :: h2o_ice                    ! H2O ice [kg.m-2]
-real(dp),    dimension(:,:), allocatable :: co2_ice                    ! CO2 ice [kg.m-2]
+real(dp),    dimension(:,:), allocatable :: h2o_ice                    ! H2O ice [kg/m2]
+real(dp),    dimension(:,:), allocatable :: co2_ice                    ! CO2 ice [kg/m2]
 real(dp)                                 :: h2oice_sublim_coverage_ini ! Initial surface area of sublimating H2O ice [m2]
 real(dp)                                 :: co2ice_sublim_coverage_ini ! Initial surface area of sublimating CO2 ice [m2]
@@ -73,4 +73,5 @@
 real(dp), dimension(:,:),   allocatable :: icetable_depth_old ! Old depth of the ice table [m]
 real(dp), dimension(:),     allocatable :: delta_icetable     ! Total mass of the H2O exchanged with the ice table [kg]
+real(dp), dimension(:,:),   allocatable :: flux_ssice_avg     ! Average of total flux exchanged with subsurface ice [kg/m2/y]
 
 ! Tracer-related:
@@ -153,4 +154,5 @@
 allocate(q_h2o_ts(ngrid,nday))
 allocate(q_co2_ts(ngrid,nday))
+allocate(flux_ssice_avg(ngrid,nslope))
 
 ps_avg(:) = 0._dp
@@ -163,4 +165,5 @@
 q_h2o_ts(:,:) = 0._dp
 q_co2_ts(:,:) = 0._dp
+flux_ssice_avg(:,:) = 0._dp
 
 END SUBROUTINE allocate_xios_state
@@ -422,4 +425,5 @@
 if (allocated(d_h2oice)) deallocate(d_h2oice)
 if (allocated(layerings_map)) deallocate(layerings_map)
+if (allocated(flux_ssice_avg)) deallocate(flux_ssice_avg)
 
 END SUBROUTINE end_planet