Changeset 4531 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

May 12, 2023, 9:27:56 PM (21 months ago)

Author:

evignon

Message:

petite reecriture pour neige soufflee dans pbl_surface

File:

• LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90

@@ -2198 +2198 @@
     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 +2212 @@
      AcoefQBS, BcoefQBS)
     ENDIF
+