Index: trunk/LMDZ.GENERIC/libf/phygeneric/physiq_mod.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phygeneric/physiq_mod.F90	(revision 4126)
+++ trunk/LMDZ.GENERIC/libf/phygeneric/physiq_mod.F90	(revision 4127)
@@ -620,4 +620,7 @@
 
          day_ini = pday
+	 DO nw=1,L_NSPECTV
+	    albedo(:,nw)=albedodat(:)
+	 ENDDO
 #endif
 
Index: trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk_rayleigh_scattering_opacity.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk_rayleigh_scattering_opacity.F90	(revision 4126)
+++ trunk/LMDZ.GENERIC/libf/phygeneric/rad_correlatedk_rayleigh_scattering_opacity.F90	(revision 4127)
@@ -48,6 +48,6 @@
       real, intent(out) :: tauray(L_LEVELS,L_NSPECTV)
       real*8 wl,wn
-      integer N,Nfine,ifine,igas,k
-      parameter(Nfine=500.0)
+      integer N,ifine,igas,k
+      integer, parameter :: Nfine = 500
       real*8 :: Fk ! King factor for the depolarization
       real*8 :: ng(L_LEVELS) ! real refractive index 
@@ -131,4 +131,5 @@
          endif
          tauvari(igas,:) = 0.
+         tauconsti(igas,:) = 0.
       enddo
       
@@ -149,8 +150,8 @@
          bwidth = (10000.0/BWNV(N)) - (10000.0/BWNV(N+1))
          do ifine=1,Nfine
-            wl=bstart+dble(ifine)*bwidth/Nfine
-            wn=BWNV(N)+dble(ifine)*(BWNV(N+1)-BWNV(N))/Nfine
-
-            tauvar(:)=0.0
+            wl=bstart+dble(ifine)*bwidth/dble(Nfine)
+            wn=BWNV(N)+dble(ifine)*(BWNV(N+1)-BWNV(N))/dble(Nfine)
+
+            tauvar(:)=0.0d0
             do igas=1,ngasmx
                if (maxval(mass_frac(igas,:)).ge.1e-2) then
@@ -319,7 +320,6 @@
                  ! pmid*scalep -> mbar to Pa
                  ! muvar/1000 -> g/mol to kg/mol
-               
                  tauvar(:)=tauvar(:)+tauconsti(igas,:)*tauvari(igas,:)
-                 
+          
                endif !greater than 0.01
 
@@ -327,5 +327,5 @@
 
             call rad_blackbody_planck_law_wavelength(dble(wl*1e-6),dble(tstellar),df)
-            df=df*bwidth/Nfine
+            df=df*bwidth/dble(Nfine)
             tauwei=tauwei+df
             tausum(:)=tausum(:)+tauvar(:)*df