Ignore:
Timestamp:
Jan 15, 2018 12:44:17 PM (3 years ago)
Author:
suehring
Message:

Bugfix, topography is rounded to the nearest discrete grid level

File:
1 edited

Legend:

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

    r2718 r2747  
    2525! -----------------
    2626! $Id$
     27! Bugfix, topography height is rounded to the nearest discrete grid level
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    806809
    807810    USE arrays_3d,                                                             &
    808         ONLY:  zw
     811        ONLY:  zu, zw
    809812
    810813    USE control_parameters,                                                    &
     
    10181021!--             Hence, in this case, de-flag the grid point and give it
    10191022!--             urban type instead.
    1020                 IF ( zw(k) - ocean_offset <= terrain_height_f%var(j,i) )  THEN
     1023                IF ( zu(k) - ocean_offset <= terrain_height_f%var(j,i) )  THEN
    10211024                    topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
    10221025                    topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 1 )
     
    10331036                                        building_id_f%var(j,i) ), DIM = 1 )
    10341037       
    1035                       IF ( zw(k) - ocean_offset <=                             &
     1038                      IF ( zu(k) - ocean_offset <=                             &
    10361039                           oro_max(nr) + buildings_f%var_2d(j,i) )  THEN
    10371040                         topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
     
    10811084          DO  j = nys, nyn
    10821085             DO  k = nzb, nzt
    1083                 IF ( zw(k) - ocean_offset <= buildings_f%var_2d(j,i) )  THEN
     1086                IF ( zu(k) - ocean_offset <= buildings_f%var_2d(j,i) )  THEN
    10841087                    topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
    10851088                    topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) !indicates terrain
Note: See TracChangeset for help on using the changeset viewer.