Changeset 4115

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:

• chemistry_model_mod.f90 (modified) (2 diffs)
• init_grid.f90 (modified) (2 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Fix faulty IF statement in decycling initialization
+! 
+! 4110 2019-07-22 17:05:21Z suehring
 ! - Decycling boundary conditions are only set at the ghost points not on the
 !   prognostic grid points 
@@ -2000 +2003 @@
           ENDIF
        ENDIF
-       IF ( .NOT. decycle_chem_lr )  THEN
+       IF ( decycle_chem_lr )  THEN
           IF ( nxl == 0  )  THEN
              DO  i = 1, nbgp   

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in setting near-surface flag 24, inidicating wall-bounded grid points 
+! 
+! 4110 2019-07-22 17:05:21Z suehring
 ! - Separate initialization of advection flags for momentum and scalars.
 ! - 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 )
 !
-!-- Set outer array for scalars to mask near-surface grid points in 
-!-- production_e
+!-- Set outer array for scalars to mask near-surface grid points. Note, on 
+!-- basis of flag 24 futher flags will be derived which are used to control
+!-- production of subgrid TKE production near walls. 
     DO i = nxl, nxr
        DO j = nys, nyn
           DO k = nzb, nzt+1
-             IF ( BTEST( wall_flags_0(k,j-1,i), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j+1,i), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j,i-1), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j-1,i-1), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j+1,i-1), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j-1,i+1), 0 )  .AND.                       &
-                  BTEST( wall_flags_0(k,j+1,i+1), 0 ) )                            &
+             IF ( BTEST( wall_flags_0(k,j-1,i), 0 )    .AND.                   &
+                  BTEST( wall_flags_0(k,j+1,i), 0 )    .AND.                   &
+                  BTEST( wall_flags_0(k,j,i-1), 0 )    .AND.                   &
+                  BTEST( wall_flags_0(k,j,i+1), 0 )    .AND.                   &
+                  BTEST( wall_flags_0(k,j-1,i-1), 0 )  .AND.                   &
+                  BTEST( wall_flags_0(k,j+1,i-1), 0 )  .AND.                   &
+                  BTEST( wall_flags_0(k,j-1,i+1), 0 )  .AND.                   &
+                  BTEST( wall_flags_0(k,j+1,i+1), 0 ) )                        &
                 wall_flags_0(k,j,i) = IBSET( wall_flags_0(k,j,i), 24 )
           ENDDO