urbanopt-cli 0.4.1 → 0.6.1

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/tests/test_enclosure.rb

@@ -0,0 +1,151 @@
+# frozen_string_literal: true
+
+require_relative '../resources/minitest_helper'
+require 'openstudio'
+require 'openstudio/measure/ShowRunnerOutput'
+require 'fileutils'
+require_relative '../measure.rb'
+require_relative '../resources/util.rb'
+require_relative 'util.rb'
+
+class HPXMLtoOpenStudioEnclosureTest < MiniTest::Test
+  def sample_files_dir
+    return File.join(File.dirname(__FILE__), '..', '..', 'workflow', 'sample_files')
+  end
+
+  def test_windows
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Check window properties
+    hpxml.windows.each do |window|
+      os_window = model.getSubSurfaces.select { |w| w.name.to_s == window.id }[0]
+      os_simple_glazing = os_window.construction.get.to_LayeredConstruction.get.getLayer(0).to_SimpleGlazing.get
+
+      assert_equal(window.shgc, os_simple_glazing.solarHeatGainCoefficient)
+      assert_in_epsilon(window.ufactor, UnitConversions.convert(os_simple_glazing.uFactor, 'W/(m^2*K)', 'Btu/(hr*ft^2*F)'), 0.001)
+    end
+  end
+
+  def test_windows_shading
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-enclosure-windows-shading.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    summer_date = OpenStudio::Date.new(OpenStudio::MonthOfYear.new('June'), 1, model.yearDescription.get.assumedYear)
+    winter_date = OpenStudio::Date.new(OpenStudio::MonthOfYear.new('January'), 1, model.yearDescription.get.assumedYear)
+
+    hpxml.windows.each do |window|
+      sf_summer = window.interior_shading_factor_summer
+      sf_winter = window.interior_shading_factor_winter
+      sf_summer *= window.exterior_shading_factor_summer unless window.exterior_shading_factor_summer.nil?
+      sf_winter *= window.exterior_shading_factor_winter unless window.exterior_shading_factor_winter.nil?
+
+      # Check shading transmittance for sky beam and sky diffuse
+      os_shading_surface = model.getShadingSurfaces.select { |ss| ss.name.to_s.start_with? window.id }[0]
+      if (sf_summer == 1) && (sf_winter == 1)
+        assert_nil(os_shading_surface) # No shading
+      else
+        refute_nil(os_shading_surface) # Shading
+        summer_transmittance = os_shading_surface.transmittanceSchedule.get.to_ScheduleRuleset.get.getDaySchedules(summer_date, summer_date).map { |ds| ds.values.sum }.sum
+        winter_transmittance = os_shading_surface.transmittanceSchedule.get.to_ScheduleRuleset.get.getDaySchedules(winter_date, winter_date).map { |ds| ds.values.sum }.sum
+        assert_equal(sf_summer, summer_transmittance)
+        assert_equal(sf_winter, winter_transmittance)
+      end
+
+      # Check subsurface view factor to ground
+      subsurface_view_factor = 0.5
+      window_actuator = model.getEnergyManagementSystemActuators.select { |w| w.actuatedComponent.get.name.to_s == window.id }[0]
+      program_values = get_ems_values(model.getEnergyManagementSystemPrograms, 'fixedwindow view factor to ground program')
+      assert_equal(subsurface_view_factor, program_values["#{window_actuator.name.to_s}"][0])
+    end
+  end
+
+  def test_skylights
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-enclosure-skylights.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Check skylight properties
+    hpxml.skylights.each do |skylight|
+      os_skylight = model.getSubSurfaces.select { |w| w.name.to_s == skylight.id }[0]
+      os_simple_glazing = os_skylight.construction.get.to_LayeredConstruction.get.getLayer(0).to_SimpleGlazing.get
+
+      assert_equal(skylight.shgc, os_simple_glazing.solarHeatGainCoefficient)
+      assert_in_epsilon(skylight.ufactor / 1.2, UnitConversions.convert(os_simple_glazing.uFactor, 'W/(m^2*K)', 'Btu/(hr*ft^2*F)'), 0.001)
+    end
+  end
+
+  def test_skylights_shading
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-enclosure-skylights-shading.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    summer_date = OpenStudio::Date.new(OpenStudio::MonthOfYear.new('June'), 1, model.yearDescription.get.assumedYear)
+    winter_date = OpenStudio::Date.new(OpenStudio::MonthOfYear.new('January'), 1, model.yearDescription.get.assumedYear)
+
+    hpxml.skylights.each do |skylight|
+      sf_summer = skylight.interior_shading_factor_summer
+      sf_winter = skylight.interior_shading_factor_winter
+      sf_summer *= skylight.exterior_shading_factor_summer unless skylight.exterior_shading_factor_summer.nil?
+      sf_winter *= skylight.exterior_shading_factor_winter unless skylight.exterior_shading_factor_winter.nil?
+
+      # Check shading transmittance for sky beam and sky diffuse
+      os_shading_surface = model.getShadingSurfaces.select { |ss| ss.name.to_s.start_with? skylight.id }[0]
+      if (sf_summer == 1) && (sf_winter == 1)
+        assert_nil(os_shading_surface) # No shading
+      else
+        refute_nil(os_shading_surface) # Shading
+        summer_transmittance = os_shading_surface.transmittanceSchedule.get.to_ScheduleRuleset.get.getDaySchedules(summer_date, summer_date).map { |ds| ds.values.sum }.sum
+        winter_transmittance = os_shading_surface.transmittanceSchedule.get.to_ScheduleRuleset.get.getDaySchedules(winter_date, winter_date).map { |ds| ds.values.sum }.sum
+        assert_equal(sf_summer, summer_transmittance)
+        assert_equal(sf_winter, winter_transmittance)
+      end
+
+      # Check subsurface view factor to ground
+      subsurface_view_factor = 0.05 # 6:12 pitch
+      skylight_actuator = model.getEnergyManagementSystemActuators.select { |w| w.actuatedComponent.get.name.to_s == skylight.id }[0]
+      program_values = get_ems_values(model.getEnergyManagementSystemPrograms, 'skylight view factor to ground program')
+      assert_equal(subsurface_view_factor, program_values["#{skylight_actuator.name.to_s}"][0])
+    end
+  end
+
+  def _test_measure(args_hash)
+    # create an instance of the measure
+    measure = HPXMLtoOpenStudio.new
+
+    runner = OpenStudio::Measure::OSRunner.new(OpenStudio::WorkflowJSON.new)
+    model = OpenStudio::Model::Model.new
+
+    # get arguments
+    args_hash['output_dir'] = 'tests'
+    arguments = measure.arguments(model)
+    argument_map = OpenStudio::Measure.convertOSArgumentVectorToMap(arguments)
+
+    # populate argument with specified hash value if specified
+    arguments.each do |arg|
+      temp_arg_var = arg.clone
+      if args_hash.has_key?(arg.name)
+        assert(temp_arg_var.setValue(args_hash[arg.name]))
+      end
+      argument_map[arg.name] = temp_arg_var
+    end
+
+    # run the measure
+    measure.run(model, runner, argument_map)
+    result = runner.result
+
+    # show the output
+    show_output(result) unless result.value.valueName == 'Success'
+
+    # assert that it ran correctly
+    assert_equal('Success', result.value.valueName)
+
+    hpxml = HPXML.new(hpxml_path: args_hash['hpxml_path'])
+
+    File.delete(File.join(File.dirname(__FILE__), 'in.xml'))
+
+    return model, hpxml
+  end
+end

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/tests/test_generator.rb

