openstudio-standards 0.1.0 → 0.1.1

data/lib/openstudio-standards/prototypes/Prototype.hvac_systems.rb

@@ -230,9 +230,10 @@ class OpenStudio::Model::Model
     if condenser_water_loop
       condenser_water_loop.addDemandBranchForComponent(chiller)
       chiller.setCondenserType('WaterCooled')
-    end
+	end
+	
 
-    #chilled water loop pipes
+	#chilled water loop pipes
     chiller_bypass_pipe = OpenStudio::Model::PipeAdiabatic.new(self)
     chilled_water_loop.addSupplyBranchForComponent(chiller_bypass_pipe)
     coil_bypass_pipe = OpenStudio::Model::PipeAdiabatic.new(self)
@@ -622,7 +623,6 @@ class OpenStudio::Model::Model
       oa_intake_controller.setEconomizerControlType("DifferentialEnthalpy")
       oa_intake_controller.setHeatRecoveryBypassControlType("BypassWhenOAFlowGreaterThanMinimum")
       oa_intake_controller.resetMaximumFractionofOutdoorAirSchedule
-      oa_intake_controller.resetMaximumFractionofOutdoorAirSchedule
       oa_intake_controller.resetEconomizerMinimumLimitDryBulbTemperature
     end
 
@@ -650,9 +650,8 @@ class OpenStudio::Model::Model
       terminal = OpenStudio::Model::AirTerminalSingleDuctVAVReheat.new(self,self.alwaysOnDiscreteSchedule,rht_coil)
       terminal.setName("#{zone.name} VAV Term")
       terminal.setZoneMinimumAirFlowMethod('Constant')
-
-      terminal.set_initial_prototype_damper_position(standard, zone.outdoor_airflow_rate_per_area)
-
+      
+      terminal.set_initial_prototype_damper_position(building_type, standard, zone.outdoor_airflow_rate_per_area)
       terminal.setMaximumFlowPerZoneFloorAreaDuringReheat(0.0)
       terminal.setMaximumFlowFractionDuringReheat(0.5)
       terminal.setMaximumReheatAirTemperature(rht_sa_temp_c)
@@ -873,7 +872,8 @@ class OpenStudio::Model::Model
               hvac_op_sch,
               oa_damper_sch,
               hot_water_loop = nil, 
-              return_plenum = nil)
+              return_plenum = nil,
+              building_type = nil)
 
     OpenStudio::logFree(OpenStudio::Info, 'openstudio.Model.Model', "Adding Packaged VAV for #{thermal_zones.size} zones.")
     thermal_zones.each do |zone|
@@ -911,7 +911,11 @@ class OpenStudio::Model::Model
     zn_dsn_htg_sa_temp_f = 122 # Design VAV box to reheat to 122F
     rht_rated_air_in_temp_f = 62 # Reheat coils designed to receive 62F
     rht_rated_air_out_temp_f = 90 # Reheat coils designed to supply 90F...but zone expects 122F...?
-    clg_sa_temp_f = 55 # Central deck clg temp operates at 55F
+    if sys_name == 'PVAV Outpatient F1'
+      clg_sa_temp_f = 52 # for AHU1 in Outpatient, SAT is 52F
+    else
+      clg_sa_temp_f = 55 # Central deck clg temp operates at 55F
+    end
 
     sys_dsn_prhtg_temp_c = OpenStudio.convert(sys_dsn_prhtg_temp_f,'F','C').get
     sys_dsn_clg_sa_temp_c = OpenStudio.convert(sys_dsn_clg_sa_temp_f,'F','C').get
@@ -985,6 +989,7 @@ class OpenStudio::Model::Model
     oa_intake.addToNode(air_loop.supplyInletNode)
     controller_mv = oa_intake_controller.controllerMechanicalVentilation
     controller_mv.setName("#{air_loop.name} Ventilation Controller")
+    controller_mv.setAvailabilitySchedule(oa_damper_sch)
 
     # Hook the VAV system to each zone
     thermal_zones.each do |zone|
