Changeset 4329 for palm/trunk/SOURCE/vdi_internal_controls.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (3 years 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 2019-08-22 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.