Changeset 2747

Timestamp:

Jan 15, 2018 12:44:17 PM (6 years ago)

Author:

suehring

Message:

Bugfix, topography is rounded to the nearest discrete grid level

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, topography height is rounded to the nearest discrete grid level
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -806 +809 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  zw
+        ONLY:  zu, zw
 
     USE control_parameters,                                                    &
@@ -1018 +1021 @@
 !--             Hence, in this case, de-flag the grid point and give it 
 !--             urban type instead. 
-                IF ( zw(k) - ocean_offset <= terrain_height_f%var(j,i) )  THEN
+                IF ( zu(k) - ocean_offset <= terrain_height_f%var(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), 1 ) 
@@ -1033 +1036 @@
                                         building_id_f%var(j,i) ), DIM = 1 )
         
-                      IF ( zw(k) - ocean_offset <=                             &
+                      IF ( zu(k) - ocean_offset <=                             &
                            oro_max(nr) + buildings_f%var_2d(j,i) )  THEN
                          topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
@@ -1081 +1084 @@
           DO  j = nys, nyn
              DO  k = nzb, nzt
-                IF ( zw(k) - ocean_offset <= buildings_f%var_2d(j,i) )  THEN
+                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