@@ -1009,9 +1014,7 @@ class OpenStudio::Model::Model
       terminal = OpenStudio::Model::AirTerminalSingleDuctVAVReheat.new(self,self.alwaysOnDiscreteSchedule,rht_coil)
       terminal.setName("#{zone.name} VAV Term")
       terminal.setZoneMinimumAirFlowMethod('Constant')
-
-      terminal.set_initial_prototype_damper_position(standard, zone.outdoor_airflow_rate_per_area)
-
+      terminal.set_initial_prototype_damper_position(building_type, standard, zone.outdoor_airflow_rate_per_area)
       air_loop.addBranchForZone(zone,terminal.to_StraightComponent)
 
       unless return_plenum.nil?
@@ -1188,7 +1191,7 @@ class OpenStudio::Model::Model
       terminal = OpenStudio::Model::AirTerminalSingleDuctVAVReheat.new(self,self.alwaysOnDiscreteSchedule,rht_coil)
       terminal.setName("#{zone.name} VAV Term")
       terminal.setZoneMinimumAirFlowMethod('Constant')
-      terminal.set_initial_prototype_damper_position(standard, zone.outdoor_airflow_rate_per_area)
+      terminal.set_initial_prototype_damper_position(building_type, standard, zone.outdoor_airflow_rate_per_area)
       terminal.setMaximumFlowPerZoneFloorAreaDuringReheat(0.0)
       terminal.setMaximumFlowFractionDuringReheat(0.5)
       terminal.setMaximumReheatAirTemperature(rht_sa_temp_c)
@@ -1743,7 +1746,15 @@ class OpenStudio::Model::Model
       oa_controller.setHeatRecoveryBypassControlType('BypassWhenOAFlowGreaterThanMinimum')
       oa_system = OpenStudio::Model::AirLoopHVACOutdoorAirSystem.new(self,oa_controller)
       oa_system.setName("#{air_loop.name} OA Sys")
-
+      econ_eff_sch = self.add_schedule('RetailStandalone PSZ_Econ_MaxOAFrac_Sch')
+      
+      case standard
+      when '90.1-2004','90.1-2007','90.1-2010','90.1-2013'
+        oa_controller.setMaximumFractionofOutdoorAirSchedule(econ_eff_sch) if building_type == "RetailStandalone" || building_type == "RetailStripmall"
+      when 'DOE Ref Pre-1980','DOE Ref 1980-2004'
+        OpenStudio::logFree(OpenStudio::Info, 'openstudio.model.Model', 'No maximum fraction outdoor air schedule in PSZ for building types except RetailStandalone')
+      end
+      
       # Add the components to the air loop
       # in order from closest to zone to furthest from zone
       supply_inlet_node = air_loop.supplyInletNode
@@ -3686,7 +3697,7 @@ class OpenStudio::Model::Model
 
   end
 
