Index: palm/trunk/SOURCE/indoor_model_mod.f90
===================================================================
--- palm/trunk/SOURCE/indoor_model_mod.f90	(revision 4686)
+++ palm/trunk/SOURCE/indoor_model_mod.f90	(revision 4687)
@@ -25,4 +25,8 @@
 ! -----------------
 ! $Id$
+! Bugfix: values of theta_m_t_prev were not stored for individual surfaces and thus re-used by all
+! surfaces and buildings, which led to excessive indoor temperatures
+!
+! 4681 2020-09-16 10:23:06Z pavelkrc
 ! Bugfix for implementation of downward surfaces
 !
@@ -126,4 +130,5 @@
 ! @author Matthias Suehring
 ! @author Sascha RiÃmann
+! @author BjÃ¶rn Maronga
 !
 !
@@ -144,5 +149,5 @@
 !> @note How to write indoor temperature output to pt array?
 !>
-!> @bug  <Enter known bugs here>
+!> @bug  Calculation of iwghf_eb and iwghf_eb_window is faulty
 !--------------------------------------------------------------------------------------------------!
  MODULE indoor_model_mod
@@ -234,4 +239,6 @@
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  t_in       !< mean building indoor temperature, height dependent
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  t_in_l     !< mean building indoor temperature on local subdomain, height dependent
+       REAL(wp), DIMENSION(:), ALLOCATABLE ::  theta_m_t_prev_h !< [degree_C] value of theta_m_t from previous time step (horizontal)
+       REAL(wp), DIMENSION(:), ALLOCATABLE ::  theta_m_t_prev_v !< [degree_C] value of theta_m_t from previous time step (vertical)
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  volume     !< total building volume, height dependent
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  vol_frac   !< fraction of local on total building volume, height dependent
@@ -319,6 +326,5 @@
     REAL(wp) ::  theta_air_set                       !< [degree_C] Setpoint_temperature for the room
     REAL(wp) ::  theta_m                             !< [degree_C} inner temperature of the RC-node
-    REAL(wp) ::  theta_m_t                           !< [degree_C] (Fictive) component temperature timestep
-    REAL(wp) ::  theta_m_t_prev                      !< [degree_C] (Fictive) component temperature previous timestep (do not change)
+    REAL(wp) ::  theta_m_t                           !< [degree_C] (Fictive) component temperature during timestep
     REAL(wp) ::  theta_op                            !< [degree_C] operative temperature
     REAL(wp) ::  theta_s                             !< [degree_C] surface temperature of the RC-node
@@ -449,5 +455,5 @@
 !--------------------------------------------------------------------------------------------------!
  SUBROUTINE im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                        &
-                                   near_facade_temperature, phi_hc_nd_dummy )
+                                   near_facade_temperature, phi_hc_nd_dummy, theta_m_t_prev )
 
     INTEGER(iwp) ::  i
@@ -458,4 +464,5 @@
     REAL(wp) ::  near_facade_temperature
     REAL(wp) ::  phi_hc_nd_dummy
+    REAL(wp), INTENT(IN) :: theta_m_t_prev
 !
 !-- Calculation of total mass specific thermal load (internal and external)
@@ -468,5 +475,5 @@
                )                                                                !< [degree_C] Eq. (C.5)
 !
