Correct bug in shadding object creation when no archetype available

This reverts commit 9adacadc2a.

hub/catalog_factories/data_models/energy_systems/non_pv_generation_system.py

@@ -25,9 +25,11 @@ class NonPvGenerationSystem(GenerationSystem):
                maximum_cooling_supply_temperature=None, minimum_cooling_supply_temperature=None, heat_output_curve=None,
                heat_fuel_consumption_curve=None, heat_efficiency_curve=None, cooling_output_curve=None,
                cooling_fuel_consumption_curve=None, cooling_efficiency_curve=None,
-               distribution_systems=None, energy_storage_systems=None, dual_supply_capability=False):
-    super().__init__(system_id=system_id, name=name, model_name=model_name, manufacturer=manufacturer, fuel_type=fuel_type,
-                     distribution_systems=distribution_systems, energy_storage_systems=energy_storage_systems)
+               distribution_systems=None, energy_storage_systems=None, domestic_hot_water=False,
+               reversible=None, simultaneous_heat_cold=None):
+    super().__init__(system_id=system_id, name=name, model_name=model_name, manufacturer=manufacturer,
+                     fuel_type=fuel_type, distribution_systems=distribution_systems,
+                     energy_storage_systems=energy_storage_systems)
     self._system_type = system_type
     self._nominal_heat_output = nominal_heat_output
     self._maximum_heat_output = maximum_heat_output
@@ -53,7 +55,9 @@ class NonPvGenerationSystem(GenerationSystem):
     self._cooling_output_curve = cooling_output_curve
     self._cooling_fuel_consumption_curve = cooling_fuel_consumption_curve
     self._cooling_efficiency_curve = cooling_efficiency_curve
-    self._dual_supply_capability = dual_supply_capability
+    self._domestic_hot_water = domestic_hot_water
+    self._reversible = reversible
+    self._simultaneous_heat_cold = simultaneous_heat_cold
 
   @property
   def system_type(self):
@@ -256,12 +260,28 @@ class NonPvGenerationSystem(GenerationSystem):
     return self._cooling_efficiency_curve
 
   @property
-  def dual_supply_capability(self):
+  def domestic_hot_water(self):
     """
-    Get dual supply capability
+    Get the ability to produce domestic hot water
     :return: bool
     """
-    return self._dual_supply_capability
+    return self._domestic_hot_water
+
+  @property
+  def reversibility(self):
+    """
+    Get the ability to produce heating and cooling
+    :return: bool
+    """
+    return self._reversible
+
+  @property
+  def simultaneous_heat_cold(self):
+    """
+    Get the ability to produce heating and cooling at the same time
+    :return: bool
+    """
+    return self._simultaneous_heat_cold
 
   def to_dictionary(self):
     """Class content to dictionary"""
@@ -269,6 +289,18 @@ class NonPvGenerationSystem(GenerationSystem):
                              self.distribution_systems] if self.distribution_systems is not None else None
     _energy_storage_systems = [_energy_storage_system.to_dictionary() for _energy_storage_system in
                                self.energy_storage_systems] if self.energy_storage_systems is not None else None
+    _heat_output_curve = self.heat_output_curve.to_dictionary() if (
+        self.heat_output_curve is not None) else None
+    _heat_fuel_consumption_curve = self.heat_fuel_consumption_curve.to_dictionary() if (
+        self.heat_fuel_consumption_curve is not None) else None
+    _heat_efficiency_curve = self.heat_efficiency_curve.to_dictionary() if (
+        self.heat_efficiency_curve is not None) else None
+    _cooling_output_curve = self.cooling_output_curve.to_dictionary() if (
+        self.cooling_output_curve is not None) else None
+    _cooling_fuel_consumption_curve = self.cooling_fuel_consumption_curve.to_dictionary() if (
+        self.cooling_fuel_consumption_curve is not None) else None
+    _cooling_efficiency_curve = self.cooling_efficiency_curve.to_dictionary() if (
+        self.cooling_efficiency_curve is not None) else None
 
     content = {
       'Energy Generation component':
@@ -298,13 +330,15 @@ class NonPvGenerationSystem(GenerationSystem):
         'minimum cooling supply temperature [Celsius]': self.minimum_cooling_supply_temperature,
         'heat output curve': self.heat_output_curve,
         'heat fuel consumption curve': self.heat_fuel_consumption_curve,
-        'heat efficiency curve': self.heat_efficiency_curve,
+        'heat efficiency curve': _heat_efficiency_curve,
         'cooling output curve': self.cooling_output_curve,
         'cooling fuel consumption curve': self.cooling_fuel_consumption_curve,
         'cooling efficiency curve': self.cooling_efficiency_curve,
         'distribution systems connected': _distribution_systems,
         'storage systems connected': _energy_storage_systems,
-        'dual supply capability': self.dual_supply_capability
+        'domestic hot water production capability': self.domestic_hot_water,
+        'reversible cycle': self.reversibility,
+        'simultaneous heat and cooling production': self.simultaneous_heat_cold
       }
     }
     return content

hub/catalog_factories/data_models/energy_systems/performance_curves.py

@@ -24,13 +24,13 @@ class PerformanceCurves:
   def curve_type(self):
     """
     The type of the fit function from the following
-    Linear =>>> y = a*x + b
+    Linear =>>> y = a + b*x
     Exponential =>>> y = a*(b**x)
-    Polynomial =>>> y = a*(x**2) + b*x + c
+    Second degree polynomial =>>> y = a + b*x + c*(x**2)
     Power =>>> y = a*(x**b)
-    Second degree multivariable =>>> y = a*(x**2) + b*x + c*x*z + d*z + e*(z**2) + f
+    Bi-Quadratic =>>> y = a + b*x + c*(x**2) + d*z + e*(z**2) + f*x*z
 
-    Get the type of function from ['linear', 'exponential', 'polynomial', 'power', 'second degree multivariable']
+    Get the type of function from ['linear', 'exponential', 'second degree polynomial', 'power', 'bi-quadratic']
     :return: string
     """
     return self._curve_type

hub/catalog_factories/data_models/energy_systems/thermal_storage_system.py

@@ -17,7 +17,7 @@ class ThermalStorageSystem(EnergyStorageSystem):
 
   def __init__(self, storage_id, type_energy_stored=None, model_name=None, manufacturer=None, storage_type=None,
                nominal_capacity=None, losses_ratio=None, volume=None, height=None, layers=None,
-               maximum_operating_temperature=None, storage_medium=None):
+               maximum_operating_temperature=None, storage_medium=None, heating_coil_capacity=None):
 
     super().__init__(storage_id, model_name, manufacturer, nominal_capacity, losses_ratio)
     self._type_energy_stored = type_energy_stored
@@ -27,6 +27,7 @@ class ThermalStorageSystem(EnergyStorageSystem):
     self._layers = layers
     self._maximum_operating_temperature = maximum_operating_temperature
     self._storage_medium = storage_medium
+    self._heating_coil_capacity = heating_coil_capacity
 
   @property
   def type_energy_stored(self):
