urbanopt-cli 0.9.2 → 0.10.0

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/resources/output.rb

@@ -20,8 +20,9 @@ end
 class EUT
   # End Use Types
   Heating = 'Heating'
-  HeatingHeatPumpBackup = 'Heating Heat Pump Backup'
   HeatingFanPump = 'Heating Fans/Pumps'
+  HeatingHeatPumpBackup = 'Heating Heat Pump Backup'
+  HeatingHeatPumpBackupFanPump = 'Heating Heat Pump Backup Fans/Pumps'
   Cooling = 'Cooling'
   CoolingFanPump = 'Cooling Fans/Pumps'
   HotWater = 'Hot Water'
@@ -69,6 +70,7 @@ end
 class LT
   # Load Types
   Heating = 'Heating: Delivered'
+  HeatingHeatPumpBackup = 'Heating: Heat Pump Backup' # Needed for ERI calculation for dual-fuel heat pumps
   Cooling = 'Cooling: Delivered'
   HotWaterDelivered = 'Hot Water: Delivered'
   HotWaterTankLosses = 'Hot Water: Tank Losses'
@@ -84,8 +86,10 @@ class CLT
   RimJoists = 'Rim Joists'
   FoundationWalls = 'Foundation Walls'
   Doors = 'Doors'
-  Windows = 'Windows'
-  Skylights = 'Skylights'
+  WindowsConduction = 'Windows Conduction'
+  WindowsSolar = 'Windows Solar'
+  SkylightsConduction = 'Skylights Conduction'
+  SkylightsSolar = 'Skylights Solar'
   Floors = 'Floors'
   Slabs = 'Slabs'
   InternalMass = 'Internal Mass'
@@ -95,6 +99,7 @@ class CLT
   WholeHouseFan = 'Whole House Fan'
   Ducts = 'Ducts'
   InternalGains = 'Internal Gains'
+  Lighting = 'Lighting'
 end
 
 class UHT
@@ -103,12 +108,6 @@ class UHT
   Cooling = 'Cooling'
 end
 
-class ILT
-  # Ideal Load Types
-  Heating = 'Heating'
-  Cooling = 'Cooling'
-end
-
 class PLT
   # Peak Load Types
   Heating = 'Heating: Delivered'
@@ -138,95 +137,3 @@ class WT
   DiffuseSolar = 'Diffuse Solar Radiation'
   DirectSolar = 'Direct Solar Radiation'
 end