-!-- Calculation of component temperature at factual timestep
+!-- Calculation of component temperature at current timestep
     theta_m_t = ( ( theta_m_t_prev                                                                 &
                     * ( ( c_m / 3600.0_wp ) - 0.5_wp * ( h_t_3 + h_t_wm ) )                        &
@@ -476,9 +483,9 @@
                 )                                                               !< [degree_C] Eq. (C.4)
 !
-!-- Calculation of mean inner temperature for the RC-node in actual timestep
+!-- Calculation of mean inner temperature for the RC-node in current timestep
     theta_m = ( theta_m_t + theta_m_t_prev ) * 0.5_wp                           !< [degree_C] Eq. (C.9)
 
 !
-!-- Calculation of mean surface temperature of the RC-node in actual timestep
+!-- Calculation of mean surface temperature of the RC-node in current timestep
     theta_s = ( (   h_t_ms * theta_m + phi_st + h_t_es * pt(k,j,i)                                 &
                   + h_t_1  * ( near_facade_temperature                                             &
@@ -490,6 +497,9 @@
 !
 !-- Calculation of the air temperature of the RC-node
+
+
     theta_air = ( h_t_is * theta_s + h_v * near_facade_temperature + phi_ia + phi_hc_nd_dummy ) /  &
                 ( h_t_is + h_v )                                                !< [degree_C] Eq. (C.11)
+
 
  END SUBROUTINE im_calc_temperatures
@@ -925,4 +935,7 @@
           ALLOCATE( buildings(nb)%l_v(1:buildings(nb)%num_facades_per_building_v_l) )
           ALLOCATE( buildings(nb)%m_v(1:buildings(nb)%num_facades_per_building_v_l) )
+
+          ALLOCATE( buildings(nb)%theta_m_t_prev_h(1:buildings(nb)%num_facades_per_building_h_l) )
+          ALLOCATE( buildings(nb)%theta_m_t_prev_v(1:buildings(nb)%num_facades_per_building_v_l) )
        ENDIF
 
@@ -1125,11 +1138,15 @@
     ENDDO
 !
-!-- Initial room temperature [K]
-!-- (after first loop, use theta_m_t as theta_m_t_prev)
-    theta_m_t_prev = initial_indoor_temperature
-!
 !-- Initialize indoor temperature. Actually only for output at initial state.
     DO  nb = 1, num_build
-       IF ( buildings(nb)%on_pe )  buildings(nb)%t_in(:) = initial_indoor_temperature
+       IF ( buildings(nb)%on_pe )  THEN
+          buildings(nb)%t_in(:) = initial_indoor_temperature
+
+!
+!--       (after first loop, use theta_m_t as theta_m_t_prev)
+          buildings(nb)%theta_m_t_prev_h(:) = initial_indoor_temperature
+          buildings(nb)%theta_m_t_prev_v(:) = initial_indoor_temperature
+
+       ENDIF
     ENDDO
 
@@ -1337,5 +1354,5 @@
 
              CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                 &
-                                          near_facade_temperature, phi_hc_nd )
+                                          near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
 !
 !--          If air temperature between border temperatures of heating and cooling, assign output
@@ -1366,5 +1383,5 @@
 
                 CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
-                                             near_facade_temperature, phi_hc_nd )
+                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
                 theta_air_10 = theta_air                                                !< temperature with 10 W/m2 of heating
                 phi_hc_nd_un = phi_hc_nd_10 * (theta_air_set - theta_air_0)                        &
@@ -1395,10 +1412,11 @@
 !--             Calculate the temperature with phi_hc_nd_ac (new)
                 CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
-                                             near_facade_temperature, phi_hc_nd )
+                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_h(fa) )
                 theta_air_ac = theta_air
              ENDIF
 !
 !--          Update theta_m_t_prev
-             theta_m_t_prev = theta_m_t
+             buildings(nb)%theta_m_t_prev_h(fa) = theta_m_t
+
 
              q_vent = h_v * ( theta_air - near_facade_temperature )
@@ -1408,4 +1426,5 @@
 !--          Will be used for thermal comfort calculations.
              theta_op     = 0.3_wp * theta_air_ac + 0.7_wp * theta_s          !< [degree_C] operative Temperature Eq. (C.12)
+
 !              surf_usm_h(l)%t_indoor(m) = theta_op                               !< not integrated now
 !
@@ -1440,5 +1459,5 @@
                             ) / facade_element_area                                             !< [W/m2] , observe the directional
                                                                                                 !< convention in PALM!
-             surf_usm_h(l)%waste_heat(m) = q_waste_heat
+             surf_usm_h(l)%waste_heat(m) = 0.0_wp !q_waste_heat
           ENDDO !< Horizontal surfaces loop
 !
@@ -1524,7 +1543,8 @@
              k = surf_usm_v(l)%k(m)
              near_facade_temperature = surf_usm_v(l)%pt_10cm(m)
-             indoor_wall_window_temperature =                                                      &
-                                    surf_usm_v(l)%frac(m,ind_veg_wall) * t_wall_v(l)%val(nzt_wall,m) &
+             indoor_wall_window_temperature =                                                          &
+                                    surf_usm_v(l)%frac(m,ind_veg_wall) * t_wall_v(l)%val(nzt_wall,m)   &
                                   + surf_usm_v(l)%frac(m,ind_wat_win)  * t_window_v(l)%val(nzt_wall,m)
