Changeset 4102 for palm/trunk/SOURCE/vertical_nesting_mod.f90

Timestamp:

Jul 17, 2019 4:00:03 PM (5 years ago)

Author:

suehring

Message:

Bugfix, set Neumann boundary conditions for the subgrid TKE at vertical walls instead of implicit Dirichlet conditions that always act as a sink term for the subgrid TKE. Therefore, add new data structure for vertical surfaces and revise the setting of the boundary grid point index space. Moreover, slightly revise setting of boundary conditions at upward- and downward facing surfaces. Finally, setting of boundary conditions for subgrid TKE and dissipation (in RANS mode) is now modularized. Update test case results.

File:

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

palm/trunk/SOURCE/vertical_nesting_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! - Slightly revise setting of boundary conditions at horizontal walls, use 
+!   data-structure offset index instead of pre-calculate it for each facing
+! 
+! 4101 2019-07-17 15:14:26Z gronemeier
 ! remove old_dt
 ! 
@@ -2711 +2715 @@
         INTEGER(iwp)                          ::  i
         INTEGER(iwp)                          ::  j
-        INTEGER(iwp)                          ::  k
-        INTEGER(iwp)                          ::  kb !< variable to set respective boundary value, depends on facing. 
+        INTEGER(iwp)                          ::  k 
         INTEGER(iwp)                          ::  l  !< running index boundary type, for up- and downward-facing walls
         INTEGER(iwp)                          ::  m  !< running index surface elements
@@ -2860 +2863 @@
             IF ( ibc_pt_b == 1 )  THEN 
                DO  l = 0, 1 
-!
-!--       Set kb, for upward-facing surfaces value at topography top (k-1) is set,
-!--       for downward-facing surfaces at topography bottom (k+1). 
-                    kb = MERGE( -1, 1, l == 0 )
-                    DO  m = 1, bc_h(l)%ns
-                       i = bc_h(l)%i(m)     
-                       j = bc_h(l)%j(m)
-                       k = bc_h(l)%k(m)
-                       pt(k+kb,j,i) = pt(k,j,i)
-                    ENDDO
-                 ENDDO
+                  DO  m = 1, bc_h(l)%ns
+                     i = bc_h(l)%i(m)     
+                     j = bc_h(l)%j(m)
+                     k = bc_h(l)%k(m)
+                     pt(k+bc_h(l)%koff,j,i) = pt(k,j,i)
+                  ENDDO
+               ENDDO
             ENDIF