openstudio-standards 0.1.4 → 0.1.5

data/lib/openstudio-standards/hvac_sizing/Siz.CoolingTowerTwoSpeed.rb

@@ -13,15 +13,71 @@ class OpenStudio::Model::CoolingTowerTwoSpeed
   def applySizingValues
 
     OpenStudio::logFree(OpenStudio::Warn, "openstudio.sizing.CoolingTowerTwoSpeed", ".applySizingValues not yet implemented for #{self.iddObject.type.valueDescription}.")
-        
+    
+  end
+
+  # returns the autosized design water flow rate as an optional double
+  def autosizedDesignWaterFlowRate
+
+    return self.model.getAutosizedValue(self, 'Design Water Flow Rate', 'm3/s')
+
   end
+  
+  # returns the autosized fan power at high speed as an optional double
+  def autosizedHighFanSpeedFanPower
 
-  # returns the autosized design supply air flow rate as an optional double
-  # def autosizedDesignSupplyAirFlowRate
+    return self.model.getAutosizedValue(self, 'Fan Power at High Fan Speed', 'W')
+    
+  end
+  
+  # returns the autosized reference design air flow rate as an optional double
+  def autosizedHighFanSpeedAirFlowRate
 
-    # return self.model.getAutosizedValue(self, 'Design Supply Air Flow Rate', 'm3/s')
+    return self.model.getAutosizedValue(self, 'Air Flow Rate at High Fan Speed', 'm3/s')
     
-  # end
+  end
+  
+  # returns the autosized u-factor times area value at high speed as an optional double
+  def autosizedHighFanSpeedUFactorTimesAreaValue
+
+    return self.model.getAutosizedValue(self, 'U-Factor Times Area Value at High Fan Speed', 'W/C')
+
+  end
   
+  # returns the autosized reference low speed air flow rate as an optional double
+  def autosizedLowFanSpeedAirFlowRate
+
+    return self.model.getAutosizedValue(self, 'Low Fan Speed Air Flow Rate', 'm3/s')
+    
+  end  
+  
+  # returns the autosized fan power at low speed as an optional double
+  def autosizedLowFanSpeedFanPower
+
+    return self.model.getAutosizedValue(self, 'Fan Power at Low Fan Speed', 'W')
+    
+  end  
+  
+  # returns the autosized u-factor times area value at design air flow rate as an optional double
+  def autosizedLowFanSpeedUFactorTimesAreaValue
+
+    return self.model.getAutosizedValue(self, 'U-Factor Times Area Value at Low Fan Speed', 'W/K')
+
+  end  
+  
+  # returns the autosized air flow rate in free convection regime as an optional double
+  def autosizedFreeConvectionRegimeAirFlowRate
+
+    return self.model.getAutosizedValue(self, 'Free Convection Regime Air Flow Rate', 'm3/s')
+
+  end
+  
+  # returns the autosized u-factor times area value in free convection as an optional double
+  def autosizedFreeConvectionRegimeUFactorTimesAreaValue
+
+    return self.model.getAutosizedValue(self, 'Free Convection U-Factor Times Area Value', 'W/K')
+    
+  end
+    
   
 end

data/lib/openstudio-standards/hvac_sizing/Siz.CoolingTowerVariableSpeed.rb

@@ -4,7 +4,11 @@ class OpenStudio::Model::CoolingTowerVariableSpeed
 
   # Sets all auto-sizeable fields to autosize
   def autosize
-    OpenStudio::logFree(OpenStudio::Warn, "openstudio.sizing.CoolingTowerVariableSpeed", ".autosize not yet implemented for #{self.iddObject.type.valueDescription}.")
+    
+    self.autosizeDesignWaterFlowRate
+    self.autosizeDesignAirFlowRate
+    self.autosizeDesignFanPower
+
   end
 
   # Takes the values calculated by the EnergyPlus sizing routines
@@ -12,16 +16,42 @@ class OpenStudio::Model::CoolingTowerVariableSpeed
   # Must have previously completed a run with sql output for this to work.
   def applySizingValues
 