-  def add_swh_end_uses_by_space(building_type, building_vintage, climate_zone, swh_loop, space_type_name, space_name, space_multiplier = nil)
+  def add_swh_end_uses_by_space(building_type, building_vintage, climate_zone, swh_loop, space_type_name, space_name, space_multiplier = nil, is_flow_per_area = true)
 
     # find the specific space_type properties from standard.json
     search_criteria = {
@@ -3708,7 +3719,11 @@ class OpenStudio::Model::Model
     # Water fixture definition
     water_fixture_def = OpenStudio::Model::WaterUseEquipmentDefinition.new(self)
     rated_flow_rate_per_area = data['service_water_heating_peak_flow_per_area'].to_f   # gal/h.ft2
-    rated_flow_rate_gal_per_hour = rated_flow_rate_per_area * space_area * space_multiplier   # gal/h
+    if is_flow_per_area
+      rated_flow_rate_gal_per_hour = rated_flow_rate_per_area * space_area * space_multiplier   # gal/h
+    else
+      rated_flow_rate_gal_per_hour = data['service_water_heating_peak_flow_rate'].to_f
+    end
     rated_flow_rate_gal_per_min = rated_flow_rate_gal_per_hour/60  # gal/h to gal/min
     rated_flow_rate_m3_per_s = OpenStudio.convert(rated_flow_rate_gal_per_min,'gal/min','m^3/s').get
     # water_use_sensible_frac_sch = OpenStudio::Model::ScheduleConstant.new(self)
@@ -3875,7 +3890,7 @@ class OpenStudio::Model::Model
     fan.setFanEfficiency(0.58175)
     fan.setPressureRise(622.5) #Pa
     if fan_max_flow_rate != nil
-      fan.setMaximumFlowRate(fan_max_flow_rate)
+      fan.setMaximumFlowRate(OpenStudio.convert(fan_max_flow_rate, 'cfm', 'm^3/s').get)  # unit of fan_max_flow_rate is cfm
     else
       fan.autosizeMaximumFlowRate
     end
@@ -4143,6 +4158,72 @@ class OpenStudio::Model::Model
 
   end
 
+  # Adds a zone ventilation design flow rate to each zone.
+  #
+  # @param availability_sch_name [String] the name of the fan availability schedule
+  # @param flow_rate [Double] the ventilation design flow rate in m^3/s for Exhaust/Natural or
+  # Flow Rate per Zone Floor Area in m^3/s-m^2 for Intake
+  # @param ventilation_type [String] the zone ventilation type either Exhaust, Natural, or Intake
+  # @param thermal_zones [Array<OpenStudio::Model::ThermalZone>] an array of thermal zones
+  # @return [Array<OpenStudio::Model::ZoneVentilationDesignFlowRate>] an array of zone ventilation objects created
+  def add_zone_ventilation(availability_sch_name,
+                           flow_rate,
+                           ventilation_type,
+                           thermal_zones)
+
+    # Make an exhaust fan for each zone
+    zone_ventilations = []
+    thermal_zones.each do |zone|
+      ventilation = OpenStudio::Model::ZoneVentilationDesignFlowRate.new(self)
+      ventilation.setName("#{zone.name} Ventilation")
+      ventilation.setSchedule(self.add_schedule(availability_sch_name))
+      ventilation.setVentilationType(ventilation_type)
+
+      ventilation.setAirChangesperHour(0)
+      ventilation.setTemperatureTermCoefficient(0)
+
+      if ventilation_type == 'Exhaust'
+        ventilation.setDesignFlowRateCalculationMethod("Flow/Zone")
+        ventilation.setDesignFlowRate(flow_rate)
+        ventilation.setFanPressureRise(31.1361206455786)
+        ventilation.setFanTotalEfficiency(0.51)
+        ventilation.setConstantTermCoefficient(1)
+        ventilation.setVelocityTermCoefficient(0)
+        ventilation.setMinimumIndoorTemperature(29.4444452244559)
+        ventilation.setMaximumIndoorTemperature(100)
+        ventilation.setDeltaTemperature(-100)
+      elsif ventilation_type == 'Natural'
+        ventilation.setDesignFlowRateCalculationMethod("Flow/Zone")
+        ventilation.setDesignFlowRate(flow_rate)
+        ventilation.setFanPressureRise(0)
+        ventilation.setFanTotalEfficiency(1)
+        ventilation.setConstantTermCoefficient(0)
+        ventilation.setVelocityTermCoefficient(0.224)
+        ventilation.setMinimumIndoorTemperature(-73.3333352760033)
+        ventilation.setMaximumIndoorTemperature(29.4444452244559)
+        ventilation.setDeltaTemperature(-100)
+      elsif ventilation_type == 'Intake'
+        ventilation.setDesignFlowRateCalculationMethod("Flow/Area")
+        ventilation.setFlowRateperZoneFloorArea(flow_rate)
+        ventilation.setFanPressureRise(49.8)
+        ventilation.setFanTotalEfficiency(0.53625)
+        ventilation.setConstantTermCoefficient(1)
+        ventilation.setVelocityTermCoefficient(0)
+        ventilation.setMinimumIndoorTemperature(7.5)
+        ventilation.setMaximumIndoorTemperature(35)
+        ventilation.setDeltaTemperature(-27.5)
+        ventilation.setMinimumOutdoorTemperature(-30.0)
+        ventilation.setMaximumOutdoorTemperature(50.0)
+        ventilation.setMaximumWindSpeed(6.0)
+      end      
+      ventilation.addToThermalZone(zone)
+      zone_ventilations << ventilation
+    end
+
+    return zone_ventilations
+
+  end
+
   # Adds a single refrigerated case connected to a rack composed
   # of a single compressor and a single air-cooled condenser.
   #

data/lib/openstudio-standards/prototypes/Prototype.outpatient.rb

@@ -15,16 +15,17 @@ class OpenStudio::Model::Model
       'CleanWork' => ['Floor 1 Clean', 'Floor 1 Clean Work', ],
       'Conference' => ['Floor 2 Conference'],
       'DressingRoom' => ['Floor 1 Dressing Room', 'Floor 3 Dressing Room'],
-      'Elec/MechRoom' => ['Floor 1 Electrical Room', 'Floor 3 Mechanical', 'NW Elevator'],
+      'Elec/MechRoom' => ['Floor 1 Electrical Room', 'Floor 3 Mechanical'],
       'ElevatorPumpRoom' => ['Floor 1 Elevator Pump Room'],
       # 'Floor 3 Treatment' same as 'Exam'
       'Exam' => ['Floor 2 Exam 1', 'Floor 2 Exam 2', 'Floor 2 Exam 3', 'Floor 2 Exam 4', 'Floor 2 Exam 5', 'Floor 2 Exam 6', 'Floor 2 Exam 7', 
         'Floor 2 Exam 8', 'Floor 2 Exam 9', 'Floor 3 Treatment'],
       # 'Floor 1 Scrub', 'Floor 1 Sub-Sterile', 'Floor 1 Vestibule' same as 'Hall'
       'Hall' => ['Floor 1 IT Hall', 'Floor 1 Lobby Hall', 'Floor 1 Locker Room Hall', 'Floor 1 MRI Hall', 'Floor 1 Nurse Hall', 'Floor 1 Pre-Op Hall',
-        'Floor 1 Reception Hall', 'Floor 1 Sterile Hall', 'Floor 1 Utility Hall', 'Floor 2 Exam Hall 1', 'Floor 2 Exam Hall 2', 'Floor 2 Exam Hall 3',
+        'Floor 1 Reception Hall', 'Floor 1 Sterile Hall', 'Floor 2 Exam Hall 1', 'Floor 2 Exam Hall 2', 'Floor 2 Exam Hall 3',
         'Floor 2 Exam Hall 4', 'Floor 2 Exam Hall 5', 'Floor 2 Exam Hall 6', 'Floor 2 Office Hall', 'Floor 2 Reception Hall', 'Floor 2 Work Hall', 
-        'Floor 3 Elevator Hall', 'Floor 3 Mechanical Hall', 'Floor 3 Office Hall', 'Floor 1 Scrub', 'Floor 1 Sub-Sterile', 'Floor 1 Vestibule'],
+        'Floor 3 Mechanical Hall', 'Floor 1 Scrub'],
+      'Hall_infil' => ['Floor 1 Utility Hall', 'Floor 1 Sub-Sterile', 'Floor 1 Vestibule', 'Floor 3 Elevator Hall', 'Floor 3 Office Hall'],
       'IT_Room' => ['Floor 1 IT Room'],
       # ['Floor 1 Sterile Storage', 'Floor 1 Storage', 'Floor 1 Utility Room', 'Floor 2 Storage 1', 'Floor 2 Storage 2', 'Floor 2 Storage 3', ...
       # 'Floor 2 Utility', 'Floor 3 Storage 1', 'Floor 3 Storage 2', 'Floor 3 Utility'] same as 'Janitor'
