Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (23 months 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:
1 edited

Legend:

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

    r4329 r4346  
    2525! -----------------
    2626! $Id$
     27! Introduction of wall_flags_total_0, which currently sets bits based on static
     28! topography information used in wall_flags_static_0
     29!
     30! 4329 2019-12-10 15:46:36Z motisi
    2731! Renamed wall_flags_0 to wall_flags_static_0
    2832!
     
    8387               nzb,                         &
    8488               nzt,                         &
    85                wall_flags_static_0
     89               wall_flags_total_0
    8690
    8791    USE kinds
     
    205209!
    206210!--    If the grid point lies in a building, a new point is defined
    207        DO WHILE ( .NOT. BTEST( wall_flags_static_0(cgp_k,cgp_j,cgp_i), 1 ) )
     211       DO WHILE ( .NOT. BTEST( wall_flags_total_0(cgp_k,cgp_j,cgp_i), 1 ) )
    208212          CALL RANDOM_NUMBER( random )
    209213          cgp_k = cgp_k + FLOOR( ( nzt - cgp_k ) * random )   !< Random number upon cgp_k
     
    454458       DO  j = nys, nyn
    455459          DO  k = nzb+1, nzt+1
    456              flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) )
     460             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), quant_type ) )
    457461             quant_av_k_l(k) = quant_av_k_l(k) + quant(k,j,i) * flag
    458462             count_2d_l(k)   = count_2d_l(k) + INT( flag, KIND=iwp )
     
    482486                                + ( quant(k,j,i) - quant_av_k(k) )**2 &
    483487                                * MERGE( 1.0_wp, 0.0_wp,              &
    484                                          BTEST( wall_flags_static_0(k,j,i), quant_type ) )
     488                                         BTEST( wall_flags_total_0(k,j,i), quant_type ) )
    485489          ENDDO
    486490       ENDDO
     
    722726      DO  j = nys, nyn
    723727         DO  k = nzb, nzt+1
    724             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) )
     728            flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i), quant_type ) )
    725729            quant_av_l = quant_av_l + quant(k,j,i) * flag
    726730            average_count_l = average_count_l + INT( flag, KIND=iwp )
     
    774778                              + u(k,j,i) * dzw(k) * dy                  &
    775779                              * MERGE( 1.0_wp, 0.0_wp,                  &
    776                                        BTEST( wall_flags_static_0(k,j,i), 1 )  &
     780                                 BTEST( wall_flags_total_0(k,j,i), 1 )  &
    777781                                     )
    778782          ENDDO
     
    788792                              - u(k,j,i) * dzw(k) * dy                  &
    789793                              * MERGE( 1.0_wp, 0.0_wp,                  &
    790                                        BTEST( wall_flags_static_0(k,j,i), 1 )  &
     794                                 BTEST( wall_flags_total_0(k,j,i), 1 )  &
    791795                                     )
    792796          ENDDO
     
    803807                              + v(k,j,i) * dzw(k) * dx                  &
    804808                              * MERGE( 1.0_wp, 0.0_wp,                  &
    805                                        BTEST( wall_flags_static_0(k,j,i), 2 )  &
     809                                 BTEST( wall_flags_total_0(k,j,i), 2 )  &
    806810                                     )
    807811          ENDDO
     
    817821                              - v(k,j,i) * dzw(k) * dx                  &
    818822                              * MERGE( 1.0_wp, 0.0_wp,                  &
    819                                        BTEST( wall_flags_static_0(k,j,i), 2 )  &
     823                                 BTEST( wall_flags_total_0(k,j,i), 2 )  &
    820824                                     )
    821825          ENDDO
Note: See TracChangeset for help on using the changeset viewer.