Changeset 4346 for palm/trunk/SOURCE/data_output_mask.f90
- Timestamp:
- Dec 18, 2019 11:55:56 AM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_mask.f90
r4331 r4346 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Introduction of wall_flags_total_0, which currently sets bits based on static 28 ! topography information used in wall_flags_static_0 29 ! 30 ! 4331 2019-12-10 18:25:02Z suehring 27 31 ! Formatting adjustment 28 32 ! … … 97 101 98 102 USE indices, & 99 ONLY: nbgp, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_ static_0103 ONLY: nbgp, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_total_0 100 104 101 105 USE kinds … … 267 271 im = mask_i(mid,i) 268 272 jm = mask_j(mid,j) 269 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_ static_0(:,jm,im), 5 )),&273 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )),& 270 274 DIM = 1 ) - 1 271 275 DO k = 1, mask_size_l(mid,3) 272 276 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 273 277 !-- Set value if not in building 274 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN278 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 275 279 local_pf(i,j,k) = fill_value 276 280 ELSE … … 336 340 im = mask_i(mid,i) 337 341 jm = mask_j(mid,j) 338 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), DIM = 1 ) - 1 342 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), & 343 DIM = 1 ) - 1 339 344 DO k = 1, mask_size_l(mid,3) 340 345 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 341 346 !-- Set value if not in building 342 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN347 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 343 348 local_pf(i,j,k) = fill_value 344 349 ELSE … … 379 384 im = mask_i(mid,i) 380 385 jm = mask_j(mid,j) 381 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), DIM = 1 ) - 1 386 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), & 387 DIM = 1 ) - 1 382 388 DO k = 1, mask_size_l(mid,3) 383 389 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 384 390 !-- Set value if not in building 385 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN391 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 386 392 local_pf(i,j,k) = fill_value 387 393 ELSE … … 473 479 im = mask_i(mid,i) 474 480 jm = mask_j(mid,j) 475 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), DIM = 1 ) - 1 481 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), & 482 DIM = 1 ) - 1 476 483 DO k = 1, mask_size_l(mid,3) 477 484 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 478 485 !-- Set value if not in building 479 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN486 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 480 487 local_pf(i,j,k) = fill_value 481 488 ELSE … … 512 519 im = mask_i(mid,i) 513 520 jm = mask_j(mid,j) 514 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), DIM = 1 ) - 1 521 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), & 522 DIM = 1 ) - 1 515 523 DO k = 1, mask_size_l(mid,3) 516 524 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 517 525 !-- Set value if not in building 518 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN526 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 519 527 local_pf(i,j,k) = fill_value 520 528 ELSE … … 654 662 im = mask_i(mid,i) 655 663 jm = mask_j(mid,j) 656 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_static_0(:,jm,im), 5 )), DIM = 1 ) - 1 664 ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,jm,im), 5 )), & 665 DIM = 1 ) - 1 657 666 DO k = 1, mask_size_l(mid,3) 658 667 kk = MIN( ktt+mask_k(mid,k), nzt+1 ) 659 668 !-- Set value if not in building 660 IF ( BTEST( wall_flags_ static_0(kk,jm,im), 6 ) ) THEN669 IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) ) THEN 661 670 local_pf(i,j,k) = fill_value 662 671 ELSE
Note: See TracChangeset
for help on using the changeset viewer.