Changeset 2823

Timestamp:

Feb 20, 2018 3:31:45 PM (7 years ago)

Author:

Giersch

Message:

Set boundary condtions for 3D topography in case of non-cyclic boundary conditions

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Set boundary conditions for 3D topography in case of non-cyclic boundary 
+! conditions
+! 
+! 2796 2018-02-08 12:25:39Z suehring
 ! Bugfix in 3D building initialization
 ! 
@@ -1589 +1593 @@
 
           CALL exchange_horiz_int( topo, nys, nyn, nxl, nxr, nzt, nbgp )
+!
+!--       Set boundary conditions also for flags. Can be interpreted as Neumann
+!--       boundary conditions for topography. 
+          IF ( .NOT. bc_ns_cyc )  THEN
+             IF ( nys == 0  )  THEN
+                DO  i = 1, nbgp     
+                   topo(:,nys-i,:)   = topo(:,nys,:)
+                ENDDO
+             ENDIF
+             IF ( nyn == ny )  THEN
+                DO  i = 1, nbgp  
+                   topo(:,nyn+i,:) = topo(:,nyn,:)
+                ENDDO
+             ENDIF
+          ENDIF
+          IF ( .NOT. bc_lr_cyc )  THEN
+             IF ( nxl == 0  )  THEN
+                DO  i = 1, nbgp   
+                   topo(:,:,nxl-i)   = topo(:,:,nxl)
+                ENDDO
+             ENDIF
+             IF ( nxr == nx )  THEN 
+                DO  i = 1, nbgp   
+                   topo(:,:,nxr+i) = topo(:,:,nxr)     
+                ENDDO
+             ENDIF      
+          ENDIF
 
        CASE ( 'single_street_canyon' )
@@ -1673 +1704 @@
 
           CALL exchange_horiz_int( topo, nys, nyn, nxl, nxr, nzt, nbgp )
+!
+!--       Set boundary conditions also for flags. Can be interpreted as Neumann
+!--       boundary conditions for topography. 
+          IF ( .NOT. bc_ns_cyc )  THEN
+             IF ( nys == 0  )  THEN
+                DO  i = 1, nbgp     
+                   topo(:,nys-i,:)   = topo(:,nys,:)
+                ENDDO
+             ENDIF
+             IF ( nyn == ny )  THEN
+                DO  i = 1, nbgp  
+                   topo(:,nyn+i,:) = topo(:,nyn,:)
+                ENDDO
+             ENDIF
+          ENDIF
+          IF ( .NOT. bc_lr_cyc )  THEN
+             IF ( nxl == 0  )  THEN
+                DO  i = 1, nbgp   
+                   topo(:,:,nxl-i)   = topo(:,:,nxl)
+                ENDDO
+             ENDIF
+             IF ( nxr == nx )  THEN 
+                DO  i = 1, nbgp   
+                   topo(:,:,nxr+i) = topo(:,:,nxr)     
+                ENDDO
+             ENDIF      
+          ENDIF
 
        CASE ( 'tunnel' )
@@ -1812 +1870 @@
 
           CALL exchange_horiz_int( topo, nys, nyn, nxl, nxr, nzt, nbgp )
+!
+!--       Set boundary conditions also for flags. Can be interpreted as Neumann
+!--       boundary conditions for topography. 
+          IF ( .NOT. bc_ns_cyc )  THEN
+             IF ( nys == 0  )  THEN
+                DO  i = 1, nbgp     
+                   topo(:,nys-i,:)   = topo(:,nys,:)
+                ENDDO
+             ENDIF
+             IF ( nyn == ny )  THEN
+                DO  i = 1, nbgp  
+                   topo(:,nyn+i,:) = topo(:,nyn,:)
+                ENDDO
+             ENDIF
+          ENDIF
+          IF ( .NOT. bc_lr_cyc )  THEN
+             IF ( nxl == 0  )  THEN
+                DO  i = 1, nbgp   
+                   topo(:,:,nxl-i)   = topo(:,:,nxl)
+                ENDDO
+             ENDIF
+             IF ( nxr == nx )  THEN 
+                DO  i = 1, nbgp   
+                   topo(:,:,nxr+i) = topo(:,:,nxr)     
+                ENDDO
+             ENDIF      
+          ENDIF
 
        CASE ( 'read_from_file' )