Ignore:
Timestamp:
Nov 20, 2018 10:55:41 AM (3 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/init_grid.f90

    r3294 r3538  
    2525! -----------------
    2626! $Id$
     27! Comment added
     28!
     29! 3294 2018-10-01 02:37:10Z raasch
    2730! ocean renamed ocean_mode
    2831!
     
    16471650          DO  j = nys, nyn
    16481651             DO  k = nzb, nzt
     1652!
     1653!--             Flag topography for all grid points which are below
     1654!--             the local topography height.
     1655!--             Note, each topography is flagged as building.
    16491656                IF ( zu(k) - ocean_offset <= buildings_f%var_2d(j,i) )  THEN
    16501657                    topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
    1651                     topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates terrain
     1658                    topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates building
    16521659                ENDIF
    16531660             ENDDO
Note: See TracChangeset for help on using the changeset viewer.