Changeset 4745 for palm/trunk/SOURCE


Ignore:
Timestamp:
Oct 15, 2020 4:37:13 PM (3 years ago)
Author:
suehring
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/pmc_interface_mod.f90

    r4734 r4745  
    2525! -----------------
    2626! $Id$
     27! Adjustement of face-area calculation to 3D-topographies
     28!
     29! 4734 2020-10-09 13:21:06Z suehring
    2730! Bugfix for RANS-RANS nesting and minor revision in the TKE treatment (only logical control)
    2831!
     
    23322335    INTEGER(iwp) ::  j       !< Running index in the y-direction
    23332336    INTEGER(iwp) ::  k       !< Running index in the z-direction
    2334     INTEGER(iwp) ::  k_wall  !< Local topography top k-index
    23352337    INTEGER(iwp) ::  n       !< Running index over boundary faces
    23362338
     
    23462348       DO  j = nys, nyn
    23472349          sub_sum = 0.0_wp
    2348           k_wall = topo_top_ind(j,i,1)
    2349           DO   k = k_wall + 1, nzt
    2350              sub_sum = sub_sum + dzw(k)
     2350          DO  k = nzb + 1, nzt
     2351             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 1 ) )
    23512352          ENDDO
    23522353          face_area_local =  face_area_local + dy * sub_sum
     
    23682369       DO  j = nys, nyn
    23692370          sub_sum = 0.0_wp
    2370           k_wall = topo_top_ind(j,i,1)
    2371           DO   k = k_wall + 1, nzt
    2372              sub_sum = sub_sum + dzw(k)
     2371          DO  k = nzb + 1, nzt
     2372             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 1 ) )
    23732373          ENDDO
    23742374          face_area_local =  face_area_local + dy * sub_sum
     
    23902390       DO  i = nxl, nxr
    23912391          sub_sum = 0.0_wp
    2392           k_wall = topo_top_ind(j,i,2)
    2393           DO  k = k_wall + 1, nzt
    2394              sub_sum = sub_sum + dzw(k)
     2392          DO  k = nzb + 1, nzt
     2393             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 2 ) )
    23952394          ENDDO
    23962395          face_area_local = face_area_local + dx * sub_sum
     
    24122411       DO  i = nxl, nxr
    24132412          sub_sum = 0.0_wp
    2414           k_wall = topo_top_ind(j,i,2)
    2415           DO  k = k_wall + 1, nzt
    2416              sub_sum = sub_sum + dzw(k)
     2413          DO  k = nzb + 1, nzt
     2414             sub_sum = sub_sum + dzw(k) * MERGE( 1.0_wp, 0.0_wp, BTEST(wall_flags_total_0(k,j,i), 2 ) )
    24172415          ENDDO
    24182416          face_area_local = face_area_local + dx * sub_sum
Note: See TracChangeset for help on using the changeset viewer.