Ignore:
Timestamp:
Sep 21, 2015 8:15:16 AM (6 years ago)
Author:
gronemeier
Message:

Bugfix: Definition of topography height in case of grid stretching

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/header.f90

    r1592 r1660  
    2020! Current revisions:
    2121! -----------------
     22! Bugfix: Definition of building/street canyon height if vertical grid stretching
     23!         starts below the maximum topography height.
    2224!
    2325!
     
    196198
    197199    USE arrays_3d,                                                             &
    198         ONLY:  pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu
     200        ONLY:  pt_init, qsws, q_init, sa_init, shf, ug, vg, w_subs, zu, zw
    199201       
    200202    USE control_parameters
     
    754756          blx = INT( building_length_x / dx )
    755757          bly = INT( building_length_y / dy )
    756           bh  = INT( building_height / dz )
     758          DO  k = 0, nzt
     759             IF  ( building_height .GT. zw(k) ) THEN
     760                bh = k+1
     761             ENDIF
     762          ENDDO
    757763
    758764          IF ( building_wall_left == 9999999.9_wp )  THEN
     
    772778
    773779       CASE ( 'single_street_canyon' )
    774           ch  = NINT( canyon_height / dz )
     780          DO  k = 0, nzt
     781             IF  ( canyon_height .GT. zw(k) ) THEN
     782                ch = k+1
     783             ENDIF
     784          ENDDO
    775785          IF ( canyon_width_x /= 9999999.9_wp )  THEN
    776786!
Note: See TracChangeset for help on using the changeset viewer.