Index: trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F	(revision 2900)
+++ trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F	(revision 2902)
@@ -170,5 +170,6 @@
       REAL masseq(nlayer),wq(nlayer+1)
       INTEGER ifils,iq2
-
+c LL: flux dependant albedo
+      real :: totflux(ngrid)
 c----------------------------------------------------------------------
 
@@ -176,4 +177,5 @@
 c   --------------
 c
+       totflux(:) = fluxsurf_sw(:,1)+fluxsurf_sw(:,2)
       ! AS: firstcall OK absolute
       IF (firstcall) THEN
@@ -500,4 +502,16 @@
 !     ----------------------------------------
       CALL albedocaps(zls,ngrid,piceco2,psolaralb,emisref)
+      DO ig = 1,ngrid
+        if(piceco2(ig).gt.0.) then
+          psolaralb(ig,1) = 0.267+0.00059*totflux(ig)
+          psolaralb(ig,2) = psolaralb(ig,1)
+          if(ig.eq.ngrid) then
+           psolaralb(ig,1) = 0.532+8.72e-4*totflux(ig)
+           psolaralb(ig,2) = psolaralb(ig,1)
+          endif
+        endif
+      ENDDO
+
+
 
 ! set pemisurf() to emissiv when there is bare surface (needed for co2snow)
@@ -739,8 +753,8 @@
 !             write(*,*) "co2condens: South pole: latitude(ngrid)=",
 !     &                                           latitude(ngrid)
-             ztcondsol(ngrid)=
-     &          1./(bcond-acond*log(.01*vmr_co2(ngrid,1)*
-     &                    (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
-             pdtsrfc(ngrid)=(ztcondsol(ngrid)-ztsrf(ngrid))/ptimestep
+!             ztcondsol(ngrid)=
+!     &          1./(bcond-acond*log(.01*vmr_co2(ngrid,1)*
+!     &                    (pplev(ngrid,1)+pdpsrf(ngrid)*ptimestep)))
+!             pdtsrfc(ngrid)=(ztcondsol(ngrid)-ztsrf(ngrid))/ptimestep
            endif
          endif