+
 !
 !--          Solar thermal gains. If net_sw_in larger than sun-protection
@@ -1563,5 +1583,5 @@
              phi_hc_nd = 0.0_wp
              CALL im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,                  &
-                                         near_facade_temperature, phi_hc_nd )
+                                         near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
 !
 !--          If air temperature between border temperatures of heating and cooling, assign output
@@ -1592,5 +1612,5 @@
 
                 CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature,              &
-                                             near_facade_temperature, phi_hc_nd )
+                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
 
                 theta_air_10 = theta_air !< Note the temperature with 10 W/m2 of heating
@@ -1623,10 +1643,11 @@
 !--             Calculate the temperature with phi_hc_nd_ac (new)
                 CALL  im_calc_temperatures ( i, j, k, indoor_wall_window_temperature, &
-                                             near_facade_temperature, phi_hc_nd )
+                                             near_facade_temperature, phi_hc_nd, buildings(nb)%theta_m_t_prev_v(fa) )
                 theta_air_ac = theta_air
              ENDIF
 !
 !--          Update theta_m_t_prev
-             theta_m_t_prev = theta_m_t
+             buildings(nb)%theta_m_t_prev_v(fa) = theta_m_t
+
 
              q_vent = h_v * ( theta_air - near_facade_temperature )
@@ -1635,4 +1656,5 @@
 !--          Will be used for thermal comfort calculations.
              theta_op     = 0.3_wp * theta_air_ac + 0.7_wp * theta_s
+
 !              surf_usm_v(l)%t_indoor(m) = theta_op                  !< not integrated yet
 !
@@ -1665,5 +1687,5 @@
                                                     ) / facade_element_area  !< [W/m2] , observe the directional convention in
                                                                              !< PALM!
-             surf_usm_v(l)%waste_heat(m) = q_waste_heat
+             surf_usm_v(l)%waste_heat(m) = 0.0_wp !q_waste_heat
           ENDDO !< Vertical surfaces loop
        ENDIF !< buildings(nb)%on_pe
Index: palm/trunk/SOURCE/time_integration_spinup.f90
===================================================================
--- palm/trunk/SOURCE/time_integration_spinup.f90	(revision 4686)
+++ palm/trunk/SOURCE/time_integration_spinup.f90	(revision 4687)
@@ -25,4 +25,7 @@
 ! -----------------
 ! $Id$
+! Indoor model is now available during spinup
+!
+! 4671 2020-09-09 20:27:58Z pavelkrc
 ! Implementation of downward facing USM and LSM surfaces
 !
@@ -107,4 +110,5 @@
                constant_flux_layer,                                                                &
                coupling_start_time,                                                                &
+               current_timestep_number,                                                            &
                data_output_during_spinup,                                                          &
                debug_output_timestep,                                                              &
@@ -120,4 +124,5 @@
                dt_3d,                                                                              &
                humidity,                                                                           &
+               indoor_model,                                                                       &
                intermediate_timestep_count,                                                        &
                intermediate_timestep_count_max,                                                    &
@@ -161,4 +166,7 @@
                nxrg
 
+    USE indoor_model_mod,                                                                          &
+        ONLY:  dt_indoor, im_main_heatcool, time_indoor
+
     USE land_surface_model_mod,                                                                    &
         ONLY:  lsm_energy_balance,                                                                 &
@@ -477,4 +485,23 @@
 
 !
+!--    If required, calculate indoor temperature, waste heat, heat flux
+!--    through wall, etc.
+!--    dt_indoor steers the frequency of the indoor model calculations.
+!--    Note, at first timestep indoor model is called, in order to provide
+!--    a waste heat flux.
+       IF ( indoor_model )  THEN
+
+          time_indoor = time_indoor + dt_3d
+
+          IF ( time_indoor >= dt_indoor  .OR.  current_timestep_number == 0 )  THEN
+
+             time_indoor = time_indoor - dt_indoor
+
+             CALL im_main_heatcool
+
+          ENDIF
+       ENDIF
+
+!
 !--    Increase simulation time and output times
        current_timestep_number_spinup = current_timestep_number_spinup + 1
