Changeset 2796

Timestamp:

Feb 8, 2018 12:25:39 PM (7 years ago)

Author:

suehring

Message:

Bugfix in 3D building initialization

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in 3D building initialization
+! 
+! 2747 2018-01-15 12:44:17Z suehring
 ! Bugfix, topography height is rounded to the nearest discrete grid level
 ! 
@@ -1058 +1061 @@
                    k2 = nzb+1
                    DO k = k_surf + 1, nzt + 1
-                      IF ( k2 <= buildings_f%nz )  THEN
+                      IF ( k2 <= buildings_f%nz-1 )  THEN
                          IF ( buildings_f%var_3d(k2,j,i) == 1 )  THEN
                             topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )