openstudio-standards 0.2.12.rc4 → 0.2.13.rc3

data/lib/openstudio-standards/standards/Standards.Space.rb

@@ -569,6 +569,7 @@ class Standard
     construction_name_to_vt_map = {}
     space.surfaces.sort.each do |surface|
       next unless surface.outsideBoundaryCondition == 'Outdoors' && surface.surfaceType == 'Wall'
+
       surface.subSurfaces.sort.each do |sub_surface|
         next unless sub_surface.outsideBoundaryCondition == 'Outdoors' && (sub_surface.subSurfaceType == 'FixedWindow' || sub_surface.subSurfaceType == 'OperableWindow' || sub_surface.subSurfaceType == 'GlassDoor')
 
@@ -684,6 +685,7 @@ class Standard
     construction_name_to_vt_map = {}
     space.surfaces.sort.each do |surface|
       next unless surface.outsideBoundaryCondition == 'Outdoors' && surface.surfaceType == 'RoofCeiling'
+
       surface.subSurfaces.sort.each do |sub_surface|
         next unless sub_surface.outsideBoundaryCondition == 'Outdoors' && sub_surface.subSurfaceType == 'Skylight'
 
@@ -821,8 +823,10 @@ class Standard
     ext_fen_area_m2 = 0
     space.surfaces.each do |surface|
       next unless surface.outsideBoundaryCondition == 'Outdoors'
+
       surface.subSurfaces.each do |sub_surface|
         next unless sub_surface.subSurfaceType == 'FixedWindow' || sub_surface.subSurfaceType == 'OperableWindow' || sub_surface.subSurfaceType == 'Skylight' || sub_surface.subSurfaceType == 'GlassDoor'
+
         ext_fen_area_m2 += sub_surface.netArea
       end
     end
@@ -1332,6 +1336,7 @@ class Standard
       next unless surface.outsideBoundaryCondition == 'Outdoors'
       # Skip non-walls
       next unless surface.surfaceType == 'Wall'
+
       # This surface
       area_m2 += surface.netArea
       # Subsurfaces in this surface
@@ -1356,6 +1361,7 @@ class Standard
       next unless surface.outsideBoundaryCondition == 'Outdoors'
       # Skip non-walls
       next unless surface.surfaceType == 'Wall' || surface.surfaceType == 'RoofCeiling'
+
       # This surface
       area_m2 += surface.netArea
       # Subsurfaces in this surface
@@ -1430,6 +1436,7 @@ class Standard
       largest_surface = nil
       space.surfaces.each do |surface|
         next unless surface.surfaceType == 'Floor' && surface.outsideBoundaryCondition == 'Surface'
+
         if surface.grossArea > largest_floor_area
           largest_floor_area = surface.grossArea
           largest_surface = surface
@@ -1591,6 +1598,7 @@ class Standard
       unless adj_surface.empty?
         space.model.getSpaces.sort.each do |other_space|
           next if other_space == space
+
           other_space.surfaces.each do |surf|
             if surf == adj_surface.get
               spaces << other_space
@@ -1602,8 +1610,10 @@ class Standard
     # If looking for only spaces adjacent on the same floor.
     if same_floor == true
       raise "Cannot get adjacent spaces of space #{space.name} since space not set to BuildingStory" if space.buildingStory.empty?
+
       spaces.each do |other_space|
         raise "One or more adjecent spaces to space #{space.name} is not assigned to a BuildingStory. Ensure all spaces are assigned." if space.buildingStory.empty?
+
         if other_space.buildingStory.get == space.buildingStory.get
           same_floor_spaces << other_space
         end
@@ -1615,6 +1625,7 @@ class Standard
     area_index = []
     array_hash = {}
     return nil if spaces.size.zero?
+
     # iterate through each surface in the space
     space.surfaces.each do |surface|
       # get the adjacent surface in another space.
@@ -1623,6 +1634,7 @@ class Standard
         # go through each of the adjeacent spaces to find the matching  surface/space.
         spaces.each_with_index do |other_space, index|
           next if other_space == space
+
           other_space.surfaces.each do |surf|
             if surf == adj_surface.get
               # initialize array index to zero for first time so += will work.
