Changeset 4745 for palm/trunk/SOURCE

Timestamp:

Oct 15, 2020 4:37:13 PM (4 years ago)

Author:

suehring

Message:

mass-flux conservation in nesting: adjustement of face-area calculation to 3D-topographies

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Adjustement of face-area calculation to 3D-topographies
+!
+! 4734 2020-10-09 13:21:06Z suehring
 ! Bugfix for RANS-RANS nesting and minor revision in the TKE treatment (only logical control)
 !
@@ -2332 +2335 @@
     INTEGER(iwp) ::  j       !< Running index in the y-direction
     INTEGER(iwp) ::  k       !< Running index in the z-direction
-    INTEGER(iwp) ::  k_wall  !< Local topography top k-index
     INTEGER(iwp) ::  n       !< Running index over boundary faces
 
@@ -2346 +2348 @@
        DO  j = nys, nyn
           sub_sum = 0.0_wp
-          k_wall = topo_top_ind(j,i,1)
-          DO   k = k_wall + 1, nzt
-             sub_sum = sub_sum + dzw(k)
+          DO  k = nzb + 1, nzt
+             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 1 ) )
           ENDDO
           face_area_local =  face_area_local + dy * sub_sum
@@ -2368 +2369 @@
        DO  j = nys, nyn
           sub_sum = 0.0_wp
-          k_wall = topo_top_ind(j,i,1)
-          DO   k = k_wall + 1, nzt
-             sub_sum = sub_sum + dzw(k)
+          DO  k = nzb + 1, nzt
+             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 1 ) )
           ENDDO
           face_area_local =  face_area_local + dy * sub_sum
@@ -2390 +2390 @@
        DO  i = nxl, nxr
           sub_sum = 0.0_wp
-          k_wall = topo_top_ind(j,i,2)
-          DO  k = k_wall + 1, nzt
-             sub_sum = sub_sum + dzw(k)
+          DO  k = nzb + 1, nzt
+             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 2 ) )
           ENDDO
           face_area_local = face_area_local + dx * sub_sum
@@ -2412 +2411 @@
        DO  i = nxl, nxr
           sub_sum = 0.0_wp
-          k_wall = topo_top_ind(j,i,2)
-          DO  k = k_wall + 1, nzt
-             sub_sum = sub_sum + dzw(k)
+          DO  k = nzb + 1, nzt
+             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 2 ) )
           ENDDO
           face_area_local = face_area_local + dx * sub_sum