openstudio-standards 0.2.12.rc4 → 0.2.13.rc3

data/lib/openstudio-standards/standards/necb/ECMS/nv.rb

@@ -0,0 +1,189 @@
+class ECMS
+
+  def apply_nv(model:, nv_type:, nv_opening_fraction:, nv_Tout_min:, nv_Delta_Tin_Tout:)
+
+    ##### If any of users' inputs are nil/false/none, do nothing.
+    ##### If users' input for 'nv_type' is 'NECB_Default', do nothing.
+    ##### If any of users' inputs for nv_opening_fraction/nv_Tout_min/nv_Delta_Tin_Tout is 'NECB_Default', use default values as defined here.
+    return if nv_type.nil? || nv_type == false || nv_type == 'none' || nv_type == 'NECB_Default'
+    return if nv_opening_fraction.nil? || nv_opening_fraction == false || nv_opening_fraction == 'none'
+    return if nv_Tout_min.nil? || nv_Tout_min == false || nv_Tout_min == 'none'
+    return if nv_Delta_Tin_Tout.nil? || nv_Delta_Tin_Tout == false || nv_Delta_Tin_Tout == 'none'
+
+    ##### Convert a string to a float (except for nv_type)
+    if nv_opening_fraction.instance_of?(String) && nv_opening_fraction != 'NECB_Default'
+      nv_opening_fraction = nv_opening_fraction.to_f
+    end
+    if nv_Tout_min.instance_of?(String) && nv_Tout_min != 'NECB_Default'
+      nv_Tout_min = nv_Tout_min.to_f
+    end
+    if nv_Delta_Tin_Tout.instance_of?(String) && nv_Delta_Tin_Tout != 'NECB_Default'
+      nv_Delta_Tin_Tout = nv_Delta_Tin_Tout.to_f
+    end
+
+    ##### Set default nv_opening_fraction as 0.1
+    if nv_opening_fraction == 'NECB_Default'
+      nv_opening_fraction = 0.1
+    end
+    ##### Set default nv_Tout_min as 13.0
+    if nv_Tout_min == 'NECB_Default'
+      nv_Tout_min = 13.0 #Note: 13.0 is based on inputs from Michel Tardif re a real school in QC
+    end
+    ##### Set default nv_Delta_Tin_Tout as 1.0
+    if nv_Delta_Tin_Tout == 'NECB_Default'
+      nv_Delta_Tin_Tout = 1.0 #Note: 1.0 is based on inputs from Michel Tardif re a real school in QC
+    end
+
+    setpoint_adjustment_for_nv = 2.0  #This is to adjust heating and cooling setpoint temperature as min and max indoor temperature to have NV
+
+    model.getZoneHVACEquipmentLists.sort.each do |zone_hvac_equipment_list|
+
+      thermal_zone = zone_hvac_equipment_list.thermalZone
+
+      thermal_zone.spaces.sort.each do |space|
+        number_of_windows = 0.0
+
+        ##### Gather OA per person and floor area of the space from the osm file
+        outdoor_air = space.designSpecificationOutdoorAir.get
+        outdoor_air_flow_per_person = outdoor_air.outdoorAirFlowperPerson
+        outdoor_air_flow_per_floor_area = outdoor_air.outdoorAirFlowperFloorArea
+
+        ##### Get heating/cooling setpoint temperature schedules from the osm file
+        # These schedules are used for min/max Tin schedules under the objects of "ZoneVentilation:DesignFlowRate" and "ZoneVentilation:WindandStackOpenArea".
+        # Note: as per E+ I/O Ref.: "If the user enters a valid schedule name, the minimum/maximum temperature values specified in this schedule will override the constant value specified in the Minimum/Maximum Indoor Temperature field." under the objects of "ZoneVentilation:DesignFlowRate" and "ZoneVentilation:WindandStackOpenArea".
+        if thermal_zone.thermostat.is_initialized
+          if thermal_zone.thermostat.get.to_ThermostatSetpointDualSetpoint.is_initialized
+            if thermal_zone.thermostat.get.to_ThermostatSetpointDualSetpoint.get.heatingSetpointTemperatureSchedule.is_initialized ||
+                thermal_zone.thermostat.get.to_ThermostatSetpointDualSetpoint.get.coolingSetpointTemperatureSchedule.is_initialized
+              zone_thermostat = thermal_zone.thermostatSetpointDualSetpoint.get
+              zone_clg_thermostat_sch = zone_thermostat.coolingSetpointTemperatureSchedule.get
+              zone_htg_thermostat_sch = zone_thermostat.heatingSetpointTemperatureSchedule.get
+
+              ##### Create schedule for max Tin to have NV on the basis of cooling setpoint temperature for default day
+              zone_clg_thermostat_sch_name = zone_clg_thermostat_sch.name
+              zone_clg_sp_schedule = zone_clg_thermostat_sch.to_ScheduleRuleset.get
+
+              max_Tin_schedule = zone_clg_sp_schedule.clone(model).to_ScheduleRuleset.get
+              max_Tin_schedule.setName('natural_ventilation_max_Tin_schedule')
+              ## default days/weekdays
+              max_Tin_schedule_defaultDay = max_Tin_schedule.defaultDaySchedule
+              max_Tin_schedule_defaultDay.setName('natural_ventilation_max_Tin_schedule_defaultDay')
+              max_Tin_schedule_defaultDay_times = max_Tin_schedule_defaultDay.times
+              max_Tin_schedule_defaultDay_values = max_Tin_schedule_defaultDay.values
+              max_Tin_schedule_defaultDay_values_adjusted = max_Tin_schedule_defaultDay_values.map { |i| i + setpoint_adjustment_for_nv }
+              i = 0.0
+              max_Tin_schedule_defaultDay_times.each do |time|
+                max_Tin_schedule_defaultDay.addValue(time, max_Tin_schedule_defaultDay_values_adjusted[i])
+                i += 1.0
+              end
+
+              ##### Create schedule for min Tin to have NV on the basis of cooling setpoint temperature for default day
+              zone_htg_thermostat_sch_name = zone_htg_thermostat_sch.name
+              zone_htg_sp_schedule = zone_htg_thermostat_sch.to_ScheduleRuleset.get
+              min_Tin_schedule = zone_htg_sp_schedule.clone(model).to_ScheduleRuleset.get
+              min_Tin_schedule.setName('natural_ventilation_min_Tin_schedule')
+              ## default days/weekdays
+              min_Tin_schedule_defaultDay = min_Tin_schedule.defaultDaySchedule
+              min_Tin_schedule_defaultDay.setName('natural_ventilation_min_Tin_schedule_defaultDay')
+              min_Tin_schedule_defaultDay_times = min_Tin_schedule_defaultDay.times
+              min_Tin_schedule_defaultDay_values = min_Tin_schedule_defaultDay.values
+              min_Tin_schedule_defaultDay_values_adjusted = min_Tin_schedule_defaultDay_values.map { |i| i - setpoint_adjustment_for_nv }
+              i = 0.0
+              min_Tin_schedule_defaultDay_times.each do |time|
+                min_Tin_schedule_defaultDay.addValue(time, min_Tin_schedule_defaultDay_values_adjusted[i])
+                i += 1.0
+              end
+
+
+            end
+          end
+        end
+
+        ##### Calculate how many windows a space has.
+        # The total number of windows is used to divide OA/person and OA/FloorArea of the space by it (i.e. number of windows).
+        # In this way, NV-driven OA in each space would be avoided to be more than required.
+        space.surfaces.sort.each do |surface|
+          surface.subSurfaces.sort.each do |subsurface|
+            if (subsurface.subSurfaceType == 'OperableWindow' || subsurface.subSurfaceType == 'FixedWindow') && subsurface.outsideBoundaryCondition == 'Outdoors'
+              number_of_windows += 1.0
+            end
+          end
+        end
+        oa_per_person_normalized_by_number_of_windows = outdoor_air_flow_per_person/number_of_windows
+        oa_per_floor_area_normalized_by_number_of_windows = outdoor_air_flow_per_floor_area/number_of_windows
+
+        ##### Add NV in each space that has window(s) using two objects: "ZoneVentilation:DesignFlowRate" and "ZoneVentilation:WindandStackOpenArea"
+        space.surfaces.sort.each do |surface|
+          surface.subSurfaces.sort.each do |subsurface|
+            if (subsurface.subSurfaceType == 'OperableWindow' || subsurface.subSurfaceType == 'FixedWindow') && subsurface.outsideBoundaryCondition == 'Outdoors'
+              window_azimuth_deg = OpenStudio::convert(subsurface.azimuth,"rad","deg").get
+              window_area = subsurface.netArea
+
+              ##### Define a constant schedule for operable windows
+              operable_window_schedule = OpenStudio::Model::ScheduleConstant.new(model)
+              operable_window_schedule.setName('operable_window_schedule_constant')
+              operable_window_schedule.setScheduleTypeLimits(BTAP::Resources::Schedules::StandardScheduleTypeLimits::get_on_off(model))
+
+              ##### Add a "ZoneVentilation:DesignFlowRate" object for NV to set OA per person.
+              zn_vent_design_flow_rate_1 = OpenStudio::Model::ZoneVentilationDesignFlowRate.new(model)
+              zn_vent_design_flow_rate_1.setDesignFlowRateCalculationMethod('Flow/Person')
+              zn_vent_design_flow_rate_1.setFlowRateperPerson(oa_per_person_normalized_by_number_of_windows)
+              zn_vent_design_flow_rate_1.setVentilationType('Natural')
+              zn_vent_design_flow_rate_1.setMinimumIndoorTemperatureSchedule(min_Tin_schedule)
+              zn_vent_design_flow_rate_1.setMaximumIndoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_design_flow_rate_1.setMinimumOutdoorTemperature(nv_Tout_min)
+              zn_vent_design_flow_rate_1.setMaximumOutdoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_design_flow_rate_1.setDeltaTemperature(nv_Delta_Tin_Tout) #E+ I/O Ref.: "This is the temperature difference between the indoor and outdoor air dry-bulb temperatures below which ventilation is shutoff."
+              zone_hvac_equipment_list.addEquipment(zn_vent_design_flow_rate_1)
+
+              ##### Add another "ZoneVentilation:DesignFlowRate" object for NV to set OA per floor area.
+              zn_vent_design_flow_rate_2 = OpenStudio::Model::ZoneVentilationDesignFlowRate.new(model)
+              zn_vent_design_flow_rate_2.setDesignFlowRateCalculationMethod('Flow/Area')
+              zn_vent_design_flow_rate_2.setFlowRateperZoneFloorArea(oa_per_floor_area_normalized_by_number_of_windows)
+              zn_vent_design_flow_rate_2.setVentilationType('Natural')
+              zn_vent_design_flow_rate_2.setMinimumIndoorTemperatureSchedule(min_Tin_schedule)
+              zn_vent_design_flow_rate_2.setMaximumIndoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_design_flow_rate_2.setMinimumOutdoorTemperature(nv_Tout_min)
+              zn_vent_design_flow_rate_2.setMaximumOutdoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_design_flow_rate_2.setDeltaTemperature(nv_Delta_Tin_Tout)
+              zone_hvac_equipment_list.addEquipment(zn_vent_design_flow_rate_2)
+
+              ##### Add the "ZoneVentilation:WindandStackOpenArea" for NV.
+              # Note: it has been assumed that 'Opening Effectiveness' and 'Discharge Coefficient for Opening' are autocalculated (which are the default assumptions).
+              zn_vent_wind_and_stack = OpenStudio::Model::ZoneVentilationWindandStackOpenArea.new(model)
+              zn_vent_wind_and_stack.setOpeningArea(window_area * nv_opening_fraction)
+              zn_vent_wind_and_stack.setOpeningAreaFractionSchedule(operable_window_schedule)
+              # (Ref: E+ I/O) The Effective Angle value "is used to calculate the angle between the wind direction and the opening outward normal to determine the opening effectiveness values when the input field Opening Effectiveness = Autocalculate."
+              # (Ref: E+ I/O) "Effective Angle is the angle in degrees counting from the North clockwise to the opening outward normal."
+              zn_vent_wind_and_stack.setEffectiveAngle(window_azimuth_deg)
+              zn_vent_wind_and_stack.setMinimumIndoorTemperatureSchedule(min_Tin_schedule)
+              zn_vent_wind_and_stack.setMaximumIndoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_wind_and_stack.setMinimumOutdoorTemperature(nv_Tout_min)
+              zn_vent_wind_and_stack.setMaximumOutdoorTemperatureSchedule(max_Tin_schedule)
+              zn_vent_wind_and_stack.setDeltaTemperature(nv_Delta_Tin_Tout)
+              zone_hvac_equipment_list.addEquipment(zn_vent_wind_and_stack)
+
+            end #if (subsurface.subSurfaceType == 'OperableWindow' || subsurface.subSurfaceType == 'FixedWindow') && subsurface.outsideBoundaryCondition == 'Outdoors'
+          end #surface.subSurfaces.sort.each do |subsurface|
+        end #space.surfaces.sort.each do |surface|
+
+
+      end #thermal_zone.spaces.sort.each do |space|
+
+    end #model.getZoneHVACEquipmentLists.sort.each do |zone_hvac_equipment_list|
+
+    ##### Add AvailabilityManagerHybridVentilation to "prevents simultaneous natural ventilation and HVAC system operation" (Ref: E+ I/O)
+    model.getAirLoopHVACs.sort.each do |air_loop|
+      air_loop.availabilityManagers.sort.each do |avail_mgr|
+        if avail_mgr.to_AvailabilityManagerHybridVentilation.empty?
+          avail_mgr_hybr_vent = OpenStudio::Model::AvailabilityManagerHybridVentilation.new(model)
+          avail_mgr_hybr_vent.setMinimumOutdoorTemperature(nv_Tout_min) #Note: since "Ventilation Control Mode" is by default set to "Temperature (i.e. 1)", only min and max Tout are needed. (see E+ I/O Ref.)  #Note: Tout_min is to avoid overcooling (see E+ I/O Ref).
+          avail_mgr_hybr_vent.setMaximumOutdoorTemperature(30.0) #Note: the AvailabilityManagerHybridVentilation obj does not have a schedule field for Tout, so it has been set to a fixed value of 30C.
+          air_loop.addAvailabilityManager(avail_mgr_hybr_vent)
+        end
+      end
+    end
+
+  end
+
+end

