Changeset 2939

Timestamp:

Mar 29, 2018 6:20:00 PM (7 years ago)

Author:

suehring

Message:

Set lateral boundary conditions for divergence in multigrid solver

Location:

palm/trunk/SOURCE

Files:

• poismg_mod.f90 (modified) (3 diffs)
• poismg_noopt_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/poismg_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set lateral boundary conditions for divergence
+! 
+! 2937 2018-03-27 14:58:33Z suehring
 ! Corrected "Former revisions" section
 ! 
@@ -149 +152 @@
 
        USE control_parameters,                                                 &
-           ONLY:  gathered_size, grid_level, grid_level_count,                 &
+           ONLY:  bc_lr_cyc, bc_ns_cyc, gathered_size, grid_level,             &
+                  grid_level_count, ibc_p_t,                                   &
                   maximum_grid_level, message_string, mgcycles, mg_cycles,     &
                   mg_switch_to_pe0_level, residual_limit, subdomain_size
@@ -195 +199 @@
        grid_level = maximum_grid_level
        CALL exchange_horiz( d, 1)
+!
+!--    Set bottom and top boundary conditions
        d(nzb,:,:) = d(nzb+1,:,:)
-
+       IF ( ibc_p_t == 1 )  d(nzt+1,:,: ) = d(nzt,:,:)
+!
+!--    Set lateral boundary conditions in non-cyclic case
+       IF ( .NOT. bc_lr_cyc )  THEN
+          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  &
+             d(:,:,nxl-1) = d(:,:,nxl)
+          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r .OR. force_bound_r )  &
+             d(:,:,nxr+1) = d(:,:,nxr)
+       ENDIF
+       IF ( .NOT. bc_ns_cyc )  THEN
+          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n .OR. force_bound_n )  &
+             d(:,nyn+1,:) = d(:,nyn,:)
+          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s .OR. force_bound_s )  &
+             d(:,nys-1,:) = d(:,nys,:)
+       ENDIF
 !
 !--    Initiation of the multigrid scheme. Does n cycles until the

palm/trunk/SOURCE/poismg_noopt_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set lateral boundary conditions for divergence
+! 
+! 2793 2018-02-07 10:54:33Z suehring
 ! Corrected "Former revisions" section
 ! 
@@ -195 +198 @@
 
        USE control_parameters,                                                 &
-           ONLY:  gathered_size, grid_level_count,                             &
+           ONLY:  bc_lr_cyc, bc_ns_cyc, gathered_size, grid_level_count,       &
+                  ibc_p_t, maximum_grid_level, message_string, mgcycles,       &
+                  mg_cycles, mg_switch_to_pe0_level, residual_limit,           &
+                  subdomain_size
+
+       USE control_parameters,                                                 &
+           ONLY:  bc_lr_cyc, bc_ns_cyc, gathered_size, grid_level,             &
+                  grid_level_count, ibc_p_t,                                   &
                   maximum_grid_level, message_string, mgcycles, mg_cycles,     &
                   mg_switch_to_pe0_level, residual_limit, subdomain_size
@@ -244 +254 @@
        grid_level = maximum_grid_level
        CALL exchange_horiz( d, 1)
+!
+!--    Set bottom and top boundary conditions
        d(nzb,:,:) = d(nzb+1,:,:)
+       IF ( ibc_p_t == 1 )  d(nzt+1,:,: ) = d(nzt,:,:)
+!
+!--    Set lateral boundary conditions in non-cyclic case
+       IF ( .NOT. bc_lr_cyc )  THEN
+          IF ( inflow_l .OR. outflow_l .OR. nest_bound_l .OR. force_bound_l )  &
+             d(:,:,nxl-1) = d(:,:,nxl)
+          IF ( inflow_r .OR. outflow_r .OR. nest_bound_r .OR. force_bound_r )  &
+             d(:,:,nxr+1) = d(:,:,nxr)
+       ENDIF
+       IF ( .NOT. bc_ns_cyc )  THEN
+          IF ( inflow_n .OR. outflow_n .OR. nest_bound_n .OR. force_bound_n )  &
+             d(:,nyn+1,:) = d(:,nyn,:)
+          IF ( inflow_s .OR. outflow_s .OR. nest_bound_s .OR. force_bound_s )  &
+             d(:,nys-1,:) = d(:,nys,:)
+       ENDIF
 
 !