@@ -84,6 +85,14 @@ class ThermalStorageSystem(EnergyStorageSystem):
     """
     return self._storage_medium
 
+  @property
+  def heating_coil_capacity(self):
+    """
+    Get heating coil capacity in Watts
+    :return: [material
+    """
+    return self._heating_coil_capacity
+
   def to_dictionary(self):
     """Class content to dictionary"""
     _layers = None
@@ -110,7 +119,8 @@ class ThermalStorageSystem(EnergyStorageSystem):
         'height [m]': self.height,
         'layers': _layers,
         'maximum operating temperature [Celsius]': self.maximum_operating_temperature,
-        'storage_medium': self.storage_medium.to_dictionary()
+        'storage_medium': self.storage_medium.to_dictionary(),
+        'heating coil capacity [W]': self.heating_coil_capacity
       }
     }
     return content

hub/catalog_factories/energy_systems/montreal_custom_catalog.py

@@ -87,7 +87,7 @@ class MontrealCustomCatalog(Catalog):
                                                   cooling_efficiency=cooling_efficiency,
                                                   electricity_efficiency=electricity_efficiency,
                                                   energy_storage_systems=storage_systems,
-                                                  dual_supply_capability=False
+                                                  domestic_hot_water=False
                                                   )
       _equipments.append(generation_system)
 
@@ -111,10 +111,7 @@ class MontrealCustomCatalog(Catalog):
         distribution_consumption_variable_flow = float(
           equipment['distribution_consumption_variable_flow']['#text']) / 100
 
-      emission_equipment = -1
-      if 'dissipation_id' in equipment:
-        emission_equipment = equipment['dissipation_id']
-
+      emission_equipment = equipment['dissipation_id']
       _emission_equipments = None
       for equipment_archetype in self._catalog_emission_equipments:
         if int(equipment_archetype.id) == int(emission_equipment):

hub/catalog_factories/energy_systems/montreal_future_system_catalogue.py

@@ -121,13 +121,26 @@ class MontrealFutureSystemCatalogue(Catalog):
           parameters = non_pv['cooling_efficiency_curve']['parameters']
           coefficients = list(non_pv['cooling_efficiency_curve']['coefficients'].values())
           cooling_efficiency_curve = PerformanceCurves(curve_type, dependant_variable, parameters, coefficients)
-        dual_supply_capability = None
-        if non_pv['dual_supply_capability'] is not None:
-          if non_pv['dual_supply_capability'] == 'True':
-            dual_supply_capability = True
+        dhw = None
+        if non_pv['domestic_hot_water'] is not None:
+          if non_pv['domestic_hot_water'] == 'True':
+            dhw = True
           else:
-            dual_supply_capability = False
+            dhw = False
 
+        reversible = None
+        if non_pv['reversible'] is not None:
+          if non_pv['reversible'] == 'True':
+            reversible = True
+          else:
+            reversible = False
+
+        dual_supply = None
+        if non_pv['simultaneous_heat_cold'] is not None:
+          if non_pv['simultaneous_heat_cold'] == 'True':
+            dual_supply = True
+          else:
+            dual_supply = False
         non_pv_component = NonPvGenerationSystem(system_id=system_id,
                                                  name=name,
                                                  system_type=system_type,
@@ -160,7 +173,9 @@ class MontrealFutureSystemCatalogue(Catalog):
                                                  cooling_efficiency_curve=cooling_efficiency_curve,
                                                  distribution_systems=distribution_systems,
                                                  energy_storage_systems=energy_storage_systems,
-                                                 dual_supply_capability=dual_supply_capability)
+                                                 domestic_hot_water=dhw,
+                                                 reversible=reversible,
+                                                 simultaneous_heat_cold=dual_supply)
         generation_components.append(non_pv_component)
     pv_generation_components = self._archetypes['EnergySystemCatalog']['energy_generation_components'][
       'pv_generation_component']
@@ -187,7 +202,6 @@ class MontrealFutureSystemCatalogue(Catalog):
           storage_component = pv['energy_storage_systems']['storage_id']
           storage_systems = self._search_storage_equipment(self._load_storage_components(), storage_component)
           energy_storage_systems = storage_systems
-
         pv_component = PvGenerationSystem(system_id=system_id,
                                           name=name,
                                           system_type=system_type,
@@ -284,6 +298,7 @@ class MontrealFutureSystemCatalogue(Catalog):
       layers = [insulation_layer, tank_layer]
       nominal_capacity = tes['nominal_capacity']
       losses_ratio = tes['losses_ratio']
+      heating_coil_capacity = tes['heating_coil_capacity']
       storage_component = ThermalStorageSystem(storage_id=storage_id,
                                                model_name=model_name,
                                                type_energy_stored=type_energy_stored,
@@ -295,7 +310,8 @@ class MontrealFutureSystemCatalogue(Catalog):
                                                height=height,
                                                layers=layers,
                                                maximum_operating_temperature=maximum_operating_temperature,
-                                               storage_medium=medium)
+                                               storage_medium=medium,
+                                               heating_coil_capacity=heating_coil_capacity)
       storage_components.append(storage_component)
 
     for template in template_storages:
@@ -303,7 +319,7 @@ class MontrealFutureSystemCatalogue(Catalog):
       storage_type = template['storage_type']
       type_energy_stored = template['type_energy_stored']
       maximum_operating_temperature = template['maximum_operating_temperature']
-      height = template['physical_characteristics']['height']
+      height = float(template['physical_characteristics']['height'])
       materials = self._load_materials()
       insulation_material_id = template['insulation']['material_id']
       insulation_material = self._search_material(materials, insulation_material_id)
@@ -322,6 +338,7 @@ class MontrealFutureSystemCatalogue(Catalog):
       nominal_capacity = template['nominal_capacity']
       losses_ratio = template['losses_ratio']
       volume = template['physical_characteristics']['volume']
+      heating_coil_capacity = template['heating_coil_capacity']
       storage_component = ThermalStorageSystem(storage_id=storage_id,
                                                model_name=model_name,
                                                type_energy_stored=type_energy_stored,
@@ -333,7 +350,8 @@ class MontrealFutureSystemCatalogue(Catalog):
                                                height=height,
                                                layers=layers,
                                                maximum_operating_temperature=maximum_operating_temperature,
-                                               storage_medium=medium)
+                                               storage_medium=medium,
+                                               heating_coil_capacity=heating_coil_capacity)
       storage_components.append(storage_component)
     return storage_components
 

hub/city_model_structure/building.py

@@ -90,7 +90,9 @@ class Building(CityObject):
         self._interior_slabs.append(surface)
       else:
         logging.error('Building %s [%s] has an unexpected surface type %s.', self.name, self.aliases, surface.type)
-    self._heating_consumption_disaggregated = {}
+    self._domestic_hot_water_peak_load = None
+    self._fuel_consumption_breakdown = {}
+    self._pv_generation = {}
 
   @property
   def shell(self) -> Polyhedron:
@@ -449,8 +451,8 @@ class Building(CityObject):
       monthly_values = PeakLoads(self).heating_peak_loads_from_methodology
     if monthly_values is None:
       return None
-    results[cte.MONTH] = monthly_values
-    results[cte.YEAR] = [max(monthly_values)]
+    results[cte.MONTH] = [x / cte.WATTS_HOUR_TO_JULES for x in monthly_values]
+    results[cte.YEAR] = [max(monthly_values) / cte.WATTS_HOUR_TO_JULES]
     return results
 
   @property
@@ -466,8 +468,24 @@ class Building(CityObject):
       monthly_values = PeakLoads(self).cooling_peak_loads_from_methodology
     if monthly_values is None:
       return None
-    results[cte.MONTH] = [x * cte.WATTS_HOUR_TO_JULES for x in monthly_values]
-    results[cte.YEAR] = [max(monthly_values)]
+    results[cte.MONTH] = [x / cte.WATTS_HOUR_TO_JULES for x in monthly_values]
+    results[cte.YEAR] = [max(monthly_values) / cte.WATTS_HOUR_TO_JULES]
+    return results
+
+  @property
+  def domestic_hot_water_peak_load(self) -> Union[None, dict]:
+    """
+    Get cooling peak load in W
+    :return: dict{[float]}
+    """
+    results = {}
+    monthly_values = None
+    if cte.HOUR in self.domestic_hot_water_heat_demand:
+      monthly_values = PeakLoads().peak_loads_from_hourly(self.domestic_hot_water_heat_demand[cte.HOUR])
+    if monthly_values is None:
+      return None
+    results[cte.MONTH] = [x / cte.WATTS_HOUR_TO_JULES for x in monthly_values]
+    results[cte.YEAR] = [max(monthly_values) / cte.WATTS_HOUR_TO_JULES]
     return results
 
   @property
@@ -825,23 +843,6 @@ class Building(CityObject):
             self._onsite_electrical_production[_key] = _results
     return self._onsite_electrical_production
 
-  @property
-  def heating_consumption_disaggregated(self) -> dict:
-    """
-    Get energy consumed for heating from different fuels in J
-    return: dict
-    """
-    return self._heating_consumption_disaggregated
-
-  @heating_consumption_disaggregated.setter
-  def heating_consumption_disaggregated(self, value):
-    """
-    Get energy consumed for heating from different fuels in J
-    return: dict
-    """
-    self._heating_consumption_disaggregated = value
-
-
   @property
   def lower_corner(self):
     """
@@ -855,3 +856,60 @@ class Building(CityObject):
     Get building upper corner.
     """
     return [self._max_x, self._max_y, self._max_z]
+
+  @property
+  def energy_consumption_breakdown(self) -> dict:
+    """
+    Get energy consumption of different sectors
+    return: dict
+    """
+    fuel_breakdown = {cte.ELECTRICITY: {cte.LIGHTING: self.lighting_electrical_demand[cte.YEAR][0],
+                                        cte.APPLIANCES: self.appliances_electrical_demand[cte.YEAR][0]}}
+    energy_systems = self.energy_systems
+    for energy_system in energy_systems:
+      demand_types = energy_system.demand_types
+      generation_systems = energy_system.generation_systems
+      for demand_type in demand_types:
+        for generation_system in generation_systems:
+          if generation_system.system_type != cte.PHOTOVOLTAIC:
+            if generation_system.fuel_type not in fuel_breakdown:
+              fuel_breakdown[generation_system.fuel_type] = {}
+            if demand_type in generation_system.energy_consumption:
+              fuel_breakdown[f'{generation_system.fuel_type}'][f'{demand_type}'] = (
+                generation_system.energy_consumption)[f'{demand_type}'][cte.YEAR][0]
+            storage_systems = generation_system.energy_storage_systems
+            if storage_systems:
+              for storage_system in storage_systems:
+                if storage_system.type_energy_stored == 'thermal' and storage_system.heating_coil_energy_consumption:
+                  fuel_breakdown[cte.ELECTRICITY][f'{demand_type}'] += storage_system.heating_coil_energy_consumption[cte.YEAR][0]
+    #TODO: When simulation models of all energy system archetypes are created, this part can be removed
+    heating_fuels = []
+    dhw_fuels = []
+    for energy_system in self.energy_systems:
+      if cte.HEATING in energy_system.demand_types:
+        for generation_system in energy_system.generation_systems:
+          heating_fuels.append(generation_system.fuel_type)
+      if cte.DOMESTIC_HOT_WATER in energy_system.demand_types:
+        for generation_system in energy_system.generation_systems:
+          dhw_fuels.append(generation_system.fuel_type)
+    for key in fuel_breakdown:
+      if key == cte.ELECTRICITY and cte.COOLING not in fuel_breakdown[key]:
+        for energy_system in energy_systems:
+          if cte.COOLING in energy_system.demand_types and cte.COOLING not in fuel_breakdown[key]:
+            for generation_system in energy_system.generation_systems:
+              fuel_breakdown[generation_system.fuel_type][cte.COOLING] = self.cooling_consumption[cte.YEAR][0]
+      for fuel in heating_fuels:
+        if cte.HEATING not in fuel_breakdown[fuel]:
+          for energy_system in energy_systems:
+            if cte.HEATING in energy_system.demand_types:
+              for generation_system in energy_system.generation_systems:
+                fuel_breakdown[generation_system.fuel_type][cte.HEATING] = self.heating_consumption[cte.YEAR][0]
+      for fuel in dhw_fuels:
+        if cte.DOMESTIC_HOT_WATER not in fuel_breakdown[fuel]:
+          for energy_system in energy_systems:
+            if cte.DOMESTIC_HOT_WATER in energy_system.demand_types:
+              for generation_system in energy_system.generation_systems:
+                fuel_breakdown[generation_system.fuel_type][cte.DOMESTIC_HOT_WATER] = self.domestic_hot_water_consumption[cte.YEAR][0]
+    self._fuel_consumption_breakdown = fuel_breakdown
+    return self._fuel_consumption_breakdown
+

hub/city_model_structure/building_demand/internal_zone.py

@@ -132,6 +132,8 @@ class InternalZone:
       _thermal_boundary = ThermalBoundary(surface, surface.solid_polygon.area, windows_areas)
       surface.associated_thermal_boundaries = [_thermal_boundary]
       _thermal_boundaries.append(_thermal_boundary)
+    if self.thermal_archetype is None:
+      return None  # there are no archetype
     _number_of_storeys = int(self.volume / self.area / self.thermal_archetype.average_storey_height)
     _thermal_zone = ThermalZone(_thermal_boundaries, self, self.volume, self.area, _number_of_storeys)
     for thermal_boundary in _thermal_zone.thermal_boundaries:

hub/city_model_structure/building_demand/surface.py

@@ -46,6 +46,8 @@ class Surface:
     self._vegetation = None
     self._percentage_shared = None
     self._solar_collectors_area_reduction_factor = None
+    self._global_irradiance_tilted = {}
+    self._installed_solar_collector_area = None
 
   @property
   def name(self):
@@ -384,3 +386,35 @@ class Surface:
     :param value: float
     """
     self._solar_collectors_area_reduction_factor = value
+
+  @property
+  def global_irradiance_tilted(self) -> dict:
+    """
+    Get global irradiance on a tilted surface in J/m2
+    :return: dict
+    """
+    return self._global_irradiance_tilted
+
+  @global_irradiance_tilted.setter
+  def global_irradiance_tilted(self, value):
+    """
+    Set global irradiance on a tilted surface in J/m2
+    :param value: dict
+    """
+    self._global_irradiance_tilted = value
+
+  @property
+  def installed_solar_collector_area(self):
+    """
+    Get installed solar collector area in m2
+    :return: dict
+    """
+    return self._installed_solar_collector_area
+
+  @installed_solar_collector_area.setter
+  def installed_solar_collector_area(self, value):
+    """
+    Set installed solar collector area in m2
+    :return: dict
+    """
+    self._installed_solar_collector_area = value

hub/city_model_structure/city_object.py

@@ -41,9 +41,10 @@ class CityObject:
     self._ground_temperature = {}
     self._global_horizontal = {}
     self._diffuse = {}
-    self._beam = {}
+    self._direct_normal = {}
     self._sensors = []
     self._neighbours = None
+    self._beam = {}
 
   @property
   def level_of_detail(self) -> LevelOfDetail:
@@ -238,20 +239,20 @@ class CityObject:
     self._diffuse = value
 
   @property
-  def beam(self) -> dict:
+  def direct_normal(self) -> dict:
     """
     Get beam radiation surrounding the city object in J/m2
     :return: dict{dict{[float]}}
     """
-    return self._beam
+    return self._direct_normal
 
-  @beam.setter
-  def beam(self, value):
+  @direct_normal.setter
+  def direct_normal(self, value):
     """
     Set beam radiation surrounding the city object in J/m2
     :param value: dict{dict{[float]}}
     """
-    self._beam = value
+    self._direct_normal = value
 
   @property
   def lower_corner(self):
@@ -302,3 +303,19 @@ class CityObject:
     Set the list of neighbour_objects and their properties associated to the current city_object
     """
     self._neighbours = value
+
+  @property
+  def beam(self) -> dict:
+    """
+    Get beam radiation surrounding the city object in J/m2
+    :return: dict{dict{[float]}}
+    """
+    return self._beam
+
+  @beam.setter
+  def beam(self, value):
+    """
+    Set beam radiation surrounding the city object in J/m2
+    :param value: dict{dict{[float]}}
+    """
+    self._beam = value

hub/city_model_structure/energy_systems/generation_system.py

@@ -11,7 +11,8 @@ from abc import ABC
 from typing import Union, List
 
 from hub.city_model_structure.energy_systems.distribution_system import DistributionSystem
-from hub.city_model_structure.energy_systems.energy_storage_system import EnergyStorageSystem
+from hub.city_model_structure.energy_systems.thermal_storage_system import ThermalStorageSystem
+from hub.city_model_structure.energy_systems.electrical_storage_system import ElectricalStorageSystem
 
 
 class GenerationSystem(ABC):
@@ -26,6 +27,7 @@ class GenerationSystem(ABC):
     self._fuel_type = None
     self._distribution_systems = None
     self._energy_storage_systems = None
+    self._number_of_units = None
 
   @property
   def system_type(self):
@@ -124,7 +126,7 @@ class GenerationSystem(ABC):
     self._distribution_systems = value
 
   @property
-  def energy_storage_systems(self) -> Union[None, List[EnergyStorageSystem]]:
+  def energy_storage_systems(self) -> Union[None, List[ThermalStorageSystem], List[ElectricalStorageSystem]]:
     """
     Get energy storage systems connected to this generation system
     :return: [EnergyStorageSystem]
@@ -138,3 +140,19 @@ class GenerationSystem(ABC):
     :param value: [EnergyStorageSystem]
     """
     self._energy_storage_systems = value
+
+  @property
+  def number_of_units(self):
+    """
+    Get number of a specific generation unit
+    :return: int
+    """
+    return self._number_of_units
+
+  @number_of_units.setter
+  def number_of_units(self, value):
+    """
+    Set number of a specific generation unit
+    :return: int
+    """
+    self._number_of_units = value

hub/city_model_structure/energy_systems/non_pv_generation_system.py

@@ -42,7 +42,12 @@ class NonPvGenerationSystem(GenerationSystem):
     self._cooling_output_curve = None
     self._cooling_fuel_consumption_curve = None
     self._cooling_efficiency_curve = None
-    self._dual_supply_capability = None
+    self._domestic_hot_water = None
+    self._heat_supply_temperature = None
+    self._cooling_supply_temperature = None
+    self._reversible = None
+    self._simultaneous_heat_cold = None
+    self._energy_consumption = {}
 
   @property
   def nominal_heat_output(self):
@@ -429,21 +434,106 @@ class NonPvGenerationSystem(GenerationSystem):
     self._cooling_efficiency_curve = value
 
   @property
-  def dual_supply_capability(self):
+  def domestic_hot_water(self):
     """
-    Get the capability of the generation component for simultaneous heat and cold production
+    Get the capability of generating domestic hot water
 
     :return: bool
     """
-    return self._dual_supply_capability
+    return self._domestic_hot_water
 
-  @dual_supply_capability.setter
-  def dual_supply_capability(self, value):
+  @domestic_hot_water.setter
+  def domestic_hot_water(self, value):
     """
-    Set the capability of the generation component for simultaneous heat and cold production
+    Set the capability of generating domestic hot water
 
     :return: bool
     """
-    self._dual_supply_capability = value
+    self._domestic_hot_water = value
 
+  @property
+  def heat_supply_temperature(self):
+    """
+    Get the hourly heat supply temperature
+    :return: list
+    """
+    return self._heat_supply_temperature
+
+  @heat_supply_temperature.setter
+  def heat_supply_temperature(self, value):
+    """
+    set the hourly heat supply temperature
+    :param value:
+    :return: list
+    """
+    self._heat_supply_temperature = value
+
+  @property
+  def cooling_supply_temperature(self):
+    """
+    Get the hourly cooling supply temperature
+    :return: list
+    """
+    return self._heat_supply_temperature
+
+  @cooling_supply_temperature.setter
+  def cooling_supply_temperature(self, value):
+    """
+    set the hourly cooling supply temperature
+    :param value:
+    :return: list
+    """
+    self._cooling_supply_temperature = value
+
+  @property
+  def reversibility(self):
+    """
+    Get the capability of generating both heating and cooling
+
+    :return: bool
+    """
+    return self._reversible
+
+  @reversibility.setter
+  def reversibility(self, value):
+    """
+    Set the capability of generating domestic hot water
+
+    :return: bool
+    """
+    self._reversible = value
+
+  @property
+  def simultaneous_heat_cold(self):
+    """
+    Get the capability of generating both heating and cooling at the same time
+
+    :return: bool
+    """
+    return self._simultaneous_heat_cold
+
+  @simultaneous_heat_cold.setter
+  def simultaneous_heat_cold(self, value):
+    """
+    Set the capability of generating domestic hot water at the same time
+
+    :return: bool
+    """
+    self._simultaneous_heat_cold = value
+
+  @property
+  def energy_consumption(self) -> dict:
+    """
+    Get energy consumption in W
+    :return: dict{[float]}
+    """
+    return self._energy_consumption
+
+  @energy_consumption.setter
+  def energy_consumption(self, value):
+    """
+    Set energy consumption in W
+    :param value: dict{[float]}
+    """
+    self._energy_consumption = value
 

hub/city_model_structure/energy_systems/performance_curve.py

@@ -24,28 +24,29 @@ class PerformanceCurves:
   def curve_type(self):
     """
     Get the type of the fit function from the following
-    Linear =>>> y = a*x + b
+    Linear =>>> y = a + b*x
     Exponential =>>> y = a*(b**x)
-    Polynomial =>>> y = a*(x**2) + b*x + c
+    Second degree polynomial =>>> y = a + b*x + c*(x**2)
     Power =>>> y = a*(x**b)
-    Second degree multivariable =>>> y = a*(x**2) + b*x + c*x*z + d*z + e*(z**2) + f
+    Bi-Quadratic =>>> y = a + b*x + c*(x**2) + d*z + e*(z**2) + f*x*z
 
-    Get the type of function from ['linear', 'exponential', 'polynomial', 'power', 'second degree multivariable']
+    Get the type of function from ['linear', 'exponential', 'second degree polynomial', 'power', 'bi-quadratic']
     :return: string
     """
+
     return self._curve_type
 
   @curve_type.setter
   def curve_type(self, value):
     """
     Set the type of the fit function from the following
-    Linear =>>> y = a*x + b
+    Linear =>>> y = a + b*x
     Exponential =>>> y = a*(b**x)
-    Polynomial =>>> y = a*(x**2) + b*x + c
+    Second degree polynomial =>>> y = a + b*x + c*(x**2)
     Power =>>> y = a*(x**b)
-    Second degree multivariable =>>> y = a*(x**2) + b*x + c*x*z + d*z + e*(z**2) + f
+    Bi-Quadratic =>>> y = a + b*x + c*(x**2) + d*z + e*(z**2) + f*x*z
 
-    Get the type of function from ['linear', 'exponential', 'polynomial', 'power', 'second degree multivariable']
+    Get the type of function from ['linear', 'exponential', 'second degree polynomial', 'power', 'bi-quadratic']
     :return: string
     """
     self._curve_type = value

hub/city_model_structure/energy_systems/pv_generation_system.py

@@ -26,6 +26,10 @@ class PvGenerationSystem(GenerationSystem):
     self._width = None
     self._height = None
     self._electricity_power = None
+    self._tilt_angle = None
+    self._surface_azimuth = None
+    self._solar_altitude_angle = None
+    self._solar_azimuth_angle = None
 
   @property
   def nominal_electricity_output(self):
@@ -202,3 +206,35 @@ class PvGenerationSystem(GenerationSystem):
     :param value: float
     """
     self._electricity_power = value
+
+  @property
+  def tilt_angle(self):
+    """
+    Get tilt angle of PV system in degrees
+    :return: float
+    """
+    return self._tilt_angle
+
+  @tilt_angle.setter
+  def tilt_angle(self, value):
+    """
+    Set PV system tilt angle in degrees
+    :param value: float
+    """
+    self._tilt_angle = value
+
+  @property
+  def surface_azimuth(self):
+    """
+    Get surface azimuth angle of PV system in degrees. 0 is North
+    :return: float
+    """
+    return self._surface_azimuth
+
+  @surface_azimuth.setter
+  def surface_azimuth(self, value):
+    """
+    Set PV system tilt angle in degrees
+    :param value: float
+    """
+    self._surface_azimuth = value

hub/city_model_structure/energy_systems/thermal_storage_system.py

@@ -22,6 +22,9 @@ class ThermalStorageSystem(EnergyStorageSystem):
     self._height = None
     self._layers = None
     self._maximum_operating_temperature = None
+    self._heating_coil_capacity = None
+    self._temperature = None
+    self._heating_coil_energy_consumption = None
 
   @property
   def volume(self):
@@ -86,3 +89,51 @@ class ThermalStorageSystem(EnergyStorageSystem):
     :param value: float
     """
     self._maximum_operating_temperature = value
+
+  @property
+  def heating_coil_capacity(self):
+    """
+    Get heating coil capacity in Watts
+    :return: float
+    """
+    return self._heating_coil_capacity
+
+  @heating_coil_capacity.setter
+  def heating_coil_capacity(self, value):
+    """
+    Set heating coil capacity in Watts
+    :param value: float
+    """
+    self._heating_coil_capacity = value
+
+  @property
+  def temperature(self) -> dict:
+    """
+    Get fuel consumption in W, m3, or kg
+    :return: dict{[float]}
+    """
+    return self._temperature
+
+  @temperature.setter
+  def temperature(self, value):
+    """
+    Set fuel consumption in W, m3, or kg
+    :param value: dict{[float]}
+    """
+    self._temperature = value
+
+  @property
+  def heating_coil_energy_consumption(self) -> dict:
+    """
+    Get fuel consumption in W, m3, or kg
+    :return: dict{[float]}
+    """
+    return self._heating_coil_energy_consumption
+
+  @heating_coil_energy_consumption.setter
+  def heating_coil_energy_consumption(self, value):
+    """
+    Set fuel consumption in W, m3, or kg
+    :param value: dict{[float]}
+    """
+    self._heating_coil_energy_consumption = value

hub/data/energy_systems/montreal_future_systems.xml

@@ -5,11 +5,11 @@
         <medium>
 			<medium_id>1</medium_id>
 			<name>Water</name>
-			<solar_absorptance></solar_absorptance>
-			<thermal_absorptance></thermal_absorptance>
-			<visible_absorptance></visible_absorptance>
-			<no_mass></no_mass>
-			<thermal_resistance></thermal_resistance>
+			<solar_absorptance/>
+			<thermal_absorptance/>
+			<visible_absorptance/>
+			<no_mass/>
+			<thermal_resistance/>
 			<density>981.0</density>
 			<specific_heat>4180.0</specific_heat>
 			<conductivity>0.6</conductivity>
@@ -26,7 +26,7 @@
             <minimum_heat_output>4.7</minimum_heat_output>
             <maximum_heat_output>23.5</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -50,7 +50,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>2</system_id>
@@ -62,7 +65,7 @@
             <minimum_heat_output>6.15</minimum_heat_output>
             <maximum_heat_output>30.8</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -86,7 +89,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>3</system_id>
@@ -98,7 +104,7 @@
             <minimum_heat_output>8.8</minimum_heat_output>
             <maximum_heat_output>44</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -122,7 +128,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>4</system_id>
@@ -134,7 +143,7 @@
             <minimum_heat_output>12.3</minimum_heat_output>
             <maximum_heat_output>61.5</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -158,7 +167,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>5</system_id>
@@ -170,7 +182,7 @@
             <minimum_heat_output>4.0</minimum_heat_output>
             <maximum_heat_output>35.2</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -194,7 +206,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>6</system_id>
@@ -206,7 +221,7 @@
             <minimum_heat_output>4.0</minimum_heat_output>
             <maximum_heat_output>35.2</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>False</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -230,7 +245,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>7</system_id>
@@ -242,7 +260,7 @@
             <minimum_heat_output>2.5</minimum_heat_output>
             <maximum_heat_output>25.0</maximum_heat_output>
             <heat_efficiency>0.96</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -266,7 +284,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>8</system_id>
@@ -278,7 +299,7 @@
             <minimum_heat_output>3.2</minimum_heat_output>
             <maximum_heat_output>32.0</maximum_heat_output>
             <heat_efficiency>0.96</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -302,7 +323,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>9</system_id>
@@ -314,7 +338,7 @@
             <minimum_heat_output>4.5</minimum_heat_output>
             <maximum_heat_output>45.0</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -338,7 +362,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>10</system_id>
@@ -350,7 +377,7 @@
             <minimum_heat_output>3.5</minimum_heat_output>
             <maximum_heat_output>35.0</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -374,7 +401,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>11</system_id>
@@ -386,7 +416,7 @@
             <minimum_heat_output>5.3</minimum_heat_output>
             <maximum_heat_output>53.0</maximum_heat_output>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi>True</combi>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -410,7 +440,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <pv_generation_component>
             <system_id>12</system_id>
@@ -441,10 +474,10 @@
             <minimum_heat_output>0</minimum_heat_output>
             <maximum_heat_output>51.7</maximum_heat_output>
             <heat_efficiency>3.32</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>Electricity</fuel_type>
             <source_medium>Air</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -460,7 +493,7 @@
             <heat_output_curve/>
             <heat_fuel_consumption_curve/>
             <heat_efficiency_curve>
-			    <curve_type>second degree multivariable function</curve_type>
+			    <curve_type>bi-quadratic</curve_type>
 			    <dependant_variable>COP</dependant_variable>
 		        <parameters>source_temperature</parameters>
 		        <parameters>supply_temperature</parameters>
@@ -471,7 +504,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>False</dual_supply_capability>
+            <domestic_hot_water>False</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>14</system_id>
@@ -483,10 +519,10 @@
             <minimum_heat_output>0</minimum_heat_output>
             <maximum_heat_output>279.3</maximum_heat_output>
             <heat_efficiency>3.07</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>Electricity</fuel_type>
             <source_medium>Air</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -502,7 +538,7 @@
             <heat_output_curve/>
             <heat_fuel_consumption_curve/>
             <heat_efficiency_curve>
-			    <curve_type>second degree multivariable function</curve_type>
+			    <curve_type>bi-quadratic</curve_type>
 			    <dependant_variable>COP</dependant_variable>
 		        <parameters>source_temperature</parameters>
 		        <parameters>supply_temperature</parameters>
@@ -513,7 +549,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>False</dual_supply_capability>
+            <domestic_hot_water>False</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>15</system_id>
@@ -525,10 +564,10 @@
             <minimum_heat_output>0</minimum_heat_output>
             <maximum_heat_output>557</maximum_heat_output>
             <heat_efficiency>3.46</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>Electricity</fuel_type>
             <source_medium>Air</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -544,7 +583,7 @@
             <heat_output_curve/>
             <heat_fuel_consumption_curve/>
             <heat_efficiency_curve>
-			    <curve_type>second degree multivariable function</curve_type>
+			    <curve_type>bi-quadratic</curve_type>
 			    <dependant_variable>COP</dependant_variable>
 		        <parameters>source_temperature</parameters>
 		        <parameters>supply_temperature</parameters>
@@ -555,7 +594,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>False</dual_supply_capability>
+            <domestic_hot_water>False</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>16</system_id>
@@ -567,7 +609,7 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>0.90</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -593,7 +635,10 @@
             <energy_storage_systems>
                 <storage_id>6</storage_id>
             </energy_storage_systems>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>17</system_id>
@@ -605,7 +650,7 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -631,22 +676,25 @@
             <energy_storage_systems>
                 <storage_id>6</storage_id>
             </energy_storage_systems>
-            <dual_supply_capability/>
+            <domestic_hot_water>False</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>18</system_id>
-            <name>template Air-to-Water heat pump</name>
+            <name>template Air-to-Water heat pump with storage</name>
 			<system_type>heat pump</system_type>
             <model_name/>
             <manufacturer/>
             <nominal_heat_output/>
             <minimum_heat_output/>
             <maximum_heat_output/>
-            <heat_efficiency>3</heat_efficiency>
-            <combi/>
+            <heat_efficiency>2</heat_efficiency>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Air</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -661,19 +709,34 @@
             <minimum_cooling_supply_temperature/>
             <heat_output_curve/>
             <heat_fuel_consumption_curve/>
-            <heat_efficiency_curve/>
+            <heat_efficiency_curve>
+			    <curve_type>bi-quadratic</curve_type>
+			    <dependant_variable>COP</dependant_variable>
+		        <parameters>source_temperature</parameters>
+		        <parameters>supply_temperature</parameters>
+		        <coefficients  a="0.132733" b="0.012322" c="0.000032" d="-0.011109" e="-0.000125" f="-0.000123"/>
+            </heat_efficiency_curve>
             <cooling_output_curve/>
             <cooling_fuel_consumption_curve/>
-            <cooling_efficiency_curve/>
+            <cooling_efficiency_curve>
+			    <curve_type>bi-quadratic</curve_type>
+			    <dependant_variable>COP</dependant_variable>
+		        <parameters>source_temperature</parameters>
+		        <parameters>supply_temperature</parameters>
+		        <coefficients  a="0.951894" b="-0.010518" c="0.000126" d="-0.003399" e="0.000183" f="-0.000206"/>
+            </cooling_efficiency_curve>
             <distribution_systems/>
             <energy_storage_systems>
                 <storage_id>6</storage_id>
             </energy_storage_systems>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>True</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>19</system_id>
-            <name>template Groundwater-to-Water heat pump</name>
+            <name>template Groundwater-to-Water heat pump with storage</name>
 			<system_type>heat pump</system_type>
             <model_name/>
             <manufacturer/>
@@ -681,10 +744,10 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>3.5</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Ground</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -707,11 +770,14 @@
             <energy_storage_systems>
                 <storage_id>6</storage_id>
             </energy_storage_systems>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>True</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>20</system_id>
-            <name>template Water-to-Water heat pump</name>
+            <name>template Water-to-Water heat pump with storage</name>
 			<system_type>heat pump</system_type>
             <model_name/>
             <manufacturer/>
@@ -719,10 +785,10 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>3.5</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Water</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -745,7 +811,10 @@
             <energy_storage_systems>
                 <storage_id>6</storage_id>
             </energy_storage_systems>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>21</system_id>
@@ -757,7 +826,7 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>0.90</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>natural gas</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -781,7 +850,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>22</system_id>
@@ -793,7 +865,7 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>0.95</heat_efficiency>
-            <combi/>
+            <reversible>False</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium/>
             <supply_medium/>
@@ -817,7 +889,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability/>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>23</system_id>
@@ -829,14 +904,14 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>3</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Air</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
-            <cooling_efficiency/>
+            <cooling_efficiency>4.5</cooling_efficiency>
             <electricity_efficiency/>
             <source_temperature/>
             <source_mass_flow/>
@@ -847,13 +922,28 @@
             <minimum_cooling_supply_temperature/>
             <heat_output_curve/>
             <heat_fuel_consumption_curve/>
-            <heat_efficiency_curve/>
+            <heat_efficiency_curve>
+			    <curve_type>bi-quadratic</curve_type>
+			    <dependant_variable>COP</dependant_variable>
+		        <parameters>source_temperature</parameters>
+		        <parameters>supply_temperature</parameters>
+		        <coefficients  a="-0.000277" b="0.019639" c="0.000004" d="0.012190" e="-0.00010" f="-0.000277"/>
+            </heat_efficiency_curve>
             <cooling_output_curve/>
             <cooling_fuel_consumption_curve/>
-            <cooling_efficiency_curve/>
+            <cooling_efficiency_curve>
+			    <curve_type>bi-quadratic</curve_type>
+			    <dependant_variable>COP</dependant_variable>
+		        <parameters>source_temperature</parameters>
+		        <parameters>supply_temperature</parameters>
+		        <coefficients  a="0.951894" b="-0.010518" c="0.000126" d="-0.003399" e="0.000183" f="-0.000206"/>
+            </cooling_efficiency_curve>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>True</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>24</system_id>
@@ -865,10 +955,10 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>3.5</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Ground</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -889,7 +979,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>True</simultaneous_heat_cold>
         </non_pv_generation_component>
         <non_pv_generation_component>
             <system_id>25</system_id>
@@ -901,10 +994,10 @@
             <minimum_heat_output/>
             <maximum_heat_output/>
             <heat_efficiency>3.5</heat_efficiency>
-            <combi/>
+            <reversible>True</reversible>
             <fuel_type>electricity</fuel_type>
             <source_medium>Water</source_medium>
-            <supply_medium>water</supply_medium>
+            <supply_medium>Water</supply_medium>
             <nominal_cooling_output/>
             <minimum_cooling_output/>
             <maximum_cooling_output/>
@@ -925,7 +1018,10 @@
             <cooling_efficiency_curve/>
             <distribution_systems/>
             <energy_storage_systems/>
-            <dual_supply_capability>True</dual_supply_capability>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>True</simultaneous_heat_cold>
         </non_pv_generation_component>
         <pv_generation_component>
             <system_id>26</system_id>
@@ -945,7 +1041,55 @@
             <height>1.0</height>
             <distribution_systems/>
             <energy_storage_systems/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
         </pv_generation_component>
+        <non_pv_generation_component>
+            <system_id>27</system_id>
+            <name>template domestic hot water heat pump</name>
+			<system_type>heat pump</system_type>
+            <model_name/>
+            <manufacturer/>
+            <nominal_heat_output/>
+            <minimum_heat_output/>
+            <maximum_heat_output/>
+            <heat_efficiency>3.5</heat_efficiency>
+            <reversible/>
+            <fuel_type>electricity</fuel_type>
+            <source_medium>Air</source_medium>
+            <supply_medium>Water</supply_medium>
+            <nominal_cooling_output/>
+            <minimum_cooling_output/>
+            <maximum_cooling_output/>
+            <cooling_efficiency/>
+            <electricity_efficiency/>
+            <source_temperature/>
+            <source_mass_flow/>
+            <nominal_electricity_output/>
+            <maximum_heat_supply_temperature/>
+            <minimum_heat_supply_temperature/>
+            <maximum_cooling_supply_temperature/>
+            <minimum_cooling_supply_temperature/>
+            <heat_output_curve/>
+            <heat_fuel_consumption_curve/>
+            <heat_efficiency_curve>
+			    <curve_type>bi-quadratic</curve_type>
+			    <dependant_variable>COP</dependant_variable>
+		        <parameters>source_temperature</parameters>
+		        <parameters>supply_temperature</parameters>
+		        <coefficients  a="1.990668" b="0" c="0" d="-0.027252" e="0.000131" f="0"/>
+            </heat_efficiency_curve>
+            <cooling_output_curve/>
+            <cooling_fuel_consumption_curve/>
+            <cooling_efficiency_curve/>
+            <distribution_systems/>
+            <energy_storage_systems>
+                <storage_id>7</storage_id>
+            </energy_storage_systems>
+            <domestic_hot_water>True</domestic_hot_water>
+            <heat_supply_temperature/>
+            <cooling_supply_temperature/>
+            <simultaneous_heat_cold>False</simultaneous_heat_cold>
+        </non_pv_generation_component>
     </energy_generation_components>
     <energy_storage_components>
         <thermalStorages>
@@ -970,8 +1114,9 @@
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
         </thermalStorages>
         <thermalStorages>
 		    <storage_id>2</storage_id>
@@ -995,8 +1140,9 @@
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
         </thermalStorages>
         <thermalStorages>
 		    <storage_id>3</storage_id>
@@ -1020,8 +1166,9 @@
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
         </thermalStorages>
         <thermalStorages>
 		    <storage_id>4</storage_id>
@@ -1044,8 +1191,9 @@
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
         </thermalStorages>
         <thermalStorages>
 		    <storage_id>5</storage_id>
@@ -1069,15 +1217,16 @@
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
         </thermalStorages>
         <templateStorages>
 		    <storage_id>6</storage_id>
 			<name>template Hot Water Storage Tank</name>
 			<type_energy_stored>thermal</type_energy_stored>
-			<model_name>HF 200</model_name>
-			<manufacturer></manufacturer>
+			<model_name/>
+			<manufacturer/>
 			<maximum_operating_temperature>95.0</maximum_operating_temperature>
             <insulation>
 				<material_id>1</material_id>
@@ -1088,47 +1237,74 @@
 				<tankThickness>0</tankThickness>
 				<height>1.5</height>
 				<tankMaterial>Steel</tankMaterial>
-				<volume></volume>
+				<volume/>
 			</physical_characteristics>
             <storage_medium>
 				<medium_id>1</medium_id>
 			</storage_medium>
 			<storage_type>sensible</storage_type>
-			<nominal_capacity></nominal_capacity>
-			<losses_ratio></losses_ratio>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity/>
+        </templateStorages>
+        <templateStorages>
+		    <storage_id>7</storage_id>
+			<name>template Hot Water Storage Tank with Heating Coil</name>
+			<type_energy_stored>thermal</type_energy_stored>
+			<model_name/>
+			<manufacturer/>
+			<maximum_operating_temperature>95.0</maximum_operating_temperature>
+            <insulation>
+				<material_id>1</material_id>
+				<insulationThickness>90.0</insulationThickness>
+			</insulation>
+            <physical_characteristics>
+				<material_id>2</material_id>
+				<tankThickness>0</tankThickness>
+				<height>1.5</height>
+				<tankMaterial>Steel</tankMaterial>
+				<volume/>
+			</physical_characteristics>
+            <storage_medium>
+				<medium_id>1</medium_id>
+			</storage_medium>
+			<storage_type>sensible</storage_type>
+			<nominal_capacity/>
+			<losses_ratio/>
+            <heating_coil_capacity>5000</heating_coil_capacity>
         </templateStorages>
   </energy_storage_components>
   <materials>
     <material>
 		<material_id>1</material_id>
 		<name>Polyurethane</name>
-		<solar_absorptance></solar_absorptance>
-		<thermal_absorptance></thermal_absorptance>
-		<visible_absorptance></visible_absorptance>
-		<no_mass></no_mass>
-		<thermal_resistance></thermal_resistance>
-		<density></density>
-		<specific_heat></specific_heat>
+		<solar_absorptance/>
+		<thermal_absorptance/>
+		<visible_absorptance/>
+		<no_mass/>
+		<thermal_resistance/>
+		<density/>
+		<specific_heat/>
 		<conductivity>0.028</conductivity>
 	</material>
     <material>
 		<material_id>2</material_id>
 		<name>Steel</name>
-		<solar_absorptance></solar_absorptance>
-		<thermal_absorptance></thermal_absorptance>
-		<visible_absorptance></visible_absorptance>
-		<no_mass></no_mass>
-		<thermal_resistance></thermal_resistance>
-		<density></density>
-		<specific_heat></specific_heat>
+		<solar_absorptance/>
+		<thermal_absorptance/>
+		<visible_absorptance/>
+		<no_mass/>
+		<thermal_resistance/>
+		<density/>
+		<specific_heat/>
 		<conductivity>18</conductivity>
 	</material>
   </materials>
   <distribution_systems>
-    <distribution_system></distribution_system>
+    <distribution_system/>
   </distribution_systems>
   <dissipation_systems>
-   <dissipation_system></dissipation_system>
+   <dissipation_system/>
   </dissipation_systems>
   <systems>
 	<system>
@@ -1138,7 +1314,6 @@
 		<demands>
 			<demand>heating</demand>
 			<demand>cooling</demand>
-			<demand>domestic_hot_water</demand>
 		</demands>
 		<components>
 			<generation_id>21</generation_id>
@@ -1226,13 +1401,75 @@
 			<generation_id>26</generation_id>
 		</components>
 	</system>
+    <system>
+		<id>8</id>
+		<name>4 pipe system with air source heat pump storage and gas boiler</name>
+		<schema>schemas/ASHP+TES+GasBoiler.jpg</schema>
+		<demands>
+			<demand>heating</demand>
+			<demand>cooling</demand>
+		</demands>
+		<components>
+			<generation_id>23</generation_id>
+			<generation_id>16</generation_id>
+		</components>
+	</system>
+    <system>
+		<id>9</id>
+		<name>4 pipe system with air source heat pump storage and electric boiler</name>
+		<schema>schemas/ASHP+TES+GasBoiler.jpg</schema>
+		<demands>
+			<demand>heating</demand>
+			<demand>cooling</demand>
+		</demands>
+		<components>
+			<generation_id>22</generation_id>
+			<generation_id>18</generation_id>
+		</components>
+	</system>
+    <system>
+		<id>10</id>
+		<name>Domestic Hot Water Heat Pump with Coiled Storage</name>
+		<schema>schemas/ASHP+TES+GasBoiler.jpg</schema>
+		<demands>
+			<demand>domestic_hot_water</demand>
+		</demands>
+		<components>
+			<generation_id>27</generation_id>
+		</components>
+	</system>
+    <system>
+		<id>11</id>
+		<name>Central Heating System َASHP Gas-Boiler TES</name>
+		<schema>schemas/ASHP+TES+GasBoiler.jpg</schema>
+		<demands>
+			<demand>heating</demand>
+		</demands>
+		<components>
+			<generation_id>23</generation_id>
+			<generation_id>16</generation_id>
+		</components>
+	</system>
+    <system>
+		<id>12</id>
+		<name>Unitary ASHP Cooling System</name>
+		<schema>schemas/ASHP+TES+GasBoiler.jpg</schema>
+		<demands>
+			<demand>cooling</demand>
+		</demands>
+		<components>
+			<generation_id>23</generation_id>
+		</components>
+	</system>
   </systems>
+
   <system_archetypes>
 	<system_archetype id="1">
 		<name>PV+ASHP+GasBoiler+TES</name>
 		<systems>
 			<system_id>7</system_id>
 			<system_id>1</system_id>
+            <system_id>10</system_id>
 		</systems>
 	</system_archetype>
 	<system_archetype id="2">
@@ -1306,6 +1543,31 @@
 			<system_id>6</system_id>
 		</systems>
 	</system_archetype>
+	<system_archetype id="13">
+		<name>PV+4Pipe+DHW</name>
+		<systems>
+			<system_id>7</system_id>
+			<system_id>8</system_id>
+            <system_id>10</system_id>
+		</systems>
+	</system_archetype>
+	<system_archetype id="14">
+		<name>Central Heating+Unitary Cooling+Unitary DHW</name>
+		<systems>
+			<system_id>10</system_id>
+			<system_id>11</system_id>
+            <system_id>12</system_id>
+		</systems>
+	</system_archetype>
+	<system_archetype id="15">
+		<name>Central Heating+Unitary Cooling+Unitary DHW+PV</name>
+		<systems>
+			<system_id>7</system_id>
+            <system_id>10</system_id>
+			<system_id>11</system_id>
+            <system_id>12</system_id>
+		</systems>
+	</system_archetype>
   </system_archetypes>
 </EnergySystemCatalog>
 

hub/exports/building_energy/idf.py

@@ -513,7 +513,7 @@ class Idf:
     self._lod = self._city.level_of_detail.geometry
     for building in self._city.buildings:
       for internal_zone in building.internal_zones:
-        if internal_zone.thermal_zones_from_internal_zones is None or len(thermal_zone.thermal_boundaries) == 0:
+        if internal_zone.thermal_zones_from_internal_zones is None:
           self._add_shading(building)  # if the building has not archetype use it as shadow object
           continue
         for thermal_zone in internal_zone.thermal_zones_from_internal_zones:

hub/exports/building_energy/insel/insel_monthly_energy_balance.py

@@ -270,7 +270,7 @@ class InselMonthlyEnergyBalance:
         global_irradiance = surface.global_irradiance[cte.MONTH]
         for j in range(0, len(global_irradiance)):
           parameters.append(f'{j + 1} '
-                            f'{global_irradiance[j] * cte.WATTS_HOUR_TO_JULES / 24 / _NUMBER_DAYS_PER_MONTH[j]}')
+                            f'{global_irradiance[j] / 24 / _NUMBER_DAYS_PER_MONTH[j]}')
       else:
         for j in range(0, 12):
           parameters.append(f'{j + 1} 0.0')

hub/exports/formats/simplified_radiosity_algorithm.py

@@ -66,9 +66,9 @@ class SimplifiedRadiosityAlgorithm:
           else:
             i = (total_days + day - 1) * 24 + hour - 1
           representative_building = self._city.buildings[0]
-          _diffuse = representative_building.diffuse[cte.HOUR][i] * cte.WATTS_HOUR_TO_JULES
-          _beam = representative_building.beam[cte.HOUR][i] * cte.WATTS_HOUR_TO_JULES
-          content += f'{day}  {month}  {hour}  {_diffuse}  {_beam}\n'
+          _global = representative_building.diffuse[cte.HOUR][i] / cte.WATTS_HOUR_TO_JULES
+          _beam = representative_building.direct_normal[cte.HOUR][i] / cte.WATTS_HOUR_TO_JULES
+          content += f'{day}  {month}  {hour}  {_global}  {_beam}\n'
     with open(file, 'w', encoding='utf-8') as file:
       file.write(content)
 

hub/helpers/constants.py

@@ -10,11 +10,11 @@ Project Coder Pilar Monsalvete Alvarez de Uribarri pilar.monsalvete@concordia.ca
 KELVIN = 273.15
 WATER_DENSITY = 1000  # kg/m3
 WATER_HEAT_CAPACITY = 4182  # J/kgK
-
+WATER_THERMAL_CONDUCTIVITY = 0.65  # W/mK
+NATURAL_GAS_LHV = 36.6e6  # J/m3
 AIR_DENSITY = 1.293  # kg/m3
 AIR_HEAT_CAPACITY = 1005.2  # J/kgK
 
-
 # converters
 HOUR_TO_MINUTES = 60
 MINUTES_TO_SECONDS = 60
@@ -292,6 +292,7 @@ WOOD = 'Wood'
 GAS = 'Gas'
 DIESEL = 'Diesel'
 COAL = 'Coal'
+BIOMASS = 'Biomass'
 AIR = 'Air'
 WATER = 'Water'
 GEOTHERMAL = 'Geothermal'

hub/imports/energy_systems/montreal_custom_energy_system_parameters.py

@@ -86,7 +86,8 @@ class MontrealCustomEnergySystemParameters:
       if archetype_generation_system.system_type == 'Photovoltaic':
         _generation_system = PvGenerationSystem()
         _type = 'PV system'
-        _generation_system.system_type = Dictionaries().montreal_generation_system_to_hub_energy_generation_system[_type]
+        _generation_system.system_type = Dictionaries().montreal_generation_system_to_hub_energy_generation_system[
+          _type]
         _fuel_type = Dictionaries().montreal_custom_fuel_to_hub_fuel[archetype_generation_system.fuel_type]
         _generation_system.fuel_type = _fuel_type
         _generation_system.electricity_efficiency = archetype_generation_system.electricity_efficiency
@@ -98,7 +99,8 @@ class MontrealCustomEnergySystemParameters:
       else:
         _generation_system = NonPvGenerationSystem()
         _type = archetype_generation_system.system_type
-        _generation_system.system_type = Dictionaries().montreal_generation_system_to_hub_energy_generation_system[_type]
+        _generation_system.system_type = Dictionaries().montreal_generation_system_to_hub_energy_generation_system[
+          _type]
         _fuel_type = Dictionaries().montreal_custom_fuel_to_hub_fuel[archetype_generation_system.fuel_type]
         _generation_system.fuel_type = _fuel_type
         _generation_system.source_types = archetype_generation_system.source_medium
@@ -136,8 +138,12 @@ class MontrealCustomEnergySystemParameters:
       _distribution_system.heat_losses = archetype_distribution_system.heat_losses
       _emission_system = None
       if archetype_distribution_system.emission_systems is not None:
-        _emission_system = EmissionSystem()
-      _distribution_system.emission_systems = [_emission_system]
+        _emission_systems = []
+        for emission_system in archetype_distribution_system.emission_systems:
+          _emission_system = EmissionSystem()
+          _emission_system.parasitic_energy_consumption = emission_system.parasitic_energy_consumption
+          _emission_systems.append(_emission_system)
+      _distribution_system.emission_systems = _emission_systems
       _distribution_systems.append(_distribution_system)
     return _distribution_systems
 

hub/imports/energy_systems/montreal_future_energy_systems_parameters.py

@@ -82,8 +82,7 @@ class MontrealFutureEnergySystemParameters:
 
     return _generic_energy_systems
 
-  @staticmethod
-  def _create_generation_systems(archetype_system):
+  def _create_generation_systems(self, archetype_system):
     _generation_systems = []
     archetype_generation_systems = archetype_system.generation_systems
     if archetype_generation_systems is not None:
@@ -107,6 +106,7 @@ class MontrealFutureEnergySystemParameters:
           _generation_system.cell_temperature_coefficient = archetype_generation_system.cell_temperature_coefficient
           _generation_system.width = archetype_generation_system.width
           _generation_system.height = archetype_generation_system.height
+          _generation_system.tilt_angle = self._city.latitude
           _generic_storage_system = None
           if archetype_generation_system.energy_storage_systems is not None:
             _generic_storage_system = ElectricalStorageSystem()
@@ -140,12 +140,13 @@ class MontrealFutureEnergySystemParameters:
           _generation_system.cooling_output_curve = archetype_generation_system.cooling_output_curve
           _generation_system.cooling_fuel_consumption_curve = archetype_generation_system.cooling_fuel_consumption_curve
           _generation_system.cooling_efficiency_curve = archetype_generation_system.cooling_efficiency_curve
-          _generation_system.dual_supply_capability = archetype_generation_system.dual_supply_capability
+          _generation_system.domestic_hot_water = archetype_generation_system.domestic_hot_water
           _generation_system.nominal_electricity_output = archetype_generation_system.nominal_electricity_output
           _generation_system.source_medium = archetype_generation_system.source_medium
           _generation_system.heat_efficiency = archetype_generation_system.heat_efficiency
           _generation_system.cooling_efficiency = archetype_generation_system.cooling_efficiency
           _generation_system.electricity_efficiency = archetype_generation_system.electricity_efficiency
+          _generation_system.reversibility = archetype_generation_system.reversibility
           _generic_storage_system = None
           if archetype_generation_system.energy_storage_systems is not None:
             _storage_systems = []
@@ -155,11 +156,19 @@ class MontrealFutureEnergySystemParameters:
                 _generic_storage_system.type_energy_stored = 'electrical'
               else:
                 _generic_storage_system = ThermalStorageSystem()
-                _generic_storage_system.type_energy_stored = 'thermal'
+                _generic_storage_system.type_energy_stored = storage_system.type_energy_stored
+                _generic_storage_system.height = storage_system.height
+                _generic_storage_system.layers = storage_system.layers
+                _generic_storage_system.storage_medium = storage_system.storage_medium
+                _generic_storage_system.heating_coil_capacity = storage_system.heating_coil_capacity
               _storage_systems.append(_generic_storage_system)
-            _generation_system.energy_storage_systems = [_storage_systems]
-          if archetype_generation_system.dual_supply_capability:
-            _generation_system.dual_supply_capability = True
+            _generation_system.energy_storage_systems = _storage_systems
+          if archetype_generation_system.domestic_hot_water:
+            _generation_system.domestic_hot_water = True
+          if archetype_generation_system.reversibility:
+            _generation_system.reversibility = True
+          if archetype_generation_system.simultaneous_heat_cold:
+            _generation_system.simultaneous_heat_cold = True
         _generation_systems.append(_generation_system)
     return _generation_systems
 
@@ -178,8 +187,12 @@ class MontrealFutureEnergySystemParameters:
         _distribution_system.heat_losses = archetype_distribution_system.heat_losses
         _emission_system = None
         if archetype_distribution_system.emission_systems is not None:
-          _emission_system = EmissionSystem()
-        _distribution_system.emission_systems = [_emission_system]
+          _emission_systems = []
+          for emission_system in archetype_distribution_system.emission_systems:
+            _emission_system = EmissionSystem()
+            _emission_system.parasitic_energy_consumption = emission_system.parasitic_energy_consumption
+            _emission_systems.append(_emission_system)
+        _distribution_system.emission_systems = _emission_systems
         _distribution_systems.append(_distribution_system)
       return _distribution_systems
 

hub/imports/results/ep_multiple_buildings.py

@@ -60,9 +60,12 @@ class EnergyPlusMultipleBuildings:
       for building in self._city.buildings:
         building.heating_demand[cte.HOUR] = building_energy_demands[f'Building {building.name} Heating Demand (J)']
         building.cooling_demand[cte.HOUR] = building_energy_demands[f'Building {building.name} Cooling Demand (J)']
-        building.domestic_hot_water_heat_demand[cte.HOUR] = building_energy_demands[f'Building {building.name} DHW Demand (W)']
-        building.appliances_electrical_demand[cte.HOUR] = building_energy_demands[f'Building {building.name} Appliances (W)']
-        building.lighting_electrical_demand[cte.HOUR] = building_energy_demands[f'Building {building.name} Lighting (W)']
+        building.domestic_hot_water_heat_demand[cte.HOUR] = \
+            [x * cte.WATTS_HOUR_TO_JULES for x in building_energy_demands[f'Building {building.name} DHW Demand (W)']]
+        building.appliances_electrical_demand[cte.HOUR] = \
+            [x * cte.WATTS_HOUR_TO_JULES for x in building_energy_demands[f'Building {building.name} Appliances (W)']]
+        building.lighting_electrical_demand[cte.HOUR] = \
+            [x * cte.WATTS_HOUR_TO_JULES for x in building_energy_demands[f'Building {building.name} Lighting (W)']]
         building.heating_demand[cte.MONTH] = MonthlyValues.get_total_month(building.heating_demand[cte.HOUR])
         building.cooling_demand[cte.MONTH] = MonthlyValues.get_total_month(building.cooling_demand[cte.HOUR])
         building.domestic_hot_water_heat_demand[cte.MONTH] = (

hub/imports/results/simplified_radiosity_algorithm.py

@@ -34,7 +34,7 @@ class SimplifiedRadiosityAlgorithm:
     for key in self._results:
       _irradiance = {}
       header_name = key.split(':')
-      result = [x / cte.WATTS_HOUR_TO_JULES for x in self._results[key]]
+      result = [x * cte.WATTS_HOUR_TO_JULES for x in self._results[key]]
       city_object_name = header_name[1]
       building = self._city.city_object(city_object_name)
       surface_id = header_name[2]

hub/imports/weather/epw_weather_parameters.py

@@ -110,22 +110,26 @@ class EpwWeatherParameters:
 #      new_value = pd.DataFrame(self._weather_values[['dry_bulb_temperature_c']].to_numpy(), columns=['epw'])
 #      number_invalid_records = new_value[new_value.epw == 99.9].count().epw
       building.external_temperature[cte.HOUR] = self._weather_values['dry_bulb_temperature_c']
-      building.global_horizontal[cte.HOUR] = [x / cte.WATTS_HOUR_TO_JULES
+      building.global_horizontal[cte.HOUR] = [x * cte.WATTS_HOUR_TO_JULES
                                               for x in self._weather_values['global_horizontal_radiation_wh_m2']]
-      building.diffuse[cte.HOUR] = [x / cte.WATTS_HOUR_TO_JULES
+      building.diffuse[cte.HOUR] = [x * cte.WATTS_HOUR_TO_JULES
                                     for x in self._weather_values['diffuse_horizontal_radiation_wh_m2']]
-      building.beam[cte.HOUR] = [x / cte.WATTS_HOUR_TO_JULES
-                                 for x in self._weather_values['direct_normal_radiation_wh_m2']]
+      building.direct_normal[cte.HOUR] = [x * cte.WATTS_HOUR_TO_JULES
+                                          for x in self._weather_values['direct_normal_radiation_wh_m2']]
+      building.beam[cte.HOUR] = [building.global_horizontal[cte.HOUR][i] -
+                                 building.diffuse[cte.HOUR][i]
+                                 for i in range(len(building.global_horizontal[cte.HOUR]))]
       building.cold_water_temperature[cte.HOUR] = wh().cold_water_temperature(building.external_temperature[cte.HOUR])
 
+
     # create the monthly and yearly values out of the hourly
     for building in self._city.buildings:
       building.external_temperature[cte.MONTH] = \
         MonthlyValues().get_mean_values(building.external_temperature[cte.HOUR])
-      building.external_temperature[cte.YEAR] = [sum(building.external_temperature[cte.HOUR]) / 9870]
+      building.external_temperature[cte.YEAR] = [sum(building.external_temperature[cte.HOUR]) / 8760]
       building.cold_water_temperature[cte.MONTH] = \
         MonthlyValues().get_mean_values(building.cold_water_temperature[cte.HOUR])
-      building.cold_water_temperature[cte.YEAR] = [sum(building.cold_water_temperature[cte.HOUR]) / 9870]
+      building.cold_water_temperature[cte.YEAR] = [sum(building.cold_water_temperature[cte.HOUR]) / 8760]
 
       # If the usage has already being imported, the domestic hot water missing values must be calculated here that
       # the cold water temperature is finally known

hub/imports/weather/helpers/weather.py

@@ -8,7 +8,7 @@ Project Coder Pilar Monsalvete Alvarez de Uribarri pilar.monsalvete@concordia.ca
 import logging
 import math
 import hub.helpers.constants as cte
-
+from datetime import datetime, timedelta
 
 class Weather:
   """
@@ -55,25 +55,19 @@ class Weather:
     # and Craig Christensen, National Renewable Energy Laboratory
     # ambient temperatures( in °C)
     # cold water temperatures( in °C)
-    ambient_temperature_fahrenheit = []
-    average_temperature = 0
-    maximum_temperature = -1000
-    minimum_temperature = 1000
-    for temperature in ambient_temperature:
-      value = temperature * 9 / 5 + 32
-      ambient_temperature_fahrenheit.append(value)
-      average_temperature += value / 8760
-      if value > maximum_temperature:
-        maximum_temperature = value
-      if value < minimum_temperature:
-        minimum_temperature = value
-    delta_temperature = maximum_temperature - minimum_temperature
-    ratio = 0.4 + 0.01 * (average_temperature - 44)
-    lag = 35 - 1 * (average_temperature - 44)
+    t_out_fahrenheit = [1.8 * t_out + 32 for t_out in ambient_temperature]
+    t_out_average = sum(t_out_fahrenheit) / len(t_out_fahrenheit)
+    max_difference = max(t_out_fahrenheit) - min(t_out_fahrenheit)
+    ratio = 0.4 + 0.01 * (t_out_average - 44)
+    lag = 35 - (t_out_average - 35)
+    number_of_day = [a for a in range(1, 366)]
+    day_of_year = [day for day in number_of_day for _ in range(24)]
+    cold_temperature_fahrenheit = []
     cold_temperature = []
-    for temperature in ambient_temperature_fahrenheit:
-      radians = (0.986 * (temperature-15-lag) - 90) * math.pi / 180
-      cold_temperature.append((average_temperature + 6 + ratio * (delta_temperature/2) * math.sin(radians) - 32) * 5/9)
+    for i in range(len(ambient_temperature)):
+      cold_temperature_fahrenheit.append(t_out_average + 6 + ratio * (max_difference / 2) *
+                                         math.sin(math.radians(0.986 * (day_of_year[i] - 15 - lag) - 90)))
+      cold_temperature.append((cold_temperature_fahrenheit[i] - 32) / 1.8)
     return cold_temperature
 
   def epw_file(self, region_code):

tests/test_geometry_factory.py

@@ -8,7 +8,7 @@ from pathlib import Path
 from unittest import TestCase
 
 import hub.exports.exports_factory
-from hub.helpers.dictionaries import MontrealFunctionToHubFunction
+from hub.helpers.dictionaries import Dictionaries
 from hub.helpers.geometry_helper import GeometryHelper
 from hub.imports.construction_factory import ConstructionFactory
 from hub.imports.geometry_factory import GeometryFactory
@@ -36,6 +36,7 @@ class TestGeometryFactory(TestCase):
                                  height_field=height_field,
                                  year_of_construction_field=year_of_construction_field,
                                  function_field=function_field,
+                                 function_to_hub=Dictionaries().montreal_function_to_hub_function
                                  ).city
     self.assertIsNotNone(self._city, 'city is none')
     return self._city
