openstudio-standards 0.2.12.rc4 → 0.2.13.rc3

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2010/ashrae_90_1_2010.PumpVariableSpeed.rb

@@ -0,0 +1,29 @@
+class ASHRAE9012010 < ASHRAE901
+  # @!group Pump
+
+  # Determine type of pump part load control type
+  #
+  # @code_sections [90.1-2010_6.5.4.1]
+  # @param pump [OpenStudio::Model::PumpVariableSpeed] OpenStudio pump object
+  # @param plant_loop_type [String] Type of plant loop
+  # @param pump_nominal_hp [Float] Pump nominal horsepower
+  # @return [String] Pump part load control type
+  def pump_variable_speed_get_control_type(pump, plant_loop_type, pump_nominal_hp)
+    threshold = 5 # hp
+
+    # Sizing factor to take into account that pumps
+    # are typically sized to handle a ~10% pressure
+    # increase and ~10% flow increase.
+    design_sizing_factor = 1.25
+
+    if plant_loop_type == 'Heating'
+      return 'Riding Curve'
+    else
+      if pump_nominal_hp * design_sizing_factor > threshold
+        return 'VSD DP Reset' # Requirement only applies to CHW pumps
+      else
+        return 'Riding Curve'
+      end
+    end
+  end
+  end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2013/ashrae_90_1_2013.Model.rb

@@ -16,10 +16,12 @@ class ASHRAE9012013 < ASHRAE901
   # 'DifferentialDryBulbAndEnthalpy'
   def model_economizer_type(model, climate_zone)
     economizer_type = case climate_zone
-                      when 'ASHRAE 169-2006-1A',
+                      when 'ASHRAE 169-2006-0A',
+                          'ASHRAE 169-2006-1A',
                           'ASHRAE 169-2006-2A',
                           'ASHRAE 169-2006-3A',
                           'ASHRAE 169-2006-4A',
+                          'ASHRAE 169-2013-0A',
                           'ASHRAE 169-2013-1A',
                           'ASHRAE 169-2013-2A',
                           'ASHRAE 169-2013-3A',
@@ -30,4 +32,157 @@ class ASHRAE9012013 < ASHRAE901
                       end
     return economizer_type
   end
