Index: LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 4530)
+++ LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 4531)
@@ -626,5 +626,5 @@
     REAL, DIMENSION(klon,klev)         :: y_flux_t, y_flux_q, y_flux_qbs
     REAL, DIMENSION(klon,klev)         :: y_flux_u, y_flux_v
-    REAL, DIMENSION(klon,klev)         :: ycoefh, ycoefm,ycoefq, ycoefqbs
+    REAL, DIMENSION(klon,klev)         :: ycoefh,ycoefm,ycoefq,ycoefqbs
     REAL, DIMENSION(klon)              :: ycdragh, ycdragq, ycdragm
     REAL, DIMENSION(klon,klev)         :: yu, yv
@@ -1905,11 +1905,15 @@
 !!!
 
- ! For blowing snow:
+! For blowing snow:
     IF (ok_bs) THEN
      ! following Bintanja et al 2000, part II
-     ! we assume that the eddy diffuvisity coefficient for
+     ! we assume that the eddy diffsivity coefficient for
      ! suspended particles is larger than Km by a factor zeta_bs
      ! which is equal to 3 by default
-     ycoefqbs=ycoefm*zeta_bs 
+     do k=1,klev
+        do j=1,knon
+           ycoefqbs(j,k)=ycoefm(j,k)*zeta_bs 
+        enddo
+     enddo
      CALL climb_qbs_down(knon, ycoefqbs, ypaprs, ypplay, &
      ydelp, yt, yqbs, dtime, & 
Index: LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90	(revision 4530)
+++ LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90	(revision 4531)
@@ -2198,8 +2198,12 @@
     IF (ok_bs) THEN
      ! following Bintanja et al 2000, part II
-     ! we assume that the eddy diffuvisity coefficient for
+     ! we assume that the eddy diffsivity coefficient for
      ! suspended particles is larger than Km by a factor zeta_bs
      ! which is equal to 3 by default
-     ycoefqbs=ycoefm*zeta_bs 
+     do k=1,klev
+        do j=1,knon
+           ycoefqbs(j,k)=ycoefm(j,k)*zeta_bs 
+        enddo
+     enddo
      CALL climb_qbs_down(knon, ycoefqbs, ypaprs, ypplay, &
      ydelp, yt, yqbs, dtime, & 
@@ -2208,4 +2212,5 @@
      AcoefQBS, BcoefQBS)
     ENDIF
+