Changeset 3255 for palm


Ignore:
Timestamp:
Sep 17, 2018 11:57:36 AM (6 years ago)
Author:
suehring
Message:

Missing wall flags added

File:
1 edited

Legend:

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

    r2718 r3255  
    2525! -----------------
    2626! $Id$
     27! Missing wall flags added
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    212215          IF ( ukomp > 0.0_wp )  THEN
    213216             tend(k,j,i) = tend(k,j,i) - ukomp *                               &
    214                                          ( w(k,j,i) - w(k,j,i-1) ) * ddx
     217                                         ( w(k,j,i) - w(k,j,i-1) ) * ddx       &
     218                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     219                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    215220          ELSE
    216221             tend(k,j,i) = tend(k,j,i) - ukomp *                               &
    217                                          ( w(k,j,i+1) - w(k,j,i) ) * ddx
     222                                         ( w(k,j,i+1) - w(k,j,i) ) * ddx       &
     223                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     224                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    218225          ENDIF
    219226!
     
    223230          IF ( vkomp > 0.0_wp )  THEN
    224231             tend(k,j,i) = tend(k,j,i) - vkomp *                               &
    225                                          ( w(k,j,i) - w(k,j-1,i) ) * ddy
     232                                         ( w(k,j,i) - w(k,j-1,i) ) * ddy       &
     233                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     234                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    226235          ELSE
    227236             tend(k,j,i) = tend(k,j,i) - vkomp *                               &
    228                                          ( w(k,j+1,i) - w(k,j,i) ) * ddy
     237                                         ( w(k,j+1,i) - w(k,j,i) ) * ddy       &
     238                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     239                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    229240          ENDIF
    230241!
     
    232243          IF ( w(k,j,i) > 0.0_wp )  THEN
    233244             tend(k,j,i) = tend(k,j,i) - w(k,j,i) *                            &
    234                                          ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k)
     245                                         ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k)   &
     246                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     247                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    235248          ELSE
    236249             tend(k,j,i) = tend(k,j,i) - w(k,j,i) *                            &
    237                                          ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)
     250                                         ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
     251                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     252                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    238253          ENDIF
    239254
Note: See TracChangeset for help on using the changeset viewer.