@@ -51,7 +52,7 @@ class OpenStudio::Model::Model
       'PreOp' => ['Floor 1 Pre-Op Room 1', 'Floor 1 Pre-Op Room 2'],
       'ProcedureRoom' => ['Floor 1 Procedure Room'],
       'Soil Work' => ['Floor 1 Soil', 'Floor 1 Soil Hold', 'Floor 1 Soil Work'],
-      'Stair' => ['NE Stair', 'NW Stair', 'SW Stair'],
+      'Stair' => ['NE Stair', 'NW Stair', 'SW Stair', 'NW Elevator'],
       'Toilet' => ['Floor 1 Nurse Toilet', 'Floor 1 Pre-Op Toilet', 'Floor 1 Lobby Toilet', 'Floor 1 MRI Toilet', 'Floor 2 Conference Toilet', 
         'Floor 2 Reception Toilet', 'Floor 2 Work Toilet', 'Floor 3 Lounge Toilet', 'Floor 3 Office Toilet', 'Floor 3 Physical Therapy Toilet'],
       'Xray' => ['Floor 2 X-Ray'],
@@ -140,6 +141,9 @@ class OpenStudio::Model::Model
               'Floor 3 Treatment', 'Floor 3 Work'
             ]
         }
+        # {
+          # TODO exhaust fans
+        # }
       ]      
     end
     
