Changeset 2232 for palm/trunk/SOURCE/advec_v_up.f90
- Timestamp:
- May 30, 2017 5:47:52 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_v_up.f90
r2101 r2232 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! topography representation via flags 23 23 ! 24 24 ! Former revisions: … … 95 95 96 96 USE indices, & 97 ONLY: nxl, nxr, nyn, nysv, nzb _v_inner, nzt97 ONLY: nxl, nxr, nyn, nysv, nzb, nzt, wall_flags_0 98 98 99 99 USE kinds … … 113 113 DO i = nxl, nxr 114 114 DO j = nysv, nyn 115 DO k = nzb _v_inner(j,i)+1, nzt115 DO k = nzb+1, nzt 116 116 ! 117 117 !-- x-direction … … 120 120 IF ( ukomp > 0.0_wp ) THEN 121 121 tend(k,j,i) = tend(k,j,i) - ukomp * & 122 ( v(k,j,i) - v(k,j,i-1) ) * ddx 122 ( v(k,j,i) - v(k,j,i-1) ) * ddx & 123 * MERGE( 1.0_wp, 0.0_wp, & 124 BTEST( wall_flags_0(k,j,i), 2 ) ) 123 125 ELSE 124 126 tend(k,j,i) = tend(k,j,i) - ukomp * & 125 ( v(k,j,i+1) - v(k,j,i) ) * ddx 127 ( v(k,j,i+1) - v(k,j,i) ) * ddx & 128 * MERGE( 1.0_wp, 0.0_wp, & 129 BTEST( wall_flags_0(k,j,i), 2 ) ) 126 130 ENDIF 127 131 ! … … 130 134 IF ( vkomp > 0.0_wp ) THEN 131 135 tend(k,j,i) = tend(k,j,i) - vkomp * & 132 ( v(k,j,i) - v(k,j-1,i) ) * ddy 136 ( v(k,j,i) - v(k,j-1,i) ) * ddy & 137 * MERGE( 1.0_wp, 0.0_wp, & 138 BTEST( wall_flags_0(k,j,i), 2 ) ) 133 139 ELSE 134 140 tend(k,j,i) = tend(k,j,i) - vkomp * & 135 ( v(k,j+1,i) - v(k,j,i) ) * ddy 141 ( v(k,j+1,i) - v(k,j,i) ) * ddy & 142 * MERGE( 1.0_wp, 0.0_wp, & 143 BTEST( wall_flags_0(k,j,i), 2 ) ) 136 144 ENDIF 137 145 ! … … 141 149 IF ( wkomp > 0.0_wp ) THEN 142 150 tend(k,j,i) = tend(k,j,i) - wkomp * & 143 ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k) 151 ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k) & 152 * MERGE( 1.0_wp, 0.0_wp, & 153 BTEST( wall_flags_0(k,j,i), 2 ) ) 144 154 ELSE 145 155 tend(k,j,i) = tend(k,j,i) - wkomp * & 146 ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) 156 ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) & 157 * MERGE( 1.0_wp, 0.0_wp, & 158 BTEST( wall_flags_0(k,j,i), 2 ) ) 147 159 ENDIF 148 160 … … 171 183 172 184 USE indices, & 173 ONLY: nzb _v_inner, nzt185 ONLY: nzb, nzt, wall_flags_0 174 186 175 187 USE kinds … … 187 199 188 200 189 DO k = nzb _v_inner(j,i)+1, nzt201 DO k = nzb+1, nzt 190 202 ! 191 203 !-- x-direction … … 194 206 IF ( ukomp > 0.0_wp ) THEN 195 207 tend(k,j,i) = tend(k,j,i) - ukomp * & 196 ( v(k,j,i) - v(k,j,i-1) ) * ddx 208 ( v(k,j,i) - v(k,j,i-1) ) * ddx & 209 * MERGE( 1.0_wp, 0.0_wp, & 210 BTEST( wall_flags_0(k,j,i), 2 ) ) 197 211 ELSE 198 212 tend(k,j,i) = tend(k,j,i) - ukomp * & 199 ( v(k,j,i+1) - v(k,j,i) ) * ddx 213 ( v(k,j,i+1) - v(k,j,i) ) * ddx & 214 * MERGE( 1.0_wp, 0.0_wp, & 215 BTEST( wall_flags_0(k,j,i), 2 ) ) 200 216 ENDIF 201 217 ! … … 204 220 IF ( vkomp > 0.0_wp ) THEN 205 221 tend(k,j,i) = tend(k,j,i) - vkomp * & 206 ( v(k,j,i) - v(k,j-1,i) ) * ddy 222 ( v(k,j,i) - v(k,j-1,i) ) * ddy & 223 * MERGE( 1.0_wp, 0.0_wp, & 224 BTEST( wall_flags_0(k,j,i), 2 ) ) 207 225 ELSE 208 226 tend(k,j,i) = tend(k,j,i) - vkomp * & 209 ( v(k,j+1,i) - v(k,j,i) ) * ddy 227 ( v(k,j+1,i) - v(k,j,i) ) * ddy & 228 * MERGE( 1.0_wp, 0.0_wp, & 229 BTEST( wall_flags_0(k,j,i), 2 ) ) 210 230 ENDIF 211 231 ! … … 214 234 IF ( wkomp > 0.0_wp ) THEN 215 235 tend(k,j,i) = tend(k,j,i) - wkomp * & 216 ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k) 236 ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k) & 237 * MERGE( 1.0_wp, 0.0_wp, & 238 BTEST( wall_flags_0(k,j,i), 2 ) ) 217 239 ELSE 218 240 tend(k,j,i) = tend(k,j,i) - wkomp * & 219 ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) 241 ( v(k+1,j,i) - v(k,j,i) ) * ddzu(k+1) & 242 * MERGE( 1.0_wp, 0.0_wp, & 243 BTEST( wall_flags_0(k,j,i), 2 ) ) 220 244 ENDIF 221 245
Note: See TracChangeset
for help on using the changeset viewer.