data/lib/openstudio-standards/standards/necb/ECMS/pv_ground.rb

@@ -0,0 +1,104 @@
+class ECMS
+
+  def apply_pv_ground(model:, pv_ground_type:, pv_ground_total_area_pv_panels_m2:, pv_ground_tilt_angle:, pv_ground_azimuth_angle:, pv_ground_module_description:)
+
+    ##### Remove leading or trailing whitespace in case users add them in inputs
+    if pv_ground_total_area_pv_panels_m2.instance_of?(String)
+      pv_ground_total_area_pv_panels_m2 = pv_ground_total_area_pv_panels_m2.strip
+    end
+    if pv_ground_tilt_angle.instance_of?(String)
+      pv_ground_tilt_angle = pv_ground_tilt_angle.strip
+    end
+    if pv_ground_azimuth_angle.instance_of?(String)
+      pv_ground_azimuth_angle = pv_ground_azimuth_angle.strip
+    end
+
+    ##### If any of users' inputs are nil/false do nothing.
+    return if pv_ground_type.nil? || pv_ground_type == false || pv_ground_type == 'none' || pv_ground_type == 'NECB_Default'
+    return if pv_ground_total_area_pv_panels_m2 == nil? || pv_ground_total_area_pv_panels_m2 == false || pv_ground_total_area_pv_panels_m2 == 'none'
+    return if pv_ground_tilt_angle == nil? || pv_ground_tilt_angle == false || pv_ground_tilt_angle == 'none'
+    return if pv_ground_azimuth_angle == nil? || pv_ground_azimuth_angle == false || pv_ground_azimuth_angle == 'none'
+    return if pv_ground_module_description == nil? || pv_ground_module_description == false || pv_ground_module_description == 'none'
+
+    ##### Convert a string to a float (except for pv_ground_type and pv_ground_module_description)
+    if pv_ground_total_area_pv_panels_m2.instance_of?(String) && pv_ground_total_area_pv_panels_m2 != 'NECB_Default'
+      pv_ground_total_area_pv_panels_m2 = pv_ground_total_area_pv_panels_m2.to_f
+    end
+    if pv_ground_tilt_angle.instance_of?(String) && pv_ground_tilt_angle != 'NECB_Default'
+      pv_ground_tilt_angle = pv_ground_tilt_angle.to_f
+    end
+    if pv_ground_azimuth_angle.instance_of?(String) && pv_ground_azimuth_angle != 'NECB_Default'
+      pv_ground_azimuth_angle = pv_ground_azimuth_angle.to_f
+    end
+
+    ##### Calculate footprint of the building model (this is used as default value for pv_ground_total_area_pv_panels_m2)
+    building_footprint_m2_array = []
+    lowest_floor = 10000000000.0 #dummy number as initialization to find the lowest floor among spaces #TODO: Question:it it fine that it has been assumed that the floor of all lowest spaces are at the same level?
+    model.getSpaces.sort.each do |space|
+      space.surfaces.sort.select{ |surface| (surface.surfaceType == 'Floor') &&  (surface.outsideBoundaryCondition != 'Surface') &&  (surface.outsideBoundaryCondition != 'Adiabatic')}.each do |surface|
+        floor_vertices = surface.vertices
+        floor_z = floor_vertices[0].z.round(1)
+        if floor_z <= lowest_floor
+          lowest_floor = floor_z
+          building_footprint_m2_array << surface.netArea
+        end
+      end
+    end
+    building_footprint_m2 = building_footprint_m2_array.sum
+    puts "building_footprint_m2 is #{building_footprint_m2}"
+
+    ##### Set default PV panels' total area as the building footprint
+    if pv_ground_total_area_pv_panels_m2 == 'NECB_Default'
+      pv_ground_total_area_pv_panels_m2 = building_footprint_m2
+    end
+
+    ##### Set default PV panels' tilt angle as the latitude
+    if pv_ground_tilt_angle == 'NECB_Default'
+      epw = BTAP::Environment::WeatherFile.new(model.weatherFile.get.path.get)
+      pv_ground_tilt_angle = epw.latitude.to_f
+    end
+
+    ##### Set default PV panels' azimuth angle as south-facing arrays
+    if pv_ground_azimuth_angle == 'NECB_Default'
+      pv_ground_azimuth_angle = 180 # EnergyPlus I/O Reference: "An azimuth angle of 180° is for a south-facing array, and an azimuth angle of 0° is for a north-facing array."
+    end
+
+    ##### Set default PV module type as the the below one
+    if pv_ground_module_description == 'NECB_Default'
+      pv_ground_module_description = 'HES-160-36PV 26.6  x 58.3 x 1.38' # Note: As per Mike Lubun's comment, assuming a typical panel is 5 ft x 2 ft, the closest standard type PV panel in the spreadsheet would be the 160W HES.
+    end
+
+    ##### Calculate number of PV panels
+    # Note: assuming 5 ft x 2 ft as PV panel's size since it seems to fit the racking system used for ground mounts as per Mike Lubun's comment.
+    pv_area_each_ft2 = 5.0 * 2.0
+    pv_area_each_m2 = (OpenStudio.convert(pv_area_each_ft2, 'ft^2', 'm^2').get) #convert pv_area_each_ft2 to m2
+    pv_number_panels = pv_ground_total_area_pv_panels_m2/pv_area_each_m2
+
+    ##### Get data of the PV panel from the json file
+    pv_info = @standards_data['tables']['pv']['table'].detect { |item| item['pv_module_description'] == pv_ground_module_description }
+    pv_ground_module_type = pv_info['pv_module_type']
+    pv_watt = pv_info['pv_module_wattage']
+
+    ##### Create the generator
+    # Assuming one PVWatts generator in E+ as per Mike Lubun's comment for simplification, however exact number of PVWatts generators (and inverters) are calculated for costing.
+    dc_system_capacity = pv_number_panels * pv_watt
+    generator = OpenStudio::Model::GeneratorPVWatts.new(model,dc_system_capacity)
+    generator.setModuleType(pv_ground_module_type)
+    generator.setArrayType('OneAxis')   # Note: "tilt and azimuth are fixed" for this array type (see E+ I/O Reference). This array type has been chosen as per Mike Lubun's costing spec.
+    generator.setTiltAngle(pv_ground_tilt_angle)
+    generator.setAzimuthAngle(pv_ground_azimuth_angle)
+
+    ##### Create the inverter
+    inverter = OpenStudio::Model::ElectricLoadCenterInverterPVWatts.new(model)
+    inverter.setDCToACSizeRatio(1.1) # Note: This is EnergyPlus' default value; This default value has been chosen for ground-mounted PV, assuming no storage as per Mike Lubun's costing spec.
+    inverter.setInverterEfficiency(0.96) # Note: This is EnergyPlus' default value; This default value has been chosen as per Mike Lubun's costing spec.
+
+    ##### Get distribution systems and set relevant parameters
+    model.getElectricLoadCenterDistributions.sort.each  do |elc_distribution|
+      elc_distribution.setInverter(inverter)
+      elc_distribution.setGeneratorOperationSchemeType('Baseload')  # E+ I/O Reference: "The Baseload scheme requests all generators scheduled ON (available) to operate, even if the amount of electric power generated exceeds the total facility electric power demand." This scheme type has been chosen as per Mike Lubun's costing spec.
+    end
+
+  end
+
+end

