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

Timestamp:

Oct 7, 2016 2:08:57 PM (8 years ago)

Author:

suehring

Message:

Bugfix, restore flags for nest boundaries in multigrid solver; bugfix: setting Neumann boundary conditions for topography arrays in case of non-cyclic boundary conditions

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: setting Neumann boundary conditions for topography required for 
+! topography flags in multigrid_noopt solver
 ! 
 ! Former revisions:
@@ -211 +212 @@
                io_blocks, io_group, inflow_l, inflow_n, inflow_r, inflow_s,    &
                masking_method, maximum_grid_level, message_string,             &
-               momentum_advec, nest_domain, ocean, outflow_l, outflow_n,       &
+               momentum_advec, nest_domain, nest_bound_l, nest_bound_n,        &
+               nest_bound_r, nest_bound_s, ocean, outflow_l, outflow_n,        &
                outflow_r, outflow_s, psolver, scalar_advec, topography,        &
                topography_grid_convention, use_surface_fluxes, use_top_fluxes, &
-               wall_adjustment_factor 
-        
+               wall_adjustment_factor
+         
     USE grid_variables,                                                        &
         ONLY:  ddx, ddx2, ddx2_mg, ddy, ddy2, ddy2_mg, dx, dx2, dy, dy2, fwxm, &
@@ -1343 +1345 @@
 !--          Set non-cyclic boundary conditions on respective multigrid level
              IF ( .NOT. bc_ns_cyc )  THEN
-                IF ( nys == 0  )  THEN
+                IF ( inflow_s .OR. outflow_s .OR. nest_bound_s  )  THEN
                    nzb_tmp(-2,:) = nzb_tmp(0,:)
                    nzb_tmp(-1,:) = nzb_tmp(0,:)
                 ENDIF
-                IF ( nyn == ny )  THEN
-                   nzb_tmp(ny+2,:) = nzb_tmp(ny,:)
-                   nzb_tmp(ny+1,:) = nzb_tmp(ny,:)
+                IF ( inflow_n .OR. outflow_n .OR. nest_bound_n )  THEN
+                   nzb_tmp(nyn_l+2,:) = nzb_tmp(nyn_l,:)
+                   nzb_tmp(nyn_l+1,:) = nzb_tmp(nyn_l,:)
                 ENDIF
              ENDIF
              IF ( .NOT. bc_lr_cyc )  THEN
-                IF ( nxl == 0  )  THEN
+                IF ( inflow_l .OR. outflow_l .OR. nest_bound_l  )  THEN
                    nzb_tmp(:,-2) = nzb_tmp(:,0)
                    nzb_tmp(:,-1) = nzb_tmp(:,0)
                 ENDIF
-                IF ( nxr == nx )  THEN
-                   nzb_tmp(:,nx+1) = nzb_tmp(:,nx)    
-                   nzb_tmp(:,nx+2) = nzb_tmp(:,nx)      
+                IF ( inflow_r .OR. outflow_r .OR. nest_bound_r )  THEN
+                   nzb_tmp(:,nxr_l+1) = nzb_tmp(:,nxr_l)    
+                   nzb_tmp(:,nxr_l+2) = nzb_tmp(:,nxr_l)      
                 ENDIF        
              ENDIF
@@ -1420 +1422 @@
 !
 !-- Allocate flags needed for masking walls. Even though these flags are only
-!-- required in the ws-scheme, the arrays need to be allocated as they are 
+!-- required in the ws-scheme, the arrays need to be allocated here as they are 
 !-- used in OpenACC directives. 
     ALLOCATE( wall_flags_0(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                     &
@@ -1427 +1429 @@
     wall_flags_00 = 0
 !
-!-- Init flags for ws-scheme to degrade order near walls
+!-- Init flags for ws-scheme to degrade order of the numerics near walls, i.e. 
+!-- to decrease the numerical stencil appropriately.
     IF ( momentum_advec == 'ws-scheme'  .OR.  scalar_advec == 'ws-scheme'  .OR.&
          scalar_advec   == 'ws-scheme-mono' )  THEN