-end
+
+  # Adjust model to comply with fenestration orientation requirements
+  #
+  # @code_sections [90.1-2013_5.5.4.5]
+  # @param [OpenStudio::Model::Model] OpenStudio model object
+  # @return [Boolean] Returns true if successful, false otherwise
+  def model_fenestration_orientation(model, climate_zone)
+    # Building rotation to meet the same code requirement for
+    # 90.1-2010 are kept
+    if model.getBuilding.standardsBuildingType.is_initialized
+      building_type = model.getBuilding.standardsBuildingType.get
+
+      case building_type
+        when 'Hospital'
+          # Rotate the building counter-clockwise
+          model_set_building_north_axis(model, 270.0)
+        when 'SmallHotel'
+          # Rotate the building clockwise
+          model_set_building_north_axis(model, 180)
+      end
+    end
+
+    wwr = false
+    # Section 6.2.1.2 in the ANSI/ASHRAE/IES Standard 90.1-2013 Determination
+    # of Energy Savings: Quantitative Analysis mentions that the SHGC trade-off
+    # path is most likely to be used by designers for compliance.
+    #
+    # The following adjustment are only made for models with simple glazing objects
+    non_simple_glazing = false
+    shgc_a = 0
+    model.getSpaces.each do |space|
+      # Get thermal zone multiplier
+      multiplier = space.thermalZone.get.multiplier
+
+      space.surfaces.each do |surface|
+        surface.subSurfaces.each do |subsurface|
+          # Get window subsurface type
+          subsurface_type = subsurface.subSurfaceType.to_s.downcase
+
+          # Window, glass doors
+          next unless (subsurface_type.include? 'window') || (subsurface_type.include? 'glass')
+
+          # Check if non simple glazing fenestration objects are used
+          subsurface_cons = subsurface.construction.get.to_Construction.get
+          non_simple_glazing = true unless subsurface_cons.layers[0].to_SimpleGlazing.is_initialized
+
+          if non_simple_glazing
+            OpenStudio.logFree(OpenStudio::Warn, 'openstudio.ashrae_90_1_2013.model', 'Fenestration objects in the model use non-simple glazing models, fenestration requirements are not applied')
+            return false
+          end
+
+          # Get subsurface's simple glazing object
+          subsurface_shgc = subsurface_cons.layers[0].to_SimpleGlazing.get.solarHeatGainCoefficient
+
+          # Get subsurface area
+          subsurface_area = subsurface.grossArea * subsurface.multiplier * multiplier
+
+          # SHGC * Area
+          shgc_a += subsurface_shgc * subsurface_area
+        end
+      end
+    end
+
+    # Calculate West, East and total fenestration area
+    a_w = model_get_window_area_info_for_orientation(model, 'W', wwr: wwr)
+    a_e = model_get_window_area_info_for_orientation(model, 'E', wwr: wwr)
+    a_t = a_w + a_e + model_get_window_area_info_for_orientation(model, 'N', wwr: wwr) + model_get_window_area_info_for_orientation(model, 'S', wwr: wwr)
+
+    return true if a_t == 0.0
+
+    # For prototypes SHGC_c assumed to be the building's weighted average SHGC
+    shgc_c = shgc_a / a_t
+    shgc_c = shgc_c.round(2)
+
+    # West and East facing WWR
+    wwr_w = model_get_window_area_info_for_orientation(model, 'W', wwr: true)
+    wwr_e = model_get_window_area_info_for_orientation(model, 'E', wwr: true)
+
+    # Calculate new SHGC for west and east facing fenestration;
+    # Create new simple glazing object and assign it to all
+    # West and East fenestration
+    #
+    # Exception 5 is applied when applicable
+    shgc_w = 0
+    shgc_e = 0
+    if !((a_w <= a_t / 4) && (a_e <= a_t / 4))
+      # Calculate new SHGC
+      if wwr_w > 0.2
+        shgc_w = a_t * shgc_c / (4 * a_w)
+      end
+      if wwr_e > 0.2
+        shgc_e = a_t * shgc_c / (4 * a_w)
+      end
+
+      # No SHGC adjustment needed
+      return true if shgc_w == 0 && shgc_e == 0
+
+      model.getSpaces.each do |space|
+        # Get thermal zone multiplier
+        multiplier = space.thermalZone.get.multiplier
+
+        space.surfaces.each do |surface|
+          # Proceed only for East and West facing surfaces that are required
+          # to have their SHGC adjusted
+          next unless (surface_cardinal_direction(surface) == 'W' && shgc_w > 0) ||
+                      (surface_cardinal_direction(surface) == 'E' && shgc_e > 0)
+
+          surface.subSurfaces.each do |subsurface|
+            # Get window subsurface type
+            subsurface_type = subsurface.subSurfaceType.to_s.downcase
+
+            # Window, glass doors
+            next unless (subsurface_type.include? 'window') || (subsurface_type.include? 'glass')
+
+            new_shgc = surface_cardinal_direction(surface) == 'W' ? shgc_w : shgc_e
+            new_shgc = new_shgc.round(2)
+
+            # Get construction/simple glazing associated with the subsurface
+            subsurface_org_cons = subsurface.construction.get.to_Construction.get
+            subsurface_org_cons_mat = subsurface_org_cons.layers[0].to_SimpleGlazing.get
+
+            # Only proceed if new SHGC is different than orignal one
+            next unless (new_shgc - subsurface_org_cons_mat.solarHeatGainCoefficient).abs > 0
+
+            # Clone construction/simple glazing associated with the subsurface
+            subsurface_new_cons = subsurface_org_cons.clone(model).to_Construction.get
+            subsurface_new_cons.setName("#{subsurface.name} Wind Cons U-#{OpenStudio.convert(subsurface_org_cons_mat.uFactor, 'W/m^2*K', 'Btu/ft^2*h*R').get.round(2)} SHGC #{new_shgc}")
+            subsurface_new_cons_mat = subsurface_org_cons_mat.clone(model).to_SimpleGlazing.get
+            subsurface_new_cons_mat.setName("#{subsurface.name} Wind SG Mat U-#{OpenStudio.convert(subsurface_org_cons_mat.uFactor, 'W/m^2*K', 'Btu/ft^2*h*R').get.round(2)} SHGC #{new_shgc}")
+            subsurface_new_cons_mat.setSolarHeatGainCoefficient(new_shgc)
+            new_layers = OpenStudio::Model::MaterialVector.new
+            new_layers << subsurface_new_cons_mat
+            subsurface_new_cons.setLayers(new_layers)
+
+            # Assign new construction to sub surface
+            subsurface.setConstruction(subsurface_new_cons)
+          end
+        end
+      end
+    end
+
+    return true
+  end
+
+  # Is transfer air required?
+  #
+  # @code_sections [90.1-2013_6.5.7.1.2]
+  # @param model [OpenStudio::Model::Model] OpenStudio model object
+  # @return [Boolean] true if transfer air is required, false otherwise
+  def model_transfer_air_required?(model)
+    # TODO: It actually is for kitchen but not implemented yet
+    return false
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2013/ashrae_90_1_2013.PumpVariableSpeed.rb