data/lib/openstudio-standards/standards/necb/NECB2011/autozone.rb

@@ -52,7 +52,7 @@ class NECB2011
   # Some expections are dwelling units wet zone and wild zones.  These spaces will have special considerations when autozoning a
   # building.
 
-  def apply_auto_zoning(model:, sizing_run_dir: Dir.pwd)
+  def apply_auto_zoning(model:, sizing_run_dir: Dir.pwd, lights_type: 'NECB_Default', lights_scale: 1.0)
     raise('validation of model failed.') unless validate_initial_model(model)
 
     # Check to see if model is using another vintage of spacetypes. If so overwrite the @standards for the object with the
@@ -82,9 +82,9 @@ class NECB2011
     # Remove any Thermal zones assigned again to start fresh.
     model.getThermalZones.each(&:remove)
     self.auto_zone_dwelling_units(model)
-    self.auto_zone_wet_spaces(model)
+    self.auto_zone_wet_spaces(model: model, lights_type: lights_type, lights_scale: lights_scale)
     self.auto_zone_all_other_spaces(model)
-    self.auto_zone_wild_spaces(model)
+    self.auto_zone_wild_spaces(model: model, lights_type: lights_type, lights_scale: lights_scale)
     #THis will color the spaces and zones.
     random = Random.new(1234)
     #Set ideal hvac in case we want to not implement the hvac yet and still run osm right after this function.
