Changeset 4531

Timestamp:

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

Author:

evignon

Message:

petite reecriture pour neige soufflee dans pbl_surface

Location:

LMDZ6/trunk/libf

Files:

• phylmd/pbl_surface_mod.F90 (modified) (2 diffs)
• phylmdiso/pbl_surface_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -626 +626 @@
     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 +1905 @@
 !!!
 
- ! 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, & 

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
+