Index: trunk/LMDZ.MARS/libf/phymars/albedocaps.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/albedocaps.F90	(revision 3338)
+++ trunk/LMDZ.MARS/libf/phymars/albedocaps.F90	(revision 3343)
@@ -11,5 +11,5 @@
 USE mod_phys_lmdz_transfert_para, ONLY: bcast
 USE mod_phys_lmdz_para, ONLY: is_master
-USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2
+USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2,albedo_co2_cap
 
 implicit none
@@ -105,7 +105,15 @@
       psolaralb(ig,2)=psolaralb(ig,1)
     else
-      psolaralb(ig,:)=albedice(icap)
-      if (paleoclimate .and. piceco2_peren(ig) > 0. .and. abs(piceco2(ig)) < 1.e-10) psolaralb(ig,:) = albedo_perennialco2
-    endif
+      if(paleoclimate) then
+            psolaralb(ig,1) = albedo_co2_cap
+            psolaralb(ig,2) = albedo_co2_cap
+      else
+            psolaralb(ig,1) = albedice(icap)
+            psolaralb(ig,2) = albedice(icap)
+      endif
+    endif
+  else if(paleoclimate .and. piceco2_peren(ig) > 0.) then
+     psolaralb(ig,1) = albedo_perennialco2
+     psolaralb(ig,2) = albedo_perennialco2
   else if (watercaptag(ig) .and. water) then
     ! there is a water ice cap: set the surface albedo to the water ice one
Index: trunk/LMDZ.MARS/libf/phymars/conf_phys.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/conf_phys.F	(revision 3338)
+++ trunk/LMDZ.MARS/libf/phymars/conf_phys.F	(revision 3343)
@@ -38,4 +38,5 @@
       USE mod_phys_lmdz_transfert_para, ONLY: bcast
       USE paleoclimate_mod,ONLY: paleoclimate,albedo_perennialco2,
+     &                           albedo_co2_cap,
      &                           lag_layer, include_waterbuoyancy 
       use microphys_h, only: mteta
@@ -375,4 +376,9 @@
          call getin_p("albedo_perennialco2",albedo_perennialco2)
          write(*,*)"albedo_perennialco2 = ",albedo_perennialco2
+
+         write(*,*)"Albedo for seasonal CO2 ice?"
+         albedo_co2_cap = 0.6 ! default value
+         call getin_p("albedo_co2_cap",albedo_co2_cap)
+         write(*,*)"albedo_co2_cap = ",albedo_co2_cap
 
          write(*,*)"Include water buoyancy effect??"
Index: trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90	(revision 3338)
+++ trunk/LMDZ.MARS/libf/phymars/paleoclimate_mod.F90	(revision 3343)
@@ -21,5 +21,6 @@
     logical, save                              :: lag_layer             ! Does lag layer is present?
     logical, save                              :: include_waterbuoyancy ! Include the effect of water buoyancy when computing the sublimation of water ice ?
-!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,albedo_perennialco2,lag_layer,include_waterbuoyancy)
+    real,    save                              :: albedo_co2_cap
+!$OMP THREADPRIVATE(h2o_ice_depth,lag_co2_ice,d_coef,albedo_perennialco2,lag_layer,include_waterbuoyancy,albedo_co2_cap)
 
 !=======================================================================
Index: trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 3338)
+++ trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90	(revision 3343)
@@ -842,5 +842,5 @@
         if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
             perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
-            qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
+            qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) 
         endif
     endif ! not found
@@ -852,5 +852,5 @@
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
         perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
-        qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
+        qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:)
     endif
   endif !startphy_file