@@ -1643,7 +1655,7 @@ class Standard
     return get_adjacent_spaces_with_touching_area(same_floor)[0][0]
   end
 
-  # todo - add related related to space_hours_of_operation like set_space_hours_of_operation and shift_and_expand_space_hours_of_operation
+  # TODO: - add related related to space_hours_of_operation like set_space_hours_of_operation and shift_and_expand_space_hours_of_operation
   # todo - ideally these could take in a date range, array of dates and or days of week. Hold off until need is a bit more defined.
 
   # If the model has an hours of operation schedule set in default schedule set for building that looks valid it will
@@ -1654,7 +1666,6 @@ class Standard
   # @param space [Space] takes space
   # @return [Hash] start and end of hours of operation, stat date, end date, bool for each day of the week
   def space_hours_of_operation(space)
-
     default_sch_type = OpenStudio::Model::DefaultScheduleType.new('HoursofOperationSchedule')
     hours_of_operation = space.getDefaultSchedule(default_sch_type)
     if !hours_of_operation.is_initialized
@@ -1682,7 +1693,7 @@ class Standard
     unexpected_val = false
     times = hours_of_operation.defaultDaySchedule.times
     values = hours_of_operation.defaultDaySchedule.values
-    times.each_with_index do |time,i|
+    times.each_with_index do |time, i|
       if values[i] == 0 && hoo_start.nil?
         hoo_start = time.totalHours
       elsif values[i] == 1 && hoo_end.nil?
@@ -1723,8 +1734,8 @@ class Standard
     end
     rule_hash[:hoo_hours] = hoo_hours
     days_used = []
-    indices_vector.each_with_index do |profile_index,i|
-      if profile_index == -1 then days_used << i+1 end
+    indices_vector.each_with_index do |profile_index, i|
+      if profile_index == -1 then days_used << i + 1 end
     end
     rule_hash[:days_used] = days_used
     profiles[-1] = rule_hash
@@ -1738,10 +1749,10 @@ class Standard
       unexpected_val = false
       times = rule.daySchedule.times
       values = rule.daySchedule.values
-      times.each_with_index do |time,i|
+      times.each_with_index do |time, i|
         if values[i] == 0 && hoo_start.nil?
           hoo_start = time.totalHours
-        elsif values[i] == 1  && hoo_end.nil?
+        elsif values[i] == 1 && hoo_end.nil?
           hoo_end = time.totalHours
         elsif values[i] != 1 && values[i] != 0
           unexpected_val = true
@@ -1778,37 +1789,34 @@ class Standard
       end
       rule_hash[:hoo_hours] = hoo_hours
       days_used = []
-      indices_vector.each_with_index do |profile_index,i|
-        if profile_index == rule.ruleIndex then days_used << i+1 end
+      indices_vector.each_with_index do |profile_index, i|
+        if profile_index == rule.ruleIndex then days_used << i + 1 end
       end
       rule_hash[:days_used] = days_used
 
-=begin
-      # todo - delete rule details below unless end up needing to use them
-      if rule.startDate.is_initialized
-        date = rule.startDate.get
-        rule_hash[:start_date] = "#{date.monthOfYear.value}/#{date.dayOfMonth}"
-      else
-        rule_hash[:start_date] = nil
-      end
-      if rule.endDate.is_initialized
-        date = rule.endDate.get
-        rule_hash[:end_date] = "#{date.monthOfYear.value}/#{date.dayOfMonth}"
-      else
-        rule_hash[:end_date] = nil
-      end
-      rule_hash[:mon] = rule.applyMonday
-      rule_hash[:tue] = rule.applyTuesday
-      rule_hash[:wed] = rule.applyWednesday
-      rule_hash[:thu] = rule.applyThursday
-      rule_hash[:fri] = rule.applyFriday
-      rule_hash[:sat] = rule.applySaturday
-      rule_hash[:sun] = rule.applySunday
-=end
+      #       # todo - delete rule details below unless end up needing to use them
+      #       if rule.startDate.is_initialized
+      #         date = rule.startDate.get
+      #         rule_hash[:start_date] = "#{date.monthOfYear.value}/#{date.dayOfMonth}"
+      #       else
+      #         rule_hash[:start_date] = nil
+      #       end
+      #       if rule.endDate.is_initialized
+      #         date = rule.endDate.get
+      #         rule_hash[:end_date] = "#{date.monthOfYear.value}/#{date.dayOfMonth}"
+      #       else
+      #         rule_hash[:end_date] = nil
+      #       end
+      #       rule_hash[:mon] = rule.applyMonday
+      #       rule_hash[:tue] = rule.applyTuesday
+      #       rule_hash[:wed] = rule.applyWednesday
+      #       rule_hash[:thu] = rule.applyThursday
+      #       rule_hash[:fri] = rule.applyFriday
+      #       rule_hash[:sat] = rule.applySaturday
+      #       rule_hash[:sun] = rule.applySunday
 
       # update hash
       profiles[rule.ruleIndex] = rule_hash
