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

Timestamp:

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

Author:

suehring

Message:

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

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Remove offset in terrain-following masked output
+!
+! 4887 2021-02-26 16:22:32Z banzhafs
 ! Unnecessary comments removed
 !
@@ -1696 +1699 @@
     INTEGER(iwp) ::  k               !< grid index along z-direction
     INTEGER(iwp) ::  kk              !< masked output index along z-direction
-    INTEGER(iwp) ::  ktt             !< k index of highest terrain surface
+    INTEGER(iwp) ::  ktt             !< k index of lowest non-terrain grid point
     INTEGER(iwp) ::  lsp
     INTEGER(iwp) ::  mid             !< masked output running index
@@ -1732 +1735 @@
                    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
                       im = mask_i(mid,i)
                       jm = mask_j(mid,j)
@@ -1738 +1741 @@
                                     DIM = 1 ) - 1
                       DO  k = 1, mask_size_l(mid,3)
-                         kk = MIN( ktt+mask_k(mid,k), nzt+1 )
-!
-!--                      Set value if not in building
+                         kk = MIN( ktt + mask_k(mid,k) - 1, nzt+1 )
+!
+!--                      Set value if not in building, else set fill value
                          IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
                             local_pf(i,j,k) = fill_value
@@ -1770 +1773 @@
                    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
                       im = mask_i(mid,i)
                       jm = mask_j(mid,j)
@@ -1776 +1779 @@
                                     DIM = 1 ) - 1
                       DO  k = 1, mask_size_l(mid,3)
-                         kk = MIN( ktt+mask_k(mid,k), nzt+1 )
-!
-!--                      Set value if not in building
+                         kk = MIN( ktt + mask_k(mid,k) - 1, nzt+1 )
+!
+!--                      Set value if not in building, else set fill value
                          IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
                             local_pf(i,j,k) = fill_value