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

Timestamp:

Oct 16, 2017 5:12:01 PM (7 years ago)

Author:

boeske

Message:

enable simulations with complex terrain

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set lateral boundary conditions for topography on all three ghost layers
+! 
+! 2478 2017-09-18 13:37:24Z suehring
 ! Bugfix, correct flag for use_top
 ! 
@@ -1041 +1044 @@
           
           IF ( .NOT. bc_ns_cyc )  THEN
-             IF ( nys == 0  )  nzb_local(-1,:)   = nzb_local(0,:)
-             IF ( nyn == ny )  nzb_local(ny+1,:) = nzb_local(ny,:)
+             IF ( nys == 0  )  THEN
+                nzb_local(-1,:)   = nzb_local(0,:)
+                nzb_local(-2,:)   = nzb_local(0,:)
+                nzb_local(-3,:)   = nzb_local(0,:)
+             ELSEIF ( nyn == ny )  THEN
+                nzb_local(ny+1,:) = nzb_local(ny,:)
+                nzb_local(ny+2,:) = nzb_local(ny,:)
+                nzb_local(ny+3,:) = nzb_local(ny,:)
+             ENDIF
           ENDIF
 
           IF ( .NOT. bc_lr_cyc )  THEN
-             IF ( nxl == 0  )  nzb_local(:,-1)   = nzb_local(:,0)
-             IF ( nxr == nx )  nzb_local(:,nx+1) = nzb_local(:,nx)          
+             IF ( nxl == 0  )  THEN
+                nzb_local(:,-1)   = nzb_local(:,0)
+                nzb_local(:,-2)   = nzb_local(:,0)
+                nzb_local(:,-3)   = nzb_local(:,0)
+             ELSEIF ( nxr == nx )  THEN
+                nzb_local(:,nx+1) = nzb_local(:,nx)
+                nzb_local(:,nx+2) = nzb_local(:,nx)          
+                nzb_local(:,nx+3) = nzb_local(:,nx)          
+             ENDIF          
           ENDIF