Changeset 4346 for palm/trunk/SOURCE/vdi_internal_controls.f90

Timestamp:

Dec 18, 2019 11:55:56 AM (5 years ago)

Author:

motisi

Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:

• palm/trunk/SOURCE/vdi_internal_controls.f90 (modified) (10 diffs)

palm/trunk/SOURCE/vdi_internal_controls.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Introduction of wall_flags_total_0, which currently sets bits based on static
+! topography information used in wall_flags_static_0
+! 
+! 4329 2019-12-10 15:46:36Z motisi
 ! Renamed wall_flags_0 to wall_flags_static_0
 ! 
@@ -83 +87 @@
                nzb,                         &
                nzt,                         &
-               wall_flags_static_0
+               wall_flags_total_0
 
     USE kinds
@@ -205 +209 @@
 !
 !--    If the grid point lies in a building, a new point is defined
-       DO WHILE ( .NOT. BTEST( wall_flags_static_0(cgp_k,cgp_j,cgp_i), 1 ) )
+       DO WHILE ( .NOT. BTEST( wall_flags_total_0(cgp_k,cgp_j,cgp_i), 1 ) )
           CALL RANDOM_NUMBER( random )
           cgp_k = cgp_k + FLOOR( ( nzt - cgp_k ) * random )   !< Random number upon cgp_k
@@ -454 +458 @@
        DO  j = nys, nyn
           DO  k = nzb+1, nzt+1
-             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) )
+             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), quant_type ) )
              quant_av_k_l(k) = quant_av_k_l(k) + quant(k,j,i) * flag
              count_2d_l(k)   = count_2d_l(k) + INT( flag, KIND=iwp )
@@ -482 +486 @@
                                 + ( quant(k,j,i) - quant_av_k(k) )**2 &
                                 * MERGE( 1.0_wp, 0.0_wp,              &
-                                         BTEST( wall_flags_static_0(k,j,i), quant_type ) )
+                                         BTEST( wall_flags_total_0(k,j,i), quant_type ) )
           ENDDO
        ENDDO
@@ -722 +726 @@
       DO  j = nys, nyn
          DO  k = nzb, nzt+1
-            flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) )
+            flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), quant_type ) )
             quant_av_l = quant_av_l + quant(k,j,i) * flag
             average_count_l = average_count_l + INT( flag, KIND=iwp )
@@ -774 +778 @@
                               + u(k,j,i) * dzw(k) * dy                  &
                               * MERGE( 1.0_wp, 0.0_wp,                  &
-                                       BTEST( wall_flags_static_0(k,j,i), 1 )  &
+                                 BTEST( wall_flags_total_0(k,j,i), 1 )  &
                                      )
           ENDDO
@@ -788 +792 @@
                               - u(k,j,i) * dzw(k) * dy                  &
                               * MERGE( 1.0_wp, 0.0_wp,                  &
-                                       BTEST( wall_flags_static_0(k,j,i), 1 )  &
+                                 BTEST( wall_flags_total_0(k,j,i), 1 )  &
                                      )
           ENDDO
@@ -803 +807 @@
                               + v(k,j,i) * dzw(k) * dx                  &
                               * MERGE( 1.0_wp, 0.0_wp,                  &
-                                       BTEST( wall_flags_static_0(k,j,i), 2 )  &
+                                 BTEST( wall_flags_total_0(k,j,i), 2 )  &
                                      )
           ENDDO
@@ -817 +821 @@
                               - v(k,j,i) * dzw(k) * dx                  &
                               * MERGE( 1.0_wp, 0.0_wp,                  &
-                                       BTEST( wall_flags_static_0(k,j,i), 2 )  &
+                                 BTEST( wall_flags_total_0(k,j,i), 2 )  &
                                      )
           ENDDO