Ignore:
Timestamp:
Jul 10, 2018 12:49:26 PM (3 years ago)
Author:
suehring
Message:

Separate bridges as 3D building objects from normal surface-mounted buildings in terms of correct referencing onto the terrain

File:
1 edited

Legend:

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

    r3103 r3115  
    2525! -----------------
    2626! $Id$
     27! Referencing of buildings onto top of terrain - special treatment for bridges.
     28!
     29! 3103 2018-07-04 17:30:52Z suehring
    2730! Reference lowest terrain height to zero level
    2831!
     
    12871290
    12881291    USE netcdf_data_input_mod,                                                 &
    1289         ONLY:  buildings_f, building_id_f, input_pids_static,                  &
     1292        ONLY:  buildings_f, building_id_f, building_type_f, input_pids_static, &
    12901293               terrain_height_f
    12911294
     
    15771580!--                terrain top is defined at upper bound of the grid box.
    15781581!--                Hence, check for zw in this case.
    1579                    DO k = topo_top_index + 1, nzt + 1     
    1580                       IF ( zw(k) - ocean_offset <= oro_max(nr) )  THEN
    1581                          topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
    1582                          topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 )
    1583                       ENDIF
    1584                    ENDDO       
    1585 !
    1586 !--                After surface irregularities are smoothen, determine lower
    1587 !--                start index where building starts.
    1588                    DO  k = nzb, nzt
    1589                       IF ( zw(k) - ocean_offset <= oro_max(nr) )               &
    1590                          topo_top_index = k
    1591                    ENDDO
     1582!--                Note, do this only for buildings which are surface mounted,
     1583!--                i.e. building types 1-6. Below bridges, which are represented
     1584!--                exclusively by building type 7, terrain shape should be
     1585!--                maintained.
     1586                   IF ( building_type_f%var(j,i) /= 7 )  THEN
     1587                      DO k = topo_top_index + 1, nzt + 1     
     1588                         IF ( zw(k) - ocean_offset <= oro_max(nr) )  THEN
     1589                            topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 )
     1590                            topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 )
     1591                         ENDIF
     1592                      ENDDO       
     1593!
     1594!--                   After surface irregularities are smoothen, determine lower
     1595!--                   start index where building starts.
     1596                      DO  k = nzb, nzt
     1597                         IF ( zw(k) - ocean_offset <= oro_max(nr) )            &
     1598                            topo_top_index = k
     1599                      ENDDO
     1600                   ENDIF
    15921601!
    15931602!--                Finally, map building on top.
Note: See TracChangeset for help on using the changeset viewer.