@@ -0,0 +1,101 @@
+# frozen_string_literal: true
+
+require_relative '../resources/minitest_helper'
+require 'openstudio'
+require 'openstudio/measure/ShowRunnerOutput'
+require 'fileutils'
+require_relative '../measure.rb'
+require_relative '../resources/util.rb'
+
+class HPXMLtoOpenStudioGeneratorTest < MiniTest::Test
+  def sample_files_dir
+    return File.join(File.dirname(__FILE__), '..', '..', 'workflow', 'sample_files')
+  end
+
+  def get_generator(model, name)
+    generator = nil
+    model.getGeneratorMicroTurbines.each do |g|
+      next unless g.name.to_s.start_with? "#{name} "
+
+      return g
+    end
+  end
+
+  def test_generator
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-misc-generators.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    hpxml.generators.each do |hpxml_generator|
+      generator = get_generator(model, hpxml_generator.id)
+
+      # Check object
+      assert_equal(EPlus.fuel_type(hpxml_generator.fuel_type), generator.fuelType)
+      assert_in_epsilon(57.1, generator.referenceElectricalPowerOutput, 0.01)
+      assert_in_epsilon(57.1, generator.minimumFullLoadElectricalPowerOutput, 0.01)
+      assert_in_epsilon(57.1, generator.maximumFullLoadElectricalPowerOutput, 0.01)
+      assert_in_epsilon(0.2, generator.referenceElectricalEfficiencyUsingLowerHeatingValue, 0.01)
+      assert_equal(generator.fuelHigherHeatingValue, generator.fuelLowerHeatingValue)
+      assert_equal(0, generator.standbyPower)
+      assert_equal(0, generator.ancillaryPower)
+    end
+  end
+
+  def test_generator_shared
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-generator.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    hpxml.generators.each do |hpxml_generator|
+      generator = get_generator(model, hpxml_generator.id)
+
+      # Check object
+      assert_equal(EPlus.fuel_type(hpxml_generator.fuel_type), generator.fuelType)
+      assert_in_epsilon(95.1, generator.referenceElectricalPowerOutput, 0.01)
+      assert_in_epsilon(95.1, generator.minimumFullLoadElectricalPowerOutput, 0.01)
+      assert_in_epsilon(95.1, generator.maximumFullLoadElectricalPowerOutput, 0.01)
+      assert_in_epsilon(0.2, generator.referenceElectricalEfficiencyUsingLowerHeatingValue, 0.01)
+      assert_equal(generator.fuelHigherHeatingValue, generator.fuelLowerHeatingValue)
+      assert_equal(0, generator.standbyPower)
+      assert_equal(0, generator.ancillaryPower)
+    end
+  end
+
+  def _test_measure(args_hash)
+    # create an instance of the measure
+    measure = HPXMLtoOpenStudio.new
+
+    runner = OpenStudio::Measure::OSRunner.new(OpenStudio::WorkflowJSON.new)
+    model = OpenStudio::Model::Model.new
+
+    # get arguments
+    args_hash['output_dir'] = 'tests'
+    arguments = measure.arguments(model)
+    argument_map = OpenStudio::Measure.convertOSArgumentVectorToMap(arguments)
+
+    # populate argument with specified hash value if specified
+    arguments.each do |arg|
+      temp_arg_var = arg.clone
+      if args_hash.has_key?(arg.name)
+        assert(temp_arg_var.setValue(args_hash[arg.name]))
+      end
+      argument_map[arg.name] = temp_arg_var
+    end
+
+    # run the measure
+    measure.run(model, runner, argument_map)
+    result = runner.result
+
+    # show the output
+    show_output(result) unless result.value.valueName == 'Success'
+
+    # assert that it ran correctly
+    assert_equal('Success', result.value.valueName)
+
+    hpxml = HPXML.new(hpxml_path: args_hash['hpxml_path'])
+
+    File.delete(File.join(File.dirname(__FILE__), 'in.xml'))
+
+    return model, hpxml
+  end
+end

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/tests/test_hotwater_appliance.rb

