Changeset 3538 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Nov 20, 2018 10:55:41 AM (5 years ago)

Author:

suehring

Message:

Remove unnecessary double masking of topography in advection and buoyancy terms

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Comment added
+! 
+! 3294 2018-10-01 02:37:10Z raasch
 ! ocean renamed ocean_mode
 ! 
@@ -1647 +1650 @@
           DO  j = nys, nyn
              DO  k = nzb, nzt
+!
+!--             Flag topography for all grid points which are below
+!--             the local topography height. 
+!--             Note, each topography is flagged as building. 
                 IF ( zu(k) - ocean_offset <= buildings_f%var_2d(j,i) )  THEN
                     topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
-                    topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates terrain
+                    topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates building
                 ENDIF
              ENDDO