Index: /trunk/LMDZ.COMMON/libf/evolution/ice_table_mod.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/ice_table_mod.F90	(revision 3031)
+++ /trunk/LMDZ.COMMON/libf/evolution/ice_table_mod.F90	(revision 3032)
@@ -149,4 +149,5 @@
     use comslope_mod, only: subslope_dist,def_slope_mean           
     use comconst_mod,only: pi       
+    use constants_marspem_mod,only:porosity
     implicit none 
 !   inputs
@@ -201,5 +202,5 @@
    do ig = 1,ngrid
       do islope = 1,nslope
-        delta_m_h2o(ig)  = delta_m_h2o(ig) +  rho(ig,islope)*(new_ice_table_thickness(ig,islope) - former_ice_table_thickness(ig,islope)) & ! convention > 0. <=> it condenses  
+        delta_m_h2o(ig)  = delta_m_h2o(ig) +  porosity*rho(ig,islope)*(new_ice_table_thickness(ig,islope) - former_ice_table_thickness(ig,islope)) & ! convention > 0. <=> it condenses  
                            *subslope_dist(ig,islope)/cos(def_slope_mean(islope)*pi/180.)
       enddo
Index: /trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3031)
+++ /trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3032)
@@ -65,5 +65,6 @@
                                   albedodat, zmea, zstd, zsig, zgam, zthe,     &
                                   hmons, summit, base,albedo_h2o_frost,        &
-                                  frost_albedo_threshold, emissiv, watercaptag, perenial_co2ice
+                                  frost_albedo_threshold, emissiv, watercaptag, perenial_co2ice, &
+                                  emisice, albedice
     use dimradmars_mod,     only: totcloudfrac, albedo
     use dust_param_mod,     only: tauscaling
@@ -72,4 +73,5 @@
     use planete_h,          only: aphelie, periheli, year_day, peri_day, obliquit
     use comcstfi_h,         only: r, mugaz
+    use paleoclimate_mod,   only: albedo_perenialco2
 #else
     use tracer_h,           only: noms, igcm_h2o_ice, igcm_co2 ! Tracer names
@@ -245,5 +247,5 @@
 
 ! Loop variables
-integer :: i, j, ig0, l, ig, nnq, t, islope, ig_loop, islope_loop, iloop, isoil
+integer :: i, j, ig0, l, ig, nnq, t, islope, ig_loop, islope_loop, iloop, isoil,icap
 
 ! Parallel variables
@@ -843,4 +845,13 @@
                 enddo
                 tsurf_ave(ig,islope) = ave/timelen
+                ! set the surface albedo to be the ice albedo
+                if (latitude_deg(ig) > 0) then
+                   icap = 1
+                else
+                   icap = 2
+                endif
+                albedo(ig,1,islope) = albedice(icap)
+                albedo(ig,2,islope) = albedice(icap)
+                emis(ig,islope) = emisice(icap)
             endif
         enddo
@@ -1042,4 +1053,6 @@
             perenial_co2ice(ig,islope) = 0.5*qsurf(ig,igcm_co2,islope)
             qsurf(ig,igcm_co2,islope)  = 0.5*qsurf(ig,igcm_co2,islope)
+            albedo(ig,1,islope) = albedo_perenialco2
+            albedo(ig,2,islope) = albedo_perenialco2
         endif
     enddo