Changeset 4329 for palm/trunk/SOURCE/vdi_internal_controls.f90
 Timestamp:
 Dec 10, 2019 3:46:36 PM (20 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/vdi_internal_controls.f90
r4182 r4329 25 25 !  26 26 ! $Id$ 27 ! Renamed wall_flags_0 to wall_flags_static_0 28 ! 29 ! 4182 20190822 15:20:23Z scharf 27 30 ! added "Authors" section 28 31 ! … … 80 83 nzb, & 81 84 nzt, & 82 wall_flags_ 085 wall_flags_static_0 83 86 84 87 USE kinds … … 202 205 ! 203 206 ! If the grid point lies in a building, a new point is defined 204 DO WHILE ( .NOT. BTEST( wall_flags_ 0(cgp_k,cgp_j,cgp_i), 1 ) )207 DO WHILE ( .NOT. BTEST( wall_flags_static_0(cgp_k,cgp_j,cgp_i), 1 ) ) 205 208 CALL RANDOM_NUMBER( random ) 206 209 cgp_k = cgp_k + FLOOR( ( nzt  cgp_k ) * random ) !< Random number upon cgp_k … … 451 454 DO j = nys, nyn 452 455 DO k = nzb+1, nzt+1 453 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), quant_type ) )456 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) ) 454 457 quant_av_k_l(k) = quant_av_k_l(k) + quant(k,j,i) * flag 455 458 count_2d_l(k) = count_2d_l(k) + INT( flag, KIND=iwp ) … … 479 482 + ( quant(k,j,i)  quant_av_k(k) )**2 & 480 483 * MERGE( 1.0_wp, 0.0_wp, & 481 BTEST( wall_flags_ 0(k,j,i), quant_type ) )484 BTEST( wall_flags_static_0(k,j,i), quant_type ) ) 482 485 ENDDO 483 486 ENDDO … … 719 722 DO j = nys, nyn 720 723 DO k = nzb, nzt+1 721 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), quant_type ) )724 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), quant_type ) ) 722 725 quant_av_l = quant_av_l + quant(k,j,i) * flag 723 726 average_count_l = average_count_l + INT( flag, KIND=iwp ) … … 771 774 + u(k,j,i) * dzw(k) * dy & 772 775 * MERGE( 1.0_wp, 0.0_wp, & 773 BTEST( wall_flags_ 0(k,j,i), 1 ) &776 BTEST( wall_flags_static_0(k,j,i), 1 ) & 774 777 ) 775 778 ENDDO … … 785 788  u(k,j,i) * dzw(k) * dy & 786 789 * MERGE( 1.0_wp, 0.0_wp, & 787 BTEST( wall_flags_ 0(k,j,i), 1 ) &790 BTEST( wall_flags_static_0(k,j,i), 1 ) & 788 791 ) 789 792 ENDDO … … 800 803 + v(k,j,i) * dzw(k) * dx & 801 804 * MERGE( 1.0_wp, 0.0_wp, & 802 BTEST( wall_flags_ 0(k,j,i), 2 ) &805 BTEST( wall_flags_static_0(k,j,i), 2 ) & 803 806 ) 804 807 ENDDO … … 814 817  v(k,j,i) * dzw(k) * dx & 815 818 * MERGE( 1.0_wp, 0.0_wp, & 816 BTEST( wall_flags_ 0(k,j,i), 2 ) &819 BTEST( wall_flags_static_0(k,j,i), 2 ) & 817 820 ) 818 821 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.