Changeset 1660 for palm/trunk/SOURCE/header.f90

Timestamp:

Sep 21, 2015 8:15:16 AM (9 years ago)

Author:

gronemeier

Message:

Bugfix: Definition of topography height in case of grid stretching

File:

• palm/trunk/SOURCE/header.f90 (modified) (4 diffs)

palm/trunk/SOURCE/header.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
+! Bugfix: Definition of building/street canyon height if vertical grid stretching
+!         starts below the maximum topography height.
 ! 
 ! 
@@ -196 +198 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu
+        ONLY:  pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu, zw
         
     USE control_parameters
@@ -754 +756 @@
           blx = INT( building_length_x / dx )
           bly = INT( building_length_y / dy )
-          bh  = INT( building_height / dz )
+          DO  k = 0, nzt
+             IF  ( building_height .GT. zw(k) ) THEN
+                bh = k+1
+             ENDIF
+          ENDDO
 
           IF ( building_wall_left == 9999999.9_wp )  THEN
@@ -772 +778 @@
 
        CASE ( 'single_street_canyon' )
-          ch  = NINT( canyon_height / dz )
+          DO  k = 0, nzt
+             IF  ( canyon_height .GT. zw(k) ) THEN
+                ch = k+1
+             ENDIF
+          ENDDO
           IF ( canyon_width_x /= 9999999.9_wp )  THEN
 !