-
-class OutputMethods
-  def self.get_timestamps(msgpackData, hpxml, use_timestamp_start_convention, add_dst_column = false, add_utc_column = false,
-                          use_dview_format = false, timeseries_frequency = nil)
-    return if msgpackData.nil?
-
-    if msgpackData.keys.include? 'MeterData'
-      # Get data for ReportUtilityBills measure
-      ep_timestamps = msgpackData['MeterData']['Monthly']['Rows'].map { |r| r.keys[0] }
-    else
-      # Get data for other reporting measures
-      ep_timestamps = msgpackData['Rows'].map { |r| r.keys[0] }
-    end
-
-    if add_dst_column || use_dview_format
-      dst_start_ts = Time.utc(hpxml.header.sim_calendar_year, hpxml.header.dst_begin_month, hpxml.header.dst_begin_day, 2)
-      dst_end_ts = Time.utc(hpxml.header.sim_calendar_year, hpxml.header.dst_end_month, hpxml.header.dst_end_day, 1)
-    end
-    if add_utc_column
-      utc_offset = hpxml.header.time_zone_utc_offset
-      utc_offset *= 3600 # seconds
-    end
-
-    timestamps = []
-    timestamps_dst = [] if add_dst_column || use_dview_format
-    timestamps_utc = [] if add_utc_column
-    year = hpxml.header.sim_calendar_year
-    ep_timestamps.each do |ep_timestamp|
-      month_day, hour_minute = ep_timestamp.split(' ')
-      month, day = month_day.split('/').map(&:to_i)
-      hour, minute, _ = hour_minute.split(':').map(&:to_i)
-
-      # Convert from EnergyPlus default (end-of-timestep) to start-of-timestep convention
-      if use_timestamp_start_convention
-        if timeseries_frequency == 'timestep'
-          ts_offset = hpxml.header.timestep * 60 # seconds
-        elsif timeseries_frequency == 'hourly'
-          ts_offset = 60 * 60 # seconds
-        elsif timeseries_frequency == 'daily'
-          ts_offset = 60 * 60 * 24 # seconds
-        elsif timeseries_frequency == 'monthly'
-          ts_offset = Constants.NumDaysInMonths(year)[month - 1] * 60 * 60 * 24 # seconds
-        else
-          fail 'Unexpected timeseries_frequency/'
-        end
-      end
-
-      ts = Time.utc(year, month, day, hour, minute)
-      ts -= ts_offset unless ts_offset.nil?
-
-      timestamps << ts.iso8601.delete('Z')
-
-      if add_dst_column || use_dview_format
-        if (ts >= dst_start_ts) && (ts < dst_end_ts)
-          ts_dst = ts + 3600 # 1 hr shift forward
-        else
-          ts_dst = ts
-        end
-        timestamps_dst << ts_dst.iso8601.delete('Z')
-      end
-
-      if add_utc_column
-        ts_utc = ts - utc_offset
-        timestamps_utc << ts_utc.iso8601
-      end
-    end
-
-    return timestamps, timestamps_dst, timestamps_utc
-  end
-
-  def self.msgpack_frequency_map
-    return {
-      'timestep' => 'TimeStep',
-      'hourly' => 'Hourly',
-      'daily' => 'Daily',
-      'monthly' => 'Monthly',
-    }
-  end
-
-  def self.get_dst_start_end_indexes(timestamps, timestamps_dst)
-    dst_start_ix = nil
-    dst_end_ix = nil
-    timestamps.zip(timestamps_dst).each_with_index do |ts, i|
-      dst_start_ix = i if ts[0] != ts[1] && dst_start_ix.nil?
-      dst_end_ix = i if ts[0] == ts[1] && dst_end_ix.nil? && !dst_start_ix.nil?
-    end
-
-    dst_end_ix = timestamps.size - 1 if dst_end_ix.nil? # run period ends before DST ends
-
-    return dst_start_ix, dst_end_ix
-  end
-end

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/resources/psychrometrics.rb

@@ -78,7 +78,7 @@ class Psychrometrics
     return psat
   end
 
-  def self.Tsat_fP(p)
+  def self.Tsat_fP(runner, p)
     '''
     Description:
     ------------
@@ -120,14 +120,14 @@ class Psychrometrics
 
     end
 
-    if not cvg
-      puts 'Warning: Tsat_fP failed to converge'
+    if !cvg && !runner.nil?
+      runner.registerWarning('Tsat_fP failed to converge')
     end
 
     return tsat
   end
 
-  def self.Tsat_fh_P(h, p)
+  def self.Tsat_fh_P(runner, h, p)
     '''
     Description:
     ------------
@@ -171,8 +171,8 @@ class Psychrometrics
 
     end
 
-    if not cvg
-      puts 'Warning: Tsat_fh_P failed to converge'
+    if !cvg && !runner.nil?
+      runner.registerWarning('Tsat_fh_P failed to converge')
     end
 
     return tdb
@@ -276,7 +276,7 @@ class Psychrometrics
     return pstd
   end
 
-  def self.Twb_fT_w_P(tdb, w, p)
+  def self.Twb_fT_w_P(runner, tdb, w, p)
     '''
     Description:
     ------------
@@ -300,7 +300,7 @@ class Psychrometrics
         Twb     float      wetbulb temperature (degF)
     '''
     # Initialize
-    tboil = self.Tsat_fP(p) # (degF)
+    tboil = self.Tsat_fP(runner, p) # (degF)
     twb = [[tdb, tboil - 0.1].min, 0.0].max # (degF)
 
     twb1 = twb # (degF)
@@ -332,8 +332,8 @@ class Psychrometrics
       end
     end
 
-    if not cvg
-      puts 'Warning: Twb_fT_w_P failed to converge'
+    if !cvg && !runner.nil?
+      runner.registerWarning('Twb_fT_w_P failed to converge')
     end
 
     if twb > tdb
@@ -366,7 +366,7 @@ class Psychrometrics
     return w
   end
 
