Changeset 4709 for palm


Ignore:
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:
1 edited

Legend:

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

    r4704 r4709  
    2020! Current revisions:
    2121! -----------------
    22 !
    23 !
     22! 
     23! 
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! Bugfix: avoid division by zero in case of zero window fraction (now also for vertical walls).
     28! Reactivated waste heat.
     29!
     30! 4704 2020-09-28 10:13:03Z maronga
    2731! Bugfix: avoid division by zero in case of zero window fraction
    2832!
     
    14871491                            ) / facade_element_area                                             !< [W/m2] , observe the directional
    14881492                                                                                                !< convention in PALM!
    1489              surf_usm_h(l)%waste_heat(m) = 0.0_wp !q_waste_heat
     1493             surf_usm_h(l)%waste_heat(m) = q_waste_heat
    14901494          ENDDO !< Horizontal surfaces loop
    14911495!
     
    16901694!--          Heat flux into the wall. Value needed in urban_surface_mod to
    16911695!--          calculate heat transfer through wall layers towards the facade
    1692              q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                               &
     1696             IF ( (facade_element_area - window_area_per_facade) > 0.0_wp )  THEN
     1697                q_wall = h_t_wm * ( indoor_wall_temperature - theta_m )                 &
    16931698                                    / ( facade_element_area - window_area_per_facade )
    1694              q_win = h_t_es * ( pt(k,j,i) - theta_s )                                              &
    1695                                     / ( window_area_per_facade )
     1699             ELSE
     1700                q_wall = 0.0_wp
     1701             ENDIF
     1702
     1703             IF ( window_area_per_facade > 0.0_wp )  THEN
     1704                q_win = h_t_es * ( pt(k,j,i) - theta_s ) / ( window_area_per_facade )
     1705             ELSE
     1706                q_win = 0.0_wp
     1707             ENDIF
     1708
    16961709!
    16971710!--          Transfer q_wall & q_win back to USM (innermost wall/window layer)
     
    17171730                                                    ) / facade_element_area  !< [W/m2] , observe the directional convention in
    17181731                                                                             !< PALM!
    1719              surf_usm_v(l)%waste_heat(m) = 0.0_wp !q_waste_heat
     1732             surf_usm_v(l)%waste_heat(m) = q_waste_heat
    17201733          ENDDO !< Vertical surfaces loop
    17211734       ENDIF !< buildings(nb)%on_pe
Note: See TracChangeset for help on using the changeset viewer.