-    OpenStudio::logFree(OpenStudio::Warn, "openstudio.sizing.CoolingTowerVariableSpeed", ".applySizingValues not yet implemented for #{self.iddObject.type.valueDescription}.")
-        
+    rated_water_flow_rate = self.autosizedDesignWaterFlowRate
+    if rated_water_flow_rate.is_initialized
+      self.setDesignWaterFlowRate(rated_water_flow_rate.get) 
+    end 
+
+    rated_air_flow_rate = self.autosizedDesignAirFlowRate
+    if rated_air_flow_rate.is_initialized
+      self.setDesignAirFlowRate(rated_air_flow_rate.get) 
+    end
+
+    rated_fan_power = self.autosizedDesignFanPower
+    if rated_fan_power.is_initialized
+      self.setDesignFanPower(rated_fan_power.get) 
+    end
+    
   end
 
-  # returns the autosized design supply air flow rate as an optional double
-  # def autosizedDesignSupplyAirFlowRate
+  # returns the autosized design water flow rate as an optional double
+  def autosizedDesignWaterFlowRate
 
-    # return self.model.getAutosizedValue(self, 'Design Supply Air Flow Rate', 'm3/s')
+    return self.model.getAutosizedValue(self, 'Design Water Flow Rate', 'm3/s')
     
-  # end
+  end
+
+  # returns the autosized air flow rate as an optional double
+  def autosizedDesignAirFlowRate
+
+    return self.model.getAutosizedValue(self, 'Design Air Flow Rate', 'm3/s')
+    
+  end
   
+  # returns the autosized design fan power as an optional double
+  def autosizedDesignFanPower
+
+    return self.model.getAutosizedValue(self, 'Fan Power at Design Air Flow Rate', 'W')
+    
+  end
   
 end

data/lib/openstudio-standards/hvac_sizing/Siz.DistrictCooling.rb

@@ -0,0 +1,27 @@
+
+# open the class to add methods to return sizing values
+class OpenStudio::Model::DistrictCooling
+
+  # Sets all auto-sizeable fields to autosize
+  def autosize
+    self.autosizeNominalCapacity
+  end
+
+  # Takes the values calculated by the EnergyPlus sizing routines
+  # and puts them into this object model in place of the autosized fields.
+  # Must have previously completed a run with sql output for this to work.
+  def applySizingValues
+
+    nominal_capacity = self.autosizedMaximumFlowRate
+    if nominal_capacity.is_initialized
+      self.setMaximumFlowRate(nominal_capacity.get)
+    end
+
+  end
+  
+  # returns the autosized maximum flow rate as an optional double
+  def autosizedNominalCapacity
+    return self.model.getAutosizedValue(self, 'Design Size Nominal Capacity', 'W')
+  end
+  
+end

data/lib/openstudio-standards/hvac_sizing/Siz.DistrictHeating.rb

@@ -0,0 +1,27 @@
+
+# open the class to add methods to return sizing values
+class OpenStudio::Model::DistrictHeating
+
+  # Sets all auto-sizeable fields to autosize
+  def autosize
+    self.autosizeNominalCapacity
+  end
+
+  # Takes the values calculated by the EnergyPlus sizing routines
+  # and puts them into this object model in place of the autosized fields.
+  # Must have previously completed a run with sql output for this to work.
+  def applySizingValues
+
+    nominal_capacity = self.autosizedMaximumFlowRate
+    if nominal_capacity.is_initialized
+      self.setMaximumFlowRate(nominal_capacity.get)
+    end
+
+  end
+  
+  # returns the autosized maximum flow rate as an optional double
+  def autosizedNominalCapacity
+    return self.model.getAutosizedValue(self, 'Design Size Nominal Capacity', 'W')
+  end
+  
+end

data/lib/openstudio-standards/hvac_sizing/Siz.HeaderedPumpsConstantSpeed.rb

