Changeset 2021 for palm/trunk/SOURCE/poismg_noopt.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/poismg_noopt.f90 (modified) (2 diffs)

palm/trunk/SOURCE/poismg_noopt.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: restore nest_bound_(l/r/s/n) in case of mg_switch_to_pe0
 ! 
 ! Former revisions:
@@ -1561 +1561 @@
 
 !
-!--       For non-cyclic lateral boundary conditions, restore the
-!--       in-/outflow conditions
+!--       For non-cyclic lateral boundary conditions and in case of nesting, 
+!--       restore the in-/outflow conditions. 
           inflow_l  = .FALSE.;  inflow_r  = .FALSE.
           inflow_n  = .FALSE.;  inflow_s  = .FALSE.
           outflow_l = .FALSE.;  outflow_r = .FALSE.
           outflow_n = .FALSE.;  outflow_s = .FALSE.
+!
+!--       In case of nesting, restore lateral boundary conditions
+          IF ( nest_domain )  THEN
+             nest_bound_l = .FALSE.
+             nest_bound_r = .FALSE.
+             nest_bound_s = .FALSE.
+             nest_bound_n = .FALSE.      
+          ENDIF
 
           IF ( pleft == MPI_PROC_NULL )  THEN