Changeset 4115 for palm/trunk


Ignore:
Timestamp:
Jul 24, 2019 12:50:49 PM (5 years ago)
Author:
suehring
Message:

Bugfix in setting flags inidicating wall-bounded grid-points (used for control TKE production near walls); Bugfix in chemistry decycling initialization

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4110 r4115  
    2727! -----------------
    2828! $Id$
     29! Fix faulty IF statement in decycling initialization
     30!
     31! 4110 2019-07-22 17:05:21Z suehring
    2932! - Decycling boundary conditions are only set at the ghost points not on the
    3033!   prognostic grid points
     
    20002003          ENDIF
    20012004       ENDIF
    2002        IF ( .NOT. decycle_chem_lr )  THEN
     2005       IF ( decycle_chem_lr )  THEN
    20032006          IF ( nxl == 0  )  THEN
    20042007             DO  i = 1, nbgp   
  • palm/trunk/SOURCE/init_grid.f90

    r4110 r4115  
    2525! -----------------
    2626! $Id$
     27! Bugfix in setting near-surface flag 24, inidicating wall-bounded grid points
     28!
     29! 4110 2019-07-22 17:05:21Z suehring
    2730! - Separate initialization of advection flags for momentum and scalars.
    2831! - Change subroutine interface for ws_init_flags_scalar to pass boundary flags
     
    25782581    CALL exchange_horiz_int( wall_flags_0, nys, nyn, nxl, nxr, nzt, nbgp )
    25792582!
    2580 !-- Set outer array for scalars to mask near-surface grid points in
    2581 !-- production_e
     2583!-- Set outer array for scalars to mask near-surface grid points. Note, on
     2584!-- basis of flag 24 futher flags will be derived which are used to control
     2585!-- production of subgrid TKE production near walls.
    25822586    DO i = nxl, nxr
    25832587       DO j = nys, nyn
    25842588          DO k = nzb, nzt+1
    2585              IF ( BTEST( wall_flags_0(k,j-1,i), 0 )  .AND.                       &
    2586                   BTEST( wall_flags_0(k,j+1,i), 0 )  .AND.                       &
    2587                   BTEST( wall_flags_0(k,j,i-1), 0 )  .AND.                       &
    2588                   BTEST( wall_flags_0(k,j-1,i-1), 0 )  .AND.                       &
    2589                   BTEST( wall_flags_0(k,j+1,i-1), 0 )  .AND.                       &
    2590                   BTEST( wall_flags_0(k,j-1,i+1), 0 )  .AND.                       &
    2591                   BTEST( wall_flags_0(k,j+1,i+1), 0 ) )                            &
     2589             IF ( BTEST( wall_flags_0(k,j-1,i), 0 )    .AND.                   &
     2590                  BTEST( wall_flags_0(k,j+1,i), 0 )    .AND.                   &
     2591                  BTEST( wall_flags_0(k,j,i-1), 0 )    .AND.                   &
     2592                  BTEST( wall_flags_0(k,j,i+1), 0 )    .AND.                   &
     2593                  BTEST( wall_flags_0(k,j-1,i-1), 0 )  .AND.                   &
     2594                  BTEST( wall_flags_0(k,j+1,i-1), 0 )  .AND.                   &
     2595                  BTEST( wall_flags_0(k,j-1,i+1), 0 )  .AND.                   &
     2596                  BTEST( wall_flags_0(k,j+1,i+1), 0 ) )                        &
    25922597                wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 24 )
    25932598          ENDDO
Note: See TracChangeset for help on using the changeset viewer.