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

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

File:
1 edited

Legend:

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

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