Changeset 2927 for palm

Timestamp:

Mar 23, 2018 3:13:00 PM (7 years ago)

Author:

suehring

Message:

Bugfix in setting boundary condition for topography index array in case of non-cyclic boundary conditions and Piascek-Williams advection scheme is used.

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, setting boundary conditions for topography index array.
+! 
+! 2918 2018-03-21 15:52:14Z gronemeier
 ! Moved init_mixing_length to turbulence_closure_mod.f90
 ! 
@@ -655 +658 @@
     IF ( .NOT. bc_ns_cyc )  THEN
        IF ( nys == 0  )  THEN
-          nzb_local(-1,:)   = nzb_local(0,:)
+          DO  i = 1, nbgp  
+             nzb_local(nys-i,:)   = nzb_local(nys,:)
+          ENDDO
        ELSEIF ( nyn == ny )  THEN
-          nzb_local(ny+1,:) = nzb_local(ny,:)
+          DO  i = 1, nbgp  
+             nzb_local(ny+i,:) = nzb_local(ny,:)
+          ENDDO
        ENDIF
     ENDIF
@@ -663 +670 @@
     IF ( .NOT. bc_lr_cyc )  THEN
        IF ( nxl == 0  )  THEN
-          nzb_local(:,-1)   = nzb_local(:,0)
-          nzb_local(:,-2)   = nzb_local(:,0)
-          nzb_local(:,-3)   = nzb_local(:,0)
+          DO  i = 1, nbgp  
+             nzb_local(:,nxl-i)   = nzb_local(:,nxl)
+          ENDDO
        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)          
+          DO  i = 1, nbgp  
+             nzb_local(:,nx+i) = nzb_local(:,nx)
+          ENDDO 
        ENDIF          
     ENDIF