@@ -155,6 +159,28 @@ class OpenStudio::Model::Model
 
     # add elevator for the elevator pump room (the fan&lights are already added via standard spreadsheet)
     self.add_extra_equip_elevator_pump_room(building_vintage)
+    # adjust cooling setpoint at vintages 1B,2B,3B
+    self.adjust_clg_setpoint(building_vintage,climate_zone)
+    # Get the hot water loop
+    hot_water_loop = nil
+    self.getPlantLoops.each do |loop|
+      # If it has a boiler:hotwater, it is the correct loop
+      if loop.supplyComponents('Boiler:HotWater'.to_IddObjectType).size > 0
+        hot_water_loop = loop
+      end
+    end
+    # add humidifier to AHU1 (contains operating room 1)
+    if hot_water_loop
+      self.add_humidifier(building_vintage, hot_water_loop)
+    else
+      OpenStudio::logFree(OpenStudio::Warn, 'openstudio.model.Model', 'Could not find hot water loop to attach humidifier to.')
+    end
+    # adjust infiltration for vintages 'DOE Ref Pre-1980', 'DOE Ref 1980-2004'
+    self.adjust_infiltration(building_vintage)
+    # add door infiltration for vertibule
+    self.add_door_infiltration(building_vintage,climate_zone)
+    # reset boiler sizing factor to 0.3 (default 1)
+    self.reset_boiler_sizing_factor
     
     OpenStudio::logFree(OpenStudio::Info, 'openstudio.model.Model', 'Finished adding HVAC')
     
@@ -179,12 +205,197 @@ class OpenStudio::Model::Model
     when 'DOE Ref Pre-1980', 'DOE Ref 1980-2004'
       elec_equip.setSchedule(add_schedule("OutPatientHealthCare BLDG_ELEVATORS_Pre2004"))
     end
+    return true
   end
 
-  def custom_swh_tweaks(building_type, building_vintage, climate_zone, prototype_input)
+  def adjust_clg_setpoint(building_vintage,climate_zone)
+    self.getSpaceTypes.sort.each do |space_type|
+      space_type_name = space_type.name.get
+      thermostat_name = space_type_name + ' Thermostat'
+      thermostat = self.getThermostatSetpointDualSetpointByName(thermostat_name).get
+      case building_vintage
+      when '90.1-2004', '90.1-2007', '90.1-2010'
+        case climate_zone
+        when 'ASHRAE 169-2006-2B', 'ASHRAE 169-2006-1B', 'ASHRAE 169-2006-3B'
+          thermostat.setCoolingSetpointTemperatureSchedule(add_schedule("OutPatientHealthCare CLGSETP_SCH_YES_OPTIMUM"))
+        end
+      end
+    end
+    return true
+  end
+  
+  def adjust_infiltration(building_vintage)
+    case building_vintage
+    when 'DOE Ref Pre-1980', 'DOE Ref 1980-2004'
+      self.getSpaces.sort.each do |space|
+        space_type = space.spaceType.get
+        if space.exterior_wall_and_window_area <= 0
+          next
+        elsif space_type.spaceInfiltrationDesignFlowRates.size <= 0
+          next
+        else
+          # get the infiltration information from the space type infiltration
+          infiltration_space_type = space_type.spaceInfiltrationDesignFlowRates[0]
+          infil_sch = infiltration_space_type.schedule.get
+          infil_rate = nil
+          infil_ach = nil
+          if infiltration_space_type.flowperExteriorWallArea.is_initialized
+            infil_rate = infiltration_space_type.flowperExteriorWallArea.get
+          elsif infiltration_space_type.airChangesperHour.is_initialized
+            infil_ach = infiltration_space_type.airChangesperHour.get
+          end
+          # Create an infiltration rate object for this space
+          infiltration = OpenStudio::Model::SpaceInfiltrationDesignFlowRate.new(self)
+          infiltration.setName("#{space.name} Infiltration")
+          infiltration.setFlowperExteriorSurfaceArea(infil_rate) unless infil_rate.nil? || infil_rate.to_f == 0
+          infiltration.setAirChangesperHour(infil_ach) unless infil_ach.nil? || infil_ach.to_f == 0
+          infiltration.setSchedule(infil_sch)
+          infiltration.setSpace(space)
+        end
+      end
+      self.getSpaceTypes.each do |space_type|
+        space_type.spaceInfiltrationDesignFlowRates.each do |infil|
+          infil.remove
+        end
+      end
+    else
+      return true
+    end
+  end  
 
