Analysis of the Influence Subjective Human Parameters in the Calculation of Thermal Comfort and Energy Consumption of Buildings

In the present work, we analyze the influence of the designer’s choice of values for the human metabolic index (met) and insulation by clothing (clo) that can be selected within the ISO 7730 for the calculation of the energy demand of buildings. To this aim, we first numerically modeled, using TRNSYS, two buildings in different countries and climatologies. Then, we consistently validated our simulations by predicting indoor temperatures and comparing them with measured data. After that, the energy demand of both buildings was obtained. Subsequently, the variability of the set-point temperature concerning the choice of clo and met, within limits prescribed in ISO 7730, was analyzed using a Monte Carlo method. This variability of the interior comfort conditions has been finally used in the numerical model previously validated, to calculate the changes in the energy demand of the two buildings. Therefore, this work demonstrated that the diversity of possibilities offered by ISO 7730 for the choice of clo and met results, depending on the values chosen by the designer, in significant differences in indoor comfort conditions, leading to non-negligible changes in the calculations of energy consumption, especially in the case of big buildings.


Introduction
The building envelope is the interface between the outdoor environment and the interior of buildings and its primary function is to act as a physical barrier to offer a comfortable place to develop different activities, in exchange for a certain demand for energy [1][2][3][4][5][6]. 40% of the primary energy consumption in the world and 17% of carbon dioxide emissions are caused by the building sector [7][8][9][10]. Thus, power consumption in buildings is one of the leading causes, among other adverse effects, of ozone layer depletion, global warming, and climate change [11,12]. Besides, residential buildings are often mentioned as one of the most profitable areas to decrease global power demand [13,14].
A rapid search on scientific databases shows that comfort in buildings is addressed by different areas, such as Engineering (≈ 39%) energy (≈14%), environmental (≈10%), computer (≈7.5%) and social sciences (≈7%). Most of the analysis emphasizes the fact that adequate levels of indoor comfort been developed that make dynamic simulations of building-occupant systems more tractable, such as visualization systems [52] agent-based models [53] or adaptive building simulations [54] among others. However, a significant problem when simulating the performance of a building is still the precise determination of the complete requirements of the models of the behavior of the occupants [55,56].
While most of the available simulators are focused on the performance of electrical and thermal networks, TRNSYS can similarly well be used to model other dynamic systems such as traffic flow, or even biological processes [57,58]. Studies using TRNSYS are reported in the literature as comparing cooling strategies for an office in different European climates [59], thermal comfort in an indoor swimming pool [60], or operating performance in cooling mode of a nearly zero energy building in China [61].
In this paper, we use TRNSYS for analyzing two different buildings: one in Spain and another one in Mexico. The numerical simulations have been carried out by multi-zone analysis, including technical/physical properties of the buildings, focusing on the cold (winter) and warm (summer) conditions in these countries, as there is more influence on thermal comfort. The novelty of our approach mostly lies in the combination of such BIM-based computer simulations with a sensitivity analysis (SA), to determine which input parameters have the maximum impact on the output [29].
The SA can be done by different popular methods such as Morris, regression, variance-based, or Monte Carlo, among others [62,63]. In our case, the SA based on Monte Carlo method [49] offers the possibility of prioritizing and fixating inputs when considering multiple outputs, in comparison with other methods [64,65]. Several studies can be found in the literature with SA performed on the effects of technical and physical parameters on the energy consumption of buildings and simulation programs [49,66,67]. However, to our knowledge, less attention has been paid to the effect of subjective human parameters [51] and thus, this aspect will be our main focus for the present work. Therefore, as a final step, we will show the significant impact of the designer's subjective choices on the calculation of the energy demand of buildings, depending on the election of the amounts of activity or clothing, based on the environmental ergonomic standard ISO 7730. We consider that our SA of the above-mentioned parameters using the Monte Carlo method will help the designers to define the most appropriate parameters to be considered, depending on the aim and type of the numerical modeling performed.