@@ -0,0 +1,55 @@
+
+# open the class to add methods to return sizing values
+class OpenStudio::Model::HeaderedPumpsConstantSpeed
+
+  # Sets all auto-sizeable fields to autosize
+  def autosize
+    self.autosizeTotalRatedFlowRate
+    self.autosizeRatedPowerConsumption
+  end
+
+  # Takes the values calculated by the EnergyPlus sizing routines
+  # and puts them into this object model in place of the autosized fields.
+  # Must have previously completed a run with sql output for this to work.
+  def applySizingValues
+
+    rated_flow_rate = self.autosizedTotalRatedFlowRate
+    if rated_flow_rate.is_initialized
+      self.setTotalRatedFlowRate(rated_flow_rate.get) 
+    end
+
+    rated_power_consumption = self.autosizedRatedPowerConsumption
+    if rated_power_consumption.is_initialized
+      self.setRatedPowerConsumption(rated_power_consumption.get)
+    end
+
+  end
+
+  # returns the autosized rated flow rate as an optional double
+  def autosizedTotalRatedFlowRate
+
+    # In E+ 8.5, (OS 1.10.5 onward) the column name changed
+    col_name = nil
+    if self.model.version < OpenStudio::VersionString.new('1.10.5')
+      col_name = 'Rated Flow Rate'
+    else
+      col_name = 'Design Flow Rate'
+    end
+
+    return self.model.getAutosizedValue(self, col_name, 'm3/s')
+  end
+
+  # returns the autosized rated power consumption as an optional double
+  def autosizedRatedPowerConsumption
+
+    # In E+ 8.5, (OS 1.10.5 onward) the column name changed
+    col_name = nil
+    if self.model.version < OpenStudio::VersionString.new('1.10.5')
+      col_name = 'Rated Power Consumption'
+    else
+      col_name = 'Design Power Consumption'
+    end
+
+    return self.model.getAutosizedValue(self, col_name, 'W')
+  end
+end

data/lib/openstudio-standards/hvac_sizing/Siz.HeaderedPumpsVariableSpeed.rb

@@ -0,0 +1,55 @@
+
+# open the class to add methods to return sizing values
+class OpenStudio::Model::HeaderedPumpsVariableSpeed
+
+  # Sets all auto-sizeable fields to autosize
+  def autosize
+    self.autosizeTotalRatedFlowRate
+    self.autosizeRatedPowerConsumption
+  end
+
+  # Takes the values calculated by the EnergyPlus sizing routines
+  # and puts them into this object model in place of the autosized fields.
+  # Must have previously completed a run with sql output for this to work.
+  def applySizingValues
+
+    rated_flow_rate = self.autosizedTotalRatedFlowRate
+    if rated_flow_rate.is_initialized
+      self.setTotalRatedFlowRate(rated_flow_rate.get) 
+    end
+
+    rated_power_consumption = self.autosizedRatedPowerConsumption
+    if rated_power_consumption.is_initialized
+      self.setRatedPowerConsumption(rated_power_consumption.get)
+    end
+
+  end
+
+  # returns the autosized rated flow rate as an optional double
+  def autosizedTotalRatedFlowRate
+
+    # In E+ 8.5, (OS 1.10.5 onward) the column name changed
+    col_name = nil
+    if self.model.version < OpenStudio::VersionString.new('1.10.5')
+      col_name = 'Rated Flow Rate'
+    else
+      col_name = 'Design Flow Rate'
+    end
+
+    return self.model.getAutosizedValue(self, col_name, 'm3/s')
+  end
+
+  # returns the autosized rated power consumption as an optional double
+  def autosizedRatedPowerConsumption
+
+    # In E+ 8.5, (OS 1.10.5 onward) the column name changed
+    col_name = nil
+    if self.model.version < OpenStudio::VersionString.new('1.10.5')
+      col_name = 'Rated Power Consumption'
+    else
+      col_name = 'Design Power Consumption'
+    end
+
+    return self.model.getAutosizedValue(self, col_name, 'W')
+  end
+end

data/lib/openstudio-standards/hvac_sizing/Siz.HeatingCoolingFuels.rb

@@ -11,6 +11,8 @@
 class OpenStudio::Model::Model
   
   # Get the heating fuel type of a plant loop
