Changeset 4895 for palm/trunk/SOURCE/diagnostic_output_quantities_mod.f90

Timestamp:

Mar 3, 2021 3:39:08 PM (4 years ago)

Author:

suehring

Message:

Remove offset in terrain-following masked output and allow only mask_k_over_surface >= 1

File:

• palm/trunk/SOURCE/diagnostic_output_quantities_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/diagnostic_output_quantities_mod.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Remove offset in terrain-following masked output
+!
+! 4866 2021-02-05 17:50:14Z suehring
 ! Implemented vertical passive scalar flux
 ! Implemented air temperature as diagnostic output quantity (previous commit)
@@ -1691 +1694 @@
              DO  j = 1, mask_size_l(mid,2)
 !
-!--             Get k index of the highest terraing surface
+!--             Get k index of the lowest non-terrain grid point, else set fill value
                 im = mask_i(mid,i)
                 jm = mask_j(mid,j)
                 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 ) ), DIM=1 ) - 1
                 DO  k = 1, mask_size_l(mid,3)
-                   kk = MIN( ktt+mask_k(mid,k), nzt+1 )
+                   kk = MIN( ktt + mask_k(mid,k) - 1, nzt+1 )
 !
 !--                Set value if not in building