buildingsync 0.2.0

data/lib/buildingsync/get_bcl_weather_file.rb

@@ -0,0 +1,326 @@
+# frozen_string_literal: true
+
+# *******************************************************************************
+# OpenStudio(R), Copyright (c) 2008-2020, Alliance for Sustainable Energy, LLC.
+# BuildingSync(R), Copyright (c) 2015-2020, Alliance for Sustainable Energy, LLC.
+# All rights reserved.
+#
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+#
+# (1) Redistributions of source code must retain the above copyright notice,
+# this list of conditions and the following disclaimer.
+#
+# (2) Redistributions in binary form must reproduce the above copyright notice,
+# this list of conditions and the following disclaimer in the documentation
+# and/or other materials provided with the distribution.
+#
+# (3) Neither the name of the copyright holder nor the names of any contributors
+# may be used to endorse or promote products derived from this software without
+# specific prior written permission from the respective party.
+#
+# (4) Other than as required in clauses (1) and (2), distributions in any form
+# of modifications or other derivative works may not use the "OpenStudio"
+# trademark, "OS", "os", or any other confusingly similar designation without
+# specific prior written permission from Alliance for Sustainable Energy, LLC.
+#
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDER(S) AND ANY CONTRIBUTORS
+# "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
+# THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER(S), ANY CONTRIBUTORS, THE
+# UNITED STATES GOVERNMENT, OR THE UNITED STATES DEPARTMENT OF ENERGY, NOR ANY OF
+# THEIR EMPLOYEES, BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+# EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
+# OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
+# STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+# OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+# *******************************************************************************
+require 'json'
+
+require 'buildingsync/constants'
+
+module BuildingSync
+  # GetBCLWeatherFile class to manage the process of getting weather files from BCL
+  class GetBCLWeatherFile
+    def initialize
+      # prefix for data weather path
+      @weather_file_path_prefix = WEATHER_DIR
+      @weather_file = File.join(@weather_file_path_prefix, 'weather_file.json')
+      # Above structured as {"weather_file_name":[],"city_name":[],"state_code":[],"weather_id":[]}
+
+      FileUtils.mkdir_p(@weather_file_path_prefix)
+      if File.exist?(@weather_file)
+        File.open(@weather_file, 'r') do |file|
+          @weather_json = JSON.parse(file.read, symbolize_names: true)
+        end
+      else
+        arr = []
+        @weather_json = {
+          'weather_file_name': arr,
+          'city_name': arr,
+          'state_code': arr,
+          'weather_id': arr
+        }
+        File.open(@weather_file, 'w') { |f| f.write(@weather_json.to_json) }
+      end
+    end
+
+    # download weather file from city name
+    # @param state [String]
+    # @param city [String]c
+    # @return string
+    def download_weather_file_from_city_name(state, city)
+      weather_file_name = get_weather_file_from_city(city)
+
+      if !weather_file_name.empty?
+        return File.join(@weather_file_path_prefix, weather_file_name)
+      else
+        wmo_no = 0
+        remote = OpenStudio::RemoteBCL.new
+
+        # Search for weather files
+        responses = remote.searchComponentLibrary(city, 'Weather File')
+        choices = OpenStudio::StringVector.new
+
+        filter_response = find_response_from_given_state(responses, state)
+
+        if !filter_response.nil?
+          choices << filter_response.uid
+          filter_response.attributes.each do |attribute|
+            if attribute.name == 'WMO'
+              wmo_no = attribute.valueAsDouble
+            end
+          end
+        end
+
+        if choices.count == 0
+          OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_weather_file_from_city_name',
+                             "Error, could not find uid for state #{state} and city #{city}. Initial count of weather files: #{responses.count}. Please try a different weather file.")
+          return false
+        end
+
+        epw_path = download_weather_file(remote, choices)
+        download_design_day_file(wmo_no, epw_path)
+        return epw_path
+      end
+    end
+
+    # download weather file from weather id
+    # @param weather_id [String]
+    # @return string
+    def download_weather_file_from_weather_id(weather_id)
+      weather_file_name = get_weather_file_from_weather_id(weather_id)
+
+      if !weather_file_name.empty?
+        return File.join(@weather_file_path_prefix, weather_file_name)
+      else
+        wmo_no = 0
+        remote = OpenStudio::RemoteBCL.new
+
+        # Search for weather files
+        responses = remote.searchComponentLibrary(weather_id, 'Weather File')
+
+        choices = OpenStudio::StringVector.new
+
+        responses.each do |response|
+          if response.name.include? 'TMY3'
+            choices << response.uid
+
+            response.attributes.each do |attribute|
+              if attribute.name == 'WMO'
+                wmo_no = attribute.valueAsDouble
+              end
+            end
+          end
+        end
+
+        if choices.count == 0
+          OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_weather_file_from_weather_id',
+                             "Error, could not find uid for #{weather_id}.  Please try a different weather file.")
+          return false
+        end
+
+        epw_path = download_weather_file(remote, choices)
+        download_design_day_file(wmo_no, epw_path)
+        return epw_path
+      end
+    end
+
+    # download weather file
+    # @param remote [OpenStudio::RemoteBCL]
+    # @param choices [OpenStudio::StringVector]
+    # @return string
+    def download_weather_file(remote, choices)
+      epw_path = ''
+
+      choices.each do |uid|
+        remote.downloadComponent(uid)
+        component = remote.waitForComponentDownload
+
+        if component.empty?
+          OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_weather_file',
+                             "Error, cannot find the EPW weather file with uid: #{uid}.  Please try a different weather file.")
+          return false
+        end
+
+        component = component.get
+
+        files = component.files('epw')
+
+        if files.empty?
+          OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_weather_file',
+                             "Error, cannot find the EPW weather file within the downloaded zip container with uid: #{uid}.  Please try a different weather file.")
+          return false
+        end
+
+        epw_weather_file_path = component.files('epw')[0]
+        dir_path = File.dirname(epw_weather_file_path)
+        weather_file_name = File.basename(epw_weather_file_path)
+
+        epw_path = File.expand_path(@weather_file_path_prefix.to_s, File.dirname(__FILE__))
+
+        Dir.glob("#{dir_path}/**/*.*").each do |filename|
+          FileUtils.mv(filename, epw_path)
+        end
+        epw_path = File.expand_path("#{@weather_file_path_prefix}/#{weather_file_name}", File.dirname(__FILE__))
+      end
+
+      puts "Successfully set weather file to #{epw_path}"
+      return epw_path
+    end
+
+    # download design day file
+    # @param wmo_no [String]
+    # @param epw_path [String]
+    def download_design_day_file(wmo_no, epw_path)
+      remote = OpenStudio::RemoteBCL.new
+      responses = remote.searchComponentLibrary(wmo_no.to_s[0, 6], 'Design Day')
+      choices = OpenStudio::StringVector.new
+
+      idf_path_collection = []
+
+      responses.each do |response|
+        choices << response.uid
+      end
+
+      choices.each do |uid|
+        remote.downloadComponent(uid)
+        component = remote.waitForComponentDownload
+
+        if !component.empty?
+
+          component = component.get
+
+          files = component.files('idf')
+
+          if !files.empty?
+            idf_path_collection.push(component.files('idf')[0])
+          else
+            OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_design_day_file',
+                               "Error, cannot find the design day file within the downloaded zip container with uid: #{uid}.  Please try a different weather file.")
+
+            raise "Error, cannot find the design day file within the downloaded zip container with uid: #{uid}.  Please try a different weather file."
+          end
+        else
+          OpenStudio.logFree(OpenStudio::Error, 'BuildingSync.GetBCLWeatherFile.download_design_day_file',
+                             "Error, cannot find local component for: #{uid}.  Please try a different weather file.")
+          raise "Error, cannot find local component for: #{uid}.  Please try a different weather file."
+        end
+      end
+
+      puts "Successfully downloaded ddy file to #{epw_path}"
+      create_ddy_file(idf_path_collection, epw_path)
+      puts "Successfully combined design day files to ddy file in #{epw_path}"
+    end
+
+    # create design day file (ddy)
+    # @param idf_path_collection [array<string>]
+    # @param epw_path [String]
+    # @return [Boolean]
+    def create_ddy_file(idf_path_collection, epw_path)
+      idf_file_lines = []
+
+      idf_path_collection.each do |idf_file_path|
+        idf_file = File.open(idf_file_path)
+        idf_file_lines.push(idf_file.readlines)
+      end
+
+      design_day_path = File.dirname(epw_path)
+      weather_file_name = File.basename(epw_path, '.*')
+      design_day_file = File.new("#{design_day_path}/#{weather_file_name}.ddy", 'w')
+
+      idf_file_lines.each do |line|
+        design_day_file.puts(line)
+      end
+      design_day_file.close
+    end
+
+    # get weather file from weather ID
+    # @param weather_id [String]
+    # @return [String]
+    def get_weather_file_from_weather_id(weather_id)
+      weather_file_name = ''
+      is_found, counter = find_weather_counter(weather_id)
+      weather_file_name = @weather_json[:weather_file_name][counter] if is_found
+      return weather_file_name
+    end
+
+    # check if weather ID is found in JSON data
+    # @param weather_id [String]
+    # @return [array<boolean, int>]
+    def find_weather_counter(weather_id)
+      counter = 0
+      @weather_json[:weather_id].each do |cname|
+        if cname.include? weather_id
+          return true, counter
+        end
+        counter += 1
+      end
+      return false, counter
+    end
+
+    # get weather file from city
+    # @param city [String]
+    # @return [String]
+    def get_weather_file_from_city(city)
+      weather_file_name = ''
+      is_found, counter = city_found_in_json_data(city)
+      weather_file_name = @weather_json[:weather_file_name][counter] if is_found
+      return weather_file_name
+    end
+
+    # city found in JSON data
+    # @param city [String]
+    # @return [array<boolean, int>]
+    def city_found_in_json_data(city)
+      counter = 0
+      @weather_json[:city_name].each do |cname|
+        if cname.include? city
+          return true, counter
+        end
+        counter += 1
+      end
+      return false, counter
+    end
+
+    # find response from given state
+    # @param responses [array<BCLSearchResult>]
+    # @param state [String]
+    # @return [BCLSearchResult]
+    def find_response_from_given_state(responses, state)
+      responses.each do |response|
+        if response.name.include? 'TMY3'
+          response.attributes.each do |attribute|
+            if attribute.name == 'State'
+              if attribute.valueAsString == state
+                return response
+              end
+            end
+          end
+        end
+      end
+      return nil
+    end
+  end
+end

