Index: LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/pbl_surface_mod.F90	(revision 1941)
+++ LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/pbl_surface_mod.F90	(revision 1942)
@@ -443,4 +443,5 @@
     REAL, DIMENSION(klon,klev+1)       :: ytke
     REAL, DIMENSION(klon,nsoilmx)      :: ytsoil
+    REAL, DIMENSION(klon)              :: yrmu0
     CHARACTER(len=80)                  :: abort_message
     CHARACTER(len=20)                  :: modname = 'pbl_surface'
@@ -554,4 +555,5 @@
     d_t_diss= 0.0 ;d_u = 0.0     ; d_v = 0.0        ; yqsol = 0.0    
     ytherm = 0.0  ; ytke=0.
+    yrmu0 = 0.0
     
     tke(:,:,is_ave)=0.
@@ -710,4 +712,5 @@
           ypaprs(j,klev+1) = paprs(i,klev+1)
           ywindsp(j) = SQRT(u10m(i,nsrf)**2 + v10m(i,nsrf)**2 )
+          yrmu0(j)=rmu0(i)
        END DO
 
@@ -856,5 +859,5 @@
                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
                yqsurf, ytsurf_new, y_dflux_t, y_dflux_q, &
-               y_flux_u1, y_flux_v1 )
+               y_flux_u1, y_flux_v1, yrmu0 )
                
      
Index: LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_mod.F90
===================================================================
--- LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_mod.F90	(revision 1941)
+++ LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_mod.F90	(revision 1942)
@@ -19,5 +19,5 @@
        z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
        qsurf, tsurf_new, dflux_s, dflux_l, &
-       flux_u1, flux_v1 ) 
+       flux_u1, flux_v1, rmu0) 
 
     USE dimphy
@@ -61,4 +61,5 @@
                                                          ! corresponds to previous sollwdown
     REAL, DIMENSION(klon), INTENT(IN)       :: q2m, t2m
+    REAL, DIMENSION(klon), INTENT(IN)       :: rmu0      ! cosine of solar zenith angle
 
 ! In/Output variables
@@ -132,5 +133,5 @@
             evap, fluxsens, fluxlat, &              
             tsol_rad, tsurf_new, alb1_new, alb2_new, &
-            emis_new, z0_new, qsurf)       
+            emis_new, z0_new, qsurf, rmu0)       
 
 !  
Index: LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_orchidee_mod.F90
===================================================================
--- LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_orchidee_mod.F90	(revision 1941)
+++ LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/surf_land_orchidee_mod.F90	(revision 1942)
@@ -39,5 +39,5 @@
        evap, fluxsens, fluxlat, &              
        tsol_rad, tsurf_new, alb1_new, alb2_new, &
-       emis_new, z0_new, qsurf)
+       emis_new, z0_new, qsurf, rmu0)
 
     USE mod_surf_para
@@ -124,4 +124,5 @@
     REAL, DIMENSION(klon), INTENT(IN)         :: lwdown, swnet, swdown, ps
     REAL, DIMENSION(klon), INTENT(IN)         :: q2m, t2m
+    REAL, DIMENSION(klon), INTENT(IN)         :: rmu0
 
 ! Parametres de sortie
@@ -198,4 +199,7 @@
     INTEGER :: orch_omp_rank
     INTEGER :: orch_omp_size
+
+    REAL, DIMENSION(knon)                     :: sinang
+
 !
 ! Fin definition
@@ -402,5 +406,5 @@
                evap, fluxsens, fluxlat, coastalflow, riverflow, &
                tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
-               lon_scat, lat_scat, q2m, t2m)
+               lon_scat, lat_scat, q2m, t2m, sinang=rmu0)
 #endif         
        ENDIF
@@ -426,5 +430,5 @@
             evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
             tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), &
-            lon_scat, lat_scat, q2m, t2m)
+            lon_scat, lat_scat, q2m, t2m, sinang=rmu0(1:knon))
 #endif       
     ENDIF