@@ -193,6 +193,7 @@ class NECB2011
     # set a larger tolerance for unmet hours from default 0.2 to 1.0C
     model.getOutputControlReportingTolerances.setToleranceforTimeHeatingSetpointNotMet(1.0)
     model.getOutputControlReportingTolerances.setToleranceforTimeCoolingSetpointNotMet(1.0)
+
   end
 
   # Method to store space sizing loads. This is needed because later when the zones are destroyed this information will be lost.
@@ -283,7 +284,8 @@ class NECB2011
   # Something that the code is silent on are smelly humid areas that should not be on the same system as the rest of the
   #  building.. These are the 'wet' spaces and have been defined as locker and washroom areas.. These will be put under
   # their own single system 4 system. These will be set to the dominant floor schedule.
-  def auto_zone_wet_spaces(model)
+
+  def auto_zone_wet_spaces(model:, lights_type: 'NECB_Default', lights_scale: 1.0)
     wet_zone_array = Array.new
     model.getSpaces.select {|space| is_an_necb_wet_space?(space)}.each do |space|
       #if this space was already assigned to something skip it.
@@ -302,7 +304,8 @@ class NECB2011
 
       #this method will determine if the right schedule was used for this wet & wild space if not.. it will reset the space
       # to use the correct schedule version of the wet and wild space type.
