Changeset 4704 for palm

Timestamp:

Sep 28, 2020 10:13:03 AM (4 years ago)

Author:

maronga

Message:

bugfix in indoor model regarding zero window fractions

File:

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

palm/trunk/SOURCE/indoor_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: avoid division by zero in case of zero window fraction
+!
+! 4702 2020-09-27 18:39:00Z maronga
 ! Removed unused variable indoor_wall_window_temperature
 !
@@ -1448 +1451 @@
 !--          calculate heat transfer through wall layers towards the facade
 !--          (use c_p * rho_surface to convert [W/m2] into [K m/s])
-             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)