+  # @Todo: If no heating equipment is found, check if there's a heat exchanger,
+  # or a WaterHeater:Mixed or stratified that is connected to a heating source on the demand side
   def plant_loop_heating_fuels(plant_loop)
     fuels = []
     # Get the heating fuels for all supply components
@@ -25,7 +27,7 @@ class OpenStudio::Model::Model
       when 'OS_Boiler_Steam' 
         component = component.to_BoilerHotWater.get
         fuels << component.fuelType
-      when 'OS_District_Heating'
+      when 'OS_DistrictHeating'
         fuels << 'DistrictHeating' 
       when 'OS_HeatPump_WaterToWater_EquationFit_Heating'
         fuels << 'Electricity'
@@ -38,9 +40,37 @@ class OpenStudio::Model::Model
       when 'OS_WaterHeater_HeatPump'
         fuels << 'Electricity'     
       when 'OS_WaterHeater_Mixed'
-        fuels << 'obj.fuelType'
+        component = component.to_WaterHeaterMixed.get
+
+        # Check if the heater actually has a capacity (otherwise it's simply a Storage Tank)
+        if component.heaterMaximumCapacity.empty? || component.heaterMaximumCapacity.get != 0
+          # If it does, we add the heater Fuel Type
+          fuels << component.heaterFuelType
+        end  # @Todo: not sure about whether it should be an elsif or not
+        # Check the plant loop connection on the source side
+        if component.secondaryPlantLoop.is_initialized
+          fuels += self.plant_loop_heating_fuels(component.secondaryPlantLoop.get)
+        end
       when 'OS_WaterHeater_Stratified'
-        fuels << 'obj.fuelType'
+        component = component.to_WaterHeaterStratified.get
+
+        # Check if the heater actually has a capacity (otherwise it's simply a Storage Tank)
+        if component.heaterMaximumCapacity.empty? || component.heaterMaximumCapacity.get != 0
+          # If it does, we add the heater Fuel Type
+          fuels << component.heaterFuelType
+        end # @Todo: not sure about whether it should be an elsif or not
+        # Check the plant loop connection on the source side
+        if component.secondaryPlantLoop.is_initialized
+          fuels += self.plant_loop_heating_fuels(component.secondaryPlantLoop.get)
+        end
+
+      when 'OS_HeatExchanger_FluidToFluid'
+        hx = component.to_HeatExchangerFluidToFluid.get
+        cooling_hx_control_types = ["CoolingSetpointModulated", "CoolingSetpointOnOff", "CoolingDifferentialOnOff", "CoolingSetpointOnOffWithComponentOverride"]
+        cooling_hx_control_types.each {|x| x.downcase!}
+        if !cooling_hx_control_types.include?(hx.controlType.downcase) && hx.secondaryPlantLoop.is_initialized
+          fuels += self.plant_loop_heating_fuels(hx.secondaryPlantLoop.get)
+        end
       when 'OS_Node', 'OS_Pump_ConstantSpeed', 'OS_Pump_VariableSpeed', 'OS_Connector_Splitter', 'OS_Connector_Mixer', 'OS_Pipe_Adiabatic'
         # To avoid extraneous debug messages
       else
@@ -127,13 +157,21 @@ class OpenStudio::Model::Model
       end
     when 'OS_Coil_Heating_Water_BaseboardRadiant'
       heating_coil = heating_coil.to_CoilHeatingWaterBaseboardRadiant.get
-      if heating_coil.plantLoop.is_initialized
+        if heating_coil.plantLoop.is_initialized
         fuels += self.plant_loop_heating_fuels(heating_coil.plantLoop.get)
       end  
     when 'OS_Coil_Heating_WaterToAirHeatPump_EquationFit'
       fuels << 'Electricity'
+      heating_coil = heating_coil.to_CoilHeatingWaterToAirHeatPumpEquationFit.get
+      if heating_coil.plantLoop.is_initialized
+        fuels += self.plant_loop_heating_fuels(heating_coil.plantLoop.get)
+      end
     when 'OS_Coil_Heating_WaterToAirHeatPump_VariableSpeedEquationFit'
       fuels << 'Electricity'
