Index: trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3197)
+++ trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3198)
@@ -920,5 +920,8 @@
                print*,'and the surface albedo is taken equal to the first visible spectral value'
 
+               albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
                fluxsurfabs_sw(1:ngrid) = fluxtop_dn(1:ngrid)*(1.-albedo(1:ngrid,1))
+               ! TB24:
+               fluxabs_sw(1:ngrid)=fluxsurfabs_sw(1:ngrid)
                fluxrad_sky(1:ngrid)    = fluxsurfabs_sw(1:ngrid)
                fluxtop_lw(1:ngrid)  = emis(1:ngrid)*sigma*tsurf(1:ngrid)**4
@@ -1662,5 +1665,5 @@
       if(callrad)then
 
-         !call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
+         call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
          !call writediagfi(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
          call writediagfi(ngrid,"ISR","incoming stellar rad.","W m-2",2,fluxtop_dn)
Index: trunk/LMDZ.PLUTO/libf/phypluto/surfini.F
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/surfini.F	(revision 3197)
+++ trunk/LMDZ.PLUTO/libf/phypluto/surfini.F	(revision 3198)
@@ -44,7 +44,7 @@
       ! Step 2 : We get the bare ground albedo from the start files.
       DO ig=1,ngrid
-         albedo_bareground(ig)=albedodat(ig)
+         albedo_bareground(ig)=0.1 ! TB24 albedodat(ig)
 	 DO nw=1,L_NSPECTV
-	    albedo(ig,nw)=albedo_bareground(ig)
+	    albedo(ig,nw)=0.1 !albedo_bareground(ig)
 	 ENDDO
       ENDDO
Index: trunk/LMDZ.PLUTO/libf/phypluto/turbdiff_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/turbdiff_mod.F90	(revision 3197)
+++ trunk/LMDZ.PLUTO/libf/phypluto/turbdiff_mod.F90	(revision 3198)
@@ -138,4 +138,8 @@
 
       IF (firstcall) THEN
+         ivap=1 !TB24
+         iliq=0
+         iliq_surf=0
+	 iice_surf=0 ! simply to make the code legible               	  
          sensibFlux(:)=0.
          firstcall=.false.