-
     end
 
     return profiles
@@ -1835,7 +1843,7 @@ class Standard
       end
     end
 
-    # todo - replace this with logic to get combined hours of operation for collection of spaces.
+    # TODO: - replace this with logic to get combined hours of operation for collection of spaces.
     # each hours_of_operation_array is hash with key for each profile.
     # each profile has hash with keys for hoo_start, hoo_end, hoo_hours, days_used
     # my goal is to compare profiles and days used across all profiles to create new entries as necessary
@@ -1845,11 +1853,10 @@ class Standard
 
     # OpenStudio.logFree(OpenStudio::Debug, 'openstudio.standards.Space', "Evaluating hours of operation for #{space_names.join(',')}: #{hours_of_operation_array}")
 
-    # todo - what is this getting max of, it isn't longest hours of operation, is it the most profiles?
+    # TODO: - what is this getting max of, it isn't longest hours of operation, is it the most profiles?
     hours_of_operation = hours_of_operation_array.max_by { |i| hours_of_operation_array.count(i) }
 
     return hours_of_operation
-
   end
 
   private

data/lib/openstudio-standards/standards/Standards.SpaceType.rb

@@ -1,4 +1,3 @@
-
 class Standard
   # @!group SpaceType
 
@@ -6,7 +5,6 @@ class Standard
   #
   # @return [hash] hash of internal loads for different load types
   def space_type_get_standards_data(space_type)
-
     standards_building_type = if space_type.standardsBuildingType.is_initialized
                                 space_type.standardsBuildingType.get
                               end
@@ -84,6 +82,7 @@ class Standard
     if space_type.name.get.to_s.downcase.include?('plenum')
       return false
     end
+
     if space_type.standardsSpaceType.is_initialized
       if space_type.standardsSpaceType.get.downcase.include?('plenum')
         return false
@@ -94,7 +93,7 @@ class Standard
     space_type_properties = space_type_get_standards_data(space_type)
 
     # Need to add a check, or it'll crash on space_type_properties['occupancy_per_area'].to_f below
-    if space_type_properties.nil?
+    if space_type_properties.nil? || space_type_properties.empty?
       OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} was not found in the standards data.")
       return false
     end
@@ -119,6 +118,7 @@ class Standard
       elsif instances.size > 1
         instances.each_with_index do |inst, i|
           next if i.zero?
+
           OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "Removed #{inst.name} from #{space_type.name}.")
           inst.remove
         end
@@ -211,6 +211,7 @@ class Standard
       elsif instances.size > 1
         instances.each_with_index do |inst, i|
           next if i.zero?
+
           OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "Removed #{inst.name} from #{space_type.name}.")
           inst.remove
         end
@@ -289,6 +290,7 @@ class Standard
       elsif instances.size > 1
         instances.each_with_index do |inst, i|
           next if i.zero?
+
           OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "Removed #{inst.name} from #{space_type.name}.")
           inst.remove
         end
@@ -331,6 +333,7 @@ class Standard
       elsif instances.size > 1
         instances.each_with_index do |inst, i|
           next if i.zero?
+
           OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "Removed #{inst.name} from #{space_type.name}.")
           inst.remove
         end