+      heating_coil = heating_coil.to_CoilHeatingWaterToAirHeatPumpVariableSpeedEquationFit.get
+      if heating_coil.plantLoop.is_initialized
+        fuels += self.plant_loop_heating_fuels(heating_coil.plantLoop.get)
+      end
     when 'OS_Coil_WaterHeating_AirToWaterHeatPump'
       fuels << 'Electricity'
     when 'OS_Coil_WaterHeating_Desuperheater'
@@ -269,14 +307,15 @@ class OpenStudio::Model::Model
         fuels << 'Electricity'
       when 'OS_ZoneHVAC_PackagedTerminalAirConditioner'
         equipment = equipment.to_ZoneHVACPackagedTerminalAirConditioner.get
-        fuel_coil = self.coil_heating_fuels(equipment.heatingCoil)
         fuels += self.coil_heating_fuels(equipment.heatingCoil)
       when 'OS_ZoneHVAC_PackagedTerminalHeatPump'
         fuels << 'Electricity'
       when 'OS_ZoneHVAC_TerminalUnit_VariableRefrigerantFlow'
         fuels << 'Electricity'
       when 'OS_ZoneHVAC_WaterToAirHeatPump'
-        fuels << 'Electricity'
+        # We also go check what fuel serves the loop on which the WSHP heating coil is
+        equipment = equipment.to_ZoneHVACWaterToAirHeatPump.get
+        fuels += self.coil_heating_fuels(equipment.heatingCoil)
       else
         #OpenStudio::logFree(OpenStudio::Debug, 'openstudio.sizing.Model', "No heating fuel types found for #{obj_type}")
       end
@@ -294,14 +333,17 @@ class OpenStudio::Model::Model
       # Get the object type
       obj_type = equipment.iddObjectType.valueName.to_s
       case obj_type    
-      when 'to_AirTerminal_SingleDuct_ConstantVolume_CooledBeam'
+      when 'OS_AirTerminal_SingleDuct_ConstantVolume_CooledBeam'
         equipment = equipment.to_AirTerminalSingleDuctConstantVolumeCooledBeam.get
         fuels += self.coil_cooling_fuels(equipment.coilCoolingCooledBeam)
-      when 'to_AirTerminal_SingleDuct_ConstantVolume_FourPipeInduction'      
+      when 'OS_AirTerminal_SingleDuct_ConstantVolume_FourPipeInduction'      
         equipment = equipment.to_AirTerminalSingleDuctConstantVolumeFourPipeInduction.get
         if equipment.coolingCoil.is_initialized
-          fuels += self.coil_heating_fuels(equipment.coolingCoil.get) 
+          fuels += self.coil_cooling_fuels(equipment.coolingCoil.get) 
         end
+      when 'OS_ZoneHVAC_FourPipeFanCoil'
+        equipment = equipment.to_ZoneHVACFourPipeFanCoil.get
+        fuels += self.coil_cooling_fuels(equipment.coolingCoil)
       when 'OS_Refrigeration_AirChiller'
         fuels << 'Electricity'
       when 'OS_ZoneHVAC_IdealLoadsAirSystem'

data/lib/openstudio-standards/hvac_sizing/Siz.Model.rb

@@ -24,12 +24,15 @@ class OpenStudio::Model::Model
   require_relative 'Siz.CoilHeatingGas'
   require_relative 'Siz.CoilHeatingWater'
   require_relative 'Siz.CoilHeatingDXSingleSpeed'
+  require_relative 'Siz.CoilHeatingDXMultiSpeed'
   require_relative 'Siz.CoilHeatingWaterToAirHeatPumpEquationFit'
   require_relative 'Siz.CoilCoolingWaterToAirHeatPumpEquationFit'
   require_relative 'Siz.CoilCoolingDXSingleSpeed'
   require_relative 'Siz.CoilCoolingDXTwoSpeed'
   require_relative 'Siz.CoilCoolingWater'
   require_relative 'Siz.ControllerOutdoorAir'
