Changeset 4314 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Nov 29, 2019 10:29:20 AM (4 years ago)

Author:

suehring

Message:

plant-canopy: Bugfix, plant canopy was still considered at building edges on for the u- and v-component; Relax restriction of LAD on building tops. LAD is only omitted at locations where building grid points emerged artificially by the topography filtering

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set additional topography flag 4 to mark topography grid points emerged
+! from the filtering process.
+! 
+! 4294 2019-11-13 18:34:16Z suehring
 ! Bugfix, always set bit 5 and 6 of wall_flags, indicating terrain- and 
 ! building surfaces in all  cases, in order to enable terrain-following output 
@@ -1600 +1604 @@
                       num_hole_l     = num_hole_l + 1
 !
-!--                   Clear flag 0 and set special flag ( bit 3) to indicate 
+!--                   Clear flag 0 and set special flag ( bit 4) to indicate 
 !--                   that new topography point is a result of filtering process.
                       topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
-                      topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 3 )
+                      topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 4 )
 !
 !--                   If filled grid point is occupied by a building, classify
@@ -2180 +2184 @@
           CALL exchange_horiz_int( topo, nys, nyn, nxl, nxr, nzt, nbgp )
 !
-!--       Set lateral boundary conditions for topography on all ghost layers          
+!--       Set lateral boundary conditions for topography on all ghost layers
           IF ( .NOT. bc_ns_cyc )  THEN
              IF ( nys == 0  )  THEN
@@ -2601 +2605 @@
     ENDDO
 !
+!-- Set flag 4, indicating new topography grid points due to filtering.
+    DO i = nxl, nxr
+       DO j = nys, nyn
+          DO k = nzb, nzt+1
+             IF ( BTEST( topo(k,j,i), 4 ) )                                    &
+                wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 4 )
+          ENDDO
+       ENDDO
+    ENDDO
+!
 !-- Exchange ghost points for wall flags
     CALL exchange_horiz_int( wall_flags_0, nys, nyn, nxl, nxr, nzt, nbgp )