@@ -421,6 +424,7 @@ class Standard
       elsif instances.size > 1
         instances.each_with_index do |inst, i|
           next if i.zero?
+
           OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "Removed #{inst.name} from #{space_type.name}.")
           inst.remove
         end
@@ -495,48 +499,39 @@ class Standard
     # Lights
     if set_lights
 
-      apply_lighting_schedule(space_type, space_type_properties,default_sch_set)
+      apply_lighting_schedule(space_type, space_type_properties, default_sch_set)
 
     end
 
-
-
     # Electric Equipment
     if set_electric_equipment
-
       elec_equip_sch = space_type_properties['electric_equipment_schedule']
       unless elec_equip_sch.nil?
         default_sch_set.setElectricEquipmentSchedule(model_add_schedule(space_type.model, elec_equip_sch))
         OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} set electric equipment schedule to #{elec_equip_sch}.")
       end
-
     end
 
     # Gas Equipment
     if set_gas_equipment
-
       gas_equip_sch = space_type_properties['gas_equipment_schedule']
       unless gas_equip_sch.nil?
         default_sch_set.setGasEquipmentSchedule(model_add_schedule(space_type.model, gas_equip_sch))
         OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} set gas equipment schedule to #{gas_equip_sch}.")
       end
-
     end
 
     # Infiltration
     if set_infiltration
-
       infiltration_sch = space_type_properties['infiltration_schedule']
       unless infiltration_sch.nil?
         default_sch_set.setInfiltrationSchedule(model_add_schedule(space_type.model, infiltration_sch))
         OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} set infiltration schedule to #{infiltration_sch}.")
       end
-
     end
 
     # Thermostat
     if make_thermostat
-
       thermostat = OpenStudio::Model::ThermostatSetpointDualSetpoint.new(space_type.model)
       thermostat.setName("#{space_type.name} Thermostat")
 
@@ -551,20 +546,17 @@ class Standard
         thermostat.setCoolingSetpointTemperatureSchedule(model_add_schedule(space_type.model, cooling_setpoint_sch))
         OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} set cooling setpoint schedule to #{cooling_setpoint_sch}.")
       end
-
     end
 
     return true
   end
 
-  def apply_lighting_schedule(space_type, space_type_properties,default_sch_set)
-
+  def apply_lighting_schedule(space_type, space_type_properties, default_sch_set)
     lighting_sch = space_type_properties['lighting_schedule']
     unless lighting_sch.nil?
       default_sch_set.setLightingSchedule(model_add_schedule(space_type.model, lighting_sch))
       OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.SpaceType', "#{space_type.name} set lighting schedule to #{lighting_sch}.")
     end
-
   end
 
   # Returns standards data for selected construction

data/lib/openstudio-standards/standards/Standards.SubSurface.rb

@@ -11,24 +11,24 @@ class Standard
 
     # Define the envelope component infiltration rates
     component_infil_rates_cfm_per_ft2 = {
-        'baseline' => {
-            'opaque_door' => 0.40,
-            'loading_dock_door' => 0.40,
-            'swinging_or_revolving_glass_door' => 1.0,
-            'vestibule' => 1.0,
-            'sliding_glass_door' => 0.40,
-            'window' => 0.40,
-            'skylight' => 0.40
-        },
-        'advanced' => {
-            'opaque_door' => 0.20,
-            'loading_dock_door' => 0.20,
-            'swinging_or_revolving_glass_door' => 1.0,
-            'vestibule' => 1.0,
-            'sliding_glass_door' => 0.20,
-            'window' => 0.20,
-            'skylight' => 0.20
-        }
+      'baseline' => {
+        'opaque_door' => 0.40,
+        'loading_dock_door' => 0.40,
+        'swinging_or_revolving_glass_door' => 1.0,
+        'vestibule' => 1.0,
+        'sliding_glass_door' => 0.40,
+        'window' => 0.40,
+        'skylight' => 0.40
+      },
+      'advanced' => {
+        'opaque_door' => 0.20,
+        'loading_dock_door' => 0.20,
+        'swinging_or_revolving_glass_door' => 1.0,
+        'vestibule' => 1.0,
+        'sliding_glass_door' => 0.20,
+        'window' => 0.20,
+        'skylight' => 0.20
+      }
     }
 
     boundary_condition = sub_surface.outsideBoundaryCondition