+  require_relative 'Siz.DistrictHeating'
+  require_relative 'Siz.DistrictCooling'
   require_relative 'Siz.HeatExchangerAirToAirSensibleAndLatent'
   require_relative 'Siz.PlantLoop'
   require_relative 'Siz.PumpConstantSpeed'
@@ -37,26 +40,19 @@ class OpenStudio::Model::Model
   require_relative 'Siz.BoilerHotWater'
   require_relative 'Siz.ChillerElectricEIR'
   require_relative 'Siz.CoolingTowerSingleSpeed'
+  require_relative 'Siz.CoolingTowerTwoSpeed'
+  require_relative 'Siz.CoolingTowerVariableSpeed'
   require_relative 'Siz.ControllerWaterCoil'
   require_relative 'Siz.SizingSystem'
   require_relative 'Siz.ThermalZone'
-
-  # Recently added and not fully tested
   require_relative 'Siz.ZoneHVACPackagedTerminalAirConditioner'
   require_relative 'Siz.ZoneHVACPackagedTerminalHeatPump'
   require_relative 'Siz.ZoneHVACTerminalUnitVariableRefrigerantFlow'
   require_relative 'Siz.AirConditionerVariableRefrigerantFlow'
   require_relative 'Siz.CoilCoolingDXVariableRefrigerantFlow'
   require_relative 'Siz.CoilHeatingDXVariableRefrigerantFlow'
-
-  # Methods not yet implemented
-  require_relative 'Siz.AirTermSnglDuctVAVReheat'
-  require_relative 'Siz.AirTermSnglDuctUncontrolled'
-  require_relative 'Siz.AirLoopHVAC'
-  require_relative 'Siz.AirLoopHVACUnitaryHeatPumpAirToAir'
-  require_relative 'Siz.FanConstantVolume'
-  require_relative 'Siz.FanVariableVolume'
-  require_relative 'Siz.FanOnOff'  
+  require_relative 'Siz.HeaderedPumpsConstantSpeed'
+  require_relative 'Siz.HeaderedPumpsVariableSpeed'
 
   # Heating and cooling fuel methods
   require_relative 'Siz.HeatingCoolingFuels'
@@ -65,21 +61,107 @@ class OpenStudio::Model::Model
   # autosizing back into the self.
   def runSizingRun(sizing_run_dir = "#{Dir.pwd}/SizingRun")
 
+
     # Change the simulation to only run the sizing days
     sim_control = self.getSimulationControl
     sim_control.setRunSimulationforSizingPeriods(true)
     sim_control.setRunSimulationforWeatherFileRunPeriods(false)
-    
+
+    #check that all zones have surfaces. 
+    raise ("Error: Sizing Run Failed. Thermal Zones with no surfaces exist.") unless self.do_all_zones_have_surfaces?
     # Run the sizing run
-    self.run_simulation_and_log_errors(sizing_run_dir)
+    success = self.run_simulation_and_log_errors(sizing_run_dir)
     
     # Change the model back to running the weather file
     sim_control.setRunSimulationforSizingPeriods(false)
     sim_control.setRunSimulationforWeatherFileRunPeriods(true)
     
-    return true
+    return success
 
   end