data/lib/buildingsync/helpers/Model.hvac.rb

@@ -0,0 +1,216 @@
+# frozen_string_literal: true
+
+# Model class
+class OpenStudio::Model::Model
+  # Adds the HVAC system as derived from the combinations of
+  # CBECS 2012 MAINHT and MAINCL fields.
+  # Mapping between combinations and HVAC systems per
+  # http://www.nrel.gov/docs/fy08osti/41956.pdf
+  # Table C-31
+  # @param standard [Standard]
+  # @param system_type [String]
+  # @param zones [Array<OpenStudio::Model::ThermalZone>]
+  def add_cbecs_hvac_system(standard, system_type, zones)
+    case system_type
+    when 'PTAC with hot water heat'
+      standard.model_add_hvac_system(self, 'PTAC', ht = 'NaturalGas', znht = nil, cl = 'Electricity', zones)
+
+    when 'PTAC with gas coil heat'
+      standard.model_add_hvac_system(self, 'PTAC', ht = nil, znht = 'NaturalGas', cl = 'Electricity', zones)
+
+    when 'PTAC with electric baseboard heat'
+      standard.model_add_hvac_system(self, 'PTAC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'PTAC with no heat'
+      standard.model_add_hvac_system(self, 'PTAC', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'PTAC with district hot water heat'
+      standard.model_add_hvac_system(self, 'PTAC', ht = 'DistrictHeating', znht = nil, cl = 'Electricity', zones)
+
+    when 'PTHP'
+      standard.model_add_hvac_system(self, 'PTHP', ht = 'Electricity', znht = nil, cl = 'Electricity', zones)
+
+    when 'PSZ-AC with gas coil heat'
+      standard.model_add_hvac_system(self, 'PSZ-AC', ht = 'NaturalGas', znht = nil, cl = 'Electricity', zones)
+
+    when 'PSZ-AC with electric baseboard heat'
+      standard.model_add_hvac_system(self, 'PSZ-AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'PSZ-AC with no heat'
+      standard.model_add_hvac_system(self, 'PSZ-AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'PSZ-AC with district hot water heat'
+      standard.model_add_hvac_system(self, 'PSZ-AC', ht = 'DistrictHeating', znht = nil, cl = 'Electricity', zones)
+
+    when 'PSZ-HP'
+      standard.model_add_hvac_system(self, 'PSZ-HP', ht = 'Electricity', znht = nil, cl = 'Electricity', zones)
+
+    when 'Fan coil district chilled water with no heat'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = nil, znht = nil, cl = 'DistrictCooling', zones)
+
+    when 'Fan coil district chilled water and boiler'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = 'NaturalGas', znht = nil, cl = 'DistrictCooling', zones)
+
+    when 'Fan coil district chilled water unit heaters'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = nil, znht = nil, cl = 'DistrictCooling', zones)
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Fan coil district chilled water electric baseboard heat'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = nil, znht = nil, cl = 'DistrictCooling', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'Fan coil district hot and chilled water'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = 'DistrictHeating', znht = nil, cl = 'DistrictCooling', zones)
+
+    when 'Fan coil district hot water and chiller'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = 'DistrictHeating', znht = nil, cl = 'Electricity', zones)
+
+    when 'Fan coil chiller with no heat'
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Baseboard district hot water heat'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'DistrictHeating', znht = nil, cl = nil, zones)
+
+    when 'Baseboard district hot water heat with direct evap coolers'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'DistrictHeating', znht = nil, cl = nil, zones)
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Baseboard electric heat'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'Baseboard electric heat with direct evap coolers'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Baseboard hot water heat'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Baseboard hot water heat with direct evap coolers'
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Window AC with no heat'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Window AC with forced air furnace'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Forced Air Furnace', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Window AC with district hot water baseboard heat'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'DistrictHeating', znht = nil, cl = nil, zones)
+
+    when 'Window AC with hot water baseboard heat'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Window AC with electric baseboard heat'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'Window AC with unit heaters'
+      standard.model_add_hvac_system(self, 'Window AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Direct evap coolers'
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Direct evap coolers with unit heaters'
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Unit heaters'
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Heat pump heat with no cooling'
+      standard.model_add_hvac_system(self, 'Residential Air Source Heat Pump', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Heat pump heat with direct evap cooler'
+      # standard.model_add_hvac_system(self, 'Residential Air Source Heat Pump', ht=nil, znht=nil, cl='Electricity', zones)
+      # Using PTHP to represent zone heat pump for this configuration
+      # because only one airloop may be connected to each thermal zone.
+      standard.model_add_hvac_system(self, 'PTHP', ht = 'Electricity', znht = nil, cl = 'Electricity', zones)
+      # disable the cooling coils in all the PTHPs
+      getZoneHVACPackagedTerminalHeatPumps.each do |pthp|
+        clg_coil = pthp.heatingCoil.to_CoilHeatingDXSingleSpeed.get
+        clg_coil.setAvailabilitySchedule(alwaysOffDiscreteSchedule)
+      end
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'VAV with reheat'
+      standard.model_add_hvac_system(self, 'VAV Reheat', ht = 'NaturalGas', znht = 'NaturalGas', cl = 'Electricity', zones)
+
+    when 'VAV with PFP boxes'
+      standard.model_add_hvac_system(self, 'VAV PFP Boxes', ht = 'NaturalGas', znht = 'NaturalGas', cl = 'Electricity', zones)
+
+    when 'VAV with gas reheat'
+      standard.model_add_hvac_system(self, 'VAV Gas Reheat', ht = 'NaturalGas', ht = 'NaturalGas', cl = 'Electricity', zones)
+
+    when 'VAV with zone unit heaters'
+      standard.model_add_hvac_system(self, 'VAV No Reheat', ht = 'NaturalGas', znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'VAV with electric baseboard heat'
+      standard.model_add_hvac_system(self, 'VAV No Reheat', ht = 'NaturalGas', znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    when 'VAV cool with zone heat pump heat'
+      standard.model_add_hvac_system(self, 'VAV No Reheat', ht = 'NaturalGas', znht = nil, cl = 'Electricity', zones)
+      # standard.model_add_hvac_system(self, 'Residential Air Source Heat Pump', ht=nil, znht=nil, cl='Electricity', zones)
+      # Using PTHP to represent zone heat pump for this configuration
+      # because only one airloop may be connected to each thermal zone.
+      standard.model_add_hvac_system(self, 'PTHP', ht = 'Electricity', znht = nil, cl = 'Electricity', zones)
+      # disable the cooling coils in all the PTHPs
+      getZoneHVACPackagedTerminalHeatPumps.each do |pthp|
+        clg_coil = pthp.heatingCoil.to_CoilHeatingDXSingleSpeed.get
+        clg_coil.setAvailabilitySchedule(alwaysOffDiscreteSchedule)
+      end
+
+    when 'PVAV with reheat', 'Packaged VAV Air Loop with Boiler' # second enumeration for backwards compatibility with Tenant Star project
+      standard.model_add_hvac_system(self, 'PVAV Reheat', ht = 'NaturalGas', znht = 'NaturalGas', cl = 'Electricity', zones)
+
+    when 'PVAV with PFP boxes'
+      standard.model_add_hvac_system(self, 'PVAV PFP Boxes', ht = 'Electricity', znht = 'Electricity', cl = 'Electricity', zones)
+
+    when 'Residential forced air'
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Residential forced air cooling hot water baseboard heat'
+      standard.model_add_hvac_system(self, 'Residential AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Residential forced air with district hot water'
+      standard.model_add_hvac_system(self, 'Residential AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Residential heat pump'
+      standard.model_add_hvac_system(self, 'Residential Air Source Heat Pump', ht = 'Electricity', znht = nil, cl = 'Electricity', zones)
+
+    when 'Forced air furnace'
+      standard.model_add_hvac_system(self, 'Forced Air Furnace', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+
+    when 'Forced air furnace district chilled water fan coil'
+      standard.model_add_hvac_system(self, 'Forced Air Furnace', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+      standard.model_add_hvac_system(self, 'Fan Coil', ht = nil, znht = nil, cl = 'DistrictCooling', zones)
+
+    when 'Forced air furnace direct evap cooler'
+      # standard.model_add_hvac_system(self, 'Forced Air Furnace', ht='NaturalGas', znht=nil, cl=nil, zones)
+      # Using unit heater to represent forced air furnace for this configuration
+      # because only one airloop may be connected to each thermal zone.
+      standard.model_add_hvac_system(self, 'Unit Heaters', ht = 'NaturalGas', znht = nil, cl = nil, zones)
+      standard.model_add_hvac_system(self, 'Evaporative Cooler', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Residential AC with no heat'
+      standard.model_add_hvac_system(self, 'Residential AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+
+    when 'Residential AC with electric baseboard heat'
+      standard.model_add_hvac_system(self, 'Residential AC', ht = nil, znht = nil, cl = 'Electricity', zones)
+      standard.model_add_hvac_system(self, 'Baseboards', ht = 'Electricity', znht = nil, cl = nil, zones)
+
+    else
+      puts "HVAC system type '#{system_type}' not recognized"
+
+    end
+  end
+end