@@ -231,7 +231,7 @@ class Standard
   # 'sub_surface_create_centered_subsurface_from_scaled_surface' method because it can handle concave surfaces.
   # However, it takes longer because it uses BTAP::Geometry::Surfaces.make_convex_surfaces which includes many nested
   # loops that cycle through the verticies in a surface.
-  def sub_surface_create_scaled_subsurfaces_from_surface (surface:, area_fraction:, model:, consturction:)
+  def sub_surface_create_scaled_subsurfaces_from_surface(surface:, area_fraction:, model:, consturction:)
     # Set geometry tolerences:
     geometry_tolerence = 12
     # Get rid of all existing subsurfaces.
@@ -243,7 +243,7 @@ class Standard
     surface.vertices.each do |surf_vert|
       surface.vertices.each do |surf_vert_2|
         unless surf_vert_2.z.to_f.round(geometry_tolerence) == surf_vert.z.to_f.round(geometry_tolerence)
-          OpenStudio.logFree(OpenStudio::Error, 'openstudio.model.Model', "Currently skylights can only be added to buildings with non-plenum flat roofs.")
+          OpenStudio.logFree(OpenStudio::Error, 'openstudio.model.Model', 'Currently skylights can only be added to buildings with non-plenum flat roofs.')
         end
       end
     end
@@ -257,7 +257,7 @@ class Standard
 
     # Turn everything back into OpenStudio stuff
     os_surf_points = []
-    os_surf_cents =[]
+    os_surf_cents = []
     for i in 0..(new_surfaces.length - 1)
       os_surf_point = []
       for j in 0..(new_surfaces[i].length - 1)
@@ -295,11 +295,11 @@ class Standard
       # 'skylight').  If there will be more than one subsurface then add a counter at the end.
       new_name = surface.name.to_s + '_' + new_sub_surface.subSurfaceType.to_s
       if new_surfaces.length > 1
-        new_name = surface.name.to_s + '_' + new_sub_surface.subSurfaceType.to_s + '_' + "#{index}"
+        new_name = surface.name.to_s + '_' + new_sub_surface.subSurfaceType.to_s + '_' + index.to_s
       end
-      #Set the skylight type to 'Skylight'
+      # Set the skylight type to 'Skylight'
       new_sub_surface.setSubSurfaceType('Skylight')
-      #Set the skylight construction to whatever was passed (should be the default skylight construction)
+      # Set the skylight construction to whatever was passed (should be the default skylight construction)
       new_sub_surface.setConstruction(consturction)
       new_sub_surface.setName(new_name)
       # There is now only one surface on the subsurface.  Enforce this
@@ -322,8 +322,6 @@ class Standard
   # This removes all of the subsurfaces from a surface.  Is a preparation for replaceing windows or clearing doors
   # before adding windows.
   def remove_All_Subsurfaces(surface:)
-    surface.subSurfaces.sort.each do |sub_surface|
-      sub_surface.remove
-    end
+    surface.subSurfaces.sort.each(&:remove)
   end
 end

data/lib/openstudio-standards/standards/Standards.Surface.rb

@@ -86,14 +86,14 @@ class Standard
       vertical_surface.subSurfaces.sort.each do |vertical_subsurface|
         # Need to fix this so that error show up in right place
         if vertical_subsurface.nil?
-          puts "Surface does not exist"
+          puts 'Surface does not exist'
         else
           vertical_subsurface.remove
         end
       end
       # corner_coords = vertical_surface.vertices
-      code_window_area = fdwr*vertical_surface.grossArea
-      code_window_width = code_window_area/window_height_m
+      code_window_area = fdwr * vertical_surface.grossArea
+      code_window_width = code_window_area / window_height_m
       min_z = 0
       vertical_surface.vertices.each_with_index do |vertex, index|
         if index == 0
@@ -134,57 +134,58 @@ class Standard
           # Note that this will not include foundation walls.
           next unless surface.surfaceType == 'Wall'
           next unless surface.outsideBoundaryCondition == 'Outdoors'
