energyplus 0.1

data/lib/energyplus/SchedulesCSV.rb

@@ -0,0 +1,169 @@
+######################################################################
+#  Copyright (c) 2008-2013, Alliance for Sustainable Energy.
+#  All rights reserved.
+#
+#  This library is free software; you can redistribute it and/or
+#  modify it under the terms of the GNU Lesser General Public
+#  License as published by the Free Software Foundation; either
+#  version 2.1 of the License, or (at your option) any later version.
+#
+#  This library is distributed in the hope that it will be useful,
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+#  Lesser General Public License for more details.
+#
+#  You should have received a copy of the GNU Lesser General Public
+#  License along with this library; if not, write to the Free Software
+#  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+######################################################################
+
+require 'pathname'
+require 'csv'
+require 'Formats/Formats.rb'
+
+def make_schedule_name(s)
+  schedule_name = s
+  schedule_name = schedule_name.gsub(/[A-Z][a-z]/, ' \0').lstrip
+  schedule_name = schedule_name.gsub(/[a-z][A-Z]/) {|s| "" << s[0] << " " << s[1]}
+  schedule_name = schedule_name.gsub('_', ' ').gsub('-',' ').gsub(',',' ').gsub('()', '')
+  schedule_name = schedule_name.split(' ').each{|p| p.capitalize!}.join(' ')
+  return schedule_name
+end
+
+class Schedule
+  attr_accessor :name, :type, :table, :hours_per_year
+
+  def initialize(s)
+    @lines = s.split("\n")
+    @name = make_schedule_name(@lines[0])
+    @type = @lines[1]
+
+    @table = Array.new
+    this_table = Array.new
+    duration = ""
+    this_cols = 0
+    open = false
+    multi_duration = false
+    (2..@lines.length-1).each do |idx|
+      parts = CSV.parse_line(@lines[idx])
+
+      if parts[0].include?('Through')
+        if open
+            if (not this_table[-1][0] == 'Total Hours/Year')
+              this_table.push(Array.new)
+              this_table[-1].push('Total Hours/Year') if not multi_duration
+              (1..this_cols-1).each {|i| this_table[-1].push("") }
+            end
+            multi_duration = true
+            close_table(this_table)
+        end
+        this_table = Array.new
+        duration = parts[0]
+        open = true
+
+      elsif parts[0]==('Hour')
+        this_table.push(Array.new)
+        this_cols = parts.length
+        this_table[-1].push(parts[0]) if not multi_duration
+        (1..this_cols-1).each {|i| this_table[-1].push(parts[i] + " (" + duration + ")") }
+
+      elsif parts[0].include?('Total Hours/Day')
+        this_table.push(Array.new)
+        this_table[-1].push(parts[0]) if not multi_duration
+        (1..this_cols-1).each {|i| this_table[-1].push(custom_print(parts[i].to_f, 4)) }
+
+      elsif parts[0].include?('Total Hours/Week')
+        this_table.push(Array.new)
+        this_table[-1].push(parts[0]) if not multi_duration
+        this_table[-1].push(custom_print(parts[1].to_f, 4))
+        (2..this_cols-1).each {|i| this_table[-1].push("") }
+
+      elsif parts[0].include?('Total Hours/Year')
+        @hours_per_year = parts[1].to_f
+
+        this_table.push(Array.new)
+        this_table[-1].push(parts[0]) if not multi_duration
+        this_table[-1].push(custom_print(parts[1].to_f, 4))
+        (2..this_cols-1).each {|i| this_table[-1].push("") }
+      else
+        parts.delete_at(0) if multi_duration
+        this_table.push(parts)
+      end
+    end
+    close_table(this_table)
+    #print_table if multi_duration
+  end
+
+  def close_table(this_table)
+    if @table.empty?
+      @table = this_table
+    else
+      @table.each_index do |idx|
+        @table[idx].concat(this_table[idx])
+      end
+    end
+  end
+
+  def print_table
+    puts @name
+    puts @type
+    @table.each_index do |idx|
+      puts @table[idx].join(', ')
+    end
+    puts
+  end
+
+  def to_s
+    @lines.join("\n")
+  end
+  def hash
+    to_s.hash
+  end
+  def ==(other)
+    eql?(other)
+  end
+  def eql?(other)
+    other.is_a?(Schedule) && to_s == other.to_s
+  end
+end
+
+
+class SchedulesCSV
+  attr_accessor :schedules
+
+  def initialize(path)
+  @schedules = Hash.new
+    if File.exists?(path)
+      @path = Pathname.new(path)
+      process_schedules
+    else
+      @path = nil
+    end
+  end
+
+  def process_schedules
+    File.open(@path) do |f|
+      name = ""
+      schedule_string = ""
+      while line = f.gets
+
+        if line.chomp.empty?
+          if not name.empty?
+            new_schedule = Schedule.new(schedule_string)
+            @schedules[new_schedule.name] = new_schedule
+          end
+          name = ""
+          schedule_string = ""
+          next
+        end
+
+        if name.empty?
+          name = line.chomp
+          schedule_string = line
+        else
+          schedule_string += line
+        end
+      end
+    end
+  end
+end