-      adjust_wildcard_spacetype_schedule(space, dominant_schedule)
+      adjust_wildcard_spacetype_schedule(space: space, schedule: dominant_schedule, lights_type: lights_type, lights_scale: lights_scale)
+
       #Find spacetype thermostat and assign it to the zone.
       thermostat_name = space.spaceType.get.name.get + ' Thermostat'
       thermostat = model.getThermostatSetpointDualSetpointByName(thermostat_name)
@@ -320,7 +323,7 @@ class NECB2011
       model.getSpaces.select {|s| is_an_necb_wet_space?(s)}.each do |space_target|
         if space_target.thermalZone.empty?
           if are_space_loads_similar?(space_1: space, space_2: space_target) && space.buildingStory().get == space_target.buildingStory().get # added since chris needs zones to not span floors for costing.
-            adjust_wildcard_spacetype_schedule(space_target, dominant_schedule)
+            adjust_wildcard_spacetype_schedule(space_target, dominant_schedule, lights_type, lights_scale, space_height)
             space_target.setThermalZone(zone)
           end
         end
@@ -379,7 +382,7 @@ class NECB2011
 
   # This will take all the wildcard spaces and merge them to be supported by a system 4. The control zone will be the
   # zone that has the largest heating load per area.
-  def auto_zone_wild_spaces(model)
+  def auto_zone_wild_spaces(model:, lights_type: 'NECB_Default', lights_scale: 1.0)
     other_tz_array = Array.new
     #iterate through wildcard spaces.
     model.getSpaces.select {|space| is_an_necb_wildcard_space?(space) and not is_an_necb_wet_space?(space)}.each do |space|