+  def add_door_infiltration(building_vintage,climate_zone)
+    # add extra infiltration for vestibule door
+    case building_vintage
+    when 'DOE Ref 1980-2004', 'DOE Ref Pre-1980'
+      return true
+    else
+      vestibule_space = self.getSpaceByName('Floor 1 Vestibule').get
+      infiltration_vestibule_door = OpenStudio::Model::SpaceInfiltrationDesignFlowRate.new(self)
+      infiltration_vestibule_door.setName("Vestibule door Infiltration")
+      infiltration_rate_vestibule_door = 0
+      case building_vintage
+      when '90.1-2004'
+        infiltration_rate_vestibule_door = 1.186002811
+        infiltration_vestibule_door.setSchedule(add_schedule('OutPatientHealthCare INFIL_Door_Opening_SCH_0.144'))
+      when '90.1-2007', '90.1-2010', '90.1-2013'
+        case climate_zone
+        when 'ASHRAE 169-2006-1A', 'ASHRAE 169-2006-2A', 'ASHRAE 169-2006-2B'
+          infiltration_rate_vestibule_door = 1.186002811
+          infiltration_vestibule_door.setSchedule(add_schedule('OutPatientHealthCare INFIL_Door_Opening_SCH_0.144'))
+        else
+          infiltration_rate_vestibule_door = 0.776824762
+          infiltration_vestibule_door.setSchedule(add_schedule('OutPatientHealthCare INFIL_Door_Opening_SCH_0.131'))
+        end
+      end
+      infiltration_vestibule_door.setDesignFlowRate(infiltration_rate_vestibule_door)
+      infiltration_vestibule_door.setSpace(vestibule_space)
+    end
+  end
+
+  def update_waterheater_loss_coefficient(building_vintage)
+    case building_vintage
+    when '90.1-2004', '90.1-2007', '90.1-2010', '90.1-2013'
+      self.getWaterHeaterMixeds.sort.each do |water_heater|
+        if water_heater.name.to_s.include?("Booster")
+          water_heater.setOffCycleLossCoefficienttoAmbientTemperature(1.053159296)
+          water_heater.setOnCycleLossCoefficienttoAmbientTemperature(1.053159296)
+        else
+          water_heater.setOffCycleLossCoefficienttoAmbientTemperature(9.643286505)
+          water_heater.setOnCycleLossCoefficienttoAmbientTemperature(9.643286505)
+        end
+      end
+    end      
+  end
+  
+  def custom_swh_tweaks(building_type, building_vintage, climate_zone, prototype_input)
+    
+    self.update_waterheater_loss_coefficient(building_vintage)
+  
     return true
     