+  
+  #Method to check if all zones have surfaces. This is required to run a simulation. 
+  def do_all_zones_have_surfaces?()
+    error_string = ""
+    error = false
+    #Check to see if all zones have surfaces. 
+    self.getThermalZones.each do |zone|
+      if  0 == BTAP::Geometry::Surfaces::get_surfaces_from_thermal_zones([zone]).size
+        error_string << "Error: Thermal zone #{zone.name} does not contain surfaces.\n"
+        error = true
+      end
+      if error == true
+        puts error_string
+        OpenStudio::logFree(OpenStudio::Error, 'openstudio.Siz.Model', error_string)
+        return false
+      else
+        return true
+      end
+    end
+  end
+  
+  # A helper method to run a sizing run and pull any values calculated during
+  # autosizing back into the self.
+  def runSpaceSizingRun(sizing_run_dir = "#{Dir.pwd}/SpaceSizingRun")
+    puts "*************Runing sizing space Run ***************************"
+    #Make copy of model
+    model = BTAP::FileIO::deep_copy(self, true)
+    space_load_array = []
+    
+    #Make sure the model is good to run. 
+    #1. Ensure External surfaces are set to a construction
+    ext_surfaces =  BTAP::Geometry::Surfaces::filter_by_boundary_condition(model.getSurfaces, ["Outdoors",
+        "Ground",
+        "GroundFCfactorMethod",
+        "GroundSlabPreprocessorAverage",
+        "GroundSlabPreprocessorCore",
+        "GroundSlabPreprocessorPerimeter",
+        "GroundBasementPreprocessorAverageWall",
+        "GroundBasementPreprocessorAverageFloor",
+        "GroundBasementPreprocessorUpperWall",
+        "GroundBasementPreprocessorLowerWall"])
+    fail = false
+    ext_surfaces.each do |surface|
+      if surface.construction.nil?
+        OpenStudio::logFree(OpenStudio::Error,'openstudio.model.Model', "Ext Surface #{surface.name} does not have a construction.Cannot perform sizing.") 
+        fail = true
+      end
+    end
+    puts "#{ext_surfaces.size} External Surfaces counted."
+    raise ("Can't run sizing since envelope is not set.") if fail == true
+    
+    #remove any thermal zones. 
+    model.getThermalZones.each {|zone| zone.remove}
+    
+    #assign a zone to each space. 
+    # Create a thermal zone for each space in the self
+    model.getSpaces.each do |space|
+      zone = OpenStudio::Model::ThermalZone.new(self)
+      zone.setName("#{space.name} ZN")
+      space.setThermalZone(zone)
+    end
+    # Add a thermostat
+    BTAP::Compliance::NECB2011::set_zones_thermostat_schedule_based_on_space_type_schedules(model)
+    OpenStudio::logFree(OpenStudio::Info, 'openstudio.model.Model', 'Finished creating thermal zones')
+    # Add ideal loads to every zone/space and run
+    # a sizing run to determine heating/cooling loads,
+    # which will impact HVAC systems.
+    model.getThermalZones.each do |zone|
+      ideal_loads = OpenStudio::Model::ZoneHVACIdealLoadsAirSystem.new(model)
+      ideal_loads.addToThermalZone(zone)
+    end
+    self.runSizingRun(sizing_run_dir)
+    model.getSpaces.each do |space|
+      if not space.thermalZone.empty?
+        space_load_array << {"space_name"=> space.name, "CoolingDesignLoad" => space.thermalZone.get.coolingDesignLoad, "HeatingDesignLoad" => space.thermalZone.get.heatingDesignLoad }
+      end
+    end
+    puts space_load_array
+    puts "*************Done Runing sizing space Run ***************************"
+    return model
+  end
+  
+  
 
   # Takes the values calculated by the EnergyPlus sizing routines
   # and puts them into all objects model in place of the autosized fields.
@@ -302,13 +384,13 @@ class OpenStudio::Model::Model
     name = object.name.get.upcase
 
     object_type = object.iddObject.type.valueDescription.gsub('OS:','')
-    if(object_type == 'Coil:Heating:Gas:MultiStage') then object_type = 'Coil:Heating:GasMultiStage' end
 
     # Special logic for two coil types which are inconsistently
     # uppercase in the sqlfile:
     object_type = object_type.upcase if object_type == 'Coil:Cooling:WaterToAirHeatPump:EquationFit'
     object_type = object_type.upcase if object_type == 'Coil:Heating:WaterToAirHeatPump:EquationFit'
-      
+		object_type = 'Coil:Heating:GasMultiStage' if object_type == 'Coil:Heating:Gas:MultiStage'
+
     sql = self.sqlFile
 
     if sql.is_initialized
@@ -426,7 +508,7 @@ class OpenStudio::Model::Model
     end
 
     # If a csv path is given, output
-    if !csv_path.nil?
+    if !csv_path.nil? && !table.first.nil?
       CSV.open(csv_path, "wb") do |csv|
         csv << table.first.keys # adds the attributes name on the first line
         table.each do |hash|