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

Timestamp:

Jan 8, 2019 6:22:50 PM (5 years ago)

Author:

suehring

Message:

Split loops in advec_ws in order to reduce bit queries; Introduce new parameter to better control order degradation of advection scheme at non-cyclic boundaries; Remove setting of Neumann conditions for horizontal velocity variances; Minor bugfix in divergence correction in advection scheme (only has implications at downward-facing wall surfaces)

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove setting of nzb_max to nzt at non-cyclic boundary PEs, instead, 
+! order degradation of advection scheme is handeled directly in advec_ws
+! 
+! 3655 2019-01-07 16:51:22Z knoop
 ! Comment added
 ! 
@@ -823 +827 @@
 !-- steering the degradation of order of the applied advection scheme, 
 !-- as well in the lpm. 
-!-- In case of non-cyclic lateral boundaries, the order of the advection
-!-- scheme has to be reduced up to nzt (required at the lateral boundaries).
     k_top = 0
     DO  i = nxl, nxr
@@ -840 +842 @@
     nzb_max = k_top + 1
 #endif
-    IF ( bc_dirichlet_l  .OR.  bc_radiation_l     .OR.                         &
-         bc_dirichlet_r  .OR.  bc_radiation_r     .OR.                         &
-         bc_dirichlet_n  .OR.  bc_radiation_n     .OR.                         &
-         bc_dirichlet_s  .OR.  bc_radiation_s )                                &
-         nzb_max = nzt
 !   
-!-- Finally, if topography extents up to the model top, limit nzb_max to nzt.
+!-- If topography extents up to the model top, limit nzb_max to nzt.
     nzb_max = MIN( nzb_max, nzt ) 
 !