@@ -0,0 +1,29 @@
+class ASHRAE9012013 < ASHRAE901
+  # @!group Pump
+
+  # Determine type of pump part load control type
+  #
+  # @code_sections [90.1-2013_6.5.4.2]
+  # @param pump [OpenStudio::Model::PumpVariableSpeed] OpenStudio pump object
+  # @param plant_loop_type [String] Type of plant loop
+  # @param pump_nominal_hp [Float] Pump nominal horsepower
+  # @return [String] Pump part load control type
+  def pump_variable_speed_get_control_type(pump, plant_loop_type, pump_nominal_hp)
+    threshold = 5 # hp
+
+    # Sizing factor to take into account that pumps
+    # are typically sized to handle a ~10% pressure
+    # increase and ~10% flow increase.
+    design_sizing_factor = 1.25
+
+    if plant_loop_type == 'Heating'
+      return 'Riding Curve'
+    else
+      if pump_nominal_hp * design_sizing_factor > threshold
+        return 'VSD DP Reset' # Requirement only applies to CHW pumps
+      else
+        return 'Riding Curve'
+      end
+    end
+  end
+  end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.AirTerminalSingleDuctVAVReheat.rb

@@ -0,0 +1,24 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group AirTerminalSingleDuctVAVReheat
+  # Set the initial minimum damper position based on OA rate of the space and the template.
+  # Zones with low OA per area get lower initial guesses.
+  # Final position will be adjusted upward as necessary by Standards.AirLoopHVAC.apply_minimum_vav_damper_positions
+  #
+  # @param zone_oa_per_area [Double] the zone outdoor air per area, m^3/s
+  # @return [Bool] returns true if successful, false if not
+  def air_terminal_single_duct_vav_reheat_apply_initial_prototype_damper_position(air_terminal_single_duct_vav_reheat, zone_oa_per_area)
+    min_damper_position = case air_terminal_single_duct_vav_reheat_reheat_type(air_terminal_single_duct_vav_reheat)
+                          when 'HotWater'
+                            0.2
+                          when 'Electricity', 'NaturalGas'
+                            0.3
+                          else
+                            0.2
+                          end
+
+    # Set the minimum flow fraction
+    air_terminal_single_duct_vav_reheat.setConstantMinimumAirFlowFraction(min_damper_position)
+
+    return true
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.FanConstantVolume.rb

@@ -0,0 +1,32 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group FanConstantVolume
+
+  # Determine the prototype fan pressure rise for a constant volume
+  # fan on an AirLoopHVAC based on the airflow of the system.
+  # @return [Double] the pressure rise (in H2O).  Defaults
+  # to the logic from ASHRAE 90.1-2004 prototypes.
+  def fan_constant_volume_airloop_fan_pressure_rise(fan_constant_volume)
+    # Get the max flow rate from the fan.
+    maximum_flow_rate_m3_per_s = nil
+    if fan_constant_volume.maximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_constant_volume.maximumFlowRate.get
+    elsif fan_constant_volume.autosizedMaximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_constant_volume.autosizedMaximumFlowRate.get
+    else
+      OpenStudio.logFree(OpenStudio::Warn, 'openstudio.prototype.FanConstantVolume', "For #{fan_constant_volume.name} max flow rate is not available, cannot apply prototype assumptions.")
+      return false
+    end
+
+    # Convert max flow rate to cfm
+    maximum_flow_rate_cfm = OpenStudio.convert(maximum_flow_rate_m3_per_s, 'm^3/s', 'cfm').get
+
+    # Determine the pressure rise
+    pressure_rise_in_h2o = if maximum_flow_rate_cfm < 7437
+                             2.5
+                           else # Over 7,437 cfm
+                             4.09
+                           end
+
+    return pressure_rise_in_h2o
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.FanOnOff.rb

