Changeset 4687 for palm


Ignore:
Timestamp:
Sep 21, 2020 7:40:16 PM (4 years ago)
Author:
maronga
Message:

bugfix in indoor model, code layout change in urban surface model, indoor model integrated in spinup mechanism

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

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

    r4681 r4687  
    2525! -----------------
    2626! $Id$
     27! Bugfix: values of theta_m_t_prev were not stored for individual surfaces and thus re-used by all
     28! surfaces and buildings, which led to excessive indoor temperatures
     29!
     30! 4681 2020-09-16 10:23:06Z pavelkrc
    2731! Bugfix for implementation of downward surfaces
    2832!
     
    126130! @author Matthias Suehring
    127131! @author Sascha Rißmann
     132! @author Björn Maronga
    128133!
    129134!
     
    144149!> @note How to write indoor temperature output to pt array?
    145150!>
    146 !> @bug  <Enter known bugs here>
     151!> @bug  Calculation of iwghf_eb and iwghf_eb_window is faulty
    147152!--------------------------------------------------------------------------------------------------!
    148153 MODULE indoor_model_mod
     
    234239       REAL(wp), DIMENSION(:), ALLOCATABLE ::  t_in       !< mean building indoor temperature, height dependent
    235240       REAL(wp), DIMENSION(:), ALLOCATABLE ::  t_in_l     !< mean building indoor temperature on local subdomain, height dependent
     241       REAL(wp), DIMENSION(:), ALLOCATABLE ::  theta_m_t_prev_h !< [degree_C] value of theta_m_t from previous time step (horizontal)
     242       REAL(wp), DIMENSION(:), ALLOCATABLE ::  theta_m_t_prev_v !< [degree_C] value of theta_m_t from previous time step (vertical)
    236243       REAL(wp), DIMENSION(:), ALLOCATABLE ::  volume     !< total building volume, height dependent
    237244       REAL(wp), DIMENSION(:), ALLOCATABLE ::  vol_frac   !< fraction of local on total building volume, height dependent
     
    319326    REAL(wp) ::  theta_air_set                       !< [degree_C] Setpoint_temperature for the room
    320327    REAL(wp) ::  theta_m                             !< [degree_C} inner temperature of the RC-node
    321     REAL(wp) ::  theta_m_t                           !< [degree_C] (Fictive) component temperature timestep
    322     REAL(wp) ::  theta_m_t_prev                      !< [degree_C] (Fictive) component temperature previous timestep (do not change)
     328    REAL(wp) ::  theta_m_t                           !< [degree_C] (Fictive) component temperature during timestep
    323329    REAL(wp) ::  theta_op                            !< [degree_C] operative temperature
    324330    REAL(wp) ::  theta_s                             !< [degree_C] surface temperature of the RC-node
     
    449455!--------------------------------------------------------------------------------------------------!
    450456 SUBROUTINE im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                        &
    451                                    near_facade_temperature, phi_hc_nd_dummy )
     457                                   near_facade_temperature, phi_hc_nd_dummy, theta_m_t_prev )
    452458
    453459    INTEGER(iwp) ::  i
     
    458464    REAL(wp) ::  near_facade_temperature
    459465    REAL(wp) ::  phi_hc_nd_dummy
     466    REAL(wp), INTENT(IN) :: theta_m_t_prev
    460467!
    461468!-- Calculation of total mass specific thermal load (internal and external)
     
    468475               )                                                                !< [degree_C] Eq. (C.5)
    469476!
    470 !-- Calculation of component temperature at factual timestep
     477!-- Calculation of component temperature at current timestep
    471478    theta_m_t = ( ( theta_m_t_prev                                                                 &
    472479                    * ( ( c_m / 3600.0_wp ) - 0.5_wp * ( h_t_3 + h_t_wm ) )                        &
     
    476483                )                                                               !< [degree_C] Eq. (C.4)
    477484!
    478 !-- Calculation of mean inner temperature for the RC-node in actual timestep
     485!-- Calculation of mean inner temperature for the RC-node in current timestep
    479486    theta_m = ( theta_m_t + theta_m_t_prev ) * 0.5_wp                           !< [degree_C] Eq. (C.9)
    480487
    481488!
    482 !-- Calculation of mean surface temperature of the RC-node in actual timestep
     489!-- Calculation of mean surface temperature of the RC-node in current timestep
    483490    theta_s = ( (   h_t_ms * theta_m + phi_st + h_t_es * pt(k,j,i)                                 &
    484491                  + h_t_1  * ( near_facade_temperature                                             &
     
    490497!
    491498!-- Calculation of the air temperature of the RC-node
     499
     500
    492501    theta_air = ( h_t_is * theta_s + h_v * near_facade_temperature + phi_ia + phi_hc_nd_dummy ) /  &
    493502                ( h_t_is + h_v )                                                !< [degree_C] Eq. (C.11)
     503
    494504
    495505 END SUBROUTINE im_calc_temperatures
     
    925935          ALLOCATE( buildings(nb)%l_v(1:buildings(nb)%num_facades_per_building_v_l) )
    926936          ALLOCATE( buildings(nb)%m_v(1:buildings(nb)%num_facades_per_building_v_l) )
     937
     938          ALLOCATE( buildings(nb)%theta_m_t_prev_h(1:buildings(nb)%num_facades_per_building_h_l) )
     939          ALLOCATE( buildings(nb)%theta_m_t_prev_v(1:buildings(nb)%num_facades_per_building_v_l) )
    927940       ENDIF
    928941
     
    11251138    ENDDO
    11261139!
    1127 !-- Initial room temperature [K]
    1128 !-- (after first loop, use theta_m_t as theta_m_t_prev)
    1129     theta_m_t_prev = initial_indoor_temperature
    1130 !
    11311140!-- Initialize indoor temperature. Actually only for output at initial state.
    11321141    DO  nb = 1, num_build
    1133        IF ( buildings(nb)%on_pe )  buildings(nb)%t_in(:) = initial_indoor_temperature
     1142       IF ( buildings(nb)%on_pe )  THEN
     1143          buildings(nb)%t_in(:) = initial_indoor_temperature
     1144
     1145!
     1146!--       (after first loop, use theta_m_t as theta_m_t_prev)
     1147          buildings(nb)%theta_m_t_prev_h(:) = initial_indoor_temperature
     1148          buildings(nb)%theta_m_t_prev_v(:) = initial_indoor_temperature
     1149
     1150       ENDIF
    11341151    ENDDO
    11351152
     
    13371354
    13381355             CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                 &
    1339                                           near_facade_temperature, phi_hc_nd )
     1356                                          near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
    13401357!
    13411358!--          If air temperature between border temperatures of heating and cooling, assign output
     
    13661383
    13671384                CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
    1368                                              near_facade_temperature, phi_hc_nd )
     1385                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
    13691386                theta_air_10 = theta_air                                                !< temperature with 10 W/m2 of heating
    13701387                phi_hc_nd_un = phi_hc_nd_10 * (theta_air_set - theta_air_0)                        &
     
    13951412!--             Calculate the temperature with phi_hc_nd_ac (new)
    13961413                CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
    1397                                              near_facade_temperature, phi_hc_nd )
     1414                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
    13981415                theta_air_ac = theta_air
    13991416             ENDIF
    14001417!
    14011418!--          Update theta_m_t_prev
    1402              theta_m_t_prev = theta_m_t
     1419             buildings(nb)%theta_m_t_prev_h(fa) = theta_m_t
     1420
    14031421
    14041422             q_vent = h_v * ( theta_air - near_facade_temperature )
     
    14081426!--          Will be used for thermal comfort calculations.
    14091427             theta_op     = 0.3_wp * theta_air_ac + 0.7_wp * theta_s          !< [degree_C] operative Temperature Eq. (C.12)
     1428
    14101429!              surf_usm_h(l)%t_indoor(m) = theta_op                               !< not integrated now
    14111430!
     
    14401459                            ) / facade_element_area                                             !< [W/m2] , observe the directional
    14411460                                                                                                !< convention in PALM!
    1442              surf_usm_h(l)%waste_heat(m) = q_waste_heat
     1461             surf_usm_h(l)%waste_heat(m) = 0.0_wp !q_waste_heat
    14431462          ENDDO !< Horizontal surfaces loop
    14441463!
     
    15241543             k = surf_usm_v(l)%k(m)
    15251544             near_facade_temperature = surf_usm_v(l)%pt_10cm(m)
    1526              indoor_wall_window_temperature =                                                      &
    1527                                     surf_usm_v(l)%frac(m,ind_veg_wall) * t_wall_v(l)%val(nzt_wall,m) &
     1545             indoor_wall_window_temperature =                                                          &
     1546                                    surf_usm_v(l)%frac(m,ind_veg_wall) * t_wall_v(l)%val(nzt_wall,m)   &
    15281547                                  + surf_usm_v(l)%frac(m,ind_wat_win)  * t_window_v(l)%val(nzt_wall,m)
     1548
    15291549!
    15301550!--          Solar thermal gains. If net_sw_in larger than sun-protection
     
    15631583             phi_hc_nd = 0.0_wp
    15641584             CALL im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                  &
    1565                                          near_facade_temperature, phi_hc_nd )
     1585                                         near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
    15661586!
    15671587!--          If air temperature between border temperatures of heating and cooling, assign output
     
    15921612
    15931613                CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
    1594                                              near_facade_temperature, phi_hc_nd )
     1614                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
    15951615
    15961616                theta_air_10 = theta_air !< Note the temperature with 10 W/m2 of heating
     
    16231643!--             Calculate the temperature with phi_hc_nd_ac (new)
    16241644                CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature, &
    1625                                              near_facade_temperature, phi_hc_nd )
     1645                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
    16261646                theta_air_ac = theta_air
    16271647             ENDIF
    16281648!
    16291649!--          Update theta_m_t_prev
    1630              theta_m_t_prev = theta_m_t
     1650             buildings(nb)%theta_m_t_prev_v(fa) = theta_m_t
     1651
    16311652
    16321653             q_vent = h_v * ( theta_air - near_facade_temperature )
     
    16351656!--          Will be used for thermal comfort calculations.
    16361657             theta_op     = 0.3_wp * theta_air_ac + 0.7_wp * theta_s
     1658
    16371659!              surf_usm_v(l)%t_indoor(m) = theta_op                  !< not integrated yet
    16381660!
     
    16651687                                                    ) / facade_element_area  !< [W/m2] , observe the directional convention in
    16661688                                                                             !< PALM!
    1667              surf_usm_v(l)%waste_heat(m) = q_waste_heat
     1689             surf_usm_v(l)%waste_heat(m) = 0.0_wp !q_waste_heat
    16681690          ENDDO !< Vertical surfaces loop
    16691691       ENDIF !< buildings(nb)%on_pe
  • palm/trunk/SOURCE/time_integration_spinup.f90

    r4671 r4687  
    2525! -----------------
    2626! $Id$
     27! Indoor model is now available during spinup
     28!
     29! 4671 2020-09-09 20:27:58Z pavelkrc
    2730! Implementation of downward facing USM and LSM surfaces
    2831!
     
    107110               constant_flux_layer,                                                                &
    108111               coupling_start_time,                                                                &
     112               current_timestep_number,                                                            &
    109113               data_output_during_spinup,                                                          &
    110114               debug_output_timestep,                                                              &
     
    120124               dt_3d,                                                                              &
    121125               humidity,                                                                           &
     126               indoor_model,                                                                       &
    122127               intermediate_timestep_count,                                                        &
    123128               intermediate_timestep_count_max,                                                    &
     
    161166               nxrg
    162167
     168    USE indoor_model_mod,                                                                          &
     169        ONLY:  dt_indoor, im_main_heatcool, time_indoor
     170
    163171    USE land_surface_model_mod,                                                                    &
    164172        ONLY:  lsm_energy_balance,                                                                 &
     
    477485
    478486!
     487!--    If required, calculate indoor temperature, waste heat, heat flux
     488!--    through wall, etc.
     489!--    dt_indoor steers the frequency of the indoor model calculations.
     490!--    Note, at first timestep indoor model is called, in order to provide
     491!--    a waste heat flux.
     492       IF ( indoor_model )  THEN
     493
     494          time_indoor = time_indoor + dt_3d
     495
     496          IF ( time_indoor >= dt_indoor  .OR.  current_timestep_number == 0 )  THEN
     497
     498             time_indoor = time_indoor - dt_indoor
     499
     500             CALL im_main_heatcool
     501
     502          ENDIF
     503       ENDIF
     504
     505!
    479506!--    Increase simulation time and output times
    480507       current_timestep_number_spinup = current_timestep_number_spinup + 1
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r4671 r4687  
    2727! -----------------
    2828! $Id$
     29! Optimized code structure for treatment of inner wall and window heat flux
     30!
     31! 4671 2020-09-09 20:27:58Z pavelkrc
    2932! Radiative transfer model RTM version 4.1
    3033! - Implementation of downward facing USM and LSM surfaces
     
    51735176       wtend(nzb_wall) = ( 1.0_wp / surf%rho_c_wall(nzb_wall,m) )                     &
    51745177                          * ( surf%lambda_h(nzb_wall,m) * wall_mod(nzb_wall)          &
    5175                               * ( t_wall%val(nzb_wall+1,m) - t_wall%val(nzb_wall,m) )              &
     5178                              * ( t_wall%val(nzb_wall+1,m) - t_wall%val(nzb_wall,m) ) &
    51765179                              * surf%ddz_wall(nzb_wall+1,m)                           &
    51775180                              + surf%frac(m,ind_veg_wall)                             &
     
    51835186                                  + surf%frac(m,ind_pav_green) )                      &
    51845187                              * ( surf%lambda_h_green(nzt_wall,m)                     &
    5185                               * wall_mod(nzt_wall)                                             &
     5188                              * wall_mod(nzt_wall)                                    &
    51865189                              * surf%ddz_green(nzt_wall,m)                            &
    51875190                              + surf%lambda_h(nzb_wall,m)                             &
    5188                               * wall_mod(nzb_wall)                                             &
     5191                              * wall_mod(nzb_wall)                                    &
    51895192                              * surf%ddz_wall(nzb_wall,m) )                           &
    51905193                              / ( surf%ddz_green(nzt_wall,m)                          &
    51915194                              + surf%ddz_wall(nzb_wall,m) )                           &
    5192                               * ( t_wall%val(nzb_wall,m) - t_green%val(nzt_wall,m) )               &
     5195                              * ( t_wall%val(nzb_wall,m) - t_green%val(nzt_wall,m) )  &
    51935196                            ) * surf%ddz_wall_stag(nzb_wall,m)
    5194 !
    5195 !--    If indoor model is used inner wall layer is calculated by using iwghf (indoor wall ground heat flux)
    5196        IF ( indoor_model )  THEN
     5197
     5198!
     5199!--    If no indoor model is used, calculate the inner wall heat flux according to
     5200!--    heat diffusion (innermost wall temperature as boundary conditions)
     5201       IF ( .NOT. indoor_model ) THEN
     5202          surf%iwghf_eb(m) = - surf%lambda_h(nzt_wall-1,m) * ( t_wall%val(nzt_wall,m) &
     5203                             - t_wall%val(nzt_wall-1,m) ) * surf%ddz_wall(nzt_wall,m)
     5204          ENDIF
     5205
    51975206          DO  kw = nzb_wall+1, nzt_wall-1
    5198              wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                           &
    5199                          * ( surf%lambda_h(kw,m) * wall_mod(kw)                       &
    5200                              * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                           &
    5201                              * surf%ddz_wall(kw+1,m)                                  &
    5202                              - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                 &
    5203                              * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                           &
    5204                              * surf%ddz_wall(kw,m)                                    &
     5207             wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                             &
     5208                         * ( surf%lambda_h(kw,m) * wall_mod(kw)                         &
     5209                             * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                &
     5210                             * surf%ddz_wall(kw+1,m)                                    &
     5211                             - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                   &
     5212                             * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                &
     5213                             * surf%ddz_wall(kw,m)                                      &
    52055214                           ) * surf%ddz_wall_stag(kw,m)
    52065215          ENDDO
    5207           wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                  &
    5208                             * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)   &
    5209                                 * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) )            &
    5210                                 * surf%ddz_wall(nzt_wall,m)                           &
    5211                                 + surf%iwghf_eb(m)                                    &
     5216          wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                    &
     5217                            * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)     &
     5218                                * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
     5219                                * surf%ddz_wall(nzt_wall,m)                             &
     5220                                + surf%iwghf_eb(m)                                      &
    52125221                              ) * surf%ddz_wall_stag(nzt_wall,m)
    5213        ELSE
    5214           DO  kw = nzb_wall+1, nzt_wall
    5215              wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                           &
    5216                         * ( surf%lambda_h(kw,m)  * wall_mod(kw)                       &
    5217                             * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                            &
    5218                             * surf%ddz_wall(kw+1,m)                                   &
    5219                             - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                  &
    5220                             * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                            &
    5221                             * surf%ddz_wall(kw,m)                                     &
    5222                           ) * surf%ddz_wall_stag(kw,m)
    5223           ENDDO
    5224        ENDIF
    5225 
    5226        t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d                 &
    5227                                          * ( tsc(2) * wtend(nzb_wall:nzt_wall) + tsc(3)        &
     5222
     5223       t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d      &
     5224                                         * ( tsc(2) * wtend(nzb_wall:nzt_wall) + tsc(3) &
    52285225                                             * surf%tt_wall_m(nzb_wall:nzt_wall,m) )
    52295226
    5230 !
     5227
    52315228!--    During spinup the tempeature inside window layers is not calculated to make larger timesteps possible
    52325229       IF ( .NOT. during_spinup )  THEN
    5233           win_absorp = -log( surf%transmissivity(m) ) /                                &
    5234                         surf%zw_window(nzt_wall,m)
     5230          win_absorp = - log( surf%transmissivity(m) ) / surf%zw_window(nzt_wall,m)
    52355231!
    52365232!--       Prognostic equation for ground/roof window temperature t_window takes absorption
    52375233!--       of shortwave radiation into account
    52385234          wintend(:) = 0.0_wp
    5239           wintend(nzb_wall) = ( 1.0_wp / surf%rho_c_window(nzb_wall,m) )               &
    5240                                 * ( surf%lambda_h_window(nzb_wall,m)                   &
    5241                                 * ( t_window%val(nzb_wall+1,m) - t_window%val(nzb_wall,m) )         &
    5242                                 * surf%ddz_window(nzb_wall+1,m)                        &
    5243                                 + surf%wghf_eb_window(m)                               &
    5244                                 + surf%rad_sw_in(m)                                    &
    5245                                 * ( 1.0_wp - exp( -win_absorp                                   &
    5246                                     * surf%zw_window(nzb_wall,m) ) )                   &
     5235          wintend(nzb_wall) = ( 1.0_wp / surf%rho_c_window(nzb_wall,m) )                    &
     5236                                * ( surf%lambda_h_window(nzb_wall,m)                        &
     5237                                * ( t_window%val(nzb_wall+1,m) - t_window%val(nzb_wall,m) ) &
     5238                                * surf%ddz_window(nzb_wall+1,m)                             &
     5239                                + surf%wghf_eb_window(m)                                    &
     5240                                + surf%rad_sw_in(m)                                         &
     5241                                * ( 1.0_wp - exp( -win_absorp                               &
     5242                                    * surf%zw_window(nzb_wall,m) ) )                        &
    52475243                              ) * surf%ddz_window_stag(nzb_wall,m)
    52485244
    5249           IF ( indoor_model ) THEN
    5250              DO  kw = nzb_wall+1, nzt_wall-1
    5251                 wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                     &
    5252                               * ( surf%lambda_h_window(kw,m)                           &
    5253                                   * ( t_window%val(kw+1,m) - t_window%val(kw,m) )                   &
    5254                                   * surf%ddz_window(kw+1,m)                            &
    5255                                   - surf%lambda_h_window(kw-1,m)                       &
    5256                                  * ( t_window%val(kw,m) - t_window%val(kw-1,m) )                    &
    5257                                  * surf%ddz_window(kw,m)                               &
    5258                                  + surf%rad_sw_in(m)                                   &
    5259                                  * ( exp( -win_absorp * surf%zw_window(kw-1,m) )       &
    5260                                      - exp(-win_absorp * surf%zw_window(kw,m) )        &
    5261                                    )                                                            &
    5262                                  ) * surf%ddz_window_stag(kw,m)
    5263 
    5264              ENDDO
    5265              wintend(nzt_wall) = ( 1.0_wp / surf%rho_c_window(nzt_wall,m) )            &
    5266                                  * ( -surf%lambda_h_window(nzt_wall-1,m)               &
    5267                                     * ( t_window%val(nzt_wall,m) - t_window%val(nzt_wall-1,m) )     &
    5268                                     * surf%ddz_window(nzt_wall,m)                      &
    5269                                     + surf%iwghf_eb_window(m)                          &
    5270                                     + surf%rad_sw_in(m)                                &
    5271                                     * ( exp( -win_absorp                                        &
    5272                                         * surf%zw_window(nzt_wall-1,m) )               &
    5273                                       - exp( -win_absorp                                        &
    5274                                         * surf%zw_window(nzt_wall,m) )                 &
    5275                                       )                                                         &
    5276                                     ) * surf%ddz_window_stag(nzt_wall,m)
    5277           ELSE
    5278              DO  kw = nzb_wall+1, nzt_wall
    5279                 wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                     &
    5280                               * ( surf%lambda_h_window(kw,m)                           &
    5281                                   * ( t_window%val(kw+1,m) - t_window%val(kw,m) )                   &
    5282                                   * surf%ddz_window(kw+1,m)                            &
    5283                                   - surf%lambda_h_window(kw-1,m)                       &
    5284                                   * ( t_window%val(kw,m)                                          &
    5285                                   - t_window%val(kw-1,m) )                                        &
    5286                                   * surf%ddz_window(kw,m)                              &
    5287                                   + surf%rad_sw_in(m)                                  &
    5288                                    * ( exp( -win_absorp * surf%zw_window(kw-1,m) )     &
    5289                                        - exp(-win_absorp * surf%zw_window(kw,m) )      &
    5290                                      )                                                          &
    5291                                  ) * surf%ddz_window_stag(kw,m)
    5292 
    5293              ENDDO
    5294           ENDIF
    5295 
    5296           t_window_p%val(nzb_wall:nzt_wall,m) = t_window%val(nzb_wall:nzt_wall,m) + dt_3d           &
     5245!
     5246!--       If no indoor model is used, calculate the inner window heat flux according to
     5247!--       heat diffusion (innermost window temperature as boundary conditions
     5248          IF ( .NOT. indoor_model ) THEN
     5249             surf%iwghf_eb_window(m) = - surf%lambda_h_window(nzt_wall-1,m)            &
     5250                                       * ( t_window%val(nzt_wall,m)                    &
     5251                                       - t_window%val(nzt_wall-1,m) )                  &
     5252                                       * surf%ddz_window(nzt_wall,m)
     5253          ENDIF
     5254
     5255
     5256          DO  kw = nzb_wall+1, nzt_wall-1
     5257             wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                        &
     5258                           * ( surf%lambda_h_window(kw,m)                              &
     5259                           * ( t_window%val(kw+1,m) - t_window%val(kw,m) )             &
     5260                           * surf%ddz_window(kw+1,m)                                   &
     5261                           - surf%lambda_h_window(kw-1,m)                              &
     5262                           * ( t_window%val(kw,m) - t_window%val(kw-1,m) )             &
     5263                           * surf%ddz_window(kw,m)                                     &
     5264                           + surf%rad_sw_in(m)                                         &
     5265                           * ( exp( -win_absorp * surf%zw_window(kw-1,m) )             &
     5266                           - exp(-win_absorp * surf%zw_window(kw,m) )                  &
     5267                             )                                                         &
     5268                             ) * surf%ddz_window_stag(kw,m)
     5269
     5270          ENDDO
     5271
     5272          wintend(nzt_wall) = ( 1.0_wp / surf%rho_c_window(nzt_wall,m) )                     &
     5273                                 * ( -surf%lambda_h_window(nzt_wall-1,m)                     &
     5274                                 * ( t_window%val(nzt_wall,m) - t_window%val(nzt_wall-1,m) ) &
     5275                                 * surf%ddz_window(nzt_wall,m)                               &
     5276                                 + surf%iwghf_eb_window(m)                                   &
     5277                                 + surf%rad_sw_in(m)                                         &
     5278                                 * ( exp( -win_absorp                                        &
     5279                                 * surf%zw_window(nzt_wall-1,m) )                            &
     5280                                 - exp( -win_absorp                                          &
     5281                                 * surf%zw_window(nzt_wall,m) )                              &
     5282                                   )                                                         &
     5283                                   ) * surf%ddz_window_stag(nzt_wall,m)
     5284
     5285          t_window_p%val(nzb_wall:nzt_wall,m) = t_window%val(nzb_wall:nzt_wall,m) + dt_3d       &
    52975286                                              * ( tsc(2) * wintend(nzb_wall:nzt_wall) + tsc(3)  &
    52985287                                              * surf%tt_window_m(nzb_wall:nzt_wall,m) )
    52995288
    53005289       ENDIF
    5301 
    53025290!
    53035291!--    Calculate t_wall tendencies for the next Runge-Kutta step
     
    76237611
    76247612    building_pars(:,2) = (/                                                                        &
    7625           0.73_wp,        &  !< parameter 0   - wall fraction above ground floor level
    7626           0.27_wp,        &  !< parameter 1   - window fraction above ground floor level
    7627           0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
    7628           0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
    7629           1.5_wp,         &  !< parameter 4   - LAI roof
    7630           1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
    7631           2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
    7632           103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
    7633           900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
    7634           0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
    7635           0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
    7636           0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
    7637           299.15_wp,      &  !< parameter 12  - indoor target summer temperature
    7638           293.15_wp,      &  !< parameter 13  - indoor target winter temperature
    7639           0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
    7640           0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
    7641           0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
    7642           0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
    7643           0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
    7644           0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
    7645           4.0_wp,         &  !< parameter 20  - ground floor level height
    7646           0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
    7647           0.22_wp,        &  !< parameter 22  - window fraction ground floor level
    7648           0.0_wp,         &  !< parameter 23  - green fraction ground floor level
    7649           0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
    7650           1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
    7651           2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
    7652           103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
    7653           900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
    7654           0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
    7655           0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
    7656           0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
    7657           0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
    7658           0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
    7659           0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
    7660           0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
    7661           0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
    7662           0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
    7663           27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
    7664           5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
    7665           27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
    7666           0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
    7667           0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
    7668           0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
    7669           0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
    7670           20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
    7671           23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
    7672           20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
    7673           20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
    7674           23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
    7675           10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
    7676           1.0_wp,         &  !< parameter 51  - wall fraction ground plate
    7677           0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
    7678           0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
    7679           0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
    7680           0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
    7681           2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
    7682           103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
    7683           900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
    7684           0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
    7685           0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
    7686           0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
    7687           0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
    7688           0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
    7689           0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
    7690           0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
    7691           27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
    7692           0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
    7693           0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
    7694           0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
    7695           0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
    7696           1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
    7697           1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
    7698           1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
    7699           0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
    7700           0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
    7701           0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
    7702           27.0_wp,        &  !< parameter 77  - window albedo ground floor level
    7703           5.0_wp,         &  !< parameter 78  - green albedo ground floor level
    7704           0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
    7705           0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
    7706           0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
    7707           0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
    7708           1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
    7709           1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
    7710           1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
    7711           0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
    7712           0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
    7713           0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
    7714           1.0_wp,         &  !< parameter 89  - wall fraction roof
    7715           0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
    7716           0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
    7717           0.5_wp,         &  !< parameter 92  - 3rd wall layer thickness roof
    7718           0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
    7719           2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
    7720           103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
    7721           900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
    7722           0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
    7723           0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
    7724           0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
    7725           0.93_wp,        &  !< parameter 100 - wall emissivity roof
    7726           27.0_wp,        &  !< parameter 101 - wall albedo roof
    7727           0.0_wp,         &  !< parameter 102 - window fraction roof
    7728           0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
    7729           0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
    7730           0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
    7731           0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
    7732           1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
    7733           1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
    7734           1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
    7735           0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
    7736           0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
    7737           0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
    7738           0.87_wp,        &  !< parameter 113 - window emissivity roof
    7739           0.7_wp,         &  !< parameter 114 - window transmissivity roof
    7740           27.0_wp,        &  !< parameter 115 - window albedo roof
    7741           0.86_wp,        &  !< parameter 116 - green emissivity roof
    7742           5.0_wp,         &  !< parameter 117 - green albedo roof
    7743           0.0_wp,         &  !< parameter 118 - green type roof
    7744           0.8_wp,         &  !< parameter 119 - shading factor
    7745           0.6_wp,         &  !< parameter 120 - g-value windows
    7746           3.0_wp,         &  !< parameter 121 - u-value windows
    7747           0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp for winter 0.1
    7748           2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp
    7749                              !< for winter 0.5
    7750           0.0_wp,         &  !< parameter 124 - heat recovery efficiency
    7751           2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
    7752           165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
    7753           4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
    7754           100.0_wp,       &  !< parameter 128 - maximal heating capacity
    7755           0.0_wp,         &  !< parameter 129 - maximal cooling capacity
    7756           2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
    7757           6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
    7758           3.0_wp,         &  !< parameter 132 - storey height
    7759           0.2_wp,         &  !< parameter 133 - ceiling construction height
    7760           0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
    7761           1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
     7613        0.73_wp,         &  !< parameter 0   - wall fraction above ground floor level
     7614       0.27_wp,         &  !< parameter 1   - window fraction above ground floor level
     7615       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
     7616       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
     7617       1.5_wp,         &  !< parameter 4   - LAI roof
     7618       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
     7619       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
     7620       103000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
     7621       900000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
     7622       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
     7623       0.38_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
     7624       0.04_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
     7625       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
     7626       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
     7627       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
     7628       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
     7629       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
     7630       0.7_wp,        &  !< parameter 17  - window transmissivity above ground floor level
     7631       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
     7632       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
     7633       4.0_wp,         &  !< parameter 20  - ground floor level height
     7634       0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
     7635       0.22_wp,        &  !< parameter 22  - window fraction ground floor level
     7636       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
     7637       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
     7638       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
     7639       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
     7640       1030000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
     7641       900000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     7642       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     7643       0.38_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     7644       0.04_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     7645       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     7646       0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
     7647       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     7648       0.7_wp,        &  !< parameter 35  - window transmissivity ground floor level
     7649       0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
     7650       0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
     7651       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     7652       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     7653       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     7654       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     7655       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     7656       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     7657       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     7658       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     7659       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     7660       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     7661       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     7662       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     7663       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     7664       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     7665       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     7666       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     7667       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     7668       0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     7669       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     7670       1030000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     7671       900000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
     7672       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     7673       0.38_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     7674       0.04_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     7675       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     7676       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     7677       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     7678       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     7679       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     7680       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     7681       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     7682       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     7683       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     7684       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     7685       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     7686       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     7687       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     7688       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     7689       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     7690       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     7691       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     7692       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     7693       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     7694       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     7695       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     7696       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     7697       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     7698       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     7699       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     7700       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     7701       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     7702       1.0_wp,         &  !< parameter 89  - wall fraction roof
     7703       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     7704       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     7705       0.5_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
     7706       0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     7707       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     7708       1030000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
     7709       900000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
     7710       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     7711       0.38_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     7712       0.04_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
     7713       0.93_wp,        &  !< parameter 100 - wall emissivity roof
     7714       27.0_wp,        &  !< parameter 101 - wall albedo roof
     7715       0.0_wp,         &  !< parameter 102 - window fraction roof
     7716       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     7717       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     7718       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     7719       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     7720       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     7721       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     7722       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     7723       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     7724       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     7725       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     7726       0.87_wp,        &  !< parameter 113 - window emissivity roof
     7727       0.7_wp,        &  !< parameter 114 - window transmissivity roof
     7728       27.0_wp,        &  !< parameter 115 - window albedo roof
     7729       0.86_wp,        &  !< parameter 116 - green emissivity roof
     7730       5.0_wp,         &  !< parameter 117 - green albedo roof
     7731       0.0_wp,         &  !< parameter 118 - green type roof
     7732       0.8_wp,         &  !< parameter 119 - shading factor
     7733       0.6_wp,        &  !< parameter 120 - g-value windows
     7734       3.0_wp,         &  !< parameter 121 - u-value windows
     7735       0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp, winter 0.1
     7736       2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp, winter 0.5
     7737       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     7738       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     7739       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heat storage
     7740       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     7741       100.0_wp,       &  !< parameter 128 - maximal heating capacity
     7742       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     7743       2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     7744       6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
     7745       3.0_wp,         &  !< parameter 132 - storey height
     7746       0.2_wp,         &  !< parameter 133 - ceiling construction height
     7747       0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
     7748       1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
    77627749                           /)
    77637750
Note: See TracChangeset for help on using the changeset viewer.