Ignore:
Timestamp:
Nov 20, 2018 10:55:41 AM (6 years ago)
Author:
suehring
Message:

Remove unnecessary double masking of topography in advection and buoyancy terms

File:
1 edited

Legend:

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

    r3302 r3538  
    2525! -----------------
    2626! $Id$
     27! Remove unnecessary double-masking of topography
     28!
     29! 3302 2018-10-03 02:39:40Z raasch
    2730! Stokes drift velocity added
    2831!
     
    116119
    117120       USE indices,                                                            &
    118            ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb,             &
    119                   nzt, wall_flags_0
     121           ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
    120122
    121123       USE kinds
     
    163165                                                   ( sk(k-1,j,i) - sk(k,j,i) ) &
    164166                                                 ) * dd2zu(k)                  &
    165                                       ) * MERGE( 1.0_wp, 0.0_wp,               &
    166                                                BTEST( wall_flags_0(k,j,i), 0 ) )
     167                                      )
    167168             ENDDO
    168169          ENDDO
     
    189190
    190191       USE indices,                                                            &
    191            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
     192           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt
    192193
    193194       USE kinds
     
    233234                                                  ( sk(k-1,j,i) - sk(k,j,i) )  &
    234235                                                ) * dd2zu(k)                   &
    235                                     ) * MERGE( 1.0_wp, 0.0_wp,                 &
    236                                                BTEST( wall_flags_0(k,j,i), 0 ) )
     236                                    )
    237237       ENDDO
    238238
Note: See TracChangeset for help on using the changeset viewer.