Changeset 4115 for palm/trunk/SOURCE/init_grid.f90
- Timestamp:
- Jul 24, 2019 12:50:49 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r4110 r4115 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix in setting near-surface flag 24, inidicating wall-bounded grid points 28 ! 29 ! 4110 2019-07-22 17:05:21Z suehring 27 30 ! - Separate initialization of advection flags for momentum and scalars. 28 31 ! - Change subroutine interface for ws_init_flags_scalar to pass boundary flags … … 2578 2581 CALL exchange_horiz_int( wall_flags_0, nys, nyn, nxl, nxr, nzt, nbgp ) 2579 2582 ! 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. 2582 2586 DO i = nxl, nxr 2583 2587 DO j = nys, nyn 2584 2588 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 ) ) & 2592 2597 wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 24 ) 2593 2598 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.