Changeset 4698 for palm/trunk/SOURCE


Ignore:
Timestamp:
Sep 25, 2020 8:37:55 AM (4 years ago)
Author:
maronga
Message:

bugfix in urban surface model (introduced in previous revision)

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4687 r4698  
    2525! -----------------
    2626! $Id$
     27! Fixed faulty characters
     28!
     29! 4687 2020-09-21 19:40:16Z maronga
    2730! Bugfix: values of theta_m_t_prev were not stored for individual surfaces and thus re-used by all
    2831! surfaces and buildings, which led to excessive indoor temperatures
     
    129132! @author Farah Kanani-Suehring
    130133! @author Matthias Suehring
    131 ! @author Sascha Rißmann
    132 ! @author Björn Maronga
     134! @author Sascha Rissmann
     135! @author Bjoern Maronga
    133136!
    134137!
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r4694 r4698  
    2727! -----------------
    2828! $Id$
     29! Bugfix in calculation of iwghf_eb and iwghf_eb_window (introduced in previous revisions)
     30!
     31! 4694 2020-09-23 15:09:19Z pavelkrc
    2932! Fix writing and reading of surface data to/from MPI restart file
    3033!
     
    3538! Bugfix for previous revision
    3639!
    37 ! 4687 2020-09-21 19:40:16Z maronga
     40! 4687 2020-09-21 19:40:16Z
    3841! Optimized code structure for treatment of inner wall and window heat flux
    3942!
     
    52065209
    52075210!
    5208 !--    If no indoor model is used, calculate the inner wall heat flux according to
    5209 !--    heat diffusion (innermost wall temperature as boundary conditions)
     5211!--    If indoor model is used inner wall layer is calculated by using iwghf (indoor
     5212!--    wall ground heat flux)
    52105213       IF ( .NOT. indoor_model ) THEN
    5211           surf%iwghf_eb(m) = - surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)           &
    5212                              * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) )  &
    5213                              * surf%ddz_wall(nzt_wall,m)
     5214          surf%iwghf_eb(m) = surf%lambda_h(nzt_wall,m)  * wall_mod(nzt_wall)          &
     5215                             * ( t_wall%val(nzt_wall+1,m) - t_wall%val(nzt_wall,m) )  &
     5216                             * surf%ddz_wall(nzt_wall+1,m)
    52145217       ENDIF
    52155218
    52165219       DO  kw = nzb_wall+1, nzt_wall-1
    5217           wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                         &
    5218                       * ( surf%lambda_h(kw,m) * wall_mod(kw)                     &
    5219                       * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )                &
    5220                       * surf%ddz_wall(kw+1,m)                                    &
    5221                       - surf%lambda_h(kw-1,m) * wall_mod(kw-1)                   &
    5222                       * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )                &
    5223                       * surf%ddz_wall(kw,m)                                      &
    5224                       ) * surf%ddz_wall_stag(kw,m)
     5220          wtend(kw) = ( 1.0_wp / surf%rho_c_wall(kw,m) )                   &
     5221                         * ( surf%lambda_h(kw,m) * wall_mod(kw)            &
     5222                             * ( t_wall%val(kw+1,m) - t_wall%val(kw,m) )   &
     5223                             * surf%ddz_wall(kw+1,m)                       &
     5224                             - surf%lambda_h(kw-1,m) * wall_mod(kw-1)      &
     5225                             * ( t_wall%val(kw,m) - t_wall%val(kw-1,m) )   &
     5226                             * surf%ddz_wall(kw,m)                         &
     5227                           ) * surf%ddz_wall_stag(kw,m)
    52255228       ENDDO
    5226        wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                    &
    5227                          * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)     &
    5228                              * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
    5229                               * surf%ddz_wall(nzt_wall,m) + surf%iwghf_eb(m)         &
    5230                            ) * surf%ddz_wall_stag(nzt_wall,m)
     5229
     5230       wtend(nzt_wall) = ( 1.0_wp / surf%rho_c_wall(nzt_wall,m) )                       &
     5231                            * ( -surf%lambda_h(nzt_wall-1,m) * wall_mod(nzt_wall-1)     &
     5232                                * ( t_wall%val(nzt_wall,m) - t_wall%val(nzt_wall-1,m) ) &
     5233                                * surf%ddz_wall(nzt_wall,m)                             &
     5234                                + surf%iwghf_eb(m)                                      &
     5235                              ) * surf%ddz_wall_stag(nzt_wall,m)
    52315236
    52325237       t_wall_p%val(nzb_wall:nzt_wall,m) = t_wall%val(nzb_wall:nzt_wall,m) + dt_3d      &
     
    52345239                                             * surf%tt_wall_m(nzb_wall:nzt_wall,m) )
    52355240
    5236 
     5241!
    52375242!--    During spinup the tempeature inside window layers is not calculated to make larger timesteps possible
    52385243       IF ( .NOT. during_spinup )  THEN
    5239           win_absorp = - log( surf%transmissivity(m) ) / surf%zw_window(nzt_wall,m)
     5244          win_absorp = -log( surf%transmissivity(m) ) / surf%zw_window(nzt_wall,m)
    52405245!
    52415246!--       Prognostic equation for ground/roof window temperature t_window takes absorption
     
    52525257                              ) * surf%ddz_window_stag(nzb_wall,m)
    52535258
    5254 !
    5255 !--       If no indoor model is used, calculate the inner window heat flux according to
    5256 !--       heat diffusion (innermost window temperature as boundary conditions
    5257           IF ( .NOT. indoor_model ) THEN
    5258              surf%iwghf_eb_window(m) = - surf%lambda_h_window(nzt_wall-1,m)            &
    5259                                        * ( t_window%val(nzt_wall,m)                    &
    5260                                        - t_window%val(nzt_wall-1,m) )                  &
    5261                                        * surf%ddz_window(nzt_wall,m)
    5262           ENDIF
    5263 
     5259
     5260          IF ( .NOT. indoor_model )  THEN
     5261             surf%iwghf_eb_window(m) = surf%lambda_h_window(nzt_wall,m)                    &
     5262                               * ( t_window%val(nzt_wall+1,m) - t_window%val(nzt_wall,m) ) &
     5263                               * surf%ddz_window(nzt_wall+1,m)
     5264          ENDIF
    52645265
    52655266          DO  kw = nzb_wall+1, nzt_wall-1
    52665267             wintend(kw) = ( 1.0_wp / surf%rho_c_window(kw,m) )                        &
    5267                            * ( surf%lambda_h_window(kw,m)                              &
    5268                            * ( t_window%val(kw+1,m) - t_window%val(kw,m) )             &
    5269                            * surf%ddz_window(kw+1,m)                                   &
    5270                            - surf%lambda_h_window(kw-1,m)                              &
    5271                            * ( t_window%val(kw,m) - t_window%val(kw-1,m) )             &
    5272                            * surf%ddz_window(kw,m)                                     &
    5273                            + surf%rad_sw_in(m)                                         &
    5274                            * ( exp( -win_absorp * surf%zw_window(kw-1,m) )             &
    5275                            - exp(-win_absorp * surf%zw_window(kw,m) )                  &
    5276                              )                                                         &
    5277                              ) * surf%ddz_window_stag(kw,m)
     5268                              * ( surf%lambda_h_window(kw,m)                           &
     5269                                  * ( t_window%val(kw+1,m) - t_window%val(kw,m) )      &
     5270                                  * surf%ddz_window(kw+1,m)                            &
     5271                                  - surf%lambda_h_window(kw-1,m)                       &
     5272                                 * ( t_window%val(kw,m) - t_window%val(kw-1,m) )       &
     5273                                 * surf%ddz_window(kw,m)                               &
     5274                                 + surf%rad_sw_in(m)                                   &
     5275                                 * ( exp( -win_absorp * surf%zw_window(kw-1,m) )       &
     5276                                     - exp(-win_absorp * surf%zw_window(kw,m) )        &
     5277                                   )                                                   &
     5278                                 ) * surf%ddz_window_stag(kw,m)
    52785279
    52795280          ENDDO
    52805281
    52815282          wintend(nzt_wall) = ( 1.0_wp / surf%rho_c_window(nzt_wall,m) )                     &
    5282                                  * ( -surf%lambda_h_window(nzt_wall-1,m)                     &
     5283                              * ( -surf%lambda_h_window(nzt_wall-1,m)                        &
    52835284                                 * ( t_window%val(nzt_wall,m) - t_window%val(nzt_wall-1,m) ) &
    52845285                                 * surf%ddz_window(nzt_wall,m)                               &
     
    52865287                                 + surf%rad_sw_in(m)                                         &
    52875288                                 * ( exp( -win_absorp                                        &
    5288                                  * surf%zw_window(nzt_wall-1,m) )                            &
    5289                                  - exp( -win_absorp                                          &
    5290                                  * surf%zw_window(nzt_wall,m) )                              &
     5289                                     * surf%zw_window(nzt_wall-1,m) )                        &
     5290                                   - exp( -win_absorp                                        &
     5291                                     * surf%zw_window(nzt_wall,m) )                          &
    52915292                                   )                                                         &
    5292                                    ) * surf%ddz_window_stag(nzt_wall,m)
     5293                                ) * surf%ddz_window_stag(nzt_wall,m)
    52935294
    52945295          t_window_p%val(nzb_wall:nzt_wall,m) = t_window%val(nzb_wall:nzt_wall,m) + dt_3d       &
     
    52975298
    52985299       ENDIF
     5300
    52995301!
    53005302!--    Calculate t_wall tendencies for the next Runge-Kutta step
     
    76727674       2000000.0_wp,   &  !< parameter 26  - heat capacity 1st/2nd wall layer ground floor level
    76737675       1030000.0_wp,   &  !< parameter 27  - heat capacity 3rd wall layer ground floor level
    7674        900000.0_wp,   &  !< parameter 28  - heat capacity 4th wall layer ground floor level
    7675        0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
    7676        0.38_wp,         &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
    7677        0.04_wp,         &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
    7678        0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
    7679        0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
    7680        0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
    7681        0.7_wp,        &  !< parameter 35  - window transmissivity ground floor level
    7682        0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
    7683        0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
    7684        27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
    7685        5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
    7686        27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
    7687        0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
    7688        0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
    7689        0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
    7690        0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
    7691        20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
    7692        23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
    7693        20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
    7694        20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
    7695        23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
    7696        10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
    7697        1.0_wp,         &  !< parameter 51  - wall fraction ground plate
    7698        0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
    7699        0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
    7700        0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
    7701        0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
    7702        2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
    7703        1030000.0_wp,   &  !< parameter 57  - heat capacity 3rd wall layer ground plate
    7704        900000.0_wp,   &  !< parameter 58  - heat capacity 4th wall layer ground plate
    7705        0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
    7706        0.38_wp,         &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
    7707        0.04_wp,         &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
    7708        0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
    7709        0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
    7710        0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
    7711        0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
    7712        27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
    7713        0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
    7714        0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
    7715        0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
    7716        0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
    7717        1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
    7718        1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
    7719        1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
    7720        0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
    7721        0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
    7722        0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
    7723        27.0_wp,        &  !< parameter 77  - window albedo ground floor level
    7724        5.0_wp,         &  !< parameter 78  - green albedo ground floor level
    7725        0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
    7726        0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
    7727        0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
    7728        0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
    7729        1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
    7730        1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
    7731        1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
    7732        0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
    7733        0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
    7734        0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
    7735        1.0_wp,         &  !< parameter 89  - wall fraction roof
    7736        0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
    7737        0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
    7738        0.5_wp,        &  !< parameter 92  - 3rd wall layer thickness roof
    7739        0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
    7740        2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
    7741        1030000.0_wp,   &  !< parameter 95  - heat capacity 3rd wall layer roof
    7742        900000.0_wp,   &  !< parameter 96  - heat capacity 4th wall layer roof
    7743        0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
    7744        0.38_wp,         &  !< parameter 98  - thermal conductivity 3rd wall layer roof
    7745        0.04_wp,         &  !< parameter 99  - thermal conductivity 4th wall layer roof
    7746        0.93_wp,        &  !< parameter 100 - wall emissivity roof
    7747        27.0_wp,        &  !< parameter 101 - wall albedo roof
    7748        0.0_wp,         &  !< parameter 102 - window fraction roof
    7749        0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
    7750        0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
    7751        0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
    7752        0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
    7753        1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
    7754        1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
    7755        1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
    7756        0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
    7757        0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
    7758        0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
    7759        0.87_wp,        &  !< parameter 113 - window emissivity roof
    7760        0.7_wp,        &  !< parameter 114 - window transmissivity roof
    7761        27.0_wp,        &  !< parameter 115 - window albedo roof
    7762        0.86_wp,        &  !< parameter 116 - green emissivity roof
    7763        5.0_wp,         &  !< parameter 117 - green albedo roof
    7764        0.0_wp,         &  !< parameter 118 - green type roof
    7765        0.8_wp,         &  !< parameter 119 - shading factor
    7766        0.6_wp,        &  !< parameter 120 - g-value windows
    7767        3.0_wp,         &  !< parameter 121 - u-value windows
    7768        0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp, winter 0.1
    7769        2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp, winter 0.5
    7770        0.0_wp,         &  !< parameter 124 - heat recovery efficiency
    7771        2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
    7772        370000.0_wp,    &  !< parameter 126 - dynamic parameter innner heat storage
    7773        4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
    7774        100.0_wp,       &  !< parameter 128 - maximal heating capacity
    7775        0.0_wp,         &  !< parameter 129 - maximal cooling capacity
    7776        2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
    7777        6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
    7778        3.0_wp,         &  !< parameter 132 - storey height
    7779        0.2_wp,         &  !< parameter 133 - ceiling construction height
    7780        0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
    7781        1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
     7676          900000.0_wp,    &  !< parameter 28  - heat capacity 4th wall layer ground floor level
     7677          0.35_wp,        &  !< parameter 29  - thermal conductivity 1st/2nd wall layer ground floor level
     7678          0.38_wp,        &  !< parameter 30  - thermal conductivity 3rd wall layer ground floor level
     7679          0.04_wp,        &  !< parameter 31  - thermal conductivity 4th wall layer ground floor level
     7680          0.92_wp,        &  !< parameter 32  - wall emissivity ground floor level
     7681          0.11_wp,        &  !< parameter 33  - window emissivity ground floor level
     7682          0.86_wp,        &  !< parameter 34  - green emissivity ground floor level
     7683          0.7_wp,         &  !< parameter 35  - window transmissivity ground floor level
     7684          0.001_wp,       &  !< parameter 36  - z0 roughness ground floor level
     7685          0.0001_wp,      &  !< parameter 37  - z0h/z0q roughness heat/humidity
     7686          27.0_wp,        &  !< parameter 38  - wall albedo above ground floor level
     7687          5.0_wp,         &  !< parameter 39  - green albedo above ground floor level
     7688          27.0_wp,        &  !< parameter 40  - window albedo above ground floor level
     7689          0.005_wp,       &  !< parameter 41  - 1st wall layer thickness above ground floor level
     7690          0.01_wp,        &  !< parameter 42  - 2nd wall layer thickness above ground floor level
     7691          0.31_wp,        &  !< parameter 43  - 3rd wall layer thickness above ground floor level
     7692          0.43_wp,        &  !< parameter 44  - 4th wall layer thickness above ground floor level
     7693          20000.0_wp,     &  !< parameter 45  - heat capacity wall surface
     7694          23.0_wp,        &  !< parameter 46  - thermal conductivity of wall surface
     7695          20000.0_wp,     &  !< parameter 47  - heat capacity of window surface
     7696          20000.0_wp,     &  !< parameter 48  - heat capacity of green surface
     7697          23.0_wp,        &  !< parameter 49  - thermal conductivity of window surface
     7698          10.0_wp,        &  !< parameter 50  - thermal conductivty of green surface
     7699          1.0_wp,         &  !< parameter 51  - wall fraction ground plate
     7700          0.005_wp,       &  !< parameter 52  - 1st wall layer thickness ground plate
     7701          0.01_wp,        &  !< parameter 53  - 2nd wall layer thickness ground plate
     7702          0.31_wp,        &  !< parameter 54  - 3rd wall layer thickness ground plate
     7703          0.42_wp,        &  !< parameter 55  - 4th wall layer thickness ground plate
     7704          2000000.0_wp,   &  !< parameter 56  - heat capacity 1st/2nd wall layer ground plate
     7705          103000.0_wp,    &  !< parameter 57  - heat capacity 3rd wall layer ground plate
     7706          900000.0_wp,    &  !< parameter 58  - heat capacity 4th wall layer ground plate
     7707          0.35_wp,        &  !< parameter 59  - thermal conductivity 1st/2nd wall layer ground plate
     7708          0.38_wp,        &  !< parameter 60  - thermal conductivity 3rd wall layer ground plate
     7709          0.04_wp,        &  !< parameter 61  - thermal conductivity 4th wall layer ground plate
     7710          0.005_wp,       &  !< parameter 62  - 1st wall layer thickness ground floor level
     7711          0.01_wp,        &  !< parameter 63  - 2nd wall layer thickness ground floor level
     7712          0.31_wp,        &  !< parameter 64  - 3rd wall layer thickness ground floor level
     7713          0.43_wp,        &  !< parameter 65  - 4th wall layer thickness ground floor level
     7714          27.0_wp,        &  !< parameter 66  - wall albedo ground floor level
     7715          0.003_wp,       &  !< parameter 67  - 1st window layer thickness ground floor level
     7716          0.006_wp,       &  !< parameter 68  - 2nd window layer thickness ground floor level
     7717          0.012_wp,       &  !< parameter 69  - 3rd window layer thickness ground floor level
     7718          0.018_wp,       &  !< parameter 70  - 4th window layer thickness ground floor level
     7719          1736000.0_wp,   &  !< parameter 71  - heat capacity 1st/2nd window layer ground floor level
     7720          1736000.0_wp,   &  !< parameter 72  - heat capacity 3rd window layer ground floor level
     7721          1736000.0_wp,   &  !< parameter 73  - heat capacity 4th window layer ground floor level
     7722          0.11_wp,        &  !< parameter 74  - thermal conductivity 1st/2nd window layer ground floor level
     7723          0.11_wp,        &  !< parameter 75  - thermal conductivity 3rd window layer ground floor level
     7724          0.11_wp,        &  !< parameter 76  - thermal conductivity 4th window layer ground floor level
     7725          27.0_wp,        &  !< parameter 77  - window albedo ground floor level
     7726          5.0_wp,         &  !< parameter 78  - green albedo ground floor level
     7727          0.003_wp,       &  !< parameter 79  - 1st window layer thickness above ground floor level
     7728          0.006_wp,       &  !< parameter 80  - 2nd thickness window layer above ground floor level
     7729          0.012_wp,       &  !< parameter 81  - 3rd window layer thickness above ground floor level
     7730          0.018_wp,       &  !< parameter 82  - 4th window layer thickness above ground floor level
     7731          1736000.0_wp,   &  !< parameter 83  - heat capacity 1st/2nd window layer above ground floor level
     7732          1736000.0_wp,   &  !< parameter 84  - heat capacity 3rd window layer above ground floor level
     7733          1736000.0_wp,   &  !< parameter 85  - heat capacity 4th window layer above ground floor level
     7734          0.11_wp,        &  !< parameter 86  - thermal conductivity 1st/2nd window layer above ground floor level
     7735          0.11_wp,        &  !< parameter 87  - thermal conductivity 3rd window layer above ground floor level
     7736          0.11_wp,        &  !< parameter 88  - thermal conductivity 4th window layer above ground floor level
     7737          1.0_wp,         &  !< parameter 89  - wall fraction roof
     7738          0.005_wp,       &  !< parameter 90  - 1st wall layer thickness roof
     7739          0.01_wp,        &  !< parameter 91  - 2nd wall layer thickness roof
     7740          0.5_wp,         &  !< parameter 92  - 3rd wall layer thickness roof
     7741          0.79_wp,        &  !< parameter 93  - 4th wall layer thickness roof
     7742          2000000.0_wp,   &  !< parameter 94  - heat capacity 1st/2nd wall layer roof
     7743          103000.0_wp,    &  !< parameter 95  - heat capacity 3rd wall layer roof
     7744          900000.0_wp,    &  !< parameter 96  - heat capacity 4th wall layer roof
     7745          0.35_wp,        &  !< parameter 97  - thermal conductivity 1st/2nd wall layer roof
     7746          0.38_wp,        &  !< parameter 98  - thermal conductivity 3rd wall layer roof
     7747          0.04_wp,        &  !< parameter 99  - thermal conductivity 4th wall layer roof
     7748          0.93_wp,        &  !< parameter 100 - wall emissivity roof
     7749          27.0_wp,        &  !< parameter 101 - wall albedo roof
     7750          0.0_wp,         &  !< parameter 102 - window fraction roof
     7751          0.003_wp,       &  !< parameter 103 - window 1st layer thickness roof
     7752          0.006_wp,       &  !< parameter 104 - window 2nd layer thickness roof
     7753          0.012_wp,       &  !< parameter 105 - window 3rd layer thickness roof
     7754          0.018_wp,       &  !< parameter 106 - window 4th layer thickness roof
     7755          1736000.0_wp,   &  !< parameter 107 - heat capacity 1st/2nd window layer roof
     7756          1736000.0_wp,   &  !< parameter 108 - heat capacity 3rd window layer roof
     7757          1736000.0_wp,   &  !< parameter 109 - heat capacity 4th window layer roof
     7758          0.11_wp,        &  !< parameter 110 - thermal conductivity 1st/2nd window layer roof
     7759          0.11_wp,        &  !< parameter 111 - thermal conductivity 3rd window layer roof
     7760          0.11_wp,        &  !< parameter 112 - thermal conductivity 4th window layer roof
     7761          0.87_wp,        &  !< parameter 113 - window emissivity roof
     7762          0.7_wp,         &  !< parameter 114 - window transmissivity roof
     7763          27.0_wp,        &  !< parameter 115 - window albedo roof
     7764          0.86_wp,        &  !< parameter 116 - green emissivity roof
     7765          5.0_wp,         &  !< parameter 117 - green albedo roof
     7766          0.0_wp,         &  !< parameter 118 - green type roof
     7767          0.8_wp,         &  !< parameter 119 - shading factor
     7768          0.6_wp,         &  !< parameter 120 - g-value windows
     7769          3.0_wp,         &  !< parameter 121 - u-value windows
     7770          0.5_wp,         &  !< parameter 122 - basic airflow without occupancy of the room for - summer 0.5_wp for winter 0.1
     7771          2.0_wp,         &  !< parameter 123 - additional airflow dependent on occupancy of the room for - summer 2.0_wp
     7772                             !< for winter 0.5
     7773          0.0_wp,         &  !< parameter 124 - heat recovery efficiency
     7774          2.5_wp,         &  !< parameter 125 - dynamic parameter specific effective surface
     7775          165000.0_wp,    &  !< parameter 126 - dynamic parameter innner heatstorage
     7776          4.5_wp,         &  !< parameter 127 - ratio internal surface/floor area
     7777          100.0_wp,       &  !< parameter 128 - maximal heating capacity
     7778          0.0_wp,         &  !< parameter 129 - maximal cooling capacity
     7779          2.0_wp,         &  !< parameter 130 - additional internal heat gains dependent on occupancy of the room
     7780          6.0_wp,         &  !< parameter 131 - basic internal heat gains without occupancy of the room
     7781          3.0_wp,         &  !< parameter 132 - storey height
     7782          0.2_wp,         &  !< parameter 133 - ceiling construction height
     7783          0.1_wp,         &  !< parameter 134 - anthropogenic heat output for heating
     7784          1.333_wp        &  !< parameter 135 - anthropogenic heat output for cooling
    77827785                           /)
    77837786
Note: See TracChangeset for help on using the changeset viewer.