Changeset 3538 for palm/trunk/SOURCE/advec_u_up.f90
- Timestamp:
- Nov 20, 2018 10:55:41 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_u_up.f90
r2718 r3538 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Remove unnecessary double-masking of topography 28 ! 29 ! 2718 2018-01-02 08:49:38Z maronga 27 30 ! Corrected "Former revisions" section 28 31 ! … … 104 107 105 108 USE indices, & 106 ONLY: nxlu, nxr, nyn, nys, nzb, nzt , wall_flags_0109 ONLY: nxlu, nxr, nyn, nys, nzb, nzt 107 110 108 111 USE kinds … … 128 131 IF ( ukomp > 0.0_wp ) THEN 129 132 tend(k,j,i) = tend(k,j,i) - ukomp * & 130 ( u(k,j,i) - u(k,j,i-1) ) * ddx & 131 * MERGE( 1.0_wp, 0.0_wp, & 132 BTEST( wall_flags_0(k,j,i), 1 ) ) 133 ( u(k,j,i) - u(k,j,i-1) ) * ddx 133 134 ELSE 134 135 tend(k,j,i) = tend(k,j,i) - ukomp * & 135 ( u(k,j,i+1) - u(k,j,i) ) * ddx & 136 * MERGE( 1.0_wp, 0.0_wp, & 137 BTEST( wall_flags_0(k,j,i), 1 ) ) 136 ( u(k,j,i+1) - u(k,j,i) ) * ddx 138 137 ENDIF 139 138 ! … … 143 142 IF ( vkomp > 0.0_wp ) THEN 144 143 tend(k,j,i) = tend(k,j,i) - vkomp * & 145 ( u(k,j,i) - u(k,j-1,i) ) * ddy & 146 * MERGE( 1.0_wp, 0.0_wp, & 147 BTEST( wall_flags_0(k,j,i), 1 ) ) 144 ( u(k,j,i) - u(k,j-1,i) ) * ddy 148 145 ELSE 149 146 tend(k,j,i) = tend(k,j,i) - vkomp * & 150 ( u(k,j+1,i) - u(k,j,i) ) * ddy & 151 * MERGE( 1.0_wp, 0.0_wp, & 152 BTEST( wall_flags_0(k,j,i), 1 ) ) 147 ( u(k,j+1,i) - u(k,j,i) ) * ddy 153 148 ENDIF 154 149 ! … … 158 153 IF ( wkomp > 0.0_wp ) THEN 159 154 tend(k,j,i) = tend(k,j,i) - wkomp * & 160 ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k) & 161 * MERGE( 1.0_wp, 0.0_wp, & 162 BTEST( wall_flags_0(k,j,i), 1 ) ) 155 ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k) 163 156 ELSE 164 157 tend(k,j,i) = tend(k,j,i) - wkomp * & 165 ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) & 166 * MERGE( 1.0_wp, 0.0_wp, & 167 BTEST( wall_flags_0(k,j,i), 1 ) ) 158 ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) 168 159 ENDIF 169 160 … … 192 183 193 184 USE indices, & 194 ONLY: nzb, nzt , wall_flags_0185 ONLY: nzb, nzt 195 186 196 187 USE kinds … … 214 205 IF ( ukomp > 0.0_wp ) THEN 215 206 tend(k,j,i) = tend(k,j,i) - ukomp * & 216 ( u(k,j,i) - u(k,j,i-1) ) * ddx & 217 * MERGE( 1.0_wp, 0.0_wp, & 218 BTEST( wall_flags_0(k,j,i), 1 ) ) 207 ( u(k,j,i) - u(k,j,i-1) ) * ddx 219 208 ELSE 220 209 tend(k,j,i) = tend(k,j,i) - ukomp * & 221 ( u(k,j,i+1) - u(k,j,i) ) * ddx & 222 * MERGE( 1.0_wp, 0.0_wp, & 223 BTEST( wall_flags_0(k,j,i), 1 ) ) 210 ( u(k,j,i+1) - u(k,j,i) ) * ddx 224 211 ENDIF 225 212 ! … … 229 216 IF ( vkomp > 0.0_wp ) THEN 230 217 tend(k,j,i) = tend(k,j,i) - vkomp * & 231 ( u(k,j,i) - u(k,j-1,i) ) * ddy & 232 * MERGE( 1.0_wp, 0.0_wp, & 233 BTEST( wall_flags_0(k,j,i), 1 ) ) 218 ( u(k,j,i) - u(k,j-1,i) ) * ddy 234 219 ELSE 235 220 tend(k,j,i) = tend(k,j,i) - vkomp * & 236 ( u(k,j+1,i) - u(k,j,i) ) * ddy & 237 * MERGE( 1.0_wp, 0.0_wp, & 238 BTEST( wall_flags_0(k,j,i), 1 ) ) 221 ( u(k,j+1,i) - u(k,j,i) ) * ddy 239 222 ENDIF 240 223 ! … … 243 226 IF ( wkomp > 0.0_wp ) THEN 244 227 tend(k,j,i) = tend(k,j,i) - wkomp * & 245 ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k) & 246 * MERGE( 1.0_wp, 0.0_wp, & 247 BTEST( wall_flags_0(k,j,i), 1 ) ) 228 ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k) 248 229 ELSE 249 230 tend(k,j,i) = tend(k,j,i) - wkomp * & 250 ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) & 251 * MERGE( 1.0_wp, 0.0_wp, & 252 BTEST( wall_flags_0(k,j,i), 1 ) ) 231 ( u(k+1,j,i) - u(k,j,i) ) * ddzu(k+1) 253 232 ENDIF 254 233
Note: See TracChangeset
for help on using the changeset viewer.