data/lib/energyplus/SlabGtpFile.rb

@@ -0,0 +1,100 @@
+######################################################################
+#  Copyright (c) 2008-2013, Alliance for Sustainable Energy.
+#  All rights reserved.
+#
+#  This library is free software; you can redistribute it and/or
+#  modify it under the terms of the GNU Lesser General Public
+#  License as published by the Free Software Foundation; either
+#  version 2.1 of the License, or (at your option) any later version.
+#
+#  This library is distributed in the hope that it will be useful,
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+#  Lesser General Public License for more details.
+#
+#  You should have received a copy of the GNU Lesser General Public
+#  License along with this library; if not, write to the Free Software
+#  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+######################################################################
+
+module EnergyPlus
+  require 'pathname'
+
+  class SlabGtpFile
+    attr_accessor :convergence
+    attr_accessor :months
+    attr_accessor :averageTemps
+
+    def initialize(path)
+      @path = Pathname.new(path)
+      @convergence = false
+      @months = []
+      @averageTemps = []
+      init
+    end
+
+    def getGroundTemps
+      result = "Site:GroundTemperature:BuildingSurface,\n"
+      result+= "  #{@averageTemps[0]},  !- January Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[1]},  !- February Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[2]},  !- March Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[3]},  !- April Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[4]},  !- May Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[5]},  !- June Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[6]},  !- July Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[7]},  !- August Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[8]},  !- September Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[9]},  !- October Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[10]},  !- November Ground Temperature {C}\n"
+      result+= "  #{@averageTemps[11]};  !- December Ground Temperature {C}\n"
+      return result
+    end
+
+    private
+
+    def init
+      if @path.exist?
+        File.open(@path) do |f|
+          text = f.read
+          parse(text)
+        end
+      end
+    end
+
+    def parse(text)
+      # check for convergence
+      regex = /(Convergence has been gained.).*\n/
+      if text.match(regex) then
+        @convergence = true
+      end
+
+      #parse out the monthly values
+      startRegex = /\s*Month\s*TAverage\s*TPerimeter\s*TCore\s*TInside\s*AverageFlux\s*PerimeterFlux\s*CoreFlux/
+      readingTable = false
+      numRowsRead = 0
+      text.each_line do |line|
+        if readingTable
+
+          # split each row and count as read
+          data = line.strip.split(/\s+/)
+          numRowsRead += 1
+
+          # insert into members
+          @months << data[2]
+          @averageTemps << data[1]
+
+          # only read 12 months
+          if numRowsRead == 12
+            break
+          end
+        end
+
+        # have we started reading the table yet
+        if line.match(startRegex)
+          readingTable = true
+        end
+      end
+    end
+
+  end
+end

data/lib/energyplus/SlabIdfFile.rb

