Changeset 437c545 in flexpart.git for src/verttransform_gfs.f90

Timestamp:

Mar 20, 2019, 3:32:30 PM (5 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug

Children:

7873bf7

Parents:

77783e3

Message:

BUGFIX: boundary overflow related to cloud height

File:

• src/verttransform_gfs.f90 (modified) (1 diff)

src/verttransform_gfs.f90

@@ -548 +548 @@
         if ((lsp.gt.0.01).or.(convp.gt.0.01)) then ! cloud and precipitation
 
-          do kz=nz,1,-1 !go Bottom up!
+          do kz=nz,2,-1 !go Bottom up!
             if (clw(ix,jy,kz,n).gt. 0) then ! is in cloud
               cloudsh(ix,jy,n)=cloudsh(ix,jy,n)+height(kz)-height(kz-1)
               clouds(ix,jy,kz,n)=1                               ! is a cloud
               if (lsp.ge.convp) then
-                clouds(ix,jy,kz,n)=3                            ! lsp in-cloud
+                clouds(ix,jy,kz,n)=3                             ! lsp in-cloud
               else
                 clouds(ix,jy,kz,n)=2                             ! convp in-cloud