Ignore:
Timestamp:
Dec 10, 2019 3:46:36 PM (5 years ago)
Author:
motisi
Message:

Renamed wall_flags_0 to wall_flags_static_0

File:
1 edited

Legend:

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

    r4282 r4329  
    2525! -----------------
    2626! $Id$
     27! Renamed wall_flags_0 to wall_flags_static_0
     28!
     29! 4282 2019-10-29 16:18:46Z schwenkel
    2730! Bugfix of particle timeseries in case of more than one particle group
    2831!
     
    145148               nzb_max, nzt,nbgp, ngp_2dh_outer,                               &
    146149               topo_top_ind,                                                   &
    147                wall_flags_0
     150               wall_flags_static_0
    148151
    149152    USE kinds
     
    13521355!--                               of overhanging structures.
    13531356                                  IF ( kp > nzt  .OR.                          &
    1354                                  .NOT. BTEST( wall_flags_0(kp,jp,ip), 0 ) )  THEN
     1357                                 .NOT. BTEST( wall_flags_static_0(kp,jp,ip), 0 ) )  THEN
    13551358                                     pos_x = pos_x + pdx(i)
    13561359                                     CYCLE xloop
     
    13621365                               ELSEIF ( .NOT. seed_follows_topography .AND.    &
    13631366                                         tmp_particle%z <= zw(k_surf)  .OR.    &
    1364                                         .NOT. BTEST( wall_flags_0(kp,jp,ip), 0 ) )&
     1367                                        .NOT. BTEST( wall_flags_static_0(kp,jp,ip), 0 ) )&
    13651368                               THEN
    13661369                                  pos_x = pos_x + pdx(i)
     
    15181521!
    15191522!--                Check if particle is within topography
    1520                    IF ( .NOT. BTEST( wall_flags_0(k,j,i), 0 ) )  THEN
     1523                   IF ( .NOT. BTEST( wall_flags_static_0(k,j,i), 0 ) )  THEN
    15211524                      particles(n)%particle_mask = .FALSE.
    15221525                      deleted_particles = deleted_particles + 1
     
    17611764          DO  k = nzb, nzt+1
    17621765
    1763              IF ( .NOT. BTEST( wall_flags_0(k,j,i-1), 0 )  .AND.               &
    1764                         BTEST( wall_flags_0(k,j,i), 0   )  .AND.               &
    1765                         BTEST( wall_flags_0(k,j,i+1), 0 ) )                    &
     1766             IF ( .NOT. BTEST( wall_flags_static_0(k,j,i-1), 0 )  .AND.               &
     1767                        BTEST( wall_flags_static_0(k,j,i), 0   )  .AND.               &
     1768                        BTEST( wall_flags_static_0(k,j,i+1), 0 ) )                    &
    17661769             THEN
    17671770                de_dx(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    17681771                               ( e(k,j,i+1) - e(k,j,i) ) * ddx
    1769              ELSEIF ( BTEST( wall_flags_0(k,j,i-1), 0 )  .AND.                 &
    1770                       BTEST( wall_flags_0(k,j,i), 0   )  .AND.                 &
    1771                 .NOT. BTEST( wall_flags_0(k,j,i+1), 0 ) )                      &
     1772             ELSEIF ( BTEST( wall_flags_static_0(k,j,i-1), 0 )  .AND.                 &
     1773                      BTEST( wall_flags_static_0(k,j,i), 0   )  .AND.                 &
     1774                .NOT. BTEST( wall_flags_static_0(k,j,i+1), 0 ) )                      &
    17721775             THEN
    17731776                de_dx(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    17741777                               ( e(k,j,i) - e(k,j,i-1) ) * ddx
    1775              ELSEIF ( .NOT. BTEST( wall_flags_0(k,j,i), 22   )  .AND.          &
    1776                       .NOT. BTEST( wall_flags_0(k,j,i+1), 22 ) )               &   
     1778             ELSEIF ( .NOT. BTEST( wall_flags_static_0(k,j,i), 22   )  .AND.          &
     1779                      .NOT. BTEST( wall_flags_static_0(k,j,i+1), 22 ) )               &   
    17771780             THEN
    17781781                de_dx(k,j,i) = 0.0_wp
    1779              ELSEIF ( .NOT. BTEST( wall_flags_0(k,j,i-1), 22 )  .AND.          &
    1780                       .NOT. BTEST( wall_flags_0(k,j,i), 22   ) )               &
     1782             ELSEIF ( .NOT. BTEST( wall_flags_static_0(k,j,i-1), 22 )  .AND.          &
     1783                      .NOT. BTEST( wall_flags_static_0(k,j,i), 22   ) )               &
    17811784             THEN
    17821785                de_dx(k,j,i) = 0.0_wp
     
    17851788             ENDIF
    17861789
    1787              IF ( .NOT. BTEST( wall_flags_0(k,j-1,i), 0 )  .AND.               &
    1788                         BTEST( wall_flags_0(k,j,i), 0   )  .AND.               &
    1789                         BTEST( wall_flags_0(k,j+1,i), 0 ) )                    &
     1790             IF ( .NOT. BTEST( wall_flags_static_0(k,j-1,i), 0 )  .AND.               &
     1791                        BTEST( wall_flags_static_0(k,j,i), 0   )  .AND.               &
     1792                        BTEST( wall_flags_static_0(k,j+1,i), 0 ) )                    &
    17901793             THEN
    17911794                de_dy(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    17921795                               ( e(k,j+1,i) - e(k,j,i) ) * ddy
    1793              ELSEIF ( BTEST( wall_flags_0(k,j-1,i), 0 )  .AND.                 &
    1794                       BTEST( wall_flags_0(k,j,i), 0   )  .AND.                 &
    1795                 .NOT. BTEST( wall_flags_0(k,j+1,i), 0 ) )                      &
     1796             ELSEIF ( BTEST( wall_flags_static_0(k,j-1,i), 0 )  .AND.                 &
     1797                      BTEST( wall_flags_static_0(k,j,i), 0   )  .AND.                 &
     1798                .NOT. BTEST( wall_flags_static_0(k,j+1,i), 0 ) )                      &
    17961799             THEN
    17971800                de_dy(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    17981801                               ( e(k,j,i) - e(k,j-1,i) ) * ddy
    1799              ELSEIF ( .NOT. BTEST( wall_flags_0(k,j,i), 22   )  .AND.          &
    1800                       .NOT. BTEST( wall_flags_0(k,j+1,i), 22 ) )               &   
     1802             ELSEIF ( .NOT. BTEST( wall_flags_static_0(k,j,i), 22   )  .AND.          &
     1803                      .NOT. BTEST( wall_flags_static_0(k,j+1,i), 22 ) )               &   
    18011804             THEN
    18021805                de_dy(k,j,i) = 0.0_wp
    1803              ELSEIF ( .NOT. BTEST( wall_flags_0(k,j-1,i), 22 )  .AND.          &
    1804                       .NOT. BTEST( wall_flags_0(k,j,i), 22   ) )               &
     1806             ELSEIF ( .NOT. BTEST( wall_flags_static_0(k,j-1,i), 22 )  .AND.          &
     1807                      .NOT. BTEST( wall_flags_static_0(k,j,i), 22   ) )               &
    18051808             THEN
    18061809                de_dy(k,j,i) = 0.0_wp
     
    18201823!
    18211824!--          Flag to mask topography
    1822              flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0  ) )
     1825             flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0  ) )
    18231826
    18241827             de_dz(k,j,i) = 2.0_wp * sgs_wf_part *                             &
     
    18951898!
    18961899!--             Flag indicating vicinity of wall
    1897                 flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 24 ) )
     1900                flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 24 ) )
    18981901
    18991902                sums_l(k,1,0)  = sums_l(k,1,0)  + u(k,j,i) * flag1
     
    19351938!
    19361939!--             Flag indicating vicinity of wall
    1937                 flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 24 ) )
     1940                flag1 = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 24 ) )
    19381941
    19391942                sums_l(k,8,0)  = sums_l(k,8,0)  + e(k,j,i)                       * flag1
     
    36663669             j = jp + MERGE( -1_iwp , 1_iwp, BTEST( nb, 1 ) )
    36673670             k = kp + MERGE( -1_iwp , 1_iwp, BTEST( nb, 0 ) )
    3668              IF ( .NOT. BTEST(wall_flags_0(k,  jp, ip), 0) .OR.                &
    3669                   .NOT. BTEST(wall_flags_0(kp, j,  ip), 0) .OR.                &
    3670                   .NOT. BTEST(wall_flags_0(kp, jp, i ), 0) )                   &
     3671             IF ( .NOT. BTEST(wall_flags_static_0(k,  jp, ip), 0) .OR.                &
     3672                  .NOT. BTEST(wall_flags_static_0(kp, j,  ip), 0) .OR.                &
     3673                  .NOT. BTEST(wall_flags_static_0(kp, jp, i ), 0) )                   &
    36713674             THEN
    36723675                subbox_at_wall = .TRUE.
     
    47634766                IF ( reach_x(t_index)                      .AND.               &
    47644767                     ABS( pos_x - xwall ) < eps            .AND.               &
    4765                      .NOT. BTEST(wall_flags_0(k3,j3,i3),0) .AND.               &
     4768                     .NOT. BTEST(wall_flags_static_0(k3,j3,i3),0) .AND.               &
    47664769                     .NOT. reflect_x )  THEN
    47674770!
     
    48034806                IF ( reach_y(t_index)                      .AND.               &
    48044807                     ABS( pos_y - ywall ) < eps            .AND.               &
    4805                      .NOT. BTEST(wall_flags_0(k3,j3,i3),0) .AND.               &
     4808                     .NOT. BTEST(wall_flags_static_0(k3,j3,i3),0) .AND.               &
    48064809                     .NOT. reflect_y )  THEN
    48074810!
     
    48434846                IF ( reach_z(t_index)                      .AND.               &
    48444847                     ABS( pos_z - zwall ) < eps            .AND.               &
    4845                      .NOT. BTEST(wall_flags_0(k3,j3,i3),0) .AND.               &
     4848                     .NOT. BTEST(wall_flags_static_0(k3,j3,i3),0) .AND.               &
    48464849                     .NOT. reflect_z )  THEN
    48474850!
     
    52015204!
    52025205!--          Predetermine flag to mask topography
    5203              flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
     5206             flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0 ) )
    52045207
    52055208             q(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i) * flag
     
    52315234!
    52325235!--    Predetermine flag to mask topography
    5233        flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
     5236       flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i), 0 ) )
    52345237
    52355238       q(k,j,i)  = q(k,j,i)  - ql_c(k,j,i) * flag
Note: See TracChangeset for help on using the changeset viewer.