Case Study
In the present work, we analyze two non-residential buildings localized in different continents (America and Europe), with particular characteristics in their structure, use, and surrounding climate. The first building is a recent one in Cambre, in north-western Spain (latitude 43 • 18'22.55" N, longitude 8 • 17'36.13" W) located at a height of 99 m above sea level. The second is a building situated in San Nicolás de los Garza, in north-eastern Mexico (latitude 27 • 32'02" N, longitude 99 • 58'33"), at a height of 239 m above sea level.

Climate Conditions
The Map of K'oppen-Geiger Map [68] define different regions classifications in function of the central values of temperature and humidity, the Spanish location is classified in the Cfb region (C-temperate, f-without dry season, b-warm summer), where Oceanic Climate has a soft summer with a minimum rainfall and does not exceed of 22 • C on average in the warmest month, and at least four months averaging above 10 • C. No significant precipitation variation between seasons.
With respect to Mexico, according to the same classification, San Nicolás de los Garza can be considered to be BSh (BS-semi-arid, h-hot). A semi-arid (or steppe) region receives precipitation below potential evapotranspiration, but not extremely. Semi-arid climates tend to support short or scrubby vegetation, with large zones usually dominated by either grasses or shrubs. Also, the average annual temperature is higher than 18 • C. From the databases available in TRANSYS, the web platform Equs [69], and Meteonorm [70], we obtain the values of different climatic variables such as air temperature, relative humidity, direct irradiance (solar) and diffuse horizontal irradiance.

Building Structure And Use
The base of the first building is of rectangular shape with an area of 218.14 m 2 (see Figure 1) and its main wall faces southwest (216 • N). The building is distributed on two floors (defined as ground and first), being the lower level used for storage and the upper zone for offices, meeting rooms, and kitchen, among others uses. This building is the headquarters of a private company specialized in air-conditioning systems. The working hours, considered from Monday to Friday throughout the year, are from 08:30 a.m. to 06:00 p.m. For the second case, the plant of the building has an "L" shape, covering an area of 1214 m 2 in each of its three heights that we define as: ground floor (GF), first floor (FF) and second floor (SF). The GF corresponds to a parking lot without vertical enclosures and a reception area, while the other two floors are classrooms, meeting rooms, laboratories, offices, etc. as is shown in Figure 2. The building hosts a research center on physics and mathematics (CICFIM), which is located behind of the School of Physics and Mathematics of the Universidad Autónoma de Nuevo León (Mexico). Its principal access is on the facade facing east, next to the campus stadium, to the west. This building starts the activities at 7:00 a.m., ending at 9:00 p.m. from Monday to Friday. Saturday's schedule is from 10:00 a.m. to 2:00 p.m. The building in Spain is made with three different types of envelopes: the vertical facing walls, the floor, and the roof. In Table 1 we indicate that their respective thicknesses (t) are 46.5 cm, 30.4 cm and 185 cm. The building in Mexico Table 2 has four types of envelopes: internal (9 cm) and external (35 cm) walls, as well as 10 cm-thick divisions. In this case, the same materials were used in the floor and the roof, with a total thickness of 40 cm.
In both Tables 1 and 2, the physical characteristics of each layer of material are specified, such as thermal transmittance (U − value), thermal conductivity (κ), heat capacity (C p ), density (ρ) and resistance (R) in specific cases. For both buildings, the properties of the windows are described in Table A1. In the first building (Spain) there are two types of zones, depending on whether there is air conditioning or not (see Table A2) in Appendix A. The second edifice (Mexico) is distributed in three partitions, because of the different uses and temperatures inside. Among its main characteristics are frame percentage, thermal transmittance (U) and solar gain (G − Value).
In Spain, the HVAC system consists on an inverter heat pump (2 × 1) and a split air-air unit, with a total refrigeration capacity of 17.2 kW, a heat power of 19.5 kW and an installed power of 6.5 kW. The diffusion of air is carried out by ducts and impulse/return grids. For the building in Mexico which a conventional technology, based on air-air heat pumps with 26 exterior units and indoor cassette units and a total cooling capacity of 349.4 kW. The heat output is 447.7 kW and the installed power is 166.4 kW.

Methodology
The main objective of this research is to show the consequences on the results of BIM-based simulations of subjective decisions made by the designers when using the environmental ergonomic standard ISO 7730 to calculate thermal comfort ranges and its impact on the energy demand of buildings, as well as its economic repercussion. The thermal sensation experienced by a human being is mainly due to the overall thermal balance of the body, which depends on the physical activity (met) and clothing (clo) of the subject, as well as by the environmental parameters: air temperature, radiant temperature, air speed, and air humidity [71]. We must stress that the definition of the parameters of human clothing and activity is done subjectively by the designer and therefore we refer to these quantities as Subjective Human Parameters (SHPs). The use of biometric data from wearable devices has been recently proposed as an interesting option for metabolic rate estimation in thermal comfort analysis [72].
The PMV and PPD indices, according to Fanger's method [71] are commonly used to express warm and cool discomfort for the body as a whole, although thermal dissatisfaction may also be caused by unwanted heating or cooling of one particular part of the body (local discomfort). Therefore, when the previous factors are known, the thermal sensation for the occupants can be predicted, using a seven-point scale ranging from cold (−3) to hot (3), being 0 the value for "neutral sensation" [73]. The PMV may be computed from tables or may be measured by an instrument [74]. The methodology proposed in this paper consists of obtaining first the thermal performances of the buildings and then performing the SA for the SHPs used in the thermal comfort calculation, with the purpose of identifying which quantity has more influence on the BIM-based simulation outputs. Figure 3 shows a scheme of the procedure applied, consisting of the following steps, for both cases studied: 1. Real-time tracking of the main variables using Equs web platform. 2. TRNSYS modeling, getting the indoor temperature. 3. Verification of the thermal simulation with measured temperature from Equs database. 4. Mathematical estimation of PMV and PPD indices, according to ISO 7730. 5. SA (Monte Carlo) to determine which parameter (activity and clothing) values yield significant variation on PMV and PPD outputs. 6. New set-point temperatures are established, depending on the changes in the thermal categories in the previous step. 7. TRNSYS recalculation of the energy demand for each set-point temperature from the Monte Carlo method. 8. Evaluation of the variations in the energy demand for the different options.

Monitoring by Web Platform
A web platform for energy management of buildings named "Equs", was developed with the purpose of control, maintenance and monitoring of energy consumptions and environmental variables as temperature, humidity and CO 2 levels. This tool also includes the possibility of remotely controlling any device by SCADA (Supervisory Control and Data Acquisition) system, ensuring the quality 6 of the statistical results. Equs collects in its database the information sent by the different sensors, allowing us to know the variation of the physical variables remotely, as access to these data can be obtained with any computer or other devices such as smartphones, which are connected to the Internet, making it very simple the treatment of the data, which can be done in the same device accessing to the values, in real time [69].

Trnsys Modeling of The Buildings
Through dynamic simulation tools, such as TRNSYS, we can define and evaluate the comfort and energy demand of the two buildings accurately; in both cases, we work with multi-zone building models. Notice that TRNSYS is integrated into two parts: TRNBuild, which reads and processes an input file, solving and plotting the system; and TRNSYS Studio, which is a library of components that can be used to model different systems. An advantage of using TRNSYS is the possibility to convert complex problems as multi-zone building models into smaller components, to allow the designer to establish the connections between the building and numerous other subsystems/components in the simulation environment [75].
To verify the energy demands of the structures, after analyzing the sensitivity with the Monte Carlo method of the set-point temperatures of the air-conditioning systems with the met and clo parameters, it has opted for its simulation in the program TRNSYS. To simulate the building, both in Spain and in Mexico, all the necessary inputs, dimensions, location, construction characteristics, etc. have been introduced into the program. This simulation will be verified later comparing the results that the TRNSYS throws in certain indoor temperature conditions, with the actual measurements of both the environmental conditions and the interior temperatures. Once the model is validated, we can calculate the variation of the energy demand in both buildings as indicated by the Monte Carlo SA to the parameters studied in ISO 7730.
Entering the buildings models in TRNSYS was done in two steps: first, in TRNBuild we implemented the different zones and all its elements as location, orientations, walls, layers of materials, gains, ventilation, infiltration, schedules, and comfort in TRNBuild. The second step corresponds to include in TRNSYS models as weather, inputs such as the building previously described and outputs that will show through graphics the thermal performance and energy demands of the buildings. Thus, first, we obtained the database of the buildings, which includes information such as orientations of the facades, materials for the walls and windows, the energy demand of the installed equipment, the working hours, the occupation, etc. The process of the definition of the buildings in TRNBuild includes the following information: • Location and orientations. This information was described in Section 2.2 • Structure of the building: the materials detailed in Tables 1, 2 and A1 were assigned from the software database.

•
The zones described in Table A2, were created as WAC and AC for the first case, GF-FF, SF, and SL for the second case in the TRNBuild Manager.

•
In Table A3 we show the different types of schedules created for each building, we set a value of 1 in the periods where there is activity in the building, in any other case, including the weekend, the value is 0.

•
In Table A4 we summarize the infiltrations for summer and winter in both cases. This information was calculated based on ASHRAE Standard 55.

•
For calculating the ventilation we use standard EN 13779, ASHRAE 62 R. It was estimated a value for each zone which is shown in Table A2.

•
The internal gains of the buildings are shown in Tables 3 and 4.

•
The values of clothing factor, metabolic rate, relative air velocity were introduced in TRNSYS based on the ISO Standard 7730.
For estimation of comfort, it is important to note that thermal sensations are different between persons sharing the same environment. Perceptions of human beings about thermal comfort are very different because even people who stay in similar places have different opinions, because there exist many elements that impact the perception of human beings.
Thermal comfort is the main goal and has a strong influence in the thermal behavior of a building, for this reason, the values of clothing factor, metabolic rate, relative air velocity were introduced in TRNSYS, based on the ISO Standard 7730. In Table 5, we show the numbers for both cases in summer and winter. We decided to set the same values of clothing and activity, because the characteristics of the occupants are very similar. It should be noted that these values are subjective decisions on the part of the designer. The analysis of the impact of these decisions is one of the main goals of this work.

Simulation of Building Models
In this section, it is presented the dynamic thermal simulation model in TRNSYS Studio of the building previously described in TRNBuild, for this purpose we use multi-zone building model "TYPE 56" that can be connected to a large number of other components, including weather data, HVAC systems, occupancy schedules, controllers, output functions, etc. Other models used were Type 62, Type 109-TMY2, and Type 65.
Under ISO Standard 7730, the evaluation of the levels of comfort was calculated using a mathematical MATLAB software package for the interactive numerical study of dynamical systems.
The change of PMV and PPD in base of its parameters such as metabolic rate, clothing insulation, and air temperature, is easier of visualizing in MATLAB than in TRNSYS because you can modify its parameters independently and realize studies more specific such SA.

Monte Carlo Sensitivity Analysis
The principal aim of the SA using a Monte Carlo code is to obtain the variation of the PPD and PMV concerning the possible choices of the clo and met that are indicated in the annexes of the 7730 for thermal comfort calculations. In this way, by changing the values of clo and met according to what is reported in the tables of ISO 7730, we try to find out how the environmental category varies and, consequently, the set-point temperature. The SA is extensively suitable for associating input parameters with the overall building performance (for example, the temperature or demand for heating or cooling) [64]. The individual contribution from specific input can be determined by SA to prepare for future optimization of energy, climate, and economic performance [76,77]. SA methods can typically be classified to the local SA where input parameters are varied one at a time and the global SA where all inputs are changed simultaneously.
Sampling-based methods for uncertainty and SA have become very popular. The sampling-based Monte Carlo method is a procedure that performs a model repeated times with random samples generated from input distributions. It provides approximate solutions to both uncertainty and sensitivity analyses by making statistical modeling investigations [78]. The method can manage complex black box models irrespective of the linearity and continuity, and generate a probability distribution for each output depending on input distribution types [62].
The analysis under consideration can be represented by a function y(x) where: and designates the outputs and inputs of the analysis, respectively. Successively, sensitivity in x results in a corresponding sensitivity in y(x). Furthermore, the distributions must be well-defined to characterize the uncertainty associated with the elements of x, where D i is the distribution associated with x i for i = 1, 2 . . . n X . The most common way for applying SA in building performance analysis consists on the following steps: determine input variations, create building energy models, run energy models, collect simulation results, run SA and presentation of SA results [64]. Determining the probability distributions of input parameters is the first step in any SA, which depends on the election of sensitivity method. In the present study, Monte Carlo SA is used, considering N input parameters defining the system: Suppose that the parameter P j can change within a range given by ∆P j . To analyze the influence of the parameter P j in the objective O i , we proceed first generating a uniform random number distribution r ∈ [0, 1], thus: gives us the first random value P 1 j . With the value of N randomly generated parameters P 1 j (j = 1, 2, . . . , N), we obtain a value of O 1 i , which is the first evaluation of the objective O i . We repeat this procedure up to the desired number of random evaluations. We can randomly change each one of the parameters or we can apply random variations to one of the parameters keeping constant the value of the rest.
For this work, the objectives O i are the PMV and PPD indices, which are a function of different parameters mentioned in the previous sections. To analyze the influence on the thermal performance of the buildings, the P j parameters are CLO and MET. Depending on the admissible ranges for PMV and PPD, three kinds of comfort zones or categories of thermal requirements are defined by UNE-EN ISO 7730 as: category I (or class A) (PPD < 6 %, i.e., −0.2 < PMV < 0.2), category II (or class B) (PPD < 10%, i.e., −0.5 < PMV < 0.5) and category III (or class C) (PPD < 15%, i.e., −0.7 < PMV < 0.7). The ranges of recommended air temperatures for different types of buildings depending on the previous categories are shown in Table 6 [26]. Based on the results of the SA, we can establish the set-point temperature for the cooling and heating seasons. We use the SA based on the Monte Carlo method, randomly generating a uniform distribution for clo from 0 to 2 and met from 0.8 to 4, since in this range the PMV index lies between −2 and +2, taking into account that this is the range suggested by the standard. A uniform distribution corresponds to the case of a random variable that can only have values between two extremes a and b, so that all the intervals of the same length (within (a, b)) have the equal probability.
The above definition shows that the density function must be the same value for all points within the interval (a, b) (and zero outside the range), i.e., The results of the SA on the comfort index (PMV and PPD) for the clothing and metabolic rate are shown in Figures 4 and 5, respectively. The thermal category changes according to PMV and PPD values, which in turn depend on the subjective decisions of clo and met. Figure 4 shows that for category A, the clo values vary from 0.46 to 0.66, for B of 0.34 to 0.46 clo and 0.66 to 0.87 clo. Category C is between 0.27 to 0.34 clo and 0.87 to 1.04 clo. For the case of met and PPD, in Figure 5 we show that category A lies between met values from 0.83 to 1.04. For the case of category B, it remains from 0.8 to 0.83 met and 1.04 to 1.16 clo. Finally, category C corresponds to the range from 1.16 to 1.31 met.
An attractive alternative in this type of SA, it is to use a normal distribution, which is one of the theoretical distributions best studied in statistical texts and more in practice, also called Gaussian distribution. Its importance is mainly due to the frequency with the variables other than natural and everyday phenomena, approximately, this distribution. The normal distribution is determined by two parameters, their mean and their standard deviation, usually denoted by µ y σ [79]. With this notation, the density of the distribution Gaussian is given by the equation: For the case of a normal distribution, we have randomly generated values for clo with µ = 1.0, setting the metabolic factor to 1.2 met and, varying the standard deviation σ clo = 0.025, 0.05, 0.075, 0.1 and 0.125 . For these values, we obtain the distribution of the PMV and PPD values and, by adjusting to a normal distribution, we obtain the values of µ and σ for the PMV and PPD. Similarly, for the case of the met (metabolic factor) a µ = 1.2 has been established, setting the clothing factor to 1.0 clo and varying the standard deviation σ met = 0.025, 0.05, 0.075, 0.1 and 0.125.
That is to say, we see easily as variations in the values of clo and met, they are transferred through the model to variations of PMV and PPD, and therefore, we can see the sensitivity of the model to these variations of the clo and met parameters.
That is, we have obtained the relation of the variations of the clothing (clo) and activity (met) with the widths (σ) of the PMV and PPD, with the objective of observe which of the two parameters that the designer enters in a subjective manner, produces a more significant impact on the output of the model.
In the Figure 6a, the relationship between the widths of the distributions of clo and met and the value of σ PMV is presented. For example, for a width of 5 % with respect to the central value of clo µ = 1 and (σ clo = 0.10), a value of σ PMV of 0.150 is obtained and for a width of 5 % with respect to the central value µ = 1.2 of met (σ met = 0.10), a value of σ PMV of 0.184 is obtained. Therefore, the impact of the metabolic factor (met) is higher than that of the clothing factor (clo) for the PMV case.
In the same way, for each percentage of variation in the values of clo and met we have a value of σ PPD , which is shown in the Figure 6b, for example: for a width of 12.5 % the value of σ PPD is 1.990 for a variation in clo and 2.594 if the parameter met varies. Also, the values predict that the impact of the metabolic factor (met) is greater than the clothing factor (clo) in the PPD.

Results
In this section, we discuss the results of the simulations of both buildings. First, we will show the verification of the monitoring using the data from Equs and the TRNSYS simulation of each building. Then, after a second simulation, we will get the energy demand using the calculated results for the category change, depending on the estimation of the PMV and PPD values as well as the sensitivity analysis of the CLO and MET variables. Finally, we analyze the economic impact of energy demand savings.
As noted before, this work joins BIM-based energy simulations and SA of two non-residential buildings in two seasons (summer and winter) under the international standard ISO 7730. This norm has been studied in previous papers from different perspectives: the effect of personal factors in comfort assessment [80], energy simulation [27], thermal comfort in outdoor urban spaces [81] or the validity of ISO-PMV for predicting comfort votes in everyday thermal environments [82], among others. However, to our knowledge, the effect on the energy demand of the diversity of possibilities offered by ISO 7730 for the choice of clo and met results, depending on the values chosen by the designer, was not previously calculated.

Real-Time Monitoring and Validation of The Models
Real-time tracking made online with Equs demonstrates the reliability of the measurements that validate the computational simulations. In Figure 7 we show the temperatures measured from 1 January to 31 October 2015 of the air-conditioning zone of the building in Spain, located 43 • 18' north and 8 • 17' west. In the first month (January, winter season) the lowest temperature of the period is 12.53 • C registered on Saturday 17. The highest temperature (30.48 • C), was achieved on Monday, 15 June (summer). The mean temperature was 22.43 • C. For the study of the outside temperature, we observe that the cold and hot temperatures are between 2 • C and 41 • C, approximately. The mean temperature outside the building was 15.49 • C.  The comparison between the real data and the simulations is required to have reliable simulations of the buildings. Currently, a simulation of a building's energy consumption is usually considered "calibrated" if the criteria set out by ASHRAE Guideline 14 is met [83]. Accordingly, the model can be considered "calibrated," according to current international criteria to accept BES models, provided that considerable accordance is reached between the measured and simulated data [84]. We must stress that there is no unique model that meets the previous criteria and therefore, several alternative models of the same building can be considered to be "calibrated". Moreover, it should be noted that our calibration is very accurate for the simulated environment (e.g., temperature profiles). As an example, the parallelism among experimental data and the modeling for internal air temperature for the month of June from the hour 3624 to 4344, is shown in Figure 8a. For the period of June, the maximum difference between the simulated and measured temperatures is 1.  In Figure 8b we show the temperature of a single day of June (24 h from hour 3914 to 3938 of the simulation). We can observe in more detail the difference between the calculated and the measured temperature for the period, where the minimum variation is of 0.04 • C (obtained in working hours) and the mean difference is of 0.35 • C. The maximum variation is 0.72 • C. These results were obtained within the situation described in the previous sections. In particular, during the occupancy periods, the metabolic rate (met) of the occupiers was selected as to be "seated light work", which corresponds to 70 W/m 2 or 1.2 met. The clothing factor (clo) was established to 0.7 for the summer period and 1.2 for winter. Notice that a clothing with factor 1.2 clo corresponds to: underwear, shirt, pants, thermal jacket, socks, shoes. Summer clothing of 0.7 clo corresponds to underwear, shirt, pants, socks, shoes [26]. Table 5 shows the input that was used for the simulation for the base case scenario. In Figure 9 we can observe the energy demand in the initial configuration, where the set point was established in 23.8 • C.

Evaluation of the Variations in the Energy Demand
The main goal of the application of a SA in this work is to visualize in a simpler and more detailed way than in TRNSYS, the impact of subjective decisions made by the designer on the values of insulation of clothes and activity for thermal comfort, which are defined in the TRNSYS simulation program. Based on ISO 7730, the insulation of clothes (clo) can be estimated directly from the combination of clothes or by adding the partial insulation values for each clothes. For example, typical value of clothes insulation for summer (cooling) and winter (heating) are 0.5 and 1 clo, respectively. Therefore, these quantities can be determinant for the different types of clothes that the occupants use and thus, the designers have a greater margin of error when determining by subjective decisions the values of clothing parameters.
Once we have made the variation of clo and met, we obtain variations in the PMV and PPD indices, these values of the indices allow us to define the thermal category based on the definition of the ISO 7730 standard. This category establishes the range of recommended temperatures for these input parameters of comfort. In this permissible temperature range, we can select a value of set-point temperature, to introduce it to our validated BIM model and thus obtain the new energy demand.
For instance: summer light clothing insulation of 0.5 clo includes items such as underpants, short-sleeved shirt, light pants, thin socks, or shoes. In our case, we choose for both cases of study a mean value of 0.7 clo that includes underwear, shirt, pants, socks, shoes, which are also considered light clothing for the same season in this type of buildings. Analogously, the estimation of the value of activity or metabolic rate (met), which represents a heat production depending on the activity level, also depends on the subjective decision made by the designer. According to ISO 7730 for an activity corresponding to this type of buildings such as office, school, and laboratory, the activity parameter can vary from 1.0 to 1.6 met. For this work we assumed a value of 1.2 met for the energy simulations on TRNSYS.
Thus, it is of fundamental importance to determine the impact of the mentioned subjective choices of clothing and metabolic factor, on the values of the PMV and PPD indexes, which show the level of thermal comfort and local thermal comfort criteria, reporting the environmental requirements that are considered acceptable for global thermal comfort, including the local discomfort.
We must notice that a variation in the clo values of summer clothing, from 0.3 to 0.8, causes a change in category from A to C. In Figure 10, the impact on the calculation of energy demand is clearly shown that has the change of environmental category.
According to the SA results and the results of Table 6, we established in TRNSYS the new values for clo and set-point temperature. The first configuration was set for a value of clo of 0.7 with a set-point temperature of 23.5 • C, which belongs to category A. For category B, we used a temperature of 24 • C and clo of 0.5. This configuration was initially defined for the validation of the TRNSYS simulation previously shown. Finally, a set-point temperature of 25 • C for category C and clo of 0.3 was considered. As an example, for the month of June, from Figure 10 for clo values of 0.3, 0.5 and 0.7, we obtain a total energy demand of 347.19 kW, 565.05 kW, and 885.2 kW, respectively.
In Table 7 it can be observed that changes in the energy demand due to variations in the designer's criteria in climates such as Mexico, can reach up to 22% depending on the choice of clo and up to 23% in energy demand, depending on the choice of the activity rate made by the designer. As it can be seen after the analysis of the Monte Carlo method, the BIM calculation is very sensitive to both the appropriate choice of the met and the clo choice. In the case of the clo we have studied variations of daily clothing in summer with clo 0.3 (panties, t-shirt, shorts, thin socks, sandals), clo 0.5 (underpants, t-shirt, light pants, thin socks, shoes), clo 0.7 (underwear, shirt, pants, socks, shoes).
All previous values of clo are assignable to the attire of the occupants of both buildings studied both in the case of Spain and in Mexico, in summer. However, after analyzing it by the Monte Carlo method and calculating the environmental category by the TRNSYS program, we see that in both buildings for a room 0.3 the comfort in the buildings studied is in the category A ranges; for a clo of 0.5 in category B and if the clo is 0.3, in category C (Figure 4). This variation in the category, due to the uncertainty in the choice of the clo on the part of the designer, can have, consequently, alterations in the energy demand.  As it can be appreciated from the values of Table 7 the change is more than 20%, for the building in Mexico. For the building in Spain, the variation in the calculation of the energy demand, despite having a higher sensitivity, is much lower (only 6%). This discrepancy is due to the different volumes of the two buildings. The higher the volume, the greater the energy demand and therefore the variations of the environmental categories that determine the thermal comfort inside have a much greater impact.
With respect to the met, according to ISO 7730, the cases of a technical office (main activity of the building of Spain) and a University (main activity of the building of Mexico) correspond to met 1 (sitting rest), and met 1.2 (office activity, school, home, laboratory). As it can be appreciated from the values of Table 7, the higher variations in the calculation of energy demand in buildings have been again in the case of Mexico (more than 20%). In the case of the building in Spain the oscillation is only of 6.6%.

Conclusions
In the present work, we have analyzed, within the mandatory standard ISO 7730, the influence of the designer's choice of values for human metabolic index and insulation by clothing that can be selected for the estimation of the energy demand of structures. We have demonstrated that the diversity of possibilities offered by ISO 7730 for the choice of clo and met results in significant differences in indoor comfort conditions, leading to non-negligible changes in the calculations of energy consumption.
It has been clearly outlined that the higher the cooling demands, the greater the sensitivity, generating greater differences according to the designer's choice. Since ISO 7730 is an international norm, given the great variety of existing climates, it is vital to take into account this peculiarity of the system at the time of its application and in the choice of the designer, who must be familiar with the procedure. The international norm ISO 7730 is a reference in most of the countries of the world and nevertheless requires of an experienced designer for an appropriate application. Therefore, we have shown in the present work the considerable impact of the designer's subjective choices on the calculation of the energy demand of buildings, as well as the extra costs cumulated, depending on the election of the amounts of activity or clothing, based on the environmental ergonomic standard ISO 7730. We consider that our SA of the parameters mentioned above using the Monte Carlo method will help to define the most appropriate metrics to be admitted, depending on the aim and type of the modeling performed.
We have calculated that the effect of to the different possibilities of election within the ISO 7730 can reach a quarter of the energy demand, showing that this standard is susceptible to the designer's choice of the clo and met. As the ISO 7730 is the only standard which applies to the building regulations in many parts of the world, we conclude that the norm could be revisited from the perspective of energy efficiency, maybe including recommendations from other standards such as EN 15251-2007 or ISO 17772-1-2017, designed explicitly for providing environmental input parameters for design and assessment of energy performance of buildings.