Changeset 1660 for palm/trunk/SOURCE

Timestamp:

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

Author:

gronemeier

Message:

Bugfix: Definition of topography height in case of grid stretching

Location:

palm/trunk/SOURCE

Files:

• header.f90 (modified) (4 diffs)
• init_grid.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
 !

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
+! Bugfix: Definition of topography grid levels if vertical grid stretching
+!         starts below the maximum topography height.
 ! 
 !
@@ -541 +543 @@
           blx = NINT( building_length_x / dx )
           bly = NINT( building_length_y / dy )
-          bh  = NINT( building_height / dz )
+          bh = 0
+          DO  k = 0, nzt
+             IF  ( building_height .GT. zw(k) ) THEN
+                bh = k+1
+             ENDIF
+          ENDDO
 
           IF ( building_wall_left == 9999999.9_wp )  THEN
@@ -601 +608 @@
           ENDIF
 
-          ch             = NINT( canyon_height / dz )
+          ch = 0
+          DO  k = 0, nzt
+             IF  ( canyon_height .GT. zw(k) ) THEN
+                ch = k+1
+             ENDIF
+          ENDDO
           dp_level_ind_b = ch
 !
@@ -675 +687 @@
           DO  i = 0, nx
              DO  j = 0, ny
-                nzb_local(j,i) = NINT( topo_height(j,i) / dz )
+                nzb_local(j,i) = 0
+                DO  k = 0, nzt
+                    IF  ( topo_height(j,i) .GT. zw(k) ) THEN
+                        nzb_local(j,i) = k+1
+                    ENDIF
+                ENDDO
              ENDDO
           ENDDO