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/subsidence_mod.f90

    r4329 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! 4329 2019-12-10 15:46:36Z motisi
    2731! Renamed wall_flags_0 to wall_flags_static_0
    2832!
     
    155159       USE indices,                                                            &
    156160           ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,        &
    157                   wall_flags_static_0
     161                  wall_flags_total_0
    158162
    159163       USE kinds
     
    189193                              ( var(k+1,j,i) - var(k,j,i) ) * ddzu(k+1) *      &
    190194                                        MERGE( 1.0_wp, 0.0_wp,                 &
    191                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     195                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    192196                ELSE                               ! large-scale ascent
    193197                   tmp_tend = - w_subs(k) *                                    &
    194198                              ( var(k,j,i) - var(k-1,j,i) ) * ddzu(k) *        &
    195199                                        MERGE( 1.0_wp, 0.0_wp,                 &
    196                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     200                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    197201                ENDIF
    198202
     
    203207                                 * weight_substep(intermediate_timestep_count) &
    204208                                 * MERGE( 1.0_wp, 0.0_wp,                      &
    205                                           BTEST( wall_flags_static_0(k,j,i), 0 ) )
     209                                   BTEST( wall_flags_total_0(k,j,i), 0 ) )
    206210                ENDIF
    207211             ENDDO
     
    269273
    270274       USE indices,                                                            &
    271            ONLY:  nxl, nxlg, nxrg, nyng, nys, nysg, nzb, nzt, wall_flags_static_0
     275           ONLY:  nxl, nxlg, nxrg, nyng, nys, nysg, nzb, nzt,                  &
     276                  wall_flags_total_0
    272277
    273278       USE kinds
     
    300305                                    * ddzu(k+1)                                &
    301306                                    * MERGE( 1.0_wp, 0.0_wp,                   &
    302                                              BTEST( wall_flags_static_0(k,j,i), 0 ) )
     307                                      BTEST( wall_flags_total_0(k,j,i), 0 ) )
    303308          ELSE                                 ! large-scale ascent
    304309             tmp_tend = - w_subs(k) * ( var(k,j,i) - var(k-1,j,i) ) * ddzu(k)  &
    305310                                    * MERGE( 1.0_wp, 0.0_wp,                   &
    306                                              BTEST( wall_flags_static_0(k,j,i), 0 ) )
     311                                      BTEST( wall_flags_total_0(k,j,i), 0 ) )
    307312          ENDIF
    308313
     
    313318                                  * weight_substep(intermediate_timestep_count)&
    314319                                  * MERGE( 1.0_wp, 0.0_wp,                     &
    315                                            BTEST( wall_flags_static_0(k,j,i), 0 ) )
     320                                    BTEST( wall_flags_total_0(k,j,i), 0 ) )
    316321          ENDIF
    317322       ENDDO
Note: See TracChangeset for help on using the changeset viewer.