+
           # Determine if the wall is vertical which I define as being between 89 and 91 degrees from horizontal.
           tilt_radian = surface.tilt
           tilt_degrees = OpenStudio.convert(tilt_radian, 'rad', 'deg').get
           sub_surface_info = []
-          if tilt_degrees <= max_angle and tilt_degrees >= min_angle
+          if tilt_degrees <= max_angle && tilt_degrees >= min_angle
             # If the wall is vertical determine if it is adjacent to a plenum.  If yes include it in the array of
             # plenum walls and add it to the plenum wall area counter (accounting for space multipliers).
             if space_plenum?(space)
               plenum_surfaces << surface
-              exp_plenum_area += surface.grossArea*space.multiplier
+              exp_plenum_area += surface.grossArea * space.multiplier
             else
               # If not a plenum then include it in the array of non-plenum walls and add it to the non-plenum area
               # counter (accounting for space multipliers).
               exposed_surfaces << surface
-              exp_nonplenum_area += surface.grossArea*space.multiplier
+              exp_nonplenum_area += surface.grossArea * space.multiplier
               surface.subSurfaces.sort.each do |sub_surface|
-                sub_surface_area += sub_surface.grossArea.to_f*space.multiplier
+                sub_surface_area += sub_surface.grossArea.to_f * space.multiplier
                 sub_surface_info << {
-                    "subsurface_name" => sub_surface.nameString,
-                    "subsurface_type" => sub_surface.subSurfaceType,
-                    "gross_area_m2" => sub_surface.grossArea.to_f,
-                    "construction_name" => sub_surface.construction.get.nameString
+                  'subsurface_name' => sub_surface.nameString,
+                  'subsurface_type' => sub_surface.subSurfaceType,
+                  'gross_area_m2' => sub_surface.grossArea.to_f,
+                  'construction_name' => sub_surface.construction.get.nameString
                 }
               end
               unless sub_surface_info.empty?
                 sub_surfaces_info << {
-                    "surface_name" => surface.nameString,
-                    "subsurfaces" => sub_surface_info
+                  'surface_name' => surface.nameString,
+                  'subsurfaces' => sub_surface_info
                 }
               end
             end
             # Regardless of if the wall is adjacent to a plenum or not add it to the exposed wall area adjacent to
             # conditioned spaces (accounting for space multipliers).
-            total_exp_area += surface.grossArea*space.multiplier
+            total_exp_area += surface.grossArea * space.multiplier
           end
         end
       end
     end
     fdwr = 999
     unless exp_nonplenum_area < 0.1
-      fdwr = sub_surface_area/exp_nonplenum_area
+      fdwr = sub_surface_area / exp_nonplenum_area
     end
     # Add everything into a hash and return that hash to whomever called the method.
     exp_surf_info = {
-        "total_exp_wall_area_m2" => total_exp_area,
-        "exp_plenum_wall_area_m2" => exp_plenum_area,
-        "exp_nonplenum_wall_area_m2" => exp_nonplenum_area,
-        "exp_plenum_walls" => plenum_surfaces,
-        "exp_nonplenum_walls" => exposed_surfaces,
-        "fdwr" => fdwr,
-        "sub_surfaces" => sub_surfaces_info
+      'total_exp_wall_area_m2' => total_exp_area,
+      'exp_plenum_wall_area_m2' => exp_plenum_area,
+      'exp_nonplenum_wall_area_m2' => exp_nonplenum_area,
+      'exp_plenum_walls' => plenum_surfaces,
+      'exp_nonplenum_walls' => exposed_surfaces,
+      'fdwr' => fdwr,
+      'sub_surfaces' => sub_surfaces_info
     }
     return exp_surf_info
   end
@@ -209,53 +210,54 @@ class Standard
           # Assume a roof is of type 'RoofCeiling' and has an 'Outdoors' boundary condition.
           next unless surface.surfaceType == 'RoofCeiling'
           next unless surface.outsideBoundaryCondition == 'Outdoors'