Index: palm/trunk/SOURCE/urban_surface_mod.f90
===================================================================
--- palm/trunk/SOURCE/urban_surface_mod.f90	(revision 4686)
+++ palm/trunk/SOURCE/urban_surface_mod.f90	(revision 4687)
@@ -27,4 +27,7 @@
 ! -----------------
 ! $Id$
+! Optimized code structure for treatment of inner wall and window heat flux
+!
+! 4671 2020-09-09 20:27:58Z pavelkrc
 ! Radiative transfer model RTM version 4.1
 ! - Implementation of downward facing USM and LSM surfaces
@@ -5173,5 +5176,5 @@
        wtend(nzb_wall) = ( 1.0_wp / surf%rho_c_wall(nzb_wall,m) )                     &
                           * ( surf%lambda_h(nzb_wall,m) * wall_mod(nzb_wall)          &
-                              * ( t_wall%val(nzb_wall+1,m) - t_wall%val(nzb_wall,m) )              &
+                              * ( t_wall%val(nzb_wall+1,m) - t_wall%val(nzb_wall,m) ) &
                               * surf%ddz_wall(nzb_wall+1,m)                           &
                               + surf%frac(m,ind_veg_wall)                             &
@@ -5183,121 +5186,106 @@
                                   + surf%frac(m,ind_pav_green) )                      &
                               * ( surf%lambda_h_green(nzt_wall,m)                     &
-                              * wall_mod(nzt_wall)                                             &
+                              * wall_mod(nzt_wall)                                    &
                               * surf%ddz_green(nzt_wall,m)                            &
                               + surf%lambda_h(nzb_wall,m)                             &
-                              * wall_mod(nzb_wall)                                             &
+                              * wall_mod(nzb_wall)                                    &
                               * surf%ddz_wall(nzb_wall,m) )                           &
                               / ( surf%ddz_green(nzt_wall,m)                          &
                               + surf%ddz_wall(nzb_wall,m) )                           &
-                              * ( t_wall%val(nzb_wall,m) - t_green%val(nzt_wall,m) )               &
+                              * ( t_wall%val(nzb_wall,m) - t_green%val(nzt_wall,m) )  &
                             ) * surf%ddz_wall_stag(nzb_wall,m)
-!
-!--    If indoor model is used inner wall layer is calculated by using iwghf (indoor wall ground heat flux)
-       IF ( indoor_model )  THEN
+
+!
+!--    If no indoor model is used, calculate the inner wall heat flux according to
+!--    heat diffusion (innermost wall temperature as boundary conditions)
+       IF ( .NOT. indoor_model ) THEN
+          surf%iwghf_eb(m) = - surf%lambda_h(nzt_wall-1,m) * ( t_wall%val(nzt_wall,m) &
+                             - t_wall%val(nzt_wall-1,m) ) * surf%ddz_wall(nzt_wall,m)
+          ENDIF
+
           DO  kw = nzb_wall+1, nzt_wall-1
-             wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                           &
-                         * ( surf%lambda_h(kw,m) * wall_mod(kw)                       &
-                             * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                           &
-                             * surf%ddz_wall(kw+1,m)                                  &
-                             - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                 &
-                             * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                           &
-                             * surf%ddz_wall(kw,m)                                    &
+             wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                             &
+                         * ( surf%lambda_h(kw,m) * wall_mod(kw)                         &
+                             * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                &
+                             * surf%ddz_wall(kw+1,m)                                    &
+                             - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                   &
+                             * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                &
+                             * surf%ddz_wall(kw,m)                                      &
                            ) * surf%ddz_wall_stag(kw,m)
           ENDDO
-          wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                  &
-                            * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)   &
-                                * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) )            &
-                                * surf%ddz_wall(nzt_wall,m)                           &
-                                + surf%iwghf_eb(m)                                    &
+          wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                    &
+                            * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)     &
+                                * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
+                                * surf%ddz_wall(nzt_wall,m)                             &
+                                + surf%iwghf_eb(m)                                      &
                               ) * surf%ddz_wall_stag(nzt_wall,m)
