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

Timestamp:

Jul 30, 2020 2:54:34 PM (4 years ago)

Author:

suehring

Message:

Land-surface model: bugfix in level 3 initialization of root-area-density; Avoid double classifiation of vertical walls (at surfaces that are alo covered by buildings); Land/urban surface: bugfix in resistance calculation - avoid potential divisions by zero; init_grid: in case of ASCII topography flag grid points as terrain and building to allow application of land/urban-surface model

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! In case of ASCII topography input flag grid points as terrain and building.
+!
+! 4601 2020-07-14 12:06:09Z suehring
 ! Minor formatting adjustments
 ! 
@@ -1347 +1350 @@
 !--             Flag topography for all grid points which are below
 !--             the local topography height. 
-!--             Note, each topography is flagged as building. 
+!--             Note, each topography is flagged as building (bit 2) as well as 
+!--             terrain (bit 1) in order to employ urban-surface as well as 
+!--             land-surface model.
                 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 building
+                   topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 1 )
+                   topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 )
                 ENDIF
              ENDDO