@@ -0,0 +1,78 @@
+######################################################################
+#  Copyright (c) 2008-2013, Alliance for Sustainable Energy.
+#  All rights reserved.
+#
+#  This library is free software; you can redistribute it and/or
+#  modify it under the terms of the GNU Lesser General Public
+#  License as published by the Free Software Foundation; either
+#  version 2.1 of the License, or (at your option) any later version.
+#
+#  This library is distributed in the hope that it will be useful,
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+#  Lesser General Public License for more details.
+#
+#  You should have received a copy of the GNU Lesser General Public
+#  License along with this library; if not, write to the Free Software
+#  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+######################################################################
+
+module EnergyPlus
+  require 'pathname'
+
+  class SlabIdfFile
+
+    def initialize(path)
+      @path = Pathname.new(path)
+      @initialized = false
+      @text = ""
+      init
+    end
+
+    def savefileas(path)
+       File.open(path, "w") do |f|
+         f << @text
+       end
+    end
+
+    # function to alter input based on options
+    def setSetpoints(cdd, hdd)
+      #options is an array of cdd and hdd... so have to wait the values
+      #create a new array with the weighted values
+      heatSP = 21  #in deg c
+      coolSP = 24  #in deg c
+
+      weight = []
+      for i in 1..12
+        num = cdd * coolSP + hdd * heatSP
+        den = cdd + hdd
+        weight[i] = (num) / (den)
+      end
+
+      @text.gsub!("%JAN%", sprintf("%.2f", weight[1].to_s))
+      @text.gsub!("%FEB%", sprintf("%.2f", weight[2].to_s))
+      @text.gsub!("%MAR%", sprintf("%.2f", weight[3].to_s))
+      @text.gsub!("%APR%", sprintf("%.2f", weight[4].to_s))
+      @text.gsub!("%MAY%", sprintf("%.2f", weight[5].to_s))
+      @text.gsub!("%JUN%", sprintf("%.2f", weight[6].to_s))
+      @text.gsub!("%JUL%", sprintf("%.2f", weight[7].to_s))
+      @text.gsub!("%AUG%", sprintf("%.2f", weight[8].to_s))
+      @text.gsub!("%SEP%", sprintf("%.2f", weight[9].to_s))
+      @text.gsub!("%OCT%", sprintf("%.2f", weight[10].to_s))
+      @text.gsub!("%NOV%", sprintf("%.2f", weight[11].to_s))
+      @text.gsub!("%DEC%", sprintf("%.2f", weight[12].to_s))
+    end
+
+    private
+
+    def init
+      @initialized = true
+
+      if @path.exist?
+        File.open(@path) do |f|
+          @text = f.read
+        end
+      end
+    end
+  end
+end #module

data/lib/energyplus/StatFile.rb

