Changeset 4709 for palm/trunk/SOURCE

Timestamp:

Sep 28, 2020 7:20:00 PM (4 years ago)

Author:

maronga

Message:

reactivated release of waste heat and additional fix in indoor model

File:

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

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
-!
+! 
+! 
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Bugfix: avoid division by zero in case of zero window fraction (now also for vertical walls).
+! Reactivated waste heat.
+!
+! 4704 2020-09-28 10:13:03Z maronga
 ! Bugfix: avoid division by zero in case of zero window fraction
 !
@@ -1487 +1491 @@
                             ) / facade_element_area                                             !< [W/m2] , observe the directional
                                                                                                 !< convention in PALM!
-             surf_usm_h(l)%waste_heat(m) = 0.0_wp !q_waste_heat
+             surf_usm_h(l)%waste_heat(m) = q_waste_heat
           ENDDO !< Horizontal surfaces loop
 !
@@ -1690 +1694 @@
 !--          Heat flux into the wall. Value needed in urban_surface_mod to
 !--          calculate heat transfer through wall layers towards the facade
-             q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                               &
+             IF ( (facade_element_area - window_area_per_facade) > 0.0_wp )  THEN
+                q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                 &
                                     / ( facade_element_area - window_area_per_facade )
-             q_win = h_t_es * ( pt(k,j,i) - theta_s )                                              &
-                                    / ( window_area_per_facade )
+             ELSE
+                q_wall = 0.0_wp
+             ENDIF
+
+             IF ( window_area_per_facade > 0.0_wp )  THEN
+                q_win = h_t_es * ( pt(k,j,i) - theta_s ) / ( window_area_per_facade )
+             ELSE
+                q_win = 0.0_wp
+             ENDIF
+
 !
 !--          Transfer q_wall & q_win back to USM (innermost wall/window layer)
@@ -1717 +1730 @@
                                                     ) / facade_element_area  !< [W/m2] , observe the directional convention in
                                                                              !< PALM!
-             surf_usm_v(l)%waste_heat(m) = 0.0_wp !q_waste_heat
+             surf_usm_v(l)%waste_heat(m) = q_waste_heat
           ENDDO !< Vertical surfaces loop
        ENDIF !< buildings(nb)%on_pe