- Timestamp:
- Sep 21, 2020 7:40:16 PM (4 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
indoor_model_mod.f90 (modified) (23 diffs)
-
time_integration_spinup.f90 (modified) (5 diffs)
-
urban_surface_mod.f90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/indoor_model_mod.f90
r4681 r4687 25 25 ! ----------------- 26 26 ! $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 27 31 ! Bugfix for implementation of downward surfaces 28 32 ! … … 126 130 ! @author Matthias Suehring 127 131 ! @author Sascha RiÃmann 132 ! @author Björn Maronga 128 133 ! 129 134 ! … … 144 149 !> @note How to write indoor temperature output to pt array? 145 150 !> 146 !> @bug <Enter known bugs here>151 !> @bug Calculation of iwghf_eb and iwghf_eb_window is faulty 147 152 !--------------------------------------------------------------------------------------------------! 148 153 MODULE indoor_model_mod … … 234 239 REAL(wp), DIMENSION(:), ALLOCATABLE :: t_in !< mean building indoor temperature, height dependent 235 240 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) 236 243 REAL(wp), DIMENSION(:), ALLOCATABLE :: volume !< total building volume, height dependent 237 244 REAL(wp), DIMENSION(:), ALLOCATABLE :: vol_frac !< fraction of local on total building volume, height dependent … … 319 326 REAL(wp) :: theta_air_set !< [degree_C] Setpoint_temperature for the room 320 327 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 323 329 REAL(wp) :: theta_op !< [degree_C] operative temperature 324 330 REAL(wp) :: theta_s !< [degree_C] surface temperature of the RC-node … … 449 455 !--------------------------------------------------------------------------------------------------! 450 456 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 ) 452 458 453 459 INTEGER(iwp) :: i … … 458 464 REAL(wp) :: near_facade_temperature 459 465 REAL(wp) :: phi_hc_nd_dummy 466 REAL(wp), INTENT(IN) :: theta_m_t_prev 460 467 ! 461 468 !-- Calculation of total mass specific thermal load (internal and external) … … 468 475 ) !< [degree_C] Eq. (C.5) 469 476 ! 470 !-- Calculation of component temperature at factualtimestep477 !-- Calculation of component temperature at current timestep 471 478 theta_m_t = ( ( theta_m_t_prev & 472 479 * ( ( c_m / 3600.0_wp ) - 0.5_wp * ( h_t_3 + h_t_wm ) ) & … … 476 483 ) !< [degree_C] Eq. (C.4) 477 484 ! 478 !-- Calculation of mean inner temperature for the RC-node in actualtimestep485 !-- Calculation of mean inner temperature for the RC-node in current timestep 479 486 theta_m = ( theta_m_t + theta_m_t_prev ) * 0.5_wp !< [degree_C] Eq. (C.9) 480 487 481 488 ! 482 !-- Calculation of mean surface temperature of the RC-node in actualtimestep489 !-- Calculation of mean surface temperature of the RC-node in current timestep 483 490 theta_s = ( ( h_t_ms * theta_m + phi_st + h_t_es * pt(k,j,i) & 484 491 + h_t_1 * ( near_facade_temperature & … … 490 497 ! 491 498 !-- Calculation of the air temperature of the RC-node 499 500 492 501 theta_air = ( h_t_is * theta_s + h_v * near_facade_temperature + phi_ia + phi_hc_nd_dummy ) / & 493 502 ( h_t_is + h_v ) !< [degree_C] Eq. (C.11) 503 494 504 495 505 END SUBROUTINE im_calc_temperatures … … 925 935 ALLOCATE( buildings(nb)%l_v(1:buildings(nb)%num_facades_per_building_v_l) ) 926 936 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) ) 927 940 ENDIF 928 941 … … 1125 1138 ENDDO 1126 1139 ! 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_temperature1130 !1131 1140 !-- Initialize indoor temperature. Actually only for output at initial state. 1132 1141 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 1134 1151 ENDDO 1135 1152 … … 1337 1354 1338 1355 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) ) 1340 1357 ! 1341 1358 !-- If air temperature between border temperatures of heating and cooling, assign output … … 1366 1383 1367 1384 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) ) 1369 1386 theta_air_10 = theta_air !< temperature with 10 W/m2 of heating 1370 1387 phi_hc_nd_un = phi_hc_nd_10 * (theta_air_set - theta_air_0) & … … 1395 1412 !-- Calculate the temperature with phi_hc_nd_ac (new) 1396 1413 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) ) 1398 1415 theta_air_ac = theta_air 1399 1416 ENDIF 1400 1417 ! 1401 1418 !-- 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 1403 1421 1404 1422 q_vent = h_v * ( theta_air - near_facade_temperature ) … … 1408 1426 !-- Will be used for thermal comfort calculations. 1409 1427 theta_op = 0.3_wp * theta_air_ac + 0.7_wp * theta_s !< [degree_C] operative Temperature Eq. (C.12) 1428 1410 1429 ! surf_usm_h(l)%t_indoor(m) = theta_op !< not integrated now 1411 1430 ! … … 1440 1459 ) / facade_element_area !< [W/m2] , observe the directional 1441 1460 !< convention in PALM! 1442 surf_usm_h(l)%waste_heat(m) = q_waste_heat1461 surf_usm_h(l)%waste_heat(m) = 0.0_wp !q_waste_heat 1443 1462 ENDDO !< Horizontal surfaces loop 1444 1463 ! … … 1524 1543 k = surf_usm_v(l)%k(m) 1525 1544 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) & 1528 1547 + surf_usm_v(l)%frac(m,ind_wat_win) * t_window_v(l)%val(nzt_wall,m) 1548 1529 1549 ! 1530 1550 !-- Solar thermal gains. If net_sw_in larger than sun-protection … … 1563 1583 phi_hc_nd = 0.0_wp 1564 1584 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) ) 1566 1586 ! 1567 1587 !-- If air temperature between border temperatures of heating and cooling, assign output … … 1592 1612 1593 1613 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) ) 1595 1615 1596 1616 theta_air_10 = theta_air !< Note the temperature with 10 W/m2 of heating … … 1623 1643 !-- Calculate the temperature with phi_hc_nd_ac (new) 1624 1644 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) ) 1626 1646 theta_air_ac = theta_air 1627 1647 ENDIF 1628 1648 ! 1629 1649 !-- 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 1631 1652 1632 1653 q_vent = h_v * ( theta_air - near_facade_temperature ) … … 1635 1656 !-- Will be used for thermal comfort calculations. 1636 1657 theta_op = 0.3_wp * theta_air_ac + 0.7_wp * theta_s 1658 1637 1659 ! surf_usm_v(l)%t_indoor(m) = theta_op !< not integrated yet 1638 1660 ! … … 1665 1687 ) / facade_element_area !< [W/m2] , observe the directional convention in 1666 1688 !< PALM! 1667 surf_usm_v(l)%waste_heat(m) = q_waste_heat1689 surf_usm_v(l)%waste_heat(m) = 0.0_wp !q_waste_heat 1668 1690 ENDDO !< Vertical surfaces loop 1669 1691 ENDIF !< buildings(nb)%on_pe -
palm/trunk/SOURCE/time_integration_spinup.f90
r4671 r4687 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Indoor model is now available during spinup 28 ! 29 ! 4671 2020-09-09 20:27:58Z pavelkrc 27 30 ! Implementation of downward facing USM and LSM surfaces 28 31 ! … … 107 110 constant_flux_layer, & 108 111 coupling_start_time, & 112 current_timestep_number, & 109 113 data_output_during_spinup, & 110 114 debug_output_timestep, & … … 120 124 dt_3d, & 121 125 humidity, & 126 indoor_model, & 122 127 intermediate_timestep_count, & 123 128 intermediate_timestep_count_max, & … … 161 166 nxrg 162 167 168 USE indoor_model_mod, & 169 ONLY: dt_indoor, im_main_heatcool, time_indoor 170 163 171 USE land_surface_model_mod, & 164 172 ONLY: lsm_energy_balance, & … … 477 485 478 486 ! 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 ! 479 506 !-- Increase simulation time and output times 480 507 current_timestep_number_spinup = current_timestep_number_spinup + 1 -
palm/trunk/SOURCE/urban_surface_mod.f90
r4671 r4687 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Optimized code structure for treatment of inner wall and window heat flux 30 ! 31 ! 4671 2020-09-09 20:27:58Z pavelkrc 29 32 ! Radiative transfer model RTM version 4.1 30 33 ! - Implementation of downward facing USM and LSM surfaces … … 5173 5176 wtend(nzb_wall) = ( 1.0_wp / surf%rho_c_wall(nzb_wall,m) ) & 5174 5177 * ( 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) ) & 5176 5179 * surf%ddz_wall(nzb_wall+1,m) & 5177 5180 + surf%frac(m,ind_veg_wall) & … … 5183 5186 + surf%frac(m,ind_pav_green) ) & 5184 5187 * ( surf%lambda_h_green(nzt_wall,m) & 5185 * wall_mod(nzt_wall) 5188 * wall_mod(nzt_wall) & 5186 5189 * surf%ddz_green(nzt_wall,m) & 5187 5190 + surf%lambda_h(nzb_wall,m) & 5188 * wall_mod(nzb_wall) 5191 * wall_mod(nzb_wall) & 5189 5192 * surf%ddz_wall(nzb_wall,m) ) & 5190 5193 / ( surf%ddz_green(nzt_wall,m) & 5191 5194 + 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) ) & 5193 5196 ) * 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 5197 5206 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) & 5205 5214 ) * surf%ddz_wall_stag(kw,m) 5206 5215 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) & 5212 5221 ) * 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) & 5228 5225 * surf%tt_wall_m(nzb_wall:nzt_wall,m) ) 5229 5226 5230 ! 5227 5231 5228 !-- During spinup the tempeature inside window layers is not calculated to make larger timesteps possible 5232 5229 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) 5235 5231 ! 5236 5232 !-- Prognostic equation for ground/roof window temperature t_window takes absorption 5237 5233 !-- of shortwave radiation into account 5238 5234 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) ) ) & 5247 5243 ) * surf%ddz_window_stag(nzb_wall,m) 5248 5244 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 & 5297 5286 * ( tsc(2) * wintend(nzb_wall:nzt_wall) + tsc(3) & 5298 5287 * surf%tt_window_m(nzb_wall:nzt_wall,m) ) 5299 5288 5300 5289 ENDIF 5301 5302 5290 ! 5303 5291 !-- Calculate t_wall tendencies for the next Runge-Kutta step … … 7623 7611 7624 7612 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 7762 7749 /) 7763 7750
Note: See TracChangeset
for help on using the changeset viewer.