@@ -2,8 +2,7 @@
 
 require_relative '../resources/minitest_helper'
 require 'openstudio'
-require 'openstudio/ruleset/ShowRunnerOutput'
-require 'minitest/autorun'
+require 'openstudio/measure/ShowRunnerOutput'
 require 'fileutils'
 require_relative '../measure.rb'
 require_relative '../resources/util.rb'
@@ -69,7 +68,7 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
     if fuel.empty?
       return
     elsif fuel.uniq.size != 1
-      fail 'different fuels'
+      flunk 'different fuels'
     else
       return fuel[0]
     end
@@ -238,12 +237,12 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
 
   def test_dhw_shared_water_heater_recirc
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-dhw-shared-water-heater-recirc.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-water-heater-recirc.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # water use equipment peak flows
     fixture_gpd = 44.60
-    dist_gpd = 10.2343
+    dist_gpd = 12.354
     cw_gpd = 3.7116
     dw_gpd = 2.7342
     assert_in_epsilon(cw_gpd, get_wu_gpd(model, Constants.ObjectNameClothesWasher), 0.001)
@@ -306,12 +305,12 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
 
   def test_dhw_shared_laundry
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-dhw-shared-laundry-room.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-laundry-room.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # water use equipment peak flows
     fixture_gpd = 44.60
-    dist_gpd = 10.2343
+    dist_gpd = 12.354
     cw_gpd = 3.7116
     dw_gpd = 2.7342
     assert_in_epsilon(cw_gpd, get_wu_gpd(model, Constants.ObjectNameClothesWasher), 0.001)
@@ -450,17 +449,6 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
     assert_in_epsilon(pump_kwh_yr, get_ee_kwh_per_year(model, Constants.ObjectNameHotWaterRecircPump), 0.001)
   end
 
