Ignore:
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:
1 edited

Legend:

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

    r1581 r1660  
    2020! Current revisions:
    2121! -----------------
     22! Bugfix: Definition of topography grid levels if vertical grid stretching
     23!         starts below the maximum topography height.
    2224!
    2325!
     
    541543          blx = NINT( building_length_x / dx )
    542544          bly = NINT( building_length_y / dy )
    543           bh  = NINT( building_height / dz )
     545          bh = 0
     546          DO  k = 0, nzt
     547             IF  ( building_height .GT. zw(k) ) THEN
     548                bh = k+1
     549             ENDIF
     550          ENDDO
    544551
    545552          IF ( building_wall_left == 9999999.9_wp )  THEN
     
    601608          ENDIF
    602609
    603           ch             = NINT( canyon_height / dz )
     610          ch = 0
     611          DO  k = 0, nzt
     612             IF  ( canyon_height .GT. zw(k) ) THEN
     613                ch = k+1
     614             ENDIF
     615          ENDDO
    604616          dp_level_ind_b = ch
    605617!
     
    675687          DO  i = 0, nx
    676688             DO  j = 0, ny
    677                 nzb_local(j,i) = NINT( topo_height(j,i) / dz )
     689                nzb_local(j,i) = 0
     690                DO  k = 0, nzt
     691                    IF  ( topo_height(j,i) .GT. zw(k) ) THEN
     692                        nzb_local(j,i) = k+1
     693                    ENDIF
     694                ENDDO
    678695             ENDDO
    679696          ENDDO
Note: See TracChangeset for help on using the changeset viewer.