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

Timestamp:

Oct 2, 2015 1:25:23 PM (9 years ago)

Author:

boeske

Message:

Bugfix concerning wall_flags at PE boundaries, added new subroutine exchange_horiz_int

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: Ghost points are included in wall_flags_0 and wall_flags_00
 ! 
 ! Former revisions:
@@ -717 +717 @@
 !-- In case of non-cyclic lateral boundaries, the order of the advection
 !-- scheme have to be reduced up to nzt (required at the lateral boundaries).
-    nzb_max = MAXVAL( nzb_local )
+    nzb_max = MAXVAL( nzb_local ) + 1
     IF ( inflow_l .OR. outflow_l .OR. inflow_r .OR. outflow_r .OR.             &
          inflow_n .OR. outflow_n .OR. inflow_s .OR. outflow_s )  THEN
@@ -1188 +1188 @@
 !
 !-- Allocate flags needed for masking walls.
-    ALLOCATE( wall_flags_0(nzb:nzt,nys:nyn,nxl:nxr), &
-              wall_flags_00(nzb:nzt,nys:nyn,nxl:nxr) )
+    ALLOCATE( wall_flags_0(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                     &
+              wall_flags_00(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     wall_flags_0  = 0
     wall_flags_00 = 0
@@ -1453 +1453 @@
 !
 !--             w component - z-direction
-!--             WS1 (33), WS3 (34), WS5 (35)
+!--             WS1 (33), WS3 (34)IF ( inflow_l .OR. outflow_l )  THEN, WS5 (35)
                 flag_set = .FALSE.
                 IF ( k == nzb_w_inner(j,i) .OR. k == nzb_w_inner(j,i) + 1      &
@@ -1479 +1479 @@
 
 !
+!-- Exchange 3D integer wall_flags.
+    IF ( momentum_advec == 'ws-scheme' .OR. scalar_advec == 'ws-scheme'     &
+    .OR. scalar_advec == 'ws-scheme-mono' )  THEN  
+!
+!--    Exchange ghost points for advection flags
+       CALL exchange_horiz_int( wall_flags_0,  nbgp )
+       CALL exchange_horiz_int( wall_flags_00, nbgp )
+!
+!--    Set boundary flags at inflow and outflow boundary in case of 
+!--    non-cyclic boundary conditions.
+       IF ( inflow_l .OR. outflow_l )  THEN
+          wall_flags_0(:,:,nxl-1)  = wall_flags_0(:,:,nxl)
+          wall_flags_00(:,:,nxl-1) = wall_flags_00(:,:,nxl)
+       ENDIF
+
+       IF ( inflow_r .OR. outflow_r )  THEN
+          wall_flags_0(:,:,nxr+1)  = wall_flags_0(:,:,nxr)
+          wall_flags_00(:,:,nxr+1) = wall_flags_00(:,:,nxr)
+       ENDIF
+
+       IF ( inflow_n .OR. outflow_n )  THEN
+          wall_flags_0(:,nyn+1,:)  = wall_flags_0(:,nyn,:)
+          wall_flags_00(:,nyn+1,:) = wall_flags_00(:,nyn,:)
+       ENDIF
+
+       IF ( inflow_s .OR. outflow_s )  THEN
+          wall_flags_0(:,nys-1,:)  = wall_flags_0(:,nys,:)
+          wall_flags_00(:,nys-1,:) = wall_flags_00(:,nys,:)
+       ENDIF
+
+    ENDIF
+
+!
 !-- In case of topography: limit near-wall mixing length l_wall further:
 !-- Go through all points of the subdomain one by one and look for the closest