-  def test_dhw_recirc_shared
-    args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-dhw-shared-water-heater-recirc.xml'))
-    model, hpxml = _test_measure(args_hash)
-
-    # Get HPXML values
-    hot_water_distribution = hpxml.hot_water_distributions[0]
-    pump_kwh_yr = 8760.0 / 1000.0 * hot_water_distribution.shared_recirculation_pump_power / hot_water_distribution.shared_recirculation_number_of_units_served
-    assert_in_epsilon(pump_kwh_yr, get_ee_kwh_per_year(model, Constants.ObjectNameHotWaterRecircPump), 0.001)
-  end
-
   def test_appliances_none
     args_hash = {}
     args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-appliances-none.xml'))
@@ -983,6 +971,7 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
     model = OpenStudio::Model::Model.new
 
     # get arguments
+    args_hash['output_dir'] = 'tests'
     arguments = measure.arguments(model)
     argument_map = OpenStudio::Measure.convertOSArgumentVectorToMap(arguments)
 
@@ -1007,6 +996,8 @@ class HPXMLtoOpenStudioHotWaterApplianceTest < MiniTest::Test
 
     hpxml = HPXML.new(hpxml_path: args_hash['hpxml_path'])
 
+    File.delete(File.join(File.dirname(__FILE__), 'in.xml'))
+
     return model, hpxml
   end
 end

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/tests/test_hvac.rb

@@ -1,14 +1,12 @@
 # frozen_string_literal: true
 
-# TODO: Add fan/pump checks
-
 require_relative '../resources/minitest_helper'
 require 'openstudio'
-require 'openstudio/ruleset/ShowRunnerOutput'
-require 'minitest/autorun'
+require 'openstudio/measure/ShowRunnerOutput'
 require 'fileutils'
 require_relative '../measure.rb'
 require_relative '../resources/util.rb'
+require_relative 'util.rb'
 
 class HPXMLtoOpenStudioHVACTest < MiniTest::Test
   def sample_files_dir
@@ -30,6 +28,12 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     cop = 4.0 # Expected value
     assert_in_epsilon(cop, clg_coil.ratedCOP.get, 0.01)
     assert_in_epsilon(capacity, clg_coil.ratedTotalCoolingCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_central_air_conditioner_2_speed
@@ -44,11 +48,17 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     # Check cooling coil
     assert_equal(1, model.getCoilCoolingDXMultiSpeeds.size)
     clg_coil = model.getCoilCoolingDXMultiSpeeds[0]
-    cops = [5.0, 4.63] # Expected values
+    cops = [4.77, 4.42] # Expected values
     cops.each_with_index do |cop, i|
       assert_in_epsilon(cop, clg_coil.stages[i].grossRatedCoolingCOP, 0.01)
     end
     assert_in_epsilon(capacity, clg_coil.stages[-1].grossRatedTotalCoolingCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_central_air_conditioner_var_speed
@@ -68,6 +78,12 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
       assert_in_epsilon(cop, clg_coil.stages[i].grossRatedCoolingCOP, 0.01)
     end
     assert_in_epsilon(capacity, clg_coil.stages[-1].grossRatedTotalCoolingCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_room_air_conditioner
@@ -88,6 +104,10 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     assert_in_epsilon(capacity, clg_coil.ratedTotalCoolingCapacity.get, 0.01)
   end
 
+  def test_evap_cooler
+    # TODO
+  end
+
   def test_furnace_gas
     args_hash = {}
     args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-furnace-gas-only.xml'))
@@ -104,7 +124,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     htg_coil = model.getCoilHeatingGass[0]
     assert_in_epsilon(afue, htg_coil.gasBurnerEfficiency, 0.01)
     assert_in_epsilon(capacity, htg_coil.nominalCapacity.get, 0.01)
-    assert_equal(EPlus.input_fuel_map(fuel), htg_coil.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), htg_coil.fuelType)
   end
 
   def test_furnace_electric
@@ -140,7 +160,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     assert_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01)
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
   end
 
   def test_boiler_coal
@@ -159,7 +179,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     assert_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01)
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
   end
 
   def test_boiler_electric
@@ -178,7 +198,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     assert_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01)
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
   end
 
   def test_electric_resistance
@@ -214,10 +234,10 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     htg_coil = model.getCoilHeatingGass[0]
     assert_in_epsilon(efficiency, htg_coil.gasBurnerEfficiency, 0.01)
     assert_in_epsilon(capacity, htg_coil.nominalCapacity.get, 0.01)
-    assert_equal(EPlus.input_fuel_map(fuel), htg_coil.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), htg_coil.fuelType)
   end
 
-  def test_central_air_to_air_heat_pump_1_speed
+  def test_air_to_air_heat_pump_1_speed
     args_hash = {}
     args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-air-to-air-heat-pump-1-speed.xml'))
     model, hpxml = _test_measure(args_hash)
