 Timestamp:
 Nov 22, 2019 10:13:38 AM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/indoor_model_mod.f90
r4267 r4299 26 26 !  27 27 ! $Id$ 28 ! Output of indoor temperature revised (to avoid nondefined values within 29 ! buildings) 30 ! 31 ! 4267 20191016 18:58:49Z suehring 28 32 ! Bugfix in initialization, some indices to access building_pars where wrong. 29 33 ! Introduction of seasonal parameters. … … 1669 1673 ! Average indoor temperature. Note, in case a building is completely 1670 1674 ! surrounded by higher buildings, it may have no facade elements 1671 ! at some height levels, will will lead to a divide by zero. If this 1672 ! is the case, indoor temperature will be set to 1.0. 1675 ! at some height levels, which will lead to a divide by zero. 1673 1676 DO k = buildings(nb)%kb_min, buildings(nb)%kb_max 1674 1677 IF ( buildings(nb)%num_facade_h(k) + & … … 1677 1680 REAL( buildings(nb)%num_facade_h(k) + & 1678 1681 buildings(nb)%num_facade_v(k), KIND = wp ) 1679 ELSE 1680 buildings(nb)%t_in(k) = 1.0_wp 1682 ENDIF 1683 ENDDO 1684 ! 1685 ! If indoor temperature is not defined because of missing facade 1686 ! elements, the values from the abovelying level will be taken. 1687 ! At least at the top of the buildings facades are defined, so that 1688 ! at least there an indoor temperature is defined. This information 1689 ! will propagate downwards the building. 1690 DO k = buildings(nb)%kb_max1, buildings(nb)%kb_min, 1 1691 IF ( buildings(nb)%num_facade_h(k) + & 1692 buildings(nb)%num_facade_v(k) <= 0 ) THEN 1693 buildings(nb)%t_in(k) = buildings(nb)%t_in(k+1) 1681 1694 ENDIF 1682 1695 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.