Ignore:
Timestamp:
Jul 17, 2019 4:00:03 PM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/vertical_nesting_mod.f90

    r4101 r4102  
    2626! -----------------
    2727! $Id$
     28! - Slightly revise setting of boundary conditions at horizontal walls, use
     29!   data-structure offset index instead of pre-calculate it for each facing
     30!
     31! 4101 2019-07-17 15:14:26Z gronemeier
    2832! remove old_dt
    2933!
     
    27112715        INTEGER(iwp)                          ::  i
    27122716        INTEGER(iwp)                          ::  j
    2713         INTEGER(iwp)                          ::  k
    2714         INTEGER(iwp)                          ::  kb !< variable to set respective boundary value, depends on facing.
     2717        INTEGER(iwp)                          ::  k
    27152718        INTEGER(iwp)                          ::  l  !< running index boundary type, for up- and downward-facing walls
    27162719        INTEGER(iwp)                          ::  m  !< running index surface elements
     
    28602863            IF ( ibc_pt_b == 1 )  THEN
    28612864               DO  l = 0, 1
    2862 !
    2863 !--       Set kb, for upward-facing surfaces value at topography top (k-1) is set,
    2864 !--       for downward-facing surfaces at topography bottom (k+1).
    2865                     kb = MERGE( -1, 1, l == 0 )
    2866                     DO  m = 1, bc_h(l)%ns
    2867                        i = bc_h(l)%i(m)     
    2868                        j = bc_h(l)%j(m)
    2869                        k = bc_h(l)%k(m)
    2870                        pt(k+kb,j,i) = pt(k,j,i)
    2871                     ENDDO
    2872                  ENDDO
     2865                  DO  m = 1, bc_h(l)%ns
     2866                     i = bc_h(l)%i(m)     
     2867                     j = bc_h(l)%j(m)
     2868                     k = bc_h(l)%k(m)
     2869                     pt(k+bc_h(l)%koff,j,i) = pt(k,j,i)
     2870                  ENDDO
     2871               ENDDO
    28732872            ENDIF
    28742873   
Note: See TracChangeset for help on using the changeset viewer.