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

Timestamp:

Dec 9, 2019 6:53:04 PM (4 years ago)

Author:

suehring

Message:

Minor change in nzb_max computation; commentation added; minor formatting adjustments in advec_ws

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor change in nzb_max computation. Commentation added.
+! 
+! 4314 2019-11-29 10:29:20Z suehring
 ! Set additional topography flag 4 to mark topography grid points emerged
 ! from the filtering process.
@@ -612 +615 @@
 !-- Determine the maximum level of topography. It is used for
 !-- steering the degradation of order of the applied advection scheme, 
-!-- as well in the lpm. 
+!-- as well in the lpm.
     k_top = 0
     DO  i = nxl, nxr
@@ -622 +625 @@
     ENDDO
 #if defined( __parallel )
-    CALL MPI_ALLREDUCE( k_top + 1, nzb_max, 1, MPI_INTEGER,                    & !is +1 really necessary here?
+    CALL MPI_ALLREDUCE( k_top, nzb_max, 1, MPI_INTEGER,                        &
                         MPI_MAX, comm2d, ierr )
 #else
-    nzb_max = k_top + 1
+    nzb_max = k_top
 #endif
-!   
-!-- If topography extents up to the model top, limit nzb_max to nzt.
-    nzb_max = MIN( nzb_max, nzt ) 
+!
+!-- Increment nzb_max by 1 in order to allow for proper diverengence correction.
+!-- Further, in case topography extents up to the model top, limit to nzt.
+    nzb_max = MIN( nzb_max+1, nzt ) 
 !
 !-- Determine minimum index of topography. Usually, this will be nzb. In case