+
           # Determine if the roof is adjacent to a plenum.
           sub_surface_info = []
           if space_plenum?(space)
             # If the roof is adjacent to a plenum add it to the plenum roof array and the plenum roof area counter
             # (accounting for space multipliers).
             plenum_surfaces << surface
-            exp_plenum_area += surface.grossArea*space.multiplier
+            exp_plenum_area += surface.grossArea * space.multiplier
           else
             # If the roof is not adjacent to a plenum add it to the non-plenum roof array and the non-plenum roof area
             # counter (accounting for space multipliers).
             exposed_surfaces << surface
-            exp_nonplenum_area += surface.grossArea*space.multiplier
+            exp_nonplenum_area += surface.grossArea * space.multiplier
             surface.subSurfaces.sort.each do |sub_surface|
-              sub_surface_area += sub_surface.grossArea.to_f*space.multiplier
+              sub_surface_area += sub_surface.grossArea.to_f * space.multiplier
               sub_surface_info << {
-                  "subsurface_name" => sub_surface.nameString,
-                  "subsurface_type" => sub_surface.subSurfaceType,
-                  "gross_area_m2" => sub_surface.grossArea.to_f,
-                  "construction_name" => sub_surface.construction.get.nameString
+                'subsurface_name' => sub_surface.nameString,
+                'subsurface_type' => sub_surface.subSurfaceType,
+                'gross_area_m2' => sub_surface.grossArea.to_f,
+                'construction_name' => sub_surface.construction.get.nameString
               }
             end
             unless sub_surface_info.empty?
               sub_surfaces_info << {
-                  "surface_name" => surface.nameString,
-                  "subsurfaces" => sub_surface_info
+                'surface_name' => surface.nameString,
+                'subsurfaces' => sub_surface_info
               }
             end
           end
           # Regardless of if the roof is adjacent to a plenum or not add it to the total roof area counter (accounting
           # for space multipliers).
-          total_exp_area += surface.grossArea*space.multiplier
+          total_exp_area += surface.grossArea * space.multiplier
         end
       end
     end
     srr = 999
     unless exp_nonplenum_area < 0.1
-      srr = sub_surface_area/exp_nonplenum_area
+      srr = sub_surface_area / exp_nonplenum_area
     end
     # Put the information into a hash and return it to whomever called this method.
     exp_surf_info = {
-        "total_exp_roof_area_m2" => total_exp_area,
-        "exp_plenum_roof_area_m2" => exp_plenum_area,
-        "exp_nonplenum_roof_area_m2" => exp_nonplenum_area,
-        "exp_plenum_roofs" => plenum_surfaces,
-        "exp_nonplenum_roofs" => exposed_surfaces,
-        "srr" => srr,
-        "sub_surfaces" => sub_surfaces_info
+      'total_exp_roof_area_m2' => total_exp_area,
+      'exp_plenum_roof_area_m2' => exp_plenum_area,
+      'exp_nonplenum_roof_area_m2' => exp_nonplenum_area,
+      'exp_plenum_roofs' => plenum_surfaces,
+      'exp_nonplenum_roofs' => exposed_surfaces,
+      'srr' => srr,
+      'sub_surfaces' => sub_surfaces_info
     }
     return exp_surf_info
   end
@@ -300,6 +302,7 @@ class Standard
       end
       # If no outdoor roofceiling go to the next space.
       next if outdoor_roof == false
+
       z_Origin = space.zOrigin.to_f
       ceiling_centroid = [0, 0, 0]
 
@@ -307,8 +310,8 @@ class Standard
       # centroid calculation.
       space_surfaces.each do |sp_surface|
         if max_surf.centroid.z.to_f.round(tol) == sp_surface.centroid.z.to_f.round(tol)
-          ceiling_centroid[0] += sp_surface.centroid.x.to_f*sp_surface.grossArea.to_f
-          ceiling_centroid[1] += sp_surface.centroid.y.to_f*sp_surface.grossArea.to_f
+          ceiling_centroid[0] += sp_surface.centroid.x.to_f * sp_surface.grossArea.to_f
+          ceiling_centroid[1] += sp_surface.centroid.y.to_f * sp_surface.grossArea.to_f
           ceiling_centroid[2] += sp_surface.grossArea.to_f
         end
       end