@@ -0,0 +1,192 @@
+######################################################################
+#  Copyright (c) 2008-2013, Alliance for Sustainable Energy.
+#  All rights reserved.
+#
+#  This library is free software; you can redistribute it and/or
+#  modify it under the terms of the GNU Lesser General Public
+#  License as published by the Free Software Foundation; either
+#  version 2.1 of the License, or (at your option) any later version.
+#
+#  This library is distributed in the hope that it will be useful,
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+#  Lesser General Public License for more details.
+#
+#  You should have received a copy of the GNU Lesser General Public
+#  License along with this library; if not, write to the Free Software
+#  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+######################################################################
+
+module EnergyPlus
+
+  require 'pathname'
+
+  class StatFile
+    attr_accessor :valid
+    attr_accessor :lat
+    attr_accessor :lon
+    attr_accessor :elevation
+    attr_accessor :gmt
+    attr_accessor :monthlyDB
+    attr_accessor :hdd18
+    attr_accessor :cdd18
+    attr_accessor :hdd10
+    attr_accessor :cdd10
+
+    def initialize(path)
+      @path = Pathname.new(path)
+      @valid = false
+      @lat = []
+      @lon = []
+      @gmt = []
+      @elevation = []
+      @hdd18 = []
+      @cdd18 = []
+      @hdd10 = []
+      @cdd10 = []
+      @monthlyDB = []
+      @deltaDB = []
+      init
+    end
+
+    def valid?
+      return @valid
+    end
+
+    # the mean of the mean monthly dry bulbs
+    def meanDB
+      if not @monthlyDB.empty? then
+        thisSum = 0
+        @monthlyDB.each { |db| thisSum += db }
+        thisMean = thisSum/@monthlyDB.size
+      else
+        thisMean = ""
+      end
+      return thisMean
+    end
+
+    # max - min of the mean monthly dry bulbs
+    def deltaDB
+      if not @monthlyDB.empty? then
+        deltaT = @monthlyDB.max-@monthlyDB.min
+      else
+        deltaT = ""
+      end
+
+      return deltaT
+    end
+
+    private
+
+    # initialize
+    def init
+      if @path.exist?
+        File.open(@path) do |f|
+          text = f.read
+          parse(text)
+        end
+      end
+    end
+
+    def parse(text)
+
+      # get lat, lon, gmt
+      regex = /\{(N|S)\s*([0-9]*).\s*([0-9]*)'\}\s*\{(E|W)\s*([0-9]*).\s*([0-9]*)'\}\s*\{GMT\s*(.*)\s*Hours\}/
+      matchData = text.match(regex)
+      if matchData.nil?
+        puts "Can't find lat/lon/gmt"
+        return
+      else
+
+        @lat = matchData[2].to_f + (matchData[3].to_f)/60.0
+        if matchData[1] == 'S'
+          @lat = -@lat
+        end
+
+        @lon = matchData[5].to_f + (matchData[6].to_f)/60.0
+        if matchData[4] == 'W'
+          @lon = -@lon
+        end
+
+        @gmt = matchData[7]
+
+      end
+
+      # get elevation
+      regex = /Elevation --\s*(.*)m (above|below) sea level/
+      matchData = text.match(regex)
+      if matchData.nil?
+        puts "Can't find elevation"
+        return
+      else
+        @elevation = matchData[1].to_f
+        if matchData[2] == 'below'
+          @elevation = -@elevation
+        end
+      end
+
+
+      # get heating and cooling degree days
+      cdd10Regex = /-\s*(.*) annual \(standard\) cooling degree-days \(10.C baseline\)/
+      matchData = text.match(cdd10Regex)
+      if matchData.nil?
+        puts "Can't find CDD 10"
+        return
+      else
+        @cdd10 = matchData[1].to_f
+      end
+
+      hdd10Regex = /-\s*(.*) annual \(standard\) heating degree-days \(10.C baseline\)/
+      matchData = text.match(hdd10Regex)
+      if matchData.nil?
+        puts "Can't find HDD 10"
+        return
+      else
+        @hdd10 = matchData[1].to_f
+      end
+
+      cdd18Regex = /-\s*(.*) annual \(standard\) cooling degree-days \(18.3.C baseline\)/
+      matchData = text.match(cdd18Regex)
+      if matchData.nil?
+        puts "Can't find CDD 18"
+        return
+      else
+        @cdd18 = matchData[1].to_f
+      end
+
+      hdd18Regex = /-\s*(.*) annual \(standard\) heating degree-days \(18.3.C baseline\)/
+      matchData = text.match(hdd18Regex)
+      if matchData.nil?
+        puts "Can't find HDD 18"
+        return
+      else
+        @hdd18 = matchData[1].to_f
+      end
+
+
+      #use regex to get the temperatures
+      regex = /Daily Avg(.*)\n/
+      matchData = text.match(regex)
+      if matchData.nil?
+        puts "Can't find outdoor air temps"
+        return
+      else
+        # first match is outdoor air temps
+        monthlyTemps = matchData[1].strip.split(/\s+/)
+
+        # have to be 12 months
+        if monthlyTemps.size != 12
+          puts "Can't find outdoor air temps"
+          return
+        end
+
+        # insert as numbers
+        monthlyTemps.each {|temp| @monthlyDB << temp.to_f}
+      end
+
+      # now we are valid
+      @valid = true
+    end
+
+  end
+end