Changeset 3115 for palm


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

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

Location:
palm/trunk/SOURCE
Files:
2 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.
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3091 r3115  
    2828! -----------------
    2929! $Id$
     30! Additional building type to represent bridges
     31!
     32! 3091 2018-06-28 16:20:35Z suehring
    3033! - Limit aerodynamic resistance at vertical walls.
    3134! - Add check for local roughness length not exceeding surface-layer height and
     
    428431
    429432
    430     CHARACTER(37), DIMENSION(0:6), PARAMETER :: building_type_name = (/     &
     433    CHARACTER(37), DIMENSION(0:7), PARAMETER :: building_type_name = (/     &
    431434                                   'user-defined                         ', & !  0
    432435                                   'residential - 1950                   ', & !  1
     
    435438                                   'office - 1950                        ', & !  4
    436439                                   'office 1951 - 2000                   ', & !  5
    437                                    'office 2001 -                        '  & !  6
     440                                   'office 2001 -                        ', & !  6
     441                                   'bridges                              '  & !  7
    438442                                                                     /)
    439443!
     
    460464!-- 45 - heat capacity of the wall surface, 46 - heat conductivity
    461465!-- Please note, only preleminary dummy values so far!
    462     REAL(wp), DIMENSION(0:46,1:6), PARAMETER :: building_pars = RESHAPE( (/    &
     466    REAL(wp), DIMENSION(0:46,1:7), PARAMETER :: building_pars = RESHAPE( (/    &
    463467        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, 1.0_wp,                        & !parameter 0-5
    464468        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
    465469        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
    466         0.001_wp, 0.0001_wp, 4.0_wp,                                             & !parameter 18-20
     470        0.001_wp, 0.0001_wp, 4.0_wp,                                           & !parameter 18-20
    467471        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
    468472        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     
    475479        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
    476480        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
    477         0.001_wp, 0.0001_wp, 4.0_wp,                                             & !parameter 18-20
     481        0.001_wp, 0.0001_wp, 4.0_wp,                                           & !parameter 18-20
    478482        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
    479483        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     
    486490        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
    487491        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
    488         0.001_wp, 0.0001_wp, 4.0_wp,                                             & !parameter 18-20
     492        0.001_wp, 0.0001_wp, 4.0_wp,                                           & !parameter 18-20
    489493        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
    490494        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     
    508512        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
    509513        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
    510         0.001_wp, 0.0001_wp, 4.0_wp,                                             & !parameter 18-20
     514        0.001_wp, 0.0001_wp, 4.0_wp,                                           & !parameter 18-20
    511515        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
    512516        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     
    519523        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
    520524        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
    521         0.001_wp, 0.0001_wp, 4.0_wp,                                             & !parameter 18-20
     525        0.001_wp, 0.0001_wp, 4.0_wp,                                           & !parameter 18-20
    522526        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
    523527        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     
    526530        24.0_wp, 24.0_wp, 24.0_wp,                                             & !parameter 38-40
    527531        0.0242_wp, 0.0969_wp, 0.346_wp, 1.0_wp,                                & !parameter 41-44
    528         20000.0_wp, 10.0_wp                                                    & !parameter 45-46 - end of type 6
     532        20000.0_wp, 10.0_wp,                                                   & !parameter 45-46 - end of type 6
     533        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, 1.0_wp,                        & !parameter 0-5
     534        1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp,      & !parameter 6-11
     535        296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp,                & !parameter 12-17
     536        0.001_wp, 0.0001_wp, 0.0_wp,                                           & !parameter 18-20
     537        1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp,                                & !parameter 21-25
     538        1000000.0_wp, 1000000.0_wp, 1000000.0_wp,                              & !parameter 26-28                     
     539        0.3_wp, 0.3_wp, 0.3_wp,                                                & !parameter 29-31       
     540        0.4_wp, 0.4_wp, 0.4_wp, 0.4_wp, 0.01_wp, 0.001_wp,                     & !parameter 32-37
     541        24.0_wp, 24.0_wp, 24.0_wp,                                             & !parameter 38-40
     542        0.0242_wp, 0.0969_wp, 0.346_wp, 1.0_wp,                                & !parameter 41-44
     543        20000.0_wp, 10.0_wp                                                    & !parameter 45-46 - end of type 7 (bridges)
    529544                                                                          /),  &
    530                                                                (/47, 6/) )
     545                                                               (/47, 7/) )
    531546
    532547!
Note: See TracChangeset for help on using the changeset viewer.