Ignore:
Timestamp:
Jan 8, 2019 6:22:50 PM (2 years ago)
Author:
suehring
Message:

Split loops in advec_ws in order to reduce bit queries; Introduce new parameter to better control order degradation of advection scheme at non-cyclic boundaries; Remove setting of Neumann conditions for horizontal velocity variances; Minor bugfix in divergence correction in advection scheme (only has implications at downward-facing wall surfaces)

File:
1 edited

Legend:

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

    r3655 r3661  
    2525! -----------------
    2626! $Id$
     27! Remove setting of nzb_max to nzt at non-cyclic boundary PEs, instead,
     28! order degradation of advection scheme is handeled directly in advec_ws
     29!
     30! 3655 2019-01-07 16:51:22Z knoop
    2731! Comment added
    2832!
     
    823827!-- steering the degradation of order of the applied advection scheme,
    824828!-- as well in the lpm.
    825 !-- In case of non-cyclic lateral boundaries, the order of the advection
    826 !-- scheme has to be reduced up to nzt (required at the lateral boundaries).
    827829    k_top = 0
    828830    DO  i = nxl, nxr
     
    840842    nzb_max = k_top + 1
    841843#endif
    842     IF ( bc_dirichlet_l  .OR.  bc_radiation_l     .OR.                         &
    843          bc_dirichlet_r  .OR.  bc_radiation_r     .OR.                         &
    844          bc_dirichlet_n  .OR.  bc_radiation_n     .OR.                         &
    845          bc_dirichlet_s  .OR.  bc_radiation_s )                                &
    846          nzb_max = nzt
    847844!   
    848 !-- Finally, if topography extents up to the model top, limit nzb_max to nzt.
     845!-- If topography extents up to the model top, limit nzb_max to nzt.
    849846    nzb_max = MIN( nzb_max, nzt )
    850847!
Note: See TracChangeset for help on using the changeset viewer.