@@ -396,10 +399,13 @@ class NECB2011
       #Assign space to the new zone.
       space.setThermalZone(zone)
 
-      #lets keep the wild schedules to be the same as what dominate the floor.
+      # lets keep the wild schedules to be the same as what dominate the floor.
       dominant_floor_schedule = determine_dominant_schedule(space.model.getSpaces)
 
-      adjust_wildcard_spacetype_schedule(space, dominant_floor_schedule)
+      adjust_wildcard_spacetype_schedule(space: space,
+                                         schedule: dominant_floor_schedule,
+                                         lights_type: lights_type,
+                                         lights_scale: lights_scale)
 
       # Add a thermostat
       space_type_name = space.spaceType.get.name.get
@@ -464,7 +470,7 @@ class NECB2011
         schedule_type = determine_dominant_schedule(space.buildingStory.get.spaces)
         zone = other_adjacent_spaces.first.thermalZone.get
         wild_adjacent_spaces.each do |space|
-          adjust_wildcard_spacetype_schedule(space, schedule_type)
+          adjust_wildcard_spacetype_schedule(space, schedule_type, @lights_type, @lights_scale, @space_height)
           space.setThermalZone(zone)
         end
       end
@@ -483,7 +489,7 @@ class NECB2011
       dominant_floor_schedule = determine_dominant_schedule(space.buildingStory().get.spaces)
       #this method will determine if the right schedule was used for this wet & wild space if not.. it will reset the space
       # to use the correct schedule version of the wet and wild space type.