-  end    
+  end
+
+  # add humidifier to AHU1 (contains operating room1)
+  def add_humidifier(building_vintage, hot_water_loop)
+    operatingroom1_space = self.getSpaceByName('Floor 1 Operating Room 1').get
+    operatingroom1_zone = operatingroom1_space.thermalZone.get
+    humidistat = OpenStudio::Model::ZoneControlHumidistat.new(self)
+    humidistat.setHumidifyingRelativeHumiditySetpointSchedule(self.add_schedule('OutPatientHealthCare MinRelHumSetSch'))
+    humidistat.setDehumidifyingRelativeHumiditySetpointSchedule(self.add_schedule('OutPatientHealthCare MaxRelHumSetSch'))
+    operatingroom1_zone.setZoneControlHumidistat(humidistat)
+    self.getAirLoopHVACs.each do |air_loop|
+      if air_loop.thermalZones.include? operatingroom1_zone
+        humidifier = OpenStudio::Model::HumidifierSteamElectric.new(self)
+        humidifier.setRatedCapacity(3.72E-5)
+        humidifier.setRatedPower(100000)
+        humidifier.setName("#{air_loop.name.get} Electric Steam Humidifier")
+        # get the water heating coil and add humidifier to the outlet of heating coil (right before fan)
+        htg_coil = nil
+        air_loop.supplyComponents.each do |equip|
+          if equip.to_CoilHeatingWater.is_initialized
+            htg_coil = equip.to_CoilHeatingWater.get
+          end
+        end
+        heating_coil_outlet_node = htg_coil.airOutletModelObject().get.to_Node.get
+        supply_outlet_node = air_loop.supplyOutletNode
+        humidifier.addToNode(heating_coil_outlet_node)
+        humidity_spm = OpenStudio::Model::SetpointManagerSingleZoneHumidityMinimum.new(self)
+        case building_vintage
+        when '90.1-2004', '90.1-2007', '90.1-2010', '90.1-2013'
+          extra_elec_htg_coil = OpenStudio::Model::CoilHeatingElectric.new(self,self.alwaysOnDiscreteSchedule)
+          extra_elec_htg_coil.setName("AHU1 extra Electric Htg Coil")
+          extra_water_htg_coil = OpenStudio::Model::CoilHeatingWater.new(self,self.alwaysOnDiscreteSchedule)
+          extra_water_htg_coil.setName("AHU1 extra Water Htg Coil")
+          hot_water_loop.addDemandBranchForComponent(extra_water_htg_coil)
+          extra_elec_htg_coil.addToNode(supply_outlet_node)
+          extra_water_htg_coil.addToNode(supply_outlet_node)
+        end
+        # humidity_spm.addToNode(supply_outlet_node)
+        humidity_spm.addToNode(humidifier.outletModelObject().get.to_Node.get)
+        humidity_spm.setControlZone(operatingroom1_zone)
+      end
+    end
+  end
+  
+  # for 90.1-2010 Outpatient, AHU2 set minimum outdoor air flow rate as 0
+  # AHU1 doesn't have economizer
+  def modify_OAcontroller(building_vintage)
+    self.getAirLoopHVACs.each do |air_loop|
+      oa_sys = air_loop.airLoopHVACOutdoorAirSystem.get
+      oa_control = oa_sys.getControllerOutdoorAir
+      if air_loop.name.get == 'PVAV Outpatient F1'        
+        oa_control.setEconomizerControlType('NoEconomizer')
+      elsif air_loop.name.get == 'PVAV Outpatient F2 F3'
+        case building_vintage
+        when '90.1-2010'
+          oa_control.setMinimumOutdoorAirFlowRate(0)
+        else
+          next
+        end
+      end
+    end
+  end
+  
+  # For operating room 1&2 in 2010 and 2013, VAV minimum air flow is set by schedule
+  def reset_or_room_vav_minimum_damper(prototype_input, building_vintage)
+    case building_vintage
+    when '90.1-2004', '90.1-2007'
+      return true
+    when '90.1-2010', '90.1-2013'
+      self.getAirTerminalSingleDuctVAVReheats.sort.each do |airterminal|
+        airterminal_name = airterminal.name.get
+        if airterminal_name.include? "Floor 1 Operating Room 1" or airterminal_name.include? "Floor 1 Operating Room 2"
+          airterminal.setZoneMinimumAirFlowMethod('Scheduled')
+          airterminal.setMinimumAirFlowFractionSchedule(add_schedule("OutPatientHealthCare OR_MinSA_Sched"))
+        end
+      end
+    end
+  end
+  
+  def reset_boiler_sizing_factor
+    self.getBoilerHotWaters.sort.each do |boiler|
+      boiler.setSizingFactor(0.3)
+    end
+  end  
 
 end