Changeset 4704 for palm


Ignore:
Timestamp:
Sep 28, 2020 10:13:03 AM (4 years ago)
Author:
maronga
Message:

bugfix in indoor model regarding zero window fractions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/indoor_model_mod.f90

    r4702 r4704  
    2525! -----------------
    2626! $Id$
     27! Bugfix: avoid division by zero in case of zero window fraction
     28!
     29! 4702 2020-09-27 18:39:00Z maronga
    2730! Removed unused variable indoor_wall_window_temperature
    2831!
     
    14481451!--          calculate heat transfer through wall layers towards the facade
    14491452!--          (use c_p * rho_surface to convert [W/m2] into [K m/s])
    1450              q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                               &
     1453             IF ( (facade_element_area - window_area_per_facade) > 0.0_wp )  THEN
     1454                q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                 &
    14511455                                    / ( facade_element_area - window_area_per_facade )
    1452              q_win = h_t_es * ( pt(k,j,i) - theta_s )                                              &
    1453                                     / ( window_area_per_facade )
     1456             ELSE
     1457                q_wall = 0.0_wp
     1458             ENDIF
     1459
     1460             IF ( window_area_per_facade > 0.0_wp )  THEN
     1461                q_win = h_t_es * ( pt(k,j,i) - theta_s ) / ( window_area_per_facade )
     1462             ELSE
     1463                q_win = 0.0_wp
     1464             ENDIF
    14541465!
    14551466!--          Transfer q_wall & q_win back to USM (innermost wall/window layer)
Note: See TracChangeset for help on using the changeset viewer.