Changeset 4329 for palm/trunk/SOURCE/diffusion_s.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_s.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 ! … … 83 86 84 87 USE indices, & 85 ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt, wall_flags_ 088 ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt, wall_flags_static_0 86 89 87 90 USE kinds … … 132 135 !$ACC PRIVATE(surf_e, surf_s, flag, mask_top, mask_bottom) & 133 136 !$ACC PRIVATE(mask_north, mask_south, mask_west, mask_east) & 134 !$ACC PRESENT(wall_flags_ 0, kh) &137 !$ACC PRESENT(wall_flags_static_0, kh) & 135 138 !$ACC PRESENT(s) & 136 139 !$ACC PRESENT(ddzu, ddzw, drho_air, rho_air_zw) & … … 156 159 ! 157 160 !-- Predetermine flag to mask topography and wall-bounded grid points 158 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 0 ) )161 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0 ) ) 159 162 ! 160 163 !-- Predetermine flag to mask wall-bounded grid points, equivalent to 161 164 !-- former s_outer array 162 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i-1), 0 ) )163 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i+1), 0 ) )164 mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j-1,i), 0 ) )165 mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j+1,i), 0 ) )165 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 0 ) ) 166 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 0 ) ) 167 mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 0 ) ) 168 mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 0 ) ) 166 169 167 170 tend(k,j,i) = tend(k,j,i) & … … 297 300 !-- flux at model top. 298 301 mask_bottom = MERGE( 1.0_wp, 0.0_wp, & 299 BTEST( wall_flags_ 0(k-1,j,i), 8 ) )302 BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 300 303 mask_top = MERGE( 1.0_wp, 0.0_wp, & 301 BTEST( wall_flags_ 0(k+1,j,i), 8 ) ) * &304 BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) * & 302 305 MERGE( 1.0_wp, 0.0_wp, & 303 BTEST( wall_flags_ 0(k+1,j,i), 9 ) )306 BTEST( wall_flags_static_0(k+1,j,i), 9 ) ) 304 307 flag = MERGE( 1.0_wp, 0.0_wp, & 305 BTEST( wall_flags_ 0(k,j,i), 0 ) )308 BTEST( wall_flags_static_0(k,j,i), 0 ) ) 306 309 307 310 tend(k,j,i) = tend(k,j,i) & … … 412 415 413 416 USE indices, & 414 ONLY: nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_ 0417 ONLY: nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_static_0 415 418 416 419 USE kinds … … 462 465 ! 463 466 !-- Predetermine flag to mask topography and wall-bounded grid points 464 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i), 0 ) )467 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0 ) ) 465 468 ! 466 469 !-- Predetermine flag to mask wall-bounded grid points, equivalent to 467 470 !-- former s_outer array 468 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i-1), 0 ) )469 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j,i+1), 0 ) )470 mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j-1,i), 0 ) )471 mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_ 0(k,j+1,i), 0 ) )471 mask_west = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 0 ) ) 472 mask_east = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 0 ) ) 473 mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 0 ) ) 474 mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 0 ) ) 472 475 ! 473 476 !-- Finally, determine flag to mask both topography itself as well … … 607 610 !-- flux at model top. 608 611 mask_bottom = MERGE( 1.0_wp, 0.0_wp, & 609 BTEST( wall_flags_ 0(k-1,j,i), 8 ) )612 BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 610 613 mask_top = MERGE( 1.0_wp, 0.0_wp, & 611 BTEST( wall_flags_ 0(k+1,j,i), 8 ) ) * &614 BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) * & 612 615 MERGE( 1.0_wp, 0.0_wp, & 613 BTEST( wall_flags_ 0(k+1,j,i), 9 ) )616 BTEST( wall_flags_static_0(k+1,j,i), 9 ) ) 614 617 flag = MERGE( 1.0_wp, 0.0_wp, & 615 BTEST( wall_flags_ 0(k,j,i), 0 ) )618 BTEST( wall_flags_static_0(k,j,i), 0 ) ) 616 619 617 620 tend(k,j,i) = tend(k,j,i) &
Note: See TracChangeset
for help on using the changeset viewer.