-      adjust_wildcard_spacetype_schedule(space, dominant_floor_schedule)
+      adjust_wildcard_spacetype_schedule(space, dominant_floor_schedule, @lights_type, @lights_scale, @space_height)
       #Find spacetype thermostat and assign it to the zone.
       thermostat_name = space.spaceType.get.name.get + ' Thermostat'
       thermostat = model.getThermostatSetpointDualSetpointByName(thermostat_name)
@@ -502,7 +508,7 @@ class NECB2011
         if space_target.thermalZone.empty?
           if are_space_loads_similar?(space_1: space, space_2: space_target) &&
               space.buildingStory().get == space_target.buildingStory().get # added since chris needs zones to not span floors for costing.
-            adjust_wildcard_spacetype_schedule(space_target, dominant_floor_schedule)
+            adjust_wildcard_spacetype_schedule(space_target, dominant_floor_schedule, @lights_type, @lights_scale, @space_height)
             space_target.setThermalZone(zone)
           end
         end
@@ -708,7 +714,7 @@ class NECB2011
   end
 
   # Set wildcard spactype schedule to NECB letter index.
-  def adjust_wildcard_spacetype_schedule(space, schedule)
+  def adjust_wildcard_spacetype_schedule(space:, schedule:, lights_type: 'NECB_Default', lights_scale: 1.0)
     if space.spaceType.empty?
       OpenStudio.logFree(OpenStudio::Error, 'Error: No spacetype assigned for #{space.name.get}. This must be assigned. Aborting.')
     end
@@ -733,7 +739,7 @@ class NECB2011
         new_spacetype.setStandardsBuildingType(space.spaceType.get.standardsBuildingType.get)
         new_spacetype.setStandardsSpaceType(new_spacetype_name)
         new_spacetype.setName("#{space.spaceType.get.standardsBuildingType.get} #{new_spacetype_name}")
-        space_type_apply_internal_loads(new_spacetype, true, true, true, true, true, true)
+        space_type_apply_internal_loads(space_type: new_spacetype, lights_type: lights_type, lights_scale: lights_scale)
         space_type_apply_internal_load_schedules(new_spacetype, true, true, true, true, true, true, true)
       end
       space.setSpaceType(new_spacetype)
@@ -881,7 +887,8 @@ class NECB2011
                                                                mau_type: mau_type,
                                                                mau_heating_coil_type: mau_heating_coil_type,
                                                                baseboard_type: baseboard_type,
-                                                               hw_loop: @hw_loop)
+                                                               hw_loop: @hw_loop,
+                                                               multispeed: false)
         end
       when 2
         group_similar_zones_together(zones).each do |zones|
@@ -894,11 +901,12 @@ class NECB2011
         end
       when 3
         group_similar_zones_together(zones).each do |zones|
-          add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating_single_speed(model: model,
-                                                                                             zones: zones,
-                                                                                             heating_coil_type: heating_coil_type_sys3,
-                                                                                             baseboard_type: baseboard_type,
-                                                                                             hw_loop: @hw_loop)
+          add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+                                                                                zones: zones,
+                                                                                heating_coil_type: heating_coil_type_sys3,
+                                                                                baseboard_type: baseboard_type,
+                                                                                hw_loop: @hw_loop,
+                                                                                multispeed: false)
         end
       when 4
         group_similar_zones_together(zones).each do |zones|