@@ -319,11 +322,11 @@ class Standard
 
       # Put the info into an array containing hashes of spaces with outdoor roofceilings
       spaces_info << {
-          space: space,
-          x: ceiling_centroid[0] + space.xOrigin.to_f,
-          y: ceiling_centroid[1] + space.yOrigin.to_f,
-          z: max_surf.centroid.z.to_f + z_Origin,
-          area_m2: ceiling_centroid[2]
+        space: space,
+        x: ceiling_centroid[0] + space.xOrigin.to_f,
+        y: ceiling_centroid[1] + space.yOrigin.to_f,
+        z: max_surf.centroid.z.to_f + z_Origin,
+        area_m2: ceiling_centroid[2]
       }
       # This is to determine which are the global highest outdoor roofceilings
       if max_height.round(tol) < (max_surf.centroid.z.to_f + z_Origin).round(tol)
@@ -336,8 +339,8 @@ class Standard
       # calculate the centroid
       if space_info[:z].to_f.round(tol) == max_height.round(tol)
         top_spaces << space_info
-        roof_centroid[0] += space_info[:x]*space_info[:area_m2]
-        roof_centroid[1] += space_info[:y]*space_info[:area_m2]
+        roof_centroid[0] += space_info[:x] * space_info[:area_m2]
+        roof_centroid[1] += space_info[:y] * space_info[:area_m2]
         roof_centroid[2] += space_info[:area_m2]
       end
     end
@@ -346,10 +349,52 @@ class Standard
     roof_centroid[0] /= roof_centroid[2]
     roof_centroid[1] /= roof_centroid[2]
     roof_cent = {
-        top_spaces: top_spaces,
-        roof_centroid: [roof_centroid[0], roof_centroid[1], max_height],
-        roof_area: roof_centroid[2]
+      top_spaces: top_spaces,
+      roof_centroid: [roof_centroid[0], roof_centroid[1], max_height],
+      roof_area: roof_centroid[2]
     }
     return roof_cent
   end
+
+  # Calculate a surface's absolute azimuth
+  # source: https://github.com/NREL/openstudio-extension-gem/blob/e354355054b83ffc26e3b69befa20d6baf5ef242/lib/openstudio/extension/core/os_lib_geometry.rb#L913
+  #
+  # @param surface [OpenStudio::Model:Surface] OpenStudio Surface object
+  # @return [Float] Surface absolute azimuth
+  def surface_absolute_azimuth(surface)
+    # Get associated space
+    space = surface.space.get
+
+    # Get model object
+    model = surface.model
+
+    # Calculate azimuth
+    surface_azimuth_rel_space = OpenStudio.convert(surface.azimuth, 'rad', 'deg').get
+    space_dir_rel_N = space.directionofRelativeNorth
+    building_dir_rel_N = model.getBuilding.northAxis
+    surface_abs_azimuth = surface_azimuth_rel_space + space_dir_rel_N + building_dir_rel_N
+    surface_abs_azimuth -= 360.0 until surface_abs_azimuth < 360.0
+
+    return surface_abs_azimuth
+  end
+
+  # Determine a surface absolute cardinal direction
+  #
+  # @param surface [OpenStudio::Model::Surface] OpenStudio Surface object
+  # @return [String] Surface absolute cardinal
+  def surface_cardinal_direction(surface)
+    # Get the surface's absolute azimuth
+    surface_abs_azimuth = surface_absolute_azimuth(surface)
+
+    # Determine the surface's cardinal direction
+    if (surface_abs_azimuth >= 0 && surface_abs_azimuth <= 45) || (surface_abs_azimuth > 315 && surface_abs_azimuth <= 360)
+      return 'N'
+    elsif surface_abs_azimuth > 45 && surface_abs_azimuth <= 135
+      return 'E'
+    elsif surface_abs_azimuth > 135 && surface_abs_azimuth <= 225
+      return 'S'
+    elsif surface_abs_azimuth > 225 && surface_abs_azimuth <= 315
+      return 'W'
+    end
+  end
 end