Changeset 4895 for palm/trunk/SOURCE/data_output_mask.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/data_output_mask.f90 (modified) (13 diffs)

palm/trunk/SOURCE/data_output_mask.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Remove offset in terrain-following masked output
+! 
+! 4828 2021-01-05 11:21:41Z Giersch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -154 +157 @@
     INTEGER(iwp) ::  k                       !< loop index
     INTEGER(iwp) ::  kk                      !< vertical index
-    INTEGER(iwp) ::  ktt                     !< k index of highest terrain surface
+    INTEGER(iwp) ::  ktt                     !< k index of lowest non-terrain surface
     INTEGER(iwp) ::  mid                     !< masked output running index
     INTEGER(iwp) ::  n                       !< loop index
@@ -288 +291 @@
                    DO  i = 1, mask_size_l(mid,1)
                       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)
@@ -294 +298 @@
                                        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
@@ -356 +360 @@
                    DO  i = 1, mask_size_l(mid,1)
                       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)
@@ -362 +367 @@
                                        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
@@ -399 +404 @@
                       DO  i = 1, mask_size_l(mid,1)
                          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)
@@ -405 +411 @@
                                           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
@@ -493 +500 @@
                    DO  i = 1, mask_size_l(mid,1)
                       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)
@@ -499 +507 @@
                                        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
@@ -532 +541 @@
                    DO  i = 1, mask_size_l(mid,1)
                       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)
@@ -538 +548 @@
                                        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
@@ -679 +689 @@
              DO  i = 1, mask_size_l(mid,1)
                 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)
-                   ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 ) ),           &
+                   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 )
-!--                   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