-  def self.CoilAoFactor(dBin, wBin, p, qdot, cfm, shr)
+  def self.CoilAoFactor(runner, dBin, p, qdot, cfm, shr, win)
     '''
     Description:
     ------------
@@ -381,25 +381,25 @@ class Psychrometrics
     Inputs:
     -------
         Tdb    float    Entering Dry Bulb (degF)
-        Twb    float    Entering Wet Bulb (degF)
         P      float    Barometric pressure (psi)
         Qdot   float    Total capacity of unit (kBtu/h)
         cfm    float    Volumetric flow rate of unit (CFM)
         shr    float    Sensible heat ratio
+        win    float    Entering humidity ratio
     Outputs:
     --------
         Ao    float    Coil Ao Factor
     '''
 
-    bf = self.CoilBypassFactor(dBin, wBin, p, qdot, cfm, shr)
-    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, wBin, p, cfm), 'lbm/min', 'kg/s')
+    bf = self.CoilBypassFactor(runner, dBin, p, qdot, cfm, shr, win)
+    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, p, cfm, win), 'lbm/min', 'kg/s')
 
     ntu = -1.0 * Math.log(bf)
     ao = ntu * mfr
     return ao
   end
 
-  def self.CoilBypassFactor(dBin, wBin, p, qdot, cfm, shr)
+  def self.CoilBypassFactor(runner, dBin, p, qdot, cfm, shr, win)
     '''
     Description:
     ------------
@@ -414,20 +414,19 @@ class Psychrometrics
     Inputs:
     -------
         Tdb    float    Entering Dry Bulb (degF)
-        Twb    float    Entering Wet Bulb (degF)
         P      float    Barometric pressure (psi)
         Qdot   float    Total capacity of unit (kBtu/h)
         cfm    float    Volumetric flow rate of unit (CFM)
         shr    float    Sensible heat ratio
+        win    float    Entering humidity ratio
     Outputs:
     --------
         CBF    float    Coil Bypass Factor
     '''
 
-    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, wBin, p, cfm), 'lbm/min', 'kg/s')
+    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, p, cfm, win), 'lbm/min', 'kg/s')
 
     tin = UnitConversions.convert(dBin, 'F', 'C')
-    win = w_fT_Twb_P(dBin, wBin, p)
     p = UnitConversions.convert(p, 'psi', 'kPa')
 
     dH = UnitConversions.convert(qdot, 'kBtu/hr', 'W') / mfr
@@ -437,16 +436,12 @@ class Psychrometrics
     dW = win - wout
     hout = hin - dH
     tout = self.T_fw_h_SI(wout, hout)
-    rH_out = self.R_fT_w_P_SI(tout, wout, p)
-
-    if rH_out > 1
-      puts 'Error: Conditions passed to CoilBypassFactor result in outlet RH > 100%'
-    end
+    # rH_out = self.R_fT_w_P_SI(tout, wout, p)
 
     dT = tin - tout
     m_c = dW / dT
 
-    t_ADP = self.Tdp_fP_w_SI(p, wout) # Initial guess for iteration
+    t_ADP = self.Tdp_fP_w_SI(runner, p, wout) # Initial guess for iteration
 
     if shr == 1
       w_ADP = w_fT_Twb_P_SI(t_ADP, t_ADP, p)
@@ -461,7 +456,6 @@ class Psychrometrics
     d_T_ADP = 5.0
 
     while (cnt < 100) && (tol > 0.001)
-      # for i in range(1,itmax+1):
 
       if cnt > 0
         t_ADP += d_T_ADP
@@ -491,7 +485,7 @@ class Psychrometrics
     return [bF, 0.01].max
   end
 
-  def self.CalculateMassflowRate(dBin, wBin, p, cfm)
+  def self.CalculateMassflowRate(dBin, p, cfm, win)
     '''
     Description:
     ------------
@@ -504,14 +498,13 @@ class Psychrometrics
     Inputs:
     -------
         Tdb    float    Entering Dry Bulb (degF)
-        Twb    float    Entering Wet Bulb (degF)
         P      float    Barometric pressure (psi)
         cfm    float    Volumetric flow rate of unit (CFM)
+        win    float    Entering humidity ratio
     Outputs:
     --------
         mfr    float    mass flow rate (lbm/min)
     '''
-    win = w_fT_Twb_P(dBin, wBin, p)
     rho_in = rhoD_fT_w_P(dBin, win, p)
     mfr = cfm * rho_in
     return mfr
@@ -566,7 +559,7 @@ class Psychrometrics
     return self.R_fT_w_P(UnitConversions.convert(tdb, 'C', 'F'), w, UnitConversions.convert(p, 'kPa', 'psi'))
   end
 