@@ -907,6 +915,12 @@ class NECB2011
                                                                        heating_coil_type: heating_coil_type_sys4,
                                                                        baseboard_type: baseboard_type,
                                                                        hw_loop: @hw_loop)
+#          add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+#                                                                                zones: zones,
+#                                                                                heating_coil_type: heating_coil_type_sys4,
+#                                                                                baseboard_type: baseboard_type,
+#                                                                                hw_loop: @hw_loop,
+#                                                                                multispeed: false)
         end
       when 5
         group_similar_zones_together(zones).each do |zones|
@@ -1024,7 +1038,8 @@ class NECB2011
                                                 mau_type: mau_type,
                                                 mau_heating_coil_type: mau_heating_coil_type,
                                                 baseboard_type: baseboard_type,
-                                                hw_loop: @hw_loop)
+                                                hw_loop: @hw_loop,
+                                                multispeed: false)
         else
           #Create a separate air loop for each unit.
           zones.each do |zone|
@@ -1033,27 +1048,29 @@ class NECB2011
                                                   mau_type: mau_type,
                                                   mau_heating_coil_type: mau_heating_coil_type,
                                                   baseboard_type: baseboard_type,
-                                                  hw_loop: @hw_loop)
+                                                  hw_loop: @hw_loop,
+                                                  multispeed: false)
 
           end
         end
 
       when 3
         if dwelling_shared_ahu
-          add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating_single_speed(model: model,
-                                                                                             zones: zones,
-                                                                                             heating_coil_type: heating_coil_type_sys3,
-                                                                                             baseboard_type: baseboard_type,
-                                                                                             hw_loop: @hw_loop)
+          add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+                                                                                zones: zones,
+                                                                                heating_coil_type: heating_coil_type_sys3,
+                                                                                baseboard_type: baseboard_type,
+                                                                                hw_loop: @hw_loop,
+                                                                                multispeed: false)
         else
           #Create a separate air loop for each unit.
           zones.each do |zone|
-            add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating_single_speed(model: model,
-                                                                                               zones: [zone],
-                                                                                               heating_coil_type: heating_coil_type_sys3,
-                                                                                               baseboard_type: baseboard_type,
-                                                                                               hw_loop: @hw_loop)
-
+            add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+                                                                                  zones: [zone],
+                                                                                  heating_coil_type: heating_coil_type_sys3,
+                                                                                  baseboard_type: baseboard_type,
+                                                                                  hw_loop: @hw_loop,
+                                                                                  multispeed: false)
           end
         end
       end
@@ -1079,6 +1096,12 @@ class NECB2011
                                                                    heating_coil_type: heating_coil_type_sys4,
                                                                    baseboard_type: baseboard_type,
                                                                    hw_loop: @hw_loop)
+#      add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+#                                                                            zones: wet_tz,
+#                                                                            heating_coil_type: heating_coil_type_sys4,
+#                                                                            baseboard_type: baseboard_type,
+#                                                                            hw_loop: @hw_loop,
+#                                                                            multispeed: false)
     end
   end
 
@@ -1102,6 +1125,12 @@ class NECB2011
                                                                    heating_coil_type: heating_coil_type_sys4,
                                                                    baseboard_type: baseboard_type,
                                                                    hw_loop: @hw_loop)
+#      add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+#                                                                            zones: tz,
+#                                                                            heating_coil_type: heating_coil_type_sys4,
+#                                                                            baseboard_type: baseboard_type,
+#                                                                            hw_loop: @hw_loop,
+#                                                                            multispeed: true)
     end
   end
 
@@ -1125,6 +1154,12 @@ class NECB2011
                                                                    heating_coil_type: heating_coil_type_sys4,
                                                                    baseboard_type: baseboard_type,
                                                                    hw_loop: @hw_loop)
+#      add_sys3and8_single_zone_packaged_rooftop_unit_with_baseboard_heating(model: model,
+#                                                                            zones: zones,
+#                                                                            heating_coil_type: heating_coil_type_sys4,
+#                                                                            baseboard_type: baseboard_type,
+#                                                                            hw_loop: @hw_loop,
+#                                                                            multispeed: true)
     end
   end