@@ -248,9 +268,15 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     supp_htg_coil = model.getCoilHeatingElectrics[0]
     assert_in_epsilon(backup_efficiency, supp_htg_coil.efficiency, 0.01)
     assert_in_epsilon(supp_htg_capacity, supp_htg_coil.nominalCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
-  def test_central_air_to_air_heat_pump_2_speed
+  def test_air_to_air_heat_pump_2_speed
     args_hash = {}
     args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-air-to-air-heat-pump-2-speed.xml'))
     model, hpxml = _test_measure(args_hash)
@@ -285,9 +311,15 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     supp_htg_coil = model.getCoilHeatingElectrics[0]
     assert_in_epsilon(backup_efficiency, supp_htg_coil.efficiency, 0.01)
     assert_in_epsilon(supp_htg_capacity, supp_htg_coil.nominalCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
-  def test_central_air_to_air_heat_pump_var_speed
+  def test_air_to_air_heat_pump_var_speed
     args_hash = {}
     args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-air-to-air-heat-pump-var-speed.xml'))
     model, hpxml = _test_measure(args_hash)
@@ -322,6 +354,12 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     supp_htg_coil = model.getCoilHeatingElectrics[0]
     assert_in_epsilon(backup_efficiency, supp_htg_coil.efficiency, 0.01)
     assert_in_epsilon(supp_htg_capacity, supp_htg_coil.nominalCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_mini_split_heat_pump
@@ -356,6 +394,37 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     assert_equal(1, model.getCoilHeatingElectrics.size)
     supp_htg_coil = model.getCoilHeatingElectrics[0]
     assert_in_delta(0, supp_htg_coil.nominalCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
+  end
+
+  def test_mini_split_air_conditioner
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-mini-split-air-conditioner-only-ductless.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    cooling_system = hpxml.cooling_systems[0]
+    clg_capacity = UnitConversions.convert(cooling_system.cooling_capacity, 'Btu/hr', 'W')
+
+    # Check cooling coil
+    assert_equal(1, model.getCoilCoolingDXMultiSpeeds.size)
+    clg_coil = model.getCoilCoolingDXMultiSpeeds[0]
+    cops = [5.76, 4.99, 4.19, 3.10] # Expected values
+    cops.each_with_index do |cop, i|
+      assert_in_epsilon(cop, clg_coil.stages[i].grossRatedCoolingCOP, 0.01)
+    end
+    assert_in_epsilon(clg_capacity * 1.2, clg_coil.stages[-1].grossRatedTotalCoolingCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_ground_to_air_heat_pump
@@ -389,11 +458,17 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     supp_htg_coil = model.getCoilHeatingElectrics[0]
     assert_in_epsilon(backup_efficiency, supp_htg_coil.efficiency, 0.01)
     assert_in_epsilon(supp_htg_capacity, supp_htg_coil.nominalCapacity.get, 0.01)
+
+    # Check EMS
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    assert(program_values.empty?) # Check no EMS program
   end
 
   def test_shared_chiller_baseboard
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-chiller-only-baseboard.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-chiller-only-baseboard.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -410,7 +485,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
 
   def test_shared_chiller_fan_coil
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-chiller-only-fan-coil.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-chiller-only-fan-coil.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -420,14 +495,14 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     # Check cooling coil
     assert_equal(1, model.getCoilCoolingDXSingleSpeeds.size)
     clg_coil = model.getCoilCoolingDXSingleSpeeds[0]
-    cop = 3.67 # Expected value
+    cop = 3.71 # Expected value
     assert_in_epsilon(cop, clg_coil.ratedCOP.get, 0.01)
     refute_in_epsilon(capacity, clg_coil.ratedTotalCoolingCapacity.get, 0.01) # Uses autosized capacity
   end
 
   def test_shared_chiller_water_loop_heat_pump
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-chiller-only-water-loop-heat-pump.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-chiller-only-water-loop-heat-pump.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -437,14 +512,14 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     # Check cooling coil
     assert_equal(1, model.getCoilCoolingDXSingleSpeeds.size)
     clg_coil = model.getCoilCoolingDXSingleSpeeds[0]
-    cop = 1.34 # Expected value
+    cop = 1.35 # Expected value
     assert_in_epsilon(cop, clg_coil.ratedCOP.get, 0.01)
     refute_in_epsilon(capacity, clg_coil.ratedTotalCoolingCapacity.get, 0.01) # Uses autosized capacity
   end
 
   def test_shared_cooling_tower_water_loop_heat_pump
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-cooling-tower-only-water-loop-heat-pump.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-cooling-tower-only-water-loop-heat-pump.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -454,14 +529,14 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     # Check cooling coil
     assert_equal(1, model.getCoilCoolingDXSingleSpeeds.size)
     clg_coil = model.getCoilCoolingDXSingleSpeeds[0]
-    cop = 3.94 # Expected value
+    cop = 3.98 # Expected value
     assert_in_epsilon(cop, clg_coil.ratedCOP.get, 0.01)
     refute_in_epsilon(capacity, clg_coil.ratedTotalCoolingCapacity.get, 0.01) # Uses autosized capacity
   end
 
   def test_shared_boiler_baseboard
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-boiler-only-baseboard.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-boiler-only-baseboard.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -475,12 +550,12 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     refute_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01) # Uses autosized capacity
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
   end
 
   def test_shared_boiler_fan_coil
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-boiler-only-fan-coil.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-boiler-only-fan-coil.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -494,12 +569,12 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     refute_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01) # Uses autosized capacity
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
   end
 
   def test_shared_boiler_water_loop_heat_pump
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-boiler-only-water-loop-heat-pump.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-boiler-only-water-loop-heat-pump.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -507,14 +582,15 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     afue = heating_system.heating_efficiency_afue
     capacity = UnitConversions.convert(heating_system.heating_capacity.to_f, 'Btu/hr', 'W')
     fuel = heating_system.heating_system_fuel
-    wlhp_cop = heating_system.wlhp_heating_efficiency_cop
+    heat_pump = hpxml.heat_pumps[0]
+    wlhp_cop = heat_pump.heating_efficiency_cop
 
     # Check boiler
     assert_equal(1, model.getBoilerHotWaters.size)
     boiler = model.getBoilerHotWaters[0]
     assert_in_epsilon(afue, boiler.nominalThermalEfficiency, 0.01)
     refute_in_epsilon(capacity, boiler.nominalCapacity.get, 0.01) # Uses autosized capacity
-    assert_equal(EPlus.input_fuel_map(fuel), boiler.fuelType)
+    assert_equal(EPlus.fuel_type(fuel), boiler.fuelType)
 
     # Check cooling coil
     assert_equal(0, model.getCoilCoolingDXSingleSpeeds.size)
@@ -534,7 +610,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
 
   def test_shared_ground_loop_ground_to_air_heat_pump
     args_hash = {}
-    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-shared-ground-loop-ground-to-air-heat-pump.xml'))
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-bldgtype-multifamily-shared-ground-loop-ground-to-air-heat-pump.xml'))
     model, hpxml = _test_measure(args_hash)
 
     # Get HPXML values
@@ -565,6 +641,168 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     assert_in_epsilon(supp_htg_capacity, supp_htg_coil.nominalCapacity.get, 0.01)
   end
 
+  def test_install_quality_air_to_air_heat_pump_1_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-air-to-air-heat-pump-1-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    heat_pump = hpxml.heat_pumps[0]
+    airflow_defect = heat_pump.airflow_defect_ratio
+    charge_defect = heat_pump.charge_defect_ratio
+    fan_watts_cfm = heat_pump.fan_watts_per_cfm
+
+    # model objects:
+    # Unitary system
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    cooling_cfm = UnitConversions.convert(unitary_system.supplyAirFlowRateDuringCoolingOperation.get, 'm^3/s', 'cfm')
+    heating_cfm = UnitConversions.convert(unitary_system.supplyAirFlowRateDuringHeatingOperation.get, 'm^3/s', 'cfm')
+
+    # Cooling coil
+    assert_equal(1, model.getCoilCoolingDXSingleSpeeds.size)
+    clg_coil = model.getCoilCoolingDXSingleSpeeds[0]
+    rated_airflow_cfm_clg = UnitConversions.convert(clg_coil.ratedAirFlowRate.get, 'm^3/s', 'cfm')
+
+    # Heating coil
+    assert_equal(1, model.getCoilHeatingDXSingleSpeeds.size)
+    htg_coil = model.getCoilHeatingDXSingleSpeeds[0]
+    rated_airflow_cfm_htg = UnitConversions.convert(htg_coil.ratedAirFlowRate.get, 'm^3/s', 'cfm')
+
+    # Fan
+    fanonoff = unitary_system.supplyFan.get.to_FanOnOff.get
+    assert_in_epsilon(fan_watts_cfm, fanonoff.pressureRise / fanonoff.fanEfficiency * UnitConversions.convert(1.0, 'cfm', 'm^3/s'), 0.01)
+
+    # Check installation quality EMS
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+
+    # defect ratios in EMS is calculated correctly
+    assert_in_epsilon(program_values['F_CH'].sum, charge_defect, 0.01)
+    assert_in_epsilon(program_values['FF_AF_c'].sum, cooling_cfm / rated_airflow_cfm_clg, 0.01)
+    assert_in_epsilon(program_values['FF_AF_h'].sum, heating_cfm / rated_airflow_cfm_htg, 0.01)
+  end
+
+  def test_install_quality_air_to_air_heat_pump_2_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-air-to-air-heat-pump-2-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    heat_pump = hpxml.heat_pumps[0]
+    _check_install_quality_multispeed_ratio(heat_pump, model, heat_pump)
+  end
+
+  def test_install_quality_air_to_air_heat_pump_var_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-air-to-air-heat-pump-var-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    heat_pump = hpxml.heat_pumps[0]
+    _check_install_quality_multispeed_ratio(heat_pump, model, heat_pump)
+  end
+
+  def test_install_quality_furnace_central_air_conditioner_1_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-furnace-gas-central-ac-1-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    cooling_system = hpxml.cooling_systems[0]
+    heating_system = hpxml.heating_systems[0]
+    airflow_defect_clg = cooling_system.airflow_defect_ratio
+    charge_defect = cooling_system.charge_defect_ratio
+    fan_watts_cfm = cooling_system.fan_watts_per_cfm
+    fan_watts_cfm2 = heating_system.fan_watts_per_cfm
+
+    # model objects:
+    # Unitary system
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    cooling_cfm = UnitConversions.convert(unitary_system.supplyAirFlowRateDuringCoolingOperation.get, 'm^3/s', 'cfm')
+
+    # Cooling coil
+    assert_equal(1, model.getCoilCoolingDXSingleSpeeds.size)
+    clg_coil = model.getCoilCoolingDXSingleSpeeds[0]
+    rated_airflow_cfm = UnitConversions.convert(clg_coil.ratedAirFlowRate.get, 'm^3/s', 'cfm')
+
+    # Fan
+    fanonoff = unitary_system.supplyFan.get.to_FanOnOff.get
+    assert_in_epsilon(fan_watts_cfm, fanonoff.pressureRise / fanonoff.fanEfficiency * UnitConversions.convert(1.0, 'cfm', 'm^3/s'), 0.01)
+    assert_in_epsilon(fan_watts_cfm2, fanonoff.pressureRise / fanonoff.fanEfficiency * UnitConversions.convert(1.0, 'cfm', 'm^3/s'), 0.01)
+
+    # Check installation quality EMS
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+
+    # defect ratios in EMS is calculated correctly
+    assert_in_epsilon(program_values['F_CH'].sum, charge_defect, 0.01)
+
+    # Fan air flow has already applied air flow defect ratio
+    assert_in_epsilon(program_values['FF_AF_c'].sum, cooling_cfm / rated_airflow_cfm, 0.01)
+  end
+
+  def test_install_quality_furnace_central_air_conditioner_2_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-furnace-gas-central-ac-2-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    cooling_system = hpxml.cooling_systems[0]
+    _check_install_quality_multispeed_ratio(cooling_system, model)
+  end
+
+  def test_install_quality_furnace_central_air_conditioner_var_speed_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-furnace-gas-central-ac-var-speed.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    cooling_system = hpxml.cooling_systems[0]
+    _check_install_quality_multispeed_ratio(cooling_system, model)
+  end
+
+  def test_install_quality_furnace_gas_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-furnace-gas-only.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    heating_system = hpxml.heating_systems[0]
+    airflow_defect_htg = heating_system.airflow_defect_ratio
+    fan_watts_cfm = heating_system.fan_watts_per_cfm
+
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+
+    # Fan
+    fanonoff = unitary_system.supplyFan.get.to_FanOnOff.get
+    assert_in_epsilon(fan_watts_cfm, fanonoff.pressureRise / fanonoff.fanEfficiency * UnitConversions.convert(1.0, 'cfm', 'm^3/s'), 0.01)
+  end
+
+  def test_install_quality_ground_to_air_heat_pump_ratio
+    # TODO
+  end
+
+  def test_install_quality_mini_split_air_conditioner_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-mini-split-air-conditioner-only-ducted.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    cooling_system = hpxml.cooling_systems[0]
+    _check_install_quality_multispeed_ratio(cooling_system, model)
+  end
+
+  def test_install_quality_mini_split_heat_pump_ratio
+    args_hash = {}
+    args_hash['hpxml_path'] = File.absolute_path(File.join(sample_files_dir, 'base-hvac-install-quality-all-mini-split-heat-pump-ducted.xml'))
+    model, hpxml = _test_measure(args_hash)
+
+    # Get HPXML values
+    heat_pump = hpxml.heat_pumps[0]
+    _check_install_quality_multispeed_ratio(heat_pump, model, heat_pump)
+  end
+
   def _test_measure(args_hash)
     # create an instance of the measure
     measure = HPXMLtoOpenStudio.new
@@ -573,6 +811,7 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
     model = OpenStudio::Model::Model.new
 
     # get arguments
+    args_hash['output_dir'] = 'tests'
     arguments = measure.arguments(model)
     argument_map = OpenStudio::Measure.convertOSArgumentVectorToMap(arguments)
 
@@ -597,6 +836,55 @@ class HPXMLtoOpenStudioHVACTest < MiniTest::Test
 
     hpxml = HPXML.new(hpxml_path: args_hash['hpxml_path'])
 
+    File.delete(File.join(File.dirname(__FILE__), 'in.xml'))
+
     return model, hpxml
   end
+
+  def _check_install_quality_multispeed_ratio(hpxml_clg_sys, model, hpxml_htg_sys = nil)
+    airflow_defect = hpxml_clg_sys.airflow_defect_ratio
+    charge_defect = hpxml_clg_sys.charge_defect_ratio
+    fan_watts_cfm = hpxml_clg_sys.fan_watts_per_cfm
+
+    # model objects:
+    # Unitary system
+    assert_equal(1, model.getAirLoopHVACUnitarySystems.size)
+    unitary_system = model.getAirLoopHVACUnitarySystems[0]
+    perf = unitary_system.designSpecificationMultispeedObject.get.to_UnitarySystemPerformanceMultispeed.get
+    clg_ratios = perf.supplyAirflowRatioFields.map { |field| field.coolingRatio.get }
+    cooling_cfm = UnitConversions.convert(unitary_system.supplyAirFlowRateDuringCoolingOperation.get, 'm^3/s', 'cfm')
+
+    # Cooling coil
+    assert_equal(1, model.getCoilCoolingDXMultiSpeeds.size)
+    clg_coil = model.getCoilCoolingDXMultiSpeeds[0]
+    rated_airflow_cfm_clg = []
+    clg_coil.stages.each do |stage|
+      rated_airflow_cfm_clg << UnitConversions.convert(stage.ratedAirFlowRate.get, 'm^3/s', 'cfm')
+    end
+
+    # Fan
+    fanonoff = unitary_system.supplyFan.get.to_FanOnOff.get
+    assert_in_epsilon(fan_watts_cfm, fanonoff.pressureRise / fanonoff.fanEfficiency * UnitConversions.convert(1.0, 'cfm', 'm^3/s'), 0.01)
+
+    # Check installation quality EMS
+    program_values = get_ems_values(model.getEnergyManagementSystemPrograms, "#{unitary_system.name} install quality")
+    clg_speed_cfms = clg_ratios.map { |ratio| cooling_cfm * ratio }
+    assert_in_epsilon(program_values['F_CH'].sum, charge_defect, 0.01)
+    assert_in_epsilon(program_values['FF_AF_c'].sum, clg_speed_cfms.zip(rated_airflow_cfm_clg).map { |cfm, rated_cfm| cfm / rated_cfm }.sum, 0.01)
+    if not hpxml_htg_sys.nil?
+      heating_cfm = UnitConversions.convert(unitary_system.supplyAirFlowRateDuringHeatingOperation.get, 'm^3/s', 'cfm')
+      htg_ratios = perf.supplyAirflowRatioFields.map { |field| field.heatingRatio.get }
+
+      # Heating coil
+      assert_equal(1, model.getCoilHeatingDXMultiSpeeds.size)
+      htg_coil = model.getCoilHeatingDXMultiSpeeds[0]
+      rated_airflow_cfm_htg = []
+      htg_coil.stages.each do |stage|
+        rated_airflow_cfm_htg << UnitConversions.convert(stage.ratedAirFlowRate.get, 'm^3/s', 'cfm')
+      end
+
+      htg_speed_cfms = htg_ratios.map { |ratio| heating_cfm * ratio }
+      assert_in_epsilon(program_values['FF_AF_h'].sum, htg_speed_cfms.zip(rated_airflow_cfm_htg).map { |cfm, rated_cfm| cfm / rated_cfm }.sum, 0.01)
+    end
+  end
 end