Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (4 years ago)
Author:
motisi
Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_mask.f90

    r4331 r4346  
    2525! -----------------
    2626! $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
    2731! Formatting adjustment
    2832!
     
    97101
    98102    USE indices,                                                               &
    99         ONLY:  nbgp, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_static_0
     103        ONLY:  nbgp, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_total_0
    100104
    101105    USE kinds
     
    267271                         im = mask_i(mid,i)
    268272                         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 )),&
    270274                                       DIM = 1 ) - 1
    271275                         DO  k = 1, mask_size_l(mid,3)
    272276                            kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    273277!--                         Set value if not in building
    274                             IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     278                            IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    275279                               local_pf(i,j,k) = fill_value
    276280                            ELSE
     
    336340                         im = mask_i(mid,i)
    337341                         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
    339344                         DO  k = 1, mask_size_l(mid,3)
    340345                            kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    341346!--                         Set value if not in building
    342                             IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     347                            IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    343348                               local_pf(i,j,k) = fill_value
    344349                            ELSE
     
    379384                            im = mask_i(mid,i)
    380385                            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
    382388                            DO  k = 1, mask_size_l(mid,3)
    383389                               kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    384390!--                            Set value if not in building
    385                                IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     391                               IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    386392                                  local_pf(i,j,k) = fill_value
    387393                               ELSE
     
    473479                         im = mask_i(mid,i)
    474480                         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
    476483                         DO  k = 1, mask_size_l(mid,3)
    477484                            kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    478485!--                         Set value if not in building
    479                             IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     486                            IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    480487                               local_pf(i,j,k) = fill_value
    481488                            ELSE
     
    512519                         im = mask_i(mid,i)
    513520                         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
    515523                         DO  k = 1, mask_size_l(mid,3)
    516524                            kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    517525!--                         Set value if not in building
    518                             IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     526                            IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    519527                               local_pf(i,j,k) = fill_value
    520528                            ELSE
     
    654662                   im = mask_i(mid,i)
    655663                   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
    657666                   DO  k = 1, mask_size_l(mid,3)
    658667                      kk = MIN( ktt+mask_k(mid,k), nzt+1 )
    659668!--                   Set value if not in building
    660                       IF ( BTEST( wall_flags_static_0(kk,jm,im), 6 ) )  THEN
     669                      IF ( BTEST( wall_flags_total_0(kk,jm,im), 6 ) )  THEN
    661670                         local_pf(i,j,k) = fill_value
    662671                      ELSE
Note: See TracChangeset for help on using the changeset viewer.