-       ELSE
-          DO  kw = nzb_wall+1, nzt_wall
-             wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                           &
-                        * ( surf%lambda_h(kw,m)  * wall_mod(kw)                       &
-                            * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                            &
-                            * surf%ddz_wall(kw+1,m)                                   &
-                            - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                  &
-                            * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                            &
-                            * surf%ddz_wall(kw,m)                                     &
-                          ) * surf%ddz_wall_stag(kw,m)
-          ENDDO
-       ENDIF
-
-       t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d                 &
-                                         * ( tsc(2) * wtend(nzb_wall:nzt_wall) + tsc(3)        &
+
+       t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d      &
+                                         * ( tsc(2) * wtend(nzb_wall:nzt_wall) + tsc(3) &
                                              * surf%tt_wall_m(nzb_wall:nzt_wall,m) )
 
-!
+
 !--    During spinup the tempeature inside window layers is not calculated to make larger timesteps possible
        IF ( .NOT. during_spinup )  THEN
-          win_absorp = -log( surf%transmissivity(m) ) /                                &
-                        surf%zw_window(nzt_wall,m)
+          win_absorp = - log( surf%transmissivity(m) ) / surf%zw_window(nzt_wall,m)
 !
 !--       Prognostic equation for ground/roof window temperature t_window takes absorption
 !--       of shortwave radiation into account
           wintend(:) = 0.0_wp
-          wintend(nzb_wall) = ( 1.0_wp / surf%rho_c_window(nzb_wall,m) )               &
-                                * ( surf%lambda_h_window(nzb_wall,m)                   &
-                                * ( t_window%val(nzb_wall+1,m) - t_window%val(nzb_wall,m) )         &
-                                * surf%ddz_window(nzb_wall+1,m)                        &
-                                + surf%wghf_eb_window(m)                               &
-                                + surf%rad_sw_in(m)                                    &
-                                * ( 1.0_wp - exp( -win_absorp                                   &
-                                    * surf%zw_window(nzb_wall,m) ) )                   &
+          wintend(nzb_wall) = ( 1.0_wp / surf%rho_c_window(nzb_wall,m) )                    &
+                                * ( surf%lambda_h_window(nzb_wall,m)                        &
+                                * ( t_window%val(nzb_wall+1,m) - t_window%val(nzb_wall,m) ) &
+                                * surf%ddz_window(nzb_wall+1,m)                             &
+                                + surf%wghf_eb_window(m)                                    &
+                                + surf%rad_sw_in(m)                                         &
+                                * ( 1.0_wp - exp( -win_absorp                               &
+                                    * surf%zw_window(nzb_wall,m) ) )                        &
                               ) * surf%ddz_window_stag(nzb_wall,m)
 
-          IF ( indoor_model ) THEN
-             DO  kw = nzb_wall+1, nzt_wall-1
-                wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                     &
-                              * ( surf%lambda_h_window(kw,m)                           &
-                                  * ( t_window%val(kw+1,m) - t_window%val(kw,m) )                   &
-                                  * surf%ddz_window(kw+1,m)                            &
-                                  - surf%lambda_h_window(kw-1,m)                       &
-                                 * ( t_window%val(kw,m) - t_window%val(kw-1,m) )                    &
-                                 * surf%ddz_window(kw,m)                               &
-                                 + surf%rad_sw_in(m)                                   &
-                                 * ( exp( -win_absorp * surf%zw_window(kw-1,m) )       &
-                                     - exp(-win_absorp * surf%zw_window(kw,m) )        &
-                                   )                                                            &
-                                 ) * surf%ddz_window_stag(kw,m)
-
-             ENDDO
-             wintend(nzt_wall) = ( 1.0_wp / surf%rho_c_window(nzt_wall,m) )            &
-                                 * ( -surf%lambda_h_window(nzt_wall-1,m)               &
-                                    * ( t_window%val(nzt_wall,m) - t_window%val(nzt_wall-1,m) )     &
-                                    * surf%ddz_window(nzt_wall,m)                      &
-                                    + surf%iwghf_eb_window(m)                          &
-                                    + surf%rad_sw_in(m)                                &
-                                    * ( exp( -win_absorp                                        &
-                                        * surf%zw_window(nzt_wall-1,m) )               &
-                                      - exp( -win_absorp                                        &
-                                        * surf%zw_window(nzt_wall,m) )                 &
-                                      )                                                         &
-                                    ) * surf%ddz_window_stag(nzt_wall,m)
-          ELSE
-             DO  kw = nzb_wall+1, nzt_wall
-                wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                     &
-                              * ( surf%lambda_h_window(kw,m)                           &
-                                  * ( t_window%val(kw+1,m) - t_window%val(kw,m) )                   &
-                                  * surf%ddz_window(kw+1,m)                            &
-                                  - surf%lambda_h_window(kw-1,m)                       &
-                                  * ( t_window%val(kw,m)                                          &
-                                  - t_window%val(kw-1,m) )                                        &
-                                  * surf%ddz_window(kw,m)                              &
-                                  + surf%rad_sw_in(m)                                  &
-                                   * ( exp( -win_absorp * surf%zw_window(kw-1,m) )     &
-                                       - exp(-win_absorp * surf%zw_window(kw,m) )      &
-                                     )                                                          &
-                                 ) * surf%ddz_window_stag(kw,m)
-
-             ENDDO
-          ENDIF
-
-          t_window_p%val(nzb_wall:nzt_wall,m) = t_window%val(nzb_wall:nzt_wall,m) + dt_3d           &
+!
+!--       If no indoor model is used, calculate the inner window heat flux according to
+!--       heat diffusion (innermost window temperature as boundary conditions
+          IF ( .NOT. indoor_model ) THEN
+             surf%iwghf_eb_window(m) = - surf%lambda_h_window(nzt_wall-1,m)            &
+                                       * ( t_window%val(nzt_wall,m)                    &
+                                       - t_window%val(nzt_wall-1,m) )                  &
+                                       * surf%ddz_window(nzt_wall,m)
+          ENDIF
+
+
+          DO  kw = nzb_wall+1, nzt_wall-1
+             wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                        &
+                           * ( surf%lambda_h_window(kw,m)                              &
+                           * ( t_window%val(kw+1,m) - t_window%val(kw,m) )             &
+                           * surf%ddz_window(kw+1,m)                                   &
+                           - surf%lambda_h_window(kw-1,m)                              &
+                           * ( t_window%val(kw,m) - t_window%val(kw-1,m) )             &
+                           * surf%ddz_window(kw,m)                                     &
+                           + surf%rad_sw_in(m)                                         &
+                           * ( exp( -win_absorp * surf%zw_window(kw-1,m) )             &
+                           - exp(-win_absorp * surf%zw_window(kw,m) )                  &
+                             )                                                         &
+                             ) * surf%ddz_window_stag(kw,m)
+
+          ENDDO
+
+          wintend(nzt_wall) = ( 1.0_wp / surf%rho_c_window(nzt_wall,m) )                     &
+                                 * ( -surf%lambda_h_window(nzt_wall-1,m)                     &
+                                 * ( t_window%val(nzt_wall,m) - t_window%val(nzt_wall-1,m) ) &
+                                 * surf%ddz_window(nzt_wall,m)                               &
+                                 + surf%iwghf_eb_window(m)                                   &
+                                 + surf%rad_sw_in(m)                                         &
+                                 * ( exp( -win_absorp                                        &
+                                 * surf%zw_window(nzt_wall-1,m) )                            &
+                                 - exp( -win_absorp                                          &
+                                 * surf%zw_window(nzt_wall,m) )                              &
+                                   )                                                         &
+                                   ) * surf%ddz_window_stag(nzt_wall,m)
+
+          t_window_p%val(nzb_wall:nzt_wall,m) = t_window%val(nzb_wall:nzt_wall,m) + dt_3d       &
                                               * ( tsc(2) * wintend(nzb_wall:nzt_wall) + tsc(3)  &
                                               * surf%tt_window_m(nzb_wall:nzt_wall,m) )
 
        ENDIF
-
 !
 !--    Calculate t_wall tendencies for the next Runge-Kutta step
@@ -7623,141 +7611,140 @@
 
     building_pars(:,2) = (/                                                                        &
-          0.73_wp,        &  !< parameter 0   - wall fraction above ground floor level
-          0.27_wp,        &  !< parameter 1   - window fraction above ground floor level
-          0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
-          0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
-          1.5_wp,         &  !< parameter 4   - LAI roof
-          1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
-          2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
-          103000.0_wp,    &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
-          900000.0_wp,    &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
-          0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
-          0.38_wp,        &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
-          0.04_wp,        &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
-          299.15_wp,      &  !< parameter 12  - indoor target summer temperature
-          293.15_wp,      &  !< parameter 13  - indoor target winter temperature
-          0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
-          0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
-          0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
-          0.7_wp,         &  !< parameter 17  - window transmissivity above ground floor level
-          0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
-          0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
-          4.0_wp,         &  !< parameter 20  - ground floor level height
-          0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
-          0.22_wp,        &  !< parameter 22  - window fraction ground floor level
-          0.0_wp,         &  !< parameter 23  - green fraction ground floor level
-          0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
-          1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
-          2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
-          103000.0_wp,    &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
-          900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
-          0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
-          0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
-          0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
-          0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
-          0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
-          0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
-          0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
-          0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
-          0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
-          27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
-          5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
-          27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
-          0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
-          0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
-          0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
-          0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
-          20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
-          23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
-          20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
-          20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
-          23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
-          10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
-          1.0_wp,         &  !< parameter 51  - wall fraction ground plate
-          0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
-          0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
-          0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
-          0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
-          2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
-          103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
-          900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
-          0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
-          0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
-          0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
-          0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
-          0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
-          0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
-          0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
-          27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
-          0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
-          0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
-          0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
-          0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
-          1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
-          1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
-          1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
-          0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
-          0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
-          0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
-          27.0_wp,        &  !< parameter 77  - window albedo ground floor level
-          5.0_wp,         &  !< parameter 78  - green albedo ground floor level
-          0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
-          0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
-          0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
-          0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
-          1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
-          1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
-          1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
-          0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
-          0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
-          0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
-          1.0_wp,         &  !< parameter 89  - wall fraction roof
-          0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
-          0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
-          0.5_wp,         &  !< parameter 92  - 3rd wall layer thickness roof
-          0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
-          2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
-          103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
-          900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
-          0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
-          0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
-          0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
-          0.93_wp,        &  !< parameter 100 - wall emissivity roof
-          27.0_wp,        &  !< parameter 101 - wall albedo roof
-          0.0_wp,         &  !< parameter 102 - window fraction roof
-          0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
-          0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
-          0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
-          0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
-          1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
-          1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
-          1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
-          0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
-          0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
-          0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
-          0.87_wp,        &  !< parameter 113 - window emissivity roof
-          0.7_wp,         &  !< parameter 114 - window transmissivity roof
-          27.0_wp,        &  !< parameter 115 - window albedo roof
-          0.86_wp,        &  !< parameter 116 - green emissivity roof
-          5.0_wp,         &  !< parameter 117 - green albedo roof
-          0.0_wp,         &  !< parameter 118 - green type roof
-          0.8_wp,         &  !< parameter 119 - shading factor
-          0.6_wp,         &  !< parameter 120 - g-value windows
-          3.0_wp,         &  !< parameter 121 - u-value windows
-          0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp for winter 0.1
-          2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp
-                             !< for winter 0.5
-          0.0_wp,         &  !< parameter 124 - heat recovery efficiency
-          2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
-          165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
-          4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
-          100.0_wp,       &  !< parameter 128 - maximal heating capacity
-          0.0_wp,         &  !< parameter 129 - maximal cooling capacity
-          2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
-          6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
-          3.0_wp,         &  !< parameter 132 - storey height
-          0.2_wp,         &  !< parameter 133 - ceiling construction height
-          0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
-          1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
+        0.73_wp,         &  !< parameter 0   - wall fraction above ground floor level
+       0.27_wp,         &  !< parameter 1   - window fraction above ground floor level
+       0.0_wp,         &  !< parameter 2   - green fraction above ground floor level
+       0.0_wp,         &  !< parameter 3   - green fraction roof above ground floor level
+       1.5_wp,         &  !< parameter 4   - LAI roof
+       1.5_wp,         &  !< parameter 5   - LAI on wall above ground floor level
+       2000000.0_wp,   &  !< parameter 6   - heat capacity 1st/2nd wall layer above ground floor level
+       103000.0_wp,   &  !< parameter 7   - heat capacity 3rd wall layer above ground floor level
+       900000.0_wp,   &  !< parameter 8   - heat capacity 4th wall layer above ground floor level
+       0.35_wp,        &  !< parameter 9   - thermal conductivity 1st/2nd wall layer above ground floor level
+       0.38_wp,         &  !< parameter 10  - thermal conductivity 3rd wall layer above ground floor level
+       0.04_wp,         &  !< parameter 11  - thermal conductivity 4th wall layer above ground floor level
+       299.15_wp,      &  !< parameter 12  - indoor target summer temperature
+       293.15_wp,      &  !< parameter 13  - indoor target winter temperature
+       0.92_wp,        &  !< parameter 14  - wall emissivity above ground floor level
+       0.86_wp,        &  !< parameter 15  - green emissivity above ground floor level
+       0.87_wp,        &  !< parameter 16  - window emissivity above ground floor level
+       0.7_wp,        &  !< parameter 17  - window transmissivity above ground floor level
+       0.001_wp,       &  !< parameter 18  - z0 roughness above ground floor level
+       0.0001_wp,      &  !< parameter 19  - z0h/z0g roughness heat/humidity above ground floor level
+       4.0_wp,         &  !< parameter 20  - ground floor level height
+       0.78_wp,        &  !< parameter 21  - wall fraction ground floor level
+       0.22_wp,        &  !< parameter 22  - window fraction ground floor level
+       0.0_wp,         &  !< parameter 23  - green fraction ground floor level
+       0.0_wp,         &  !< parameter 24  - green fraction roof ground floor level
+       1.5_wp,         &  !< parameter 25  - LAI on wall ground floor level
+       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
+       1030000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
+       900000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
+       0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
+       0.38_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
+       0.04_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
+       0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
+       0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
+       0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
+       0.7_wp,        &  !< parameter 35  - window transmissivity ground floor level
+       0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
+       0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
+       27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
+       5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
+       27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
+       0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
+       0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
+       0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
+       0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
+       20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
+       23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
+       20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
+       20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
+       23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
+       10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
+       1.0_wp,         &  !< parameter 51  - wall fraction ground plate
+       0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
+       0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
+       0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
+       0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
+       2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
+       1030000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
+       900000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
+       0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
+       0.38_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
+       0.04_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
+       0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
+       0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
+       0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
+       0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
+       27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
+       0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
+       0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
+       0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
+       0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
+       1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
+       1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
+       0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
+       0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
+       0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
+       27.0_wp,        &  !< parameter 77  - window albedo ground floor level
+       5.0_wp,         &  !< parameter 78  - green albedo ground floor level
+       0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
+       0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
+       0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
+       0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
+       1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
+       1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
+       0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
+       0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
+       0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
+       1.0_wp,         &  !< parameter 89  - wall fraction roof
+       0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
+       0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
+       0.5_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
+       0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
+       2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
+       1030000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
+       900000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
+       0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
+       0.38_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
+       0.04_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
+       0.93_wp,        &  !< parameter 100 - wall emissivity roof
+       27.0_wp,        &  !< parameter 101 - wall albedo roof
+       0.0_wp,         &  !< parameter 102 - window fraction roof
+       0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
+       0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
+       0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
+       0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
+       1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
+       1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
+       1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
+       0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
+       0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
+       0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
+       0.87_wp,        &  !< parameter 113 - window emissivity roof
+       0.7_wp,        &  !< parameter 114 - window transmissivity roof
+       27.0_wp,        &  !< parameter 115 - window albedo roof
+       0.86_wp,        &  !< parameter 116 - green emissivity roof
+       5.0_wp,         &  !< parameter 117 - green albedo roof
+       0.0_wp,         &  !< parameter 118 - green type roof
+       0.8_wp,         &  !< parameter 119 - shading factor
+       0.6_wp,        &  !< parameter 120 - g-value windows
+       3.0_wp,         &  !< parameter 121 - u-value windows
+       0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp, winter 0.1
+       2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp, winter 0.5
+       0.0_wp,         &  !< parameter 124 - heat recovery efficiency
+       2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
+       370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heat storage
+       4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
+       100.0_wp,       &  !< parameter 128 - maximal heating capacity
+       0.0_wp,         &  !< parameter 129 - maximal cooling capacity
+       2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
+       6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
+       3.0_wp,         &  !< parameter 132 - storey height
+       0.2_wp,         &  !< parameter 133 - ceiling construction height
+       0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
+       1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
                            /)
 