-  def self.Tdp_fP_w_SI(p, w)
+  def self.Tdp_fP_w_SI(runner, p, w)
     '''
     Description:
     ------------
@@ -592,7 +585,7 @@ class Psychrometrics
     --------
         Tdp     float      dewpoint temperature  (degC)
     '''
-    return UnitConversions.convert(self.Tdp_fP_w(UnitConversions.convert(p, 'kPa', 'psi'), w), 'F', 'C')
+    return UnitConversions.convert(self.Tdp_fP_w(runner, UnitConversions.convert(p, 'kPa', 'psi'), w), 'F', 'C')
   end
 
   def self.w_fT_Twb_P_SI(tdb, twb, p)
@@ -699,7 +692,7 @@ class Psychrometrics
     return pw
   end
 
-  def self.Tdp_fP_w(p, w)
+  def self.Tdp_fP_w(runner, p, w)
     '''
     Description:
     ------------
@@ -753,14 +746,14 @@ class Psychrometrics
         tdp = -999.0
       else
         pw = self.Pw_fP_w(p, w)
-        tdp = self.Tsat_fP(pw)
+        tdp = self.Tsat_fP(runner, pw)
       end
 
     end
     return tdp
   end
 
-  def self.CalculateSHR(dBin, wBin, p, q, cfm, ao)
+  def self.CalculateSHR(runner, dBin, p, q, cfm, ao, win)
     '''
     Description:
     ------------
@@ -774,7 +767,6 @@ class Psychrometrics
     Inputs:
     -------
         Tdb    float    Entering Dry Bulb (degF)
-        Twb    float    Entering Wet Bulb (degF)
         P      float    Barometric pressure (psi)
         Q      float    Total capacity of unit (kBtu/h)
         cfm    float    Volumetric flow rate of unit (CFM)
@@ -784,10 +776,9 @@ class Psychrometrics
         SHR    float    Sensible Heat Ratio
     '''
 
-    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, wBin, p, cfm), 'lbm/min', 'kg/s')
+    mfr = UnitConversions.convert(self.CalculateMassflowRate(dBin, p, cfm, win), 'lbm/min', 'kg/s')
     bf = Math.exp(-1.0 * ao / mfr)
 
-    win = w_fT_Twb_P(dBin, wBin, p)
     p = UnitConversions.convert(p, 'psi', 'kPa')
     tin = UnitConversions.convert(dBin, 'F', 'C')
     hin = h_fT_w_SI(tin, win)
@@ -798,7 +789,7 @@ class Psychrometrics
     # W_ADP = self.w_fT_h_SI(T_ADP, H_ADP)
 
     # Initialize
-    t_ADP = self.Tdp_fP_w_SI(p, win)
+    t_ADP = self.Tdp_fP_w_SI(runner, p, win)
     t_ADP_1 = t_ADP # (degC)
     t_ADP_2 = t_ADP # (degC)
     w_ADP = w_fT_R_P_SI(t_ADP, 1.0, p)
@@ -821,7 +812,7 @@ class Psychrometrics
     end
 
     if not cvg
-      puts 'Warning: Tsat_fh_P failed to converge'
+      runnger.registerWarning('CalculateSHR failed to converge')
     end
 
     h_Tin_Wadp = h_fT_w_SI(tin, w_ADP)
@@ -889,7 +880,7 @@ class Psychrometrics
     return w
   end
 
-  def self.Twb_fT_R_P(tdb, r, p)
+  def self.Twb_fT_R_P(runner, tdb, r, p)
     '''
     Description:
     ------------
@@ -912,7 +903,7 @@ class Psychrometrics
     '''
 
     w = w_fT_R_P(tdb, r, p)
-    twb = self.Twb_fT_w_P(tdb, w, p)
+    twb = self.Twb_fT_w_P(runner, tdb, w, p)
     return twb
   end
 end

data/example_files/resources/hpxml-measures/HPXMLtoOpenStudio/resources/pv.rb

@@ -20,7 +20,7 @@ class PV
 
       ipvwatts = OpenStudio::Model::ElectricLoadCenterInverterPVWatts.new(model)
       ipvwatts.setName('PVSystem inverter')
-      ipvwatts.setInverterEfficiency(pv_system.inverter_efficiency)
+      ipvwatts.setInverterEfficiency(pv_system.inverter.inverter_efficiency)
 
       elcd.setInverter(ipvwatts)
     else