@@ -0,0 +1,33 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group FanOnOff
+
+  # Determine the prototype fan pressure rise for an on off
+  # fan on an AirLoopHVAC or inside a unitary system
+  # based on the airflow of the system.
+  # @return [Double] the pressure rise (in H2O).  Defaults
+  # to the logic from ASHRAE 90.1-2004 prototypes.
+  def fan_on_off_airloop_or_unitary_fan_pressure_rise(fan_on_off)
+    # Get the max flow rate from the fan.
+    maximum_flow_rate_m3_per_s = nil
+    if fan_on_off.maximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_on_off.maximumFlowRate.get
+    elsif fan_on_off.autosizedMaximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_on_off.autosizedMaximumFlowRate.get
+    else
+      OpenStudio.logFree(OpenStudio::Warn, 'openstudio.prototype.FanOnOff', "For #{fan_on_off.name} max flow rate is not available, cannot apply prototype assumptions.")
+      return false
+    end
+
+    # Convert max flow rate to cfm
+    maximum_flow_rate_cfm = OpenStudio.convert(maximum_flow_rate_m3_per_s, 'm^3/s', 'cfm').get
+
+    # Determine the pressure rise
+    pressure_rise_in_h2o = if maximum_flow_rate_cfm < 7437
+                             2.5
+                           else # Over 7,437 cfm
+                             4.09
+                           end
+
+    return pressure_rise_in_h2o
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.FanVariableVolume.rb

@@ -0,0 +1,32 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group FanVariableVolume
+
+  # Determine the prototype fan pressure rise for a variable volume
+  # fan on an AirLoopHVAC based on the airflow of the system.
+  # @return [Double] the pressure rise (in H2O).  Defaults
+  # to the logic from ASHRAE 90.1-2004 prototypes.
+  def fan_variable_volume_airloop_fan_pressure_rise(fan_variable_volume)
+    # Get the max flow rate from the fan.
+    maximum_flow_rate_m3_per_s = nil
+    if fan_variable_volume.maximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_variable_volume.maximumFlowRate.get
+    elsif fan_variable_volume.autosizedMaximumFlowRate.is_initialized
+      maximum_flow_rate_m3_per_s = fan_variable_volume.autosizedMaximumFlowRate.get
+    else
+      OpenStudio.logFree(OpenStudio::Warn, 'openstudio.prototype.FanVariableVolume', "For #{fan_variable_volume.name} max flow rate is not available, cannot apply prototype assumptions.")
+      return false
+    end
+
+    # Convert max flow rate to cfm
+    maximum_flow_rate_cfm = OpenStudio.convert(maximum_flow_rate_m3_per_s, 'm^3/s', 'cfm').get
+
+    # Determine the pressure rise
+    pressure_rise_in_h2o = if maximum_flow_rate_cfm < 4648
+                             4.0
+                           else # Over 7,437 cfm
+                             5.58
+                           end
+
+    return pressure_rise_in_h2o
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.Model.elevators.rb

@@ -0,0 +1,25 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group elevators
+
+  # Determines the percentage of the elevator cab
+  # lighting that is incandescent.  The remainder
+  # is assumed to be LED.  Defaults to 0% incandescent
+  # (100% LED), representing newer elevators.
+  def model_elevator_lighting_pct_incandescent(model)
+    pct_incandescent = 0.0 # 100% LED
+    return pct_incandescent
+  end
+
+  # Determines the power of the elevator ventilation fan.
+  # 90.1-2016 has a requirement
+  # for ventilation fan efficiency.
+  # @return [Double] the ventilaton fan power (W)
+  def model_elevator_fan_pwr(model, vent_rate_cfm)
+    vent_pwr_per_flow_w_per_cfm = 0.33
+    vent_pwr_w = vent_pwr_per_flow_w_per_cfm * vent_rate_cfm
+    # addendum 90.1-2007 aj has requirement on efficiency
+    vent_pwr_w = vent_pwr_w * 0.29 / 0.70
+
+    return vent_pwr_w
+  end
+end

