Changeset 4329 for palm/trunk/SOURCE/diffusion_v.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_v.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 ! … … 70 73 71 74 USE indices, & 72 ONLY: nxl, nxr, nyn, nysv, nzb, nzt, wall_flags_ 075 ONLY: nxl, nxr, nyn, nysv, nzb, nzt, wall_flags_static_0 73 76 74 77 USE kinds … … 101 104 !$ACC PRIVATE(surf_e, surf_s, flag, kmxm, kmxp, kmzm, kmzp) & 102 105 !$ACC PRIVATE(mask_bottom, mask_east, mask_west, mask_top) & 103 !$ACC PRESENT(wall_flags_ 0, km) &106 !$ACC PRESENT(wall_flags_static_0, km) & 104 107 !$ACC PRESENT(u, v, w) & 105 108 !$ACC PRESENT(ddzu, ddzw, drho_air, rho_air_zw) & … … 118 121 !-- It is sufficient to masked only east- and west-facing surfaces, which 119 122 !-- need special treatment for the v-component. 120 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 2 ) )121 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i+1), 2 ) )122 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i-1), 2 ) )123 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 2 ) ) 124 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 2 ) ) 125 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 2 ) ) 123 126 ! 124 127 !-- Interpolate eddy diffusivities on staggered gridpoints … … 192 195 !-- momentum flux at model top. 193 196 mask_bottom = MERGE( 1.0_wp, 0.0_wp, & 194 BTEST( wall_flags_ 0(k-1,j,i), 8 ) )197 BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 195 198 mask_top = MERGE( 1.0_wp, 0.0_wp, & 196 BTEST( wall_flags_ 0(k+1,j,i), 8 ) ) * &199 BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) * & 197 200 MERGE( 1.0_wp, 0.0_wp, & 198 BTEST( wall_flags_ 0(k+1,j,i), 9 ) )201 BTEST( wall_flags_static_0(k+1,j,i), 9 ) ) 199 202 flag = MERGE( 1.0_wp, 0.0_wp, & 200 BTEST( wall_flags_ 0(k,j,i), 2 ) )203 BTEST( wall_flags_static_0(k,j,i), 2 ) ) 201 204 ! 202 205 !-- Interpolate eddy diffusivities on staggered gridpoints … … 310 313 311 314 USE indices, & 312 ONLY: nzb, nzt, wall_flags_ 0315 ONLY: nzb, nzt, wall_flags_static_0 313 316 314 317 USE kinds … … 346 349 !-- It is sufficient to masked only east- and west-facing surfaces, which 347 350 !-- need special treatment for the v-component. 348 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 2 ) )349 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i+1), 2 ) )350 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i-1), 2 ) )351 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 2 ) ) 352 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 2 ) ) 353 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 2 ) ) 351 354 ! 352 355 !-- Interpolate eddy diffusivities on staggered gridpoints … … 416 419 !-- momentum flux at model top. 417 420 mask_bottom = MERGE( 1.0_wp, 0.0_wp, & 418 BTEST( wall_flags_ 0(k-1,j,i), 8 ) )421 BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 419 422 mask_top = MERGE( 1.0_wp, 0.0_wp, & 420 BTEST( wall_flags_ 0(k+1,j,i), 8 ) ) * &423 BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) * & 421 424 MERGE( 1.0_wp, 0.0_wp, & 422 BTEST( wall_flags_ 0(k+1,j,i), 9 ) )425 BTEST( wall_flags_static_0(k+1,j,i), 9 ) ) 423 426 flag = MERGE( 1.0_wp, 0.0_wp, & 424 BTEST( wall_flags_ 0(k,j,i), 2 ) )427 BTEST( wall_flags_static_0(k,j,i), 2 ) ) 425 428 ! 426 429 !-- Interpolate eddy diffusivities on staggered gridpoints
Note: See TracChangeset
for help on using the changeset viewer.