@@ -104,6 +105,7 @@ class TestGeometryFactory(TestCase):
     """
     file = 'FZK_Haus_LoD_2.gml'
     city = self._get_city(file, 'citygml')
+    print('test')
     self.assertTrue(len(city.buildings) == 1)
     self._check_buildings(city)
     for building in city.buildings:
@@ -132,7 +134,7 @@ class TestGeometryFactory(TestCase):
                            year_of_construction_field='ANNEE_CONS',
                            aliases_field=['ID_UEV', 'CIVIQUE_DE', 'NOM_RUE'],
                            function_field='CODE_UTILI',
-                           function_to_hub=MontrealFunctionToHubFunction().dictionary).city
+                           function_to_hub=Dictionaries().montreal_function_to_hub_function).city
     hub.exports.exports_factory.ExportsFactory('obj', city, self._output_path).export()
     for building in city.building_alias('01002777'):
       self.assertEqual('1', building.name, 'Wrong building name when looking for alias')

tests/test_systems_catalog.py

@@ -38,12 +38,12 @@ class TestSystemsCatalog(TestCase):
     catalog = EnergySystemsCatalogFactory('montreal_future').catalog
 
     catalog_categories = catalog.names()
-    archetypes = catalog.names('archetypes')
-    self.assertEqual(12, len(archetypes['archetypes']))
+    archetypes = catalog.names()
+    self.assertEqual(15, len(archetypes['archetypes']))
     systems = catalog.names('systems')
-    self.assertEqual(7, len(systems['systems']))
+    self.assertEqual(12, len(systems['systems']))
     generation_equipments = catalog.names('generation_equipments')
-    self.assertEqual(26, len(generation_equipments['generation_equipments']))
+    self.assertEqual(27, len(generation_equipments['generation_equipments']))
     with self.assertRaises(ValueError):
       catalog.names('unknown')
 

tests/test_systems_factory.py

@@ -122,14 +122,13 @@ class TestSystemsFactory(TestCase):
       for energy_system in building.energy_systems:
         if cte.HEATING in energy_system.demand_types:
           _generation_system = cast(NonPvGenerationSystem, energy_system.generation_systems[0])
-          _generation_system.heat_power = building.heating_peak_load[cte.YEAR][0]
+          _generation_system.nominal_heat_output = building.heating_peak_load[cte.YEAR][0]
         if cte.COOLING in energy_system.demand_types:
           _generation_system = cast(NonPvGenerationSystem, energy_system.generation_systems[0])
-          _generation_system.cooling_power = building.cooling_peak_load[cte.YEAR][0]
+          _generation_system.nominal_cooling_output = building.cooling_peak_load[cte.YEAR][0]
 
     for building in self._city.buildings:
       self.assertLess(0, building.heating_consumption[cte.YEAR][0])
       self.assertLess(0, building.cooling_consumption[cte.YEAR][0])
       self.assertLess(0, building.domestic_hot_water_consumption[cte.YEAR][0])
-      self.assertLess(0, building.onsite_electrical_production[cte.YEAR][0])
-  print('test')
+      self.assertLess(0, building.onsite_electrical_production[cte.YEAR][0])