Changeset 4681 for palm/trunk/SOURCE/indoor_model_mod.f90

Timestamp:

Sep 16, 2020 10:23:06 AM (4 years ago)

Author:

pavelkrc

Message:

Fix indexing of horizontal surfaces in indoor model

File:

• palm/trunk/SOURCE/indoor_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix for implementation of downward surfaces
+!
+! 4671 2020-09-09 20:27:58Z pavelkrc
 ! Implementation of downward facing USM and LSM surfaces
 !
@@ -180 +183 @@
        INTEGER(iwp) ::  ventilation_int_loads             !< [-] allocation of activity in the building
 
+       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  l_h            !< index array linking surface-element orientation index
+                                                                  !< for horizontal surfaces with building
        INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  l_v            !< index array linking surface-element orientation index
                                                                   !< for vertical surfaces with building
@@ -915 +920 @@
 !--       required in case a building is located at the subdomain boundaries, where the building and
 !--       the corresponding surface elements are located on different subdomains.
+          ALLOCATE( buildings(nb)%l_h(1:buildings(nb)%num_facades_per_building_h_l) )
           ALLOCATE( buildings(nb)%m_h(1:buildings(nb)%num_facades_per_building_h_l) )
 
@@ -957 +963 @@
 
           IF ( buildings(nb)%on_pe )  THEN
+             buildings(nb)%l_h(n_fa(nb)) = l
              buildings(nb)%m_h(n_fa(nb)) = m
              n_fa(nb) = n_fa(nb) + 1
@@ -1220 +1227 @@
 !
 !--          Determine indices where corresponding surface-type information is stored.
-             l = buildings(nb)%l_v(fa)
-             m = buildings(nb)%m_v(fa)
+             l = buildings(nb)%l_h(fa)
+             m = buildings(nb)%m_h(fa)
 !
 !--          Determine building height level index.
@@ -1233 +1240 @@
                                             ( buildings(nb)%num_facades_per_building_h +           &
                                               buildings(nb)%num_facades_per_building_v )              !< [m2] area of total facade
-             window_area_per_facade       = surf_usm_h(l)%frac(m,ind_wat_win)  * facade_element_area     !< [m2] window area per facade
+             window_area_per_facade       = surf_usm_h(l)%frac(m,ind_wat_win)  * facade_element_area  !< [m2] window area per facade
                                                                                                       !< element