Changeset 4329 for palm/trunk/SOURCE/diagnostic_output_quantities_mod.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diagnostic_output_quantities_mod.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 ! Corrected "Former revisions" section 28 31 ! … … 92 95 ! 93 96 USE indices, & 94 ONLY: nxl, nxr, nyn, nys, nzb, nzt, wall_flags_ 097 ONLY: nxl, nxr, nyn, nys, nzb, nzt, wall_flags_static_0 95 98 ! 96 99 USE kinds … … 582 585 local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i), & 583 586 REAL( fill_value, KIND = wp ), & 584 BTEST( wall_flags_ 0(k,j,i), flag_nr ) )587 BTEST( wall_flags_static_0(k,j,i), flag_nr ) ) 585 588 ENDDO 586 589 ENDDO … … 685 688 local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i), & 686 689 REAL( fill_value, KIND = wp ), & 687 BTEST( wall_flags_ 0(k,j,i), flag_nr ) )690 BTEST( wall_flags_static_0(k,j,i), flag_nr ) ) 688 691 ENDDO 689 692 ENDDO … … 802 805 im = mask_i(mid,i) 803 806 jm = mask_j(mid,j) 804 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_ 0(:,jm,im), 5 )), &807 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), & 805 808 DIM = 1 ) - 1 806 809 DO k = 1, mask_size_l(mid,3) … … 808 811 ! 809 812 !-- Set value if not in building 810 IF ( BTEST( wall_flags_ 0(kk,jm,im), 6 ) ) THEN813 IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) ) THEN 811 814 local_pf(i,j,k) = fill_value 812 815 ELSE … … 933 936 - ( u(k,j+1,i) + u(k,j+1,i+1) & 934 937 - u(k,j-1,i) - u(k,j-1,i+1) ) * ddy )**2 ) & 935 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 0) )938 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0) ) 936 939 ENDDO 937 940 ENDDO … … 944 947 DO k = nzb+1, nzt 945 948 uu(k,j,i) = u(k,j,i) * u(k,j,i) & 946 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 1) )949 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 1) ) 947 950 ENDDO 948 951 ENDDO … … 955 958 DO k = nzb+1, nzt 956 959 vv(k,j,i) = v(k,j,i) * v(k,j,i) & 957 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 2) )960 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 2) ) 958 961 ENDDO 959 962 ENDDO … … 966 969 DO k = nzb+1, nzt-1 967 970 ww(k,j,i) = w(k,j,i) * w(k,j,i) & 968 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 3) )971 * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 3) ) 969 972 ENDDO 970 973 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.