data/lib/openstudio-standards/prototypes/ashrae_90_1/ashrae_90_1_2016/ashrae_90_1_2016.Model.rb

@@ -0,0 +1,366 @@
+class ASHRAE9012016 < ASHRAE901
+  # @!group Model
+
+  # Determine the prototypical economizer type for the model.
+  #
+  # @param model [OpenStudio::Model::Model] the model
+  # @param climate_zone [String] the climate zone
+  # @return [String] the economizer type.  Possible values are:
+  # 'NoEconomizer'
+  # 'FixedDryBulb'
+  # 'FixedEnthalpy'
+  # 'DifferentialDryBulb'
+  # 'DifferentialEnthalpy'
+  # 'FixedDewPointAndDryBulb'
+  # 'ElectronicEnthalpy'
+  # 'DifferentialDryBulbAndEnthalpy'
+  def model_economizer_type(model, climate_zone)
+    economizer_type = case climate_zone
+                      when 'ASHRAE 169-2006-0A',
+                          'ASHRAE 169-2006-1A',
+                          'ASHRAE 169-2006-2A',
+                          'ASHRAE 169-2006-3A',
+                          'ASHRAE 169-2006-4A',
+                          'ASHRAE 169-2013-0A',
+                          'ASHRAE 169-2013-1A',
+                          'ASHRAE 169-2013-2A',
+                          'ASHRAE 169-2013-3A',
+                          'ASHRAE 169-2013-4A'
+                        'DifferentialEnthalpy'
+                      else
+                        'DifferentialDryBulb'
+                      end
+    return economizer_type
+  end
+
+  # Adjust model to comply with fenestration orientation requirements
+  #
+  # @code_sections [90.1-2013_5.5.4.5]
+  # @param [OpenStudio::Model::Model] OpenStudio model object
+  # @return [Boolean] Returns true if successful, false otherwise
+  def model_fenestration_orientation(model, climate_zone)
+    # Building rotation to meet the same code requirement for
+    # 90.1-2010 are kept
+    if model.getBuilding.standardsBuildingType.is_initialized
+      building_type = model.getBuilding.standardsBuildingType.get
+
+      case building_type
+        when 'Hospital'
+          # Rotate the building counter-clockwise
+          model_set_building_north_axis(model, 270.0)
+        when 'SmallHotel'
+          # Rotate the building clockwise
+          model_set_building_north_axis(model, 180)
+      end
+    end
+
+    wwr = false
+    # Section 6.2.1.2 in the ANSI/ASHRAE/IES Standard 90.1-2013 Determination
+    # of Energy Savings: Quantitative Analysis mentions that the SHGC trade-off
+    # path is most likely to be used by designers for compliance.
+    #
+    # The following adjustment are only made for models with simple glazing objects
+    non_simple_glazing = false
+    shgc_a = 0
+    model.getSpaces.each do |space|
+      # Get thermal zone multiplier
+      multiplier = space.thermalZone.get.multiplier
+
+      space.surfaces.each do |surface|
+        surface.subSurfaces.each do |subsurface|
+          # Get window subsurface type
+          subsurface_type = subsurface.subSurfaceType.to_s.downcase
+
+          # Window, glass doors
+          next unless (subsurface_type.include? 'window') || (subsurface_type.include? 'glass')
+
+          # Check if non simple glazing fenestration objects are used
+          subsurface_cons = subsurface.construction.get.to_Construction.get
+          non_simple_glazing = true unless subsurface_cons.layers[0].to_SimpleGlazing.is_initialized
+
+          if non_simple_glazing
+            OpenStudio.logFree(OpenStudio::Warn, 'openstudio.ashrae_90_1_2013.model', 'Fenestration objects in the model use non-simple glazing models, fenestration requirements are not applied')
+            return false
+          end
+
+          # Get subsurface's simple glazing object
+          subsurface_shgc = subsurface_cons.layers[0].to_SimpleGlazing.get.solarHeatGainCoefficient
+
+          # Get subsurface area
+          subsurface_area = subsurface.grossArea * subsurface.multiplier * multiplier
+
+          # SHGC * Area
+          shgc_a += subsurface_shgc * subsurface_area
+        end
+      end
+    end
+
+    # Calculate West, East and total fenestration area
+    a_w = model_get_window_area_info_for_orientation(model, 'W', wwr: wwr)
+    a_e = model_get_window_area_info_for_orientation(model, 'E', wwr: wwr)
+    a_t = a_w + a_e + model_get_window_area_info_for_orientation(model, 'N', wwr: wwr) + model_get_window_area_info_for_orientation(model, 'S', wwr: wwr)
+
+    return true if a_t == 0.0
+
+    # For prototypes SHGC_c assumed to be the building's weighted average SHGC
+    shgc_c = shgc_a / a_t
+    shgc_c = shgc_c.round(2)
+
+    # West and East facing WWR
+    wwr_w = model_get_window_area_info_for_orientation(model, 'W', wwr: true)
+    wwr_e = model_get_window_area_info_for_orientation(model, 'E', wwr: true)
+
+    # Calculate new SHGC for west and east facing fenestration;
+    # Create new simple glazing object and assign it to all
+    # West and East fenestration
+    #
+    # Exception 5 is applied when applicable
+    shgc_w = 0
+    shgc_e = 0
+
+    # Determine requirement criteria
+    case climate_zone
+      when 'ASHRAE 169-2006-0A',
+           'ASHRAE 169-2006-0B',
+           'ASHRAE 169-2006-1A',
+           'ASHRAE 169-2006-1B',
+           'ASHRAE 169-2006-2A',
+           'ASHRAE 169-2006-2B',
+           'ASHRAE 169-2006-3A',
+           'ASHRAE 169-2006-3B',
+           'ASHRAE 169-2006-3C',
+           'ASHRAE 169-2013-0A',
+           'ASHRAE 169-2013-0B',
+           'ASHRAE 169-2013-1A',
+           'ASHRAE 169-2013-1B',
+           'ASHRAE 169-2013-2A',
+           'ASHRAE 169-2013-2B',
+           'ASHRAE 169-2013-3A',
+           'ASHRAE 169-2013-3B',
+           'ASHRAE 169-2013-3C'
+        criteria = 4
+      when 'ASHRAE 169-2006-4A',
+           'ASHRAE 169-2006-4B',
+           'ASHRAE 169-2006-4C',
+           'ASHRAE 169-2006-5A',
+           'ASHRAE 169-2006-5B',
+           'ASHRAE 169-2006-5C',
+           'ASHRAE 169-2006-6A',
+           'ASHRAE 169-2006-6B',
+           'ASHRAE 169-2006-7A',
+           'ASHRAE 169-2006-7B',
+           'ASHRAE 169-2006-8A',
+           'ASHRAE 169-2006-8B',
+           'ASHRAE 169-2013-4A',
+           'ASHRAE 169-2013-4B',
+           'ASHRAE 169-2013-4C',
+           'ASHRAE 169-2013-5A',
+           'ASHRAE 169-2013-5B',
+           'ASHRAE 169-2013-5C',
+           'ASHRAE 169-2013-6A',
+           'ASHRAE 169-2013-6B',
+           'ASHRAE 169-2013-7A',
+           'ASHRAE 169-2013-7B',
+           'ASHRAE 169-2013-8A',
+           'ASHRAE 169-2013-8B'
+        criteria = 5
+      else
+        return false
+    end
+
+    if !((a_w <= a_t / criteria) && (a_e <= a_t / criteria))
+      # Calculate new SHGC
+      if wwr_w > 0.2
+        shgc_w = a_t * shgc_c / (criteria * a_w)
+      end
+      if wwr_e > 0.2
+        shgc_e = a_t * shgc_c / (criteria * a_w)
+      end
+
+      # No SHGC adjustment needed
+      return true if shgc_w == 0 && shgc_e == 0
+
+      model.getSpaces.each do |space|
+        # Get thermal zone multiplier
+        multiplier = space.thermalZone.get.multiplier
+
+        space.surfaces.each do |surface|
+          # Proceed only for East and West facing surfaces that are required
+          # to have their SHGC adjusted
+          next unless (surface_cardinal_direction(surface) == 'W' && shgc_w > 0) ||
+                      (surface_cardinal_direction(surface) == 'E' && shgc_e > 0)
+
+          surface.subSurfaces.each do |subsurface|
+            # Get window subsurface type
+            subsurface_type = subsurface.subSurfaceType.to_s.downcase
+
+            # Window, glass doors
+            next unless (subsurface_type.include? 'window') || (subsurface_type.include? 'glass')
+
+            new_shgc = surface_cardinal_direction(surface) == 'W' ? shgc_w : shgc_e
+            new_shgc = new_shgc.round(2)
+
+            # Get construction/simple glazing associated with the subsurface
+            subsurface_org_cons = subsurface.construction.get.to_Construction.get
+            subsurface_org_cons_mat = subsurface_org_cons.layers[0].to_SimpleGlazing.get
+
+            # Only proceed if new SHGC is different than orignal one
+            next unless (new_shgc - subsurface_org_cons_mat.solarHeatGainCoefficient).abs > 0
+
+            # Clone construction/simple glazing associated with the subsurface
+            subsurface_new_cons = subsurface_org_cons.clone(model).to_Construction.get
+            subsurface_new_cons.setName("#{subsurface.name} Wind Cons U-#{OpenStudio.convert(subsurface_org_cons_mat.uFactor, 'W/m^2*K', 'Btu/ft^2*h*R').get.round(2)} SHGC #{new_shgc}")
+            subsurface_new_cons_mat = subsurface_org_cons_mat.clone(model).to_SimpleGlazing.get
+            subsurface_new_cons_mat.setName("#{subsurface.name} Wind SG Mat U-#{OpenStudio.convert(subsurface_org_cons_mat.uFactor, 'W/m^2*K', 'Btu/ft^2*h*R').get.round(2)} SHGC #{new_shgc}")
+            subsurface_new_cons_mat.setSolarHeatGainCoefficient(new_shgc)
+            new_layers = OpenStudio::Model::MaterialVector.new
+            new_layers << subsurface_new_cons_mat
+            subsurface_new_cons.setLayers(new_layers)
+
+            # Assign new construction to sub surface
+            subsurface.setConstruction(subsurface_new_cons)
+          end
+        end
+      end
+    end
+
+    return true
+  end
+
+  # Is transfer air required?
+  #
+  # @code_sections [90.1-2016_6.5.7.1]
+  # @param model [OpenStudio::Model::Model] OpenStudio model object
+  # @return [Boolean] true if transfer air is required, false otherwise
+  def model_transfer_air_required?(model)
+    return true
+  end
+
+  # Metal coiling door code minimum infiltration rate at 75 Pa
+  #
+  # @code_sections [90.1-2019_5.4.3.2]
+  # @param [String] Climate zone
+  # @return [Float] Minimum infiltration rate for metal coiling doors
+  def model_door_infil_flow_rate_metal_coiling_cfm_ft2(climate_zone)
+    case climate_zone
+      when 'ASHRAE 169-2006-7A',
+           'ASHRAE 169-2006-7B',
+           'ASHRAE 169-2006-8A',
+           'ASHRAE 169-2006-8B',
+           'ASHRAE 169-2013-7A',
+           'ASHRAE 169-2013-7B',
+           'ASHRAE 169-2013-8A',
+           'ASHRAE 169-2013-8B'
+        return 0.4
+      else
+        return 1.0
+    end
+  end
+
+  # Implement occupancy based lighting level threshold (0.02 W/sqft). This is only for ASHRAE 90.1 2016 onwards.
+  #
+  # @code_sections [90.1-2016_9.4.1.1.h/i]
+  # @author Xuechen (Jerry) Lei, PNNL
+  # @param model [OpenStudio::Model::Model] OpenStudio Model
+  #
+  def model_add_lights_shutoff(model)
+    zones = model.getThermalZones
+    num_zones = 0
+    business_sch_name = prototype_input['business_schedule']
+    return if business_sch_name.nil? # This is only for 10 prototypes that do not have continuous operation.
+
+    # Add business schedule
+    model_add_schedule(model, business_sch_name)
+
+    # Add EMS object for business schedule variable
+    business_sensor = OpenStudio::Model::EnergyManagementSystemSensor.new(model, 'Schedule Value')
+    business_sensor.setKeyName(business_sch_name)
+    business_sensor.setName('Business_Sensor')
+    business_sensor_name = business_sensor.name.to_s
+
+    zones.each do |zone|
+      spaces = zone.spaces
+      if spaces.length != 1
+        puts 'warning, there are more than one spaces in the zone, need to confirm the implementation'
+      end
+      space = spaces[0]
+      space_lights = space.lights
+      if space_lights.empty?
+        space_lights = space.spaceType.get.lights
+      end
+      space_people = space.people
+      if space_people.empty?
+        space_people = space.spaceType.get.people
+      end
+
+      next if space_lights.empty? # skip space with no lights
+
+      zone_name = zone.name.to_s
+      next if zone_name =~ /data\s*center/i # skip data centers
+
+      # EnergyPlus v9.4 name change for EMS actuators
+      # https://github.com/NREL/OpenStudio/pull/4104
+      if model.version < OpenStudio::VersionString.new('3.1.0')
+        light_sensor = OpenStudio::Model::EnergyManagementSystemSensor.new(model, 'Zone Lights Electric Power')
+      else
+        light_sensor = OpenStudio::Model::EnergyManagementSystemSensor.new(model, 'Zone Lights Electricity Rate')
+      end
+      light_sensor.setKeyName(zone_name)
+      light_sensor.setName("#{zone_name}_LSr".gsub(/[\s-]/, ''))
+      light_sensor_name = light_sensor.name.to_s
+
+      floor_area = OpenStudio::Model::EnergyManagementSystemInternalVariable.new(model, 'Zone Floor Area')
+      floor_area.setInternalDataIndexKeyName(zone_name)
+      floor_area.setName("#{zone_name}_Area".gsub(/[\s-]/, ''))
+      floor_area_name = floor_area.name.to_s
+
+      # account for multiple lights (also work for single light)
+      big_light = space_lights[0] # find the light with highest power (assuming specified by watts/area)
+      space_lights.each do |light_x|
+        big_light_power = big_light.definition.to_LightsDefinition.get.wattsperSpaceFloorArea.to_f
+        light_x_power = light_x.definition.to_LightsDefinition.get.wattsperSpaceFloorArea.to_f
+        if light_x_power > big_light_power
+          big_light = light_x
+        end
+      end
+
+      add_lights_prog_0 = ''
+      add_lights_prog_null = ''
+      light_id = 0
+      space_lights.each do |light_x|
+        light_id += 1
+        # EnergyPlus v9.4 name change for EMS actuators
+        # https://github.com/NREL/OpenStudio/pull/4104
+        if model.version < OpenStudio::VersionString.new('3.1.0')
+          light_x_actuator = OpenStudio::Model::EnergyManagementSystemActuator.new(light_x, 'Lights', 'Electric Power Level')
+        else
+          light_x_actuator = OpenStudio::Model::EnergyManagementSystemActuator.new(light_x, 'Lights', 'Electricity Rate')
+        end
+        light_x_actuator.setName("#{zone_name}_Light#{light_id}_Actuator".gsub(/[\s-]/, ''))
+        light_x_actuator_name = light_x_actuator.name.to_s
+        add_lights_prog_null += "\n      SET #{light_x_actuator_name} = NULL,"
+        if light_x == big_light
+          add_lights_prog_0 += "\n      SET #{light_x_actuator_name} = 0.02*#{floor_area_name}/0.09290304,"
+          next
+        end
+        add_lights_prog_0 += "\n      SET #{light_x_actuator_name} = 0,"
+      end
+
+      light_ems_prog = OpenStudio::Model::EnergyManagementSystemProgram.new(model)
+      light_ems_prog.setName("SET_#{zone_name}_Light_EMS_Program".gsub(/[\s-]/, ''))
+      light_ems_prog_body = <<-EMS
+      SET #{light_sensor_name}_IP=0.093*#{light_sensor_name}/#{floor_area_name},
+      IF (#{business_sensor_name} <= 0) && (#{light_sensor_name}_IP >= 0.02),#{add_lights_prog_0}
+      ELSE,#{add_lights_prog_null}
+      ENDIF
+      EMS
+      light_ems_prog.setBody(light_ems_prog_body)
+
+      light_ems_prog_manager = OpenStudio::Model::EnergyManagementSystemProgramCallingManager.new(model)
+      light_ems_prog_manager.setName("SET_#{zone_name}_Light_EMS_Program_Manager")
+      light_ems_prog_manager.setCallingPoint('AfterPredictorAfterHVACManagers')
+      light_ems_prog_manager.addProgram(light_ems_prog)
+    end
+    return true
+  end
+end