PVLIB_Ruby 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 7d3ae3a823ad98dd6eaad6afeea6fc3cc3ff5902
+  data.tar.gz: 9791749f4b811b64481206aec941f5289d1d543c
+SHA512:
+  metadata.gz: 8043fbd98a76d68a1deb94449201eb34fca2365d54c5055878637531f62b445d236f810c85245de32b6774a25dd34fe3aa49e02102bc4490df4c35445ee9d0d2
+  data.tar.gz: d920617b1a837ffd83d763d3b921ab0bee605aa8aa575564eb390a22f63613c5f587d9e7e628fe33636c05895f45e04c0f9eac33a7bf6f1a0f3085b146406c2e

data/.gitignore

@@ -0,0 +1,10 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+.DS_Store

data/.rspec

@@ -0,0 +1,2 @@
+--format documentation
+--color

data/.ruby-gemset

@@ -0,0 +1 @@
+PVLIB_Ruby

data/.ruby-version

@@ -0,0 +1 @@
+2.2.2

data/.travis.yml

@@ -0,0 +1,3 @@
+language: ruby
+rvm:
+  - 2.2.2

data/CODE_OF_CONDUCT.md

@@ -0,0 +1,13 @@
+# Contributor Code of Conduct
+
+As contributors and maintainers of this project, we pledge to respect all people who contribute through reporting issues, posting feature requests, updating documentation, submitting pull requests or patches, and other activities.
+
+We are committed to making participation in this project a harassment-free experience for everyone, regardless of level of experience, gender, gender identity and expression, sexual orientation, disability, personal appearance, body size, race, age, or religion.
+
+Examples of unacceptable behavior by participants include the use of sexual language or imagery, derogatory comments or personal attacks, trolling, public or private harassment, insults, or other unprofessional conduct.
+
+Project maintainers have the right and responsibility to remove, edit, or reject comments, commits, code, wiki edits, issues, and other contributions that are not aligned to this Code of Conduct. Project maintainers who do not follow the Code of Conduct may be removed from the project team.
+
+Instances of abusive, harassing, or otherwise unacceptable behavior may be reported by opening an issue or contacting one or more of the project maintainers.
+
+This Code of Conduct is adapted from the [Contributor Covenant](http:contributor-covenant.org), version 1.0.0, available at [http://contributor-covenant.org/version/1/0/0/](http://contributor-covenant.org/version/1/0/0/)

data/Gemfile

@@ -0,0 +1,8 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in PVLIB_Ruby.gemspec
+gemspec
+
+group :development, :test do
+  gem 'pry-byebug'
+end

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2015 Tadatoshi Takahashi
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/PVLIB_Ruby.gemspec

@@ -0,0 +1,35 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'pvlib_ruby/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "PVLIB_Ruby"
+  spec.version       = PvlibRuby::VERSION
+  spec.authors       = ["Tadatoshi Takahashi"]
+  spec.email         = ["tadatoshi@gmail.com"]
+
+  spec.summary       = %q{PVLIB written in Ruby.}
+  spec.description   = %q{PVLIB_MatLab translated to Ruby code in Object-Oriented way.}
+  spec.homepage      = "https://github.com/tadatoshi/PVLIB_Ruby"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org by setting 'allowed_push_host', or
+  # delete this section to allow pushing this gem to any host.
+  if spec.respond_to?(:metadata)
+    spec.metadata['allowed_push_host'] = "https://rubygems.org"
+  else
+    raise "RubyGems 2.0 or newer is required to protect against public gem pushes."
+  end
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_dependency "activemodel"
+
+  spec.add_development_dependency "bundler", "~> 1.9"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec"
+end

data/README.md

@@ -0,0 +1,41 @@
+# PVLIB_Ruby
+
+PVLIB_MatLab translated to Ruby code in Object-Oriented way.
+
+Since the naming convention is PVLIB_[language name], e.g. PVLIB_MatLab, PVLIB_Python, this gem is following the same naming convention, i.e. PVLIB_Ruby. 
+
+However, internally, this gem is following Ruby naming convention. Hence, the code is under 'lib/pvlib_ruby'. 
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'PVLIB_Ruby'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install PVLIB_Ruby
+
+## Usage
+
+See examples under examples directory. 
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `bin/console` for an interactive prompt that will allow you to experiment.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release` to create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+1. Fork it ( https://github.com/[my-github-username]/PVLIB_Ruby/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1 @@
+require "bundler/gem_tasks"

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "pvlib_ruby"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start

data/bin/setup

@@ -0,0 +1,7 @@
+#!/bin/bash
+set -euo pipefail
+IFS=$'\n\t'
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/examples/data/sandia_inverter_example.csv

@@ -0,0 +1,3 @@
+"name","Vac","Pac0","Pdc0","Vdc0","Ps0","C0","C1","C2","C3","Pnt","Vdcmax","Idcmax","MPPTLow","MPPTHi","LibraryType","LibraryName",
+"SunPower Corp (Original Mfg - PV Powered): SPR-2500 240V [CEC 2006]",240,2500,2630.0909,219.2317,41.3618,-1.265e-05,6.6747e-05,0.0017417,0.00061375,3.9,500,20,140,450,"SandiaInverter","Sandia Inverters",
+

data/examples/data/sandia_module_example.csv

@@ -0,0 +1,3 @@
+"name","vintage","material","area","AlphaIsc","AlphaImp","Isc0","Imp0","Voc0","Vmp0","BetaVoc","BetaVmp","mBetaVoc","mBetaVmp","Ns","Np","delT","fd","n","Ix0","Ixx0","a_wind","b_wind","c",,,,,,,,"a",,,,,"b",,,,,,
+"Sample Module",2006,"c-Si",1.244,0.000232,-0.00036,5.988,5.56,48.53,40.03,-0.152,-0.162,0,0,72,1,3,1,1.241,5.93,4.12,-3.62,-0.075,1.0072,-0.0072,0.32304,-3.4984,0.9966,0.0034,1.0827,-0.0827,-0.0001223,0.002416,-0.01912,0.07365,0.9259,-2.99e-09,5.35e-07,-3.4e-05,0.000862,-0.00699,1,
+

data/examples/pv_power_generation_example.rb

@@ -0,0 +1,120 @@
+# NOTE: PVLIB_Ruby gem must be installed in local machine in order to run this code. 
+#       Execute "rake install" to use the code in this gem. 
+#       Or if you use the gem downloaded from RubyGems.org, you don't have to execute "rake install" and just do "gem install PVLIB_Ruby"
+
+require 'pvlib_ruby'
+
+require 'bigdecimal'
+
+DATA_DIRECTORY = File.expand_path('../data', __FILE__)
+
+# for 360th row in PVSC40Tutorial_Master (360/30 = 12, i.e. noon, note: measurement is every two minutes)
+direct_normal_irradiance = BigDecimal('631.3100')
+global_horizontal_irradiance = BigDecimal('862.0619')
+diffuse_horizontal_irradiance = BigDecimal('405.3100')
+day_of_year = BigDecimal('294')
+utc_offset = '-07:00'
+# TODO: Get time from year and day_of_year:
+time = Time.new(2008, 10, 20, 11, 58, 12, utc_offset)
+albedo = BigDecimal('0.1500')
+array_tilt = BigDecimal('35')
+array_azimuth = BigDecimal('180')
+latitude = BigDecimal('35.0500')
+longitude = BigDecimal('-106.5400')
+altitude = BigDecimal('1660')
+location = Location.new(latitude, longitude, altitude)
+pressure = BigDecimal('62963')
+reference_solar_irradiance = BigDecimal('1000')
+wind_speed = BigDecimal('2.8786')
+air_temperature = BigDecimal('20.7700')
+
+solar_ephemeris = SolarEphemeris.new(time, location, pressure: pressure, temperature: air_temperature)
+
+sun_azimuth = solar_ephemeris.sun_azimuth
+apparent_sun_elevation = solar_ephemeris.apparent_sun_elevation
+
+puts "------------- Solar Ephemeris -------------"
+puts "  Sun Azimuth [º]: #{sun_azimuth.round(4).to_s('F')}" # Slight difference from 182.5229 for 360th row in PVSC40Tutorial_Master
+puts "  Apparent Sun Elevation [º]: #{apparent_sun_elevation.round(4).to_s('F')}" # Slight difference from 44.2308 for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+sun_zenith = BigDecimal('90') - apparent_sun_elevation
+
+plain_of_array_irradiance = PlainOfArrayIrradiance.new(direct_normal_irradiance, global_horizontal_irradiance, diffuse_horizontal_irradiance, day_of_year, albedo, array_tilt, array_azimuth, sun_zenith, sun_azimuth)    
+
+angle_of_incidence = plain_of_array_irradiance.angle_of_incidence
+
+beam_irradiance = plain_of_array_irradiance.beam_irradiance
+ground_diffuse_irradiance = plain_of_array_irradiance.ground_diffuse_irradiance
+sky_diffuse_irradiance = plain_of_array_irradiance.sky_diffuse_irradiance
+
+puts "----- Plain Of Array (POA) Irradiance -----"
+puts "  Angle of Incidence [º]: #{angle_of_incidence.round(4).to_s('F')}" # Slight difference from 10.8703 for 360th row in PVSC40Tutorial_Master
+puts ""
+puts "  Beam Irradiance [W/m^2]: #{beam_irradiance.round(4).to_s('F')}" # Slight difference from 619.9822 for 360th row in PVSC40Tutorial_Master
+puts "  Ground Diffuse Irradiance [W/m^2]: #{ground_diffuse_irradiance.round(4).to_s('F')}" # Matches to 11.6927 for 360th row in PVSC40Tutorial_Master
+puts "  Sky Diffuse Irradiance [W/m^2]: #{sky_diffuse_irradiance.round(4).to_s('F')}" # Slight difference from 464.8004 for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+solar_irradiance_incident_on_module_surface = beam_irradiance + ground_diffuse_irradiance + sky_diffuse_irradiance
+
+sandia_module_data_filepath = File.join(DATA_DIRECTORY, 'sandia_module_example.csv')
+sandia_pv_module = PvModule.create(sandia_module_data_filepath)
+
+pv_temperature = PvTemperature.new(sandia_pv_module, solar_irradiance_incident_on_module_surface, reference_solar_irradiance, wind_speed, air_temperature)
+
+puts "------------- PV temperature --------------"
+puts "  Estimated Cell Temperature [ºC]: #{pv_temperature.cell_temperature.round(4).to_s('F')}" # Slight difference from 47.7235 for 360th row in PVSC40Tutorial_Master
+puts "  Estimated Module Temperature [ºC]: #{pv_temperature.module_temperature.round(4).to_s('F')}" # Slight difference from 44.4341 for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+estimated_cell_temperature = pv_temperature.cell_temperature
+
+air_mass = AirMass.new(sun_zenith, pressure)
+
+puts "---------------- Air Mass -----------------"
+puts "  Absolute Air Mass: #{air_mass.absolute_air_mass.round(4).to_s('F')}" # Slight difference from 0.8894 for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+absolute_air_mass = air_mass.absolute_air_mass
+
+pv_performance_characterization = PvPerformanceCharacterization.new(sandia_pv_module, absolute_air_mass, angle_of_incidence, 
+                                                                    beam_irradiance, ground_diffuse_irradiance, sky_diffuse_irradiance, 
+                                                                    estimated_cell_temperature)
+
+puts "--- PV module IV curve characterization ---"
+puts "  Short Circuit Current [A]: #{pv_performance_characterization.short_circuit_current.round(4).to_s('F')}" # Slight difference from 6.4151 for 360th row in PVSC40Tutorial_Master
+puts "  Maximum Power Point Current [A]: #{pv_performance_characterization.maximum_power_point_current.round(4).to_s('F')}" # Matches to 5.8741 for 360th row in PVSC40Tutorial_Master
+puts "  Open Circuit Voltage [V]: #{pv_performance_characterization.open_circuit_voltage.round(4).to_s('F')}" # Matches to 45.2333 for 360th row in PVSC40Tutorial_Master
+puts "  Maximum Power Point Voltage [V]: #{pv_performance_characterization.maximum_power_point_voltage.round(4).to_s('F')}" # Matches to 36.3984 for 360th row in PVSC40Tutorial_Master
+puts "  Fourth Point Current [A]: #{pv_performance_characterization.fourth_point_current.round(4).to_s('F')}" # Matches to 6.3544 for 360th row in PVSC40Tutorial_Master
+puts "  Fifth Point Current [A]: #{pv_performance_characterization.fifth_point_current.round(4).to_s('F')}" # Slight difference from 4.3899 for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+parallel_string = BigDecimal('1')
+series_modules = BigDecimal('5')
+array_current = pv_performance_characterization.maximum_power_point_current * parallel_string
+array_voltage = pv_performance_characterization.maximum_power_point_voltage * series_modules
+array_power = array_current * array_voltage
+
+puts "------------- PV array output -------------"
+puts "  Array current [A]: #{array_current.round(4).to_s('F')}" # Matches to 5.8741 for 360th row in PVSC40Tutorial_Master
+puts "  Array voltage [V]: #{array_voltage.round(4).to_s('F')}" # Slight difference from 181.9919 for 360th row in PVSC40Tutorial_Master
+puts "  Array power [W]: #{array_power.round(4).to_s('F')}" # Matches to 1069.0 (with rounding) for 360th row in PVSC40Tutorial_Master
+puts "-------------------------------------------"
+puts ""
+
+sandia_inverter_data_filename = File.join(DATA_DIRECTORY, 'sandia_inverter_example.csv')
+sandia_inverter = Inverter.create(sandia_inverter_data_filename)
+
+dc_to_ac_conversion = DcToAcConversion.new(sandia_inverter)
+
+puts "--- Inverter DC to AC conversion ---"
+puts "  AC power [W]: #{dc_to_ac_conversion.ac_power(array_voltage, array_power).round(4).to_s('F')}" # Matches to 1016.3 (with rounding) for 360th row in PVSC40Tutorial_Master
+puts "------------------------------------"
+puts ""

data/lib/PVLIB_Ruby/air_mass.rb

@@ -0,0 +1,26 @@
+require 'bigdecimal'
+
+# The default model for relative air mass is 'kastenyoung1989'. Different model is realized by subclass. 
+# 'kastenyoung1989' model is based on:
+# Fritz Kasten and Andrew Young. "Revised optical air mass tables and approximation formula". Applied Optics 28:4735–4738
+class AirMass
+  include CalculationHelper
+
+  def initialize(sun_zenith, pressure)
+    @sun_zenith = sun_zenith
+    @pressure = pressure
+  end
+
+  # Relative air mass at sea level
+  def relative_air_mass
+    relative_air_mass = BigDecimal('1') / (Math.cos(degree_to_radian(@sun_zenith).to_f) + BigDecimal('0.50572') * ((BigDecimal('6.07995') + (BigDecimal('90') - @sun_zenith)).power(-1.6364)))
+    relative_air_mass = BigDecimal('0') if relative_air_mass.nan?
+    relative_air_mass
+  end
+
+  # Airmass for locations not at sea level (i.e. not at standard pressure)
+  def absolute_air_mass
+    relative_air_mass * @pressure / BigDecimal('101325')
+  end
+
+end

data/lib/PVLIB_Ruby/calculation_helper.rb

@@ -0,0 +1,43 @@
+module CalculationHelper
+
+  def degree_to_radian(degree)
+    degree * BigDecimal(Math::PI.to_s) / BigDecimal('180')
+  end
+
+  def radian_to_degree(radian)
+    radian * BigDecimal('180') / BigDecimal(Math::PI.to_s)
+  end
+
+  def bigdecimal_exp(bigdecimal_value)
+    BigDecimal(Math.exp(bigdecimal_value.to_f).to_s)
+  end
+
+  def bigdecimal_cos(angle_in_radian)
+    BigDecimal(Math.cos(angle_in_radian.to_f).to_s)
+  end 
+
+  def bigdecimal_sin(angle_in_radian)
+    BigDecimal(Math.sin(angle_in_radian.to_f).to_s)
+  end
+
+  def bigdecimal_tan(angle_in_radian)
+    BigDecimal(Math.tan(angle_in_radian.to_f).to_s)
+  end
+  
+  def bigdecimal_asin(sin_value)
+    BigDecimal(Math.asin(sin_value.to_f).to_s)
+  end
+
+  def bigdecimal_acos(cos_value)
+    BigDecimal(Math.acos(cos_value.to_f).to_s)
+  end
+
+  def bigdecimal_atan2(y, x)
+    BigDecimal(Math.atan2(y.to_f, x.to_f).to_s)
+  end
+
+  def bigdecimal_sqrt(bigdecimal_value)
+    BigDecimal(Math.sqrt(bigdecimal_value.to_f).to_s)
+  end
+
+end

data/lib/PVLIB_Ruby/dc_to_ac_conversion.rb

@@ -0,0 +1,40 @@
+require 'bigdecimal'
+
+# This class uses Sandia's Grid-Connected PV Inverter model, which is based on:
+# SAND2007-5036, "Performance Model for Grid-Connected Photovoltaic Inverters by D. King, S. Gonzalez, G. Galbraith, W. Boyson
+# Since this paper was not available any more, we referred to https://pvpmc.sandia.gov/modeling-steps/dc-to-ac-conversion/sandia-inverter-model/
+# Other inverter models are realized by subclasses. 
+class DcToAcConversion
+
+  def initialize(inverter)  
+    @inverter = inverter
+  end
+
+  def ac_power(dc_input_voltage, dc_input_power)
+    # A = Pdc0 * ( 1 + C1 * (Vdc - Vdc0) )
+    dc_power_adjustment_parameter = @inverter.dc_power_for_maximum_ac_power_rating * (1 + @inverter.coefficient_for_dc_power_for_maximum_ac_power_rating * (dc_input_voltage - @inverter.dc_power_level_for_maximum_ac_power_rating))
+
+    # B = Ps0 * ( 1 + C2 * (Vdc - Vdc0) )
+    starting_dc_power_adjustment_parameter = @inverter.starting_dc_power * (1 + @inverter.coefficient_for_starting_dc_power * (dc_input_voltage - @inverter.dc_power_level_for_maximum_ac_power_rating))
+
+    # C = C0 * ( 1 + C3 * (Vdc - Vdc0) )
+    coefficient_adjustment_parameter = @inverter.power_adjustment_coefficient * (1 + @inverter.coefficient_for_power_adjustment_coefficient * (dc_input_voltage - @inverter.dc_power_level_for_maximum_ac_power_rating))
+
+    ac_power = (@inverter.maximum_ac_power_rating / (dc_power_adjustment_parameter - starting_dc_power_adjustment_parameter) - coefficient_adjustment_parameter * (dc_power_adjustment_parameter - starting_dc_power_adjustment_parameter)) * (dc_input_power - starting_dc_power_adjustment_parameter) + 
+    coefficient_adjustment_parameter * (dc_input_power - starting_dc_power_adjustment_parameter).power(2)
+
+    adjust_ac_power(ac_power)
+  end
+
+  private
+    def adjust_ac_power(ac_power)
+      if ac_power > @inverter.maximum_ac_power_rating # Inverter clipping at maximum rated AC Power
+        maximum_ac_power_rating
+      elsif ac_power < @inverter.starting_dc_power # Inverter night tare losses
+        -(@inverter.night_tare_loss.abs)
+      else
+        ac_power
+      end
+    end 
+
+end

data/lib/PVLIB_Ruby/models/active_csv.rb

@@ -0,0 +1,60 @@
+require 'active_model'
+require 'csv'
+
+# This is a model in model-view-controller design pattern. 
+# It is analogous to ActiveRecord. Hence, it's a placeholder for data from CSV file. 
+# Borrowed the idea from "http://kyle.conarro.com/importing-from-csv-in-rails".
+class ActiveCsv
+  include ActiveModel::Model
+
+  def self.create(csv_filepath)
+    new_instance = self.new
+    new_instance.load_data(csv_filepath)
+    new_instance
+  end
+
+  def persisted?
+    false
+  end
+
+  def valid?
+    # TODO: implement this
+  end
+
+  def load_data(csv_filepath)
+    csv = CSV.new(File.new(csv_filepath), headers: true, header_converters: :symbol)
+
+    data = csv.shift
+
+    previous_header = nil;
+    index = 0;
+
+    data.headers.each do |header|
+
+      if "#{header}=" == "="
+
+        unless previous_header.blank?
+
+          assigned_data_for_previous_header = self.send(previous_header.to_sym)
+          
+          if assigned_data_for_previous_header.instance_of? Array
+            assigned_data_for_previous_header << data[index]
+            self.send("#{previous_header}=".to_sym, assigned_data_for_previous_header)
+          else
+            self.send("#{previous_header}=".to_sym, [assigned_data_for_previous_header,data[index]])
+          end
+
+        end
+
+      else
+        self.send("#{header}=".to_sym, data[header])
+      end
+
+      previous_header = header unless header.blank?
+      index = index + 1;
+
+    end
+
+  end  
+
+end

data/lib/PVLIB_Ruby/models/inverter.rb

@@ -0,0 +1,46 @@
+require 'bigdecimal'
+
+# This is a model in model-view-controller design pattern. 
+# In this case, it is expending ActiveCsv, i.e. it is analogous to ActiveRecord. Hence, it's a placeholder for data from CSV file. 
+class Inverter < ActiveCsv
+  
+  attr_accessor :name, :vac, :pac0, :pdc0, :vdc0, :ps0, :c0, :c1, :c2, :c3 
+  attr_accessor :pnt, :vdcmax, :idcmax, :mpptlow, :mppthi, :librarytype, :libraryname
+
+  def maximum_ac_power_rating
+    BigDecimal(self.pac0.to_s)
+  end
+
+  def dc_power_for_maximum_ac_power_rating
+    BigDecimal(self.pdc0.to_s)
+  end
+
+  def coefficient_for_dc_power_for_maximum_ac_power_rating
+    BigDecimal(self.c1.to_s)
+  end
+
+  def dc_power_level_for_maximum_ac_power_rating
+    BigDecimal(self.vdc0.to_s)
+  end
+
+  def starting_dc_power
+    BigDecimal(self.ps0.to_s)
+  end
+
+  def coefficient_for_starting_dc_power
+    BigDecimal(self.c2.to_s)
+  end
+
+  def power_adjustment_coefficient
+    BigDecimal(self.c0.to_s)
+  end
+
+  def coefficient_for_power_adjustment_coefficient
+    BigDecimal(self.c3.to_s)
+  end
+
+  def night_tare_loss
+    BigDecimal(self.pnt.to_s)
+  end
+
+end

data/lib/PVLIB_Ruby/models/location.rb

@@ -0,0 +1,11 @@
+class Location
+
+  attr_accessor :latitude, :longitude, :altitude
+
+  def initialize(latitude, longitude, altitude)
+    @latitude = latitude
+    @longitude = longitude
+    @altitude = altitude
+  end
+
+end

data/lib/PVLIB_Ruby/models/pv_module.rb

@@ -0,0 +1,100 @@
+require 'bigdecimal'
+
+# This is a model in model-view-controller design pattern. 
+# In this case, it is expending ActiveCsv, i.e. it is analogous to ActiveRecord. Hence, it's a placeholder for data from CSV file. 
+class PvModule < ActiveCsv
+
+  attr_accessor :name, :vintage, :material, :area, :alphaisc, :alphaimp, :isc0, :imp0, :voc0, :vmp0 
+  attr_accessor :betavoc, :betavmp, :mbetavoc, :mbetavmp, :ns, :np, :delt, :fd, :n, :ix0, :ixx0, :a_wind, :b_wind
+  attr_accessor :c, :a, :b
+
+  def short_circuit_current_0
+    BigDecimal(self.isc0.to_s)
+  end
+
+  def air_mass_variation_polynomial_coefficients
+    convert_coefficients_to_bigdecimal(self.a)
+  end
+
+  def angle_of_incidence_polynominal_coefficients
+    convert_coefficients_to_bigdecimal(self.b)
+  end
+
+  def diffuse_irradiance_factor
+    BigDecimal(self.fd)
+  end
+
+  def normalized_temperature_coefficient_for_short_circuit_current
+    BigDecimal(self.alphaisc.to_s)
+  end
+
+  def maximum_power_point_current_0
+    BigDecimal(self.imp0.to_s)
+  end
+
+  def performance_coefficients
+    convert_coefficients_to_bigdecimal(self.c)
+  end
+
+  def normalized_temperature_coefficient_for_maximum_power_point_current
+    BigDecimal(self.alphaimp.to_s)
+  end
+
+  def open_circuit_voltage_0
+    BigDecimal(self.voc0.to_s)
+  end
+
+  def number_of_cells_in_series
+    BigDecimal(self.ns.to_s)
+  end
+
+  def diode_factor
+    BigDecimal(self.n.to_s)
+  end
+
+  def temperature_coefficient_for_open_circuit_voltage
+    BigDecimal(self.betavoc.to_s)
+  end
+
+  def irradiance_dependent_temperature_coefficient_for_open_circuit_voltage
+    BigDecimal(self.mbetavoc.to_s)
+  end
+
+  def maximum_power_point_voltage_0
+    BigDecimal(self.vmp0.to_s)
+  end
+
+  def temperature_coefficient_for_maximum_power_point_voltage
+    BigDecimal(self.betavmp.to_s)
+  end
+
+  def irradiance_dependent_temperature_coefficient_for_maximum_power_point_voltage
+    BigDecimal(self.mbetavmp.to_s)
+  end
+
+  def fourth_point_current_0
+    BigDecimal(self.ix0.to_s)
+  end
+
+  def fifth_point_current_0
+    BigDecimal(self.ixx0.to_s)
+  end
+
+  def temperature_upper_limit_coefficient
+    BigDecimal(self.a_wind.to_s)
+  end
+
+  def wind_speed_decrease_rate_coefficient
+    BigDecimal(self.b_wind.to_s)
+  end
+
+  def cell_module_back_surface_temperature_difference
+    BigDecimal(self.delt.to_s)
+  end
+
+  private
+    def convert_coefficients_to_bigdecimal(coefficients)
+      coefficients.map {|coefficient| BigDecimal(coefficient.to_s)}.compact
+    end
+
+end

data/lib/PVLIB_Ruby/plain_of_array_irradiance.rb

@@ -0,0 +1,125 @@
+require 'bigdecimal'
+
+# For Sky Diffuse Irradiance, "Reindl's 1990 model" is used. 
+# Other models are realized by subclasses. 
+#
+# References
+#   [1] Loutzenhiser P.G. et. al. "Empirical validation of models to compute
+#   solar irradiance on inclined surfaces for building energy simulation"
+#   2007, Solar Energy vol. 81. pp. 254-267
+#   [2] Reindl, D.T., Beckmann, W.A., Duffie, J.A., 1990a. Diffuse fraction
+#   correlations. Solar Energy 45 (1), 1–7.
+#   [3] Reindl, D.T., Beckmann, W.A., Duffie, J.A., 1990b. Evaluation of hourly
+#   tilted surface radiation models. Solar Energy 45 (1), 9–17.
+#   [4] Solar radiation basics http://solardat.uoregon.edu/SolarRadiationBasics.html
+# 
+# TODO: Consider to pass necessary parameters to the respective methods instead of passing all of them to the constructor. 
+#       Haven't decided which way would be better. 
+class PlainOfArrayIrradiance
+  include CalculationHelper
+
+  AVEARGE_EXTRATERRESTRIAL_IRRADIANCE = BigDecimal('1367') # [W/m^2]
+  SMALL_VALUE_FOR_SKY_DIFFUSE_IRRADIANCE = BigDecimal('0.000001') # In order to make the calculated value consistent with the one by PVLIB_MatLab
+
+  def initialize(direct_normal_irradiance, global_horizontal_irradiance, diffuse_horizontal_irradiance, day_of_year, albedo, surface_tilt, surface_azimuth, sun_zenith, sun_azimuth)
+    @direct_normal_irradiance = direct_normal_irradiance           # DNI [W/m^2]
+    @global_horizontal_irradiance = global_horizontal_irradiance   # GHI [W/m^2]
+    @diffuse_horizontal_irradiance = diffuse_horizontal_irradiance # DHI [W/m^2]
+    @day_of_year = day_of_year
+    @albedo = albedo
+    @surface_tilt = surface_tilt
+    @surface_azimuth = surface_azimuth
+    @sun_zenith = sun_zenith
+    @sun_azimuth = sun_azimuth
+  end
+
+  def beam_irradiance
+    @direct_normal_irradiance * bigdecimal_cos(degree_to_radian(angle_of_incidence))
+  end
+
+  def ground_diffuse_irradiance
+    @global_horizontal_irradiance * @albedo * (1 - bigdecimal_cos(degree_to_radian(@surface_tilt))) * BigDecimal('0.5')
+  end
+
+  # Based on the equation (5) of [3]
+  # The comment in PVLIB_MatLab says it's based on the modified version of the equation (8) of [1]. 
+  # But I observe that the code is exactly same as the equation (5) of [3]. 
+  def sky_diffuse_irradiance
+    @diffuse_horizontal_irradiance * (anisotropy_index * geometric_factor + 
+                                      (1 - anisotropy_index) * ((BigDecimal('1') + bigdecimal_cos(degree_to_radian(@surface_tilt))) / BigDecimal('2')) * anisotropic_correction_factor)
+  end
+
+  # This is used only for calculating Sky Diffuse Irradiance. 
+  # But it's made public method just in case some other code needs it. 
+  # Reference: [4]
+  def extraterrestrial_irradiance
+    earth_orbit_angle_in_radian = BigDecimal('2') * BigDecimal(Math::PI.to_s) * @day_of_year / BigDecimal('365')
+    sun_earth_distance_factor_square = BigDecimal('1.00011') + BigDecimal('0.034221') * bigdecimal_cos(earth_orbit_angle_in_radian) + BigDecimal('0.001280') * bigdecimal_sin(earth_orbit_angle_in_radian) +
+                                       BigDecimal('0.000719') * bigdecimal_cos(BigDecimal('2') * earth_orbit_angle_in_radian) + BigDecimal('0.000077') * bigdecimal_sin(BigDecimal('2') * earth_orbit_angle_in_radian)
+    AVEARGE_EXTRATERRESTRIAL_IRRADIANCE * sun_earth_distance_factor_square
+  end
+
+  def angle_of_incidence
+    cosine_of_angle_of_incidence = bigdecimal_cos(degree_to_radian(@sun_zenith)) * bigdecimal_cos(degree_to_radian(@surface_tilt)) + 
+                                   bigdecimal_sin(degree_to_radian(@surface_tilt)) * bigdecimal_sin(degree_to_radian(@sun_zenith)) * bigdecimal_cos(degree_to_radian(@sun_azimuth - @surface_azimuth))
+
+    if cosine_of_angle_of_incidence > BigDecimal('1')
+      cosine_of_angle_of_incidence = BigDecimal('1')
+    elsif cosine_of_angle_of_incidence < BigDecimal('-1')
+      cosine_of_angle_of_incidence = BigDecimal('-1')
+    end
+
+    radian_to_degree(bigdecimal_acos(cosine_of_angle_of_incidence))
+  end  
+
+  private
+    # Rb 
+    # Reference:
+    #   3. Titled Surface Models
+    #   in [3]
+    def geometric_factor
+      hourly_beam_radiation_on_tilted_surface / hourly_beam_radiation_on_horizontal_surface
+    end
+
+    # I[b,T]
+    def hourly_beam_radiation_on_tilted_surface
+      hourly_bean_radiation = bigdecimal_cos(degree_to_radian(@surface_tilt)) * bigdecimal_cos(degree_to_radian(@sun_zenith)) + 
+                              bigdecimal_sin(degree_to_radian(@surface_tilt)) * bigdecimal_sin(degree_to_radian(@sun_zenith)) * bigdecimal_cos(degree_to_radian(@sun_azimuth - @surface_azimuth))
+      [hourly_bean_radiation, BigDecimal('0')].max
+    end
+
+    # Ib
+    def hourly_beam_radiation_on_horizontal_surface
+      [bigdecimal_cos(degree_to_radian(@sun_zenith)), BigDecimal('0.01745')].max
+    end
+
+    # AI
+    #   "defines the portion of the diffuse radiation to be treated as circumsolar with the remaining portion considered isotropic"[3]
+    def anisotropy_index
+      @direct_normal_irradiance / extraterrestrial_irradiance
+    end
+
+    def anisotropic_correction_factor
+      sin_cube_of_half_surface_tilt = (bigdecimal_sin(degree_to_radian(@surface_tilt/BigDecimal('2')))).power(3)
+      BigDecimal('1') + horizontal_brighting_correction_factor * sin_cube_of_half_surface_tilt
+    end
+
+    # f: Horizontal brighting correction factor
+    #   page 13 of [3]
+    def horizontal_brighting_correction_factor
+      adjusted_global_horizontal_irradiance = @global_horizontal_irradiance < SMALL_VALUE_FOR_SKY_DIFFUSE_IRRADIANCE ? SMALL_VALUE_FOR_SKY_DIFFUSE_IRRADIANCE : @global_horizontal_irradiance
+      # In PVLIB_MatLab as well as PVLIB_Python, there is a following adjustment. 
+      # However, all the negative numbers are already converted to SMALL_VALUE_FOR_SKY_DIFFUSE_IRRADIANCE in the above statement, 
+      # the following statement doesn't do anything. 
+      # I'm just writing it to remind that it is there in PVLIB_MatLab and PVLIB_Python. 
+      # adjusted_global_horizontal_irradiance = @global_horizontal_irradiance < BigDecimal('0') ? BigDecimal('0') : @global_horizontal_irradiance
+      bigdecimal_sqrt(horizontal_beam_irradiance / adjusted_global_horizontal_irradiance)
+    end
+
+    # HB
+    def horizontal_beam_irradiance
+      calculated_horizontal_beam_irradiance = @direct_normal_irradiance * bigdecimal_cos(degree_to_radian(@sun_zenith))
+      calculated_horizontal_beam_irradiance < BigDecimal('0') ? BigDecimal('0') : calculated_horizontal_beam_irradiance
+    end    
+
+end

data/lib/PVLIB_Ruby/pv_performance_characterization.rb

@@ -0,0 +1,111 @@
+require 'bigdecimal'
+require 'bigdecimal/math'
+
+# This class uses Sandia PV Array Performance Model, which is based on:
+# King, D. et al, 2004, "Sandia Photovoltaic Array Performance Model", SAND Report 3535, Sandia National Laboratories, Albuquerque, NM
+# Other models are realized by subclasses. 
+class PvPerformanceCharacterization
+
+  REFERECE_SOLAR_IRRADIANCE = BigDecimal('1000') # E0 [W/m^2]
+  REFERECE_CELL_TEMPERATURE = BigDecimal('25')   # T0 [ºC]
+  BOLTZMANN_CONSTANT = BigDecimal('1.38066E-23') # k [J/K]
+  ELEMENTARY_CHARGE = BigDecimal('1.60218E-19') # q [C(coulomb)] 
+
+  def initialize(pv_module, absolute_air_mass, angle_of_incidence, beam_irradiance, ground_irradiance, sky_diffuse_irradiance, estimated_cell_temperature)
+    @pv_module = pv_module
+    @effective_solar_irradiance = empirical_effective_solar_irradiance(absolute_air_mass, angle_of_incidence, beam_irradiance, ground_irradiance, sky_diffuse_irradiance)
+    @estimated_cell_temperature = estimated_cell_temperature
+  end
+
+  # Isc
+  def short_circuit_current
+    @pv_module.short_circuit_current_0 * @effective_solar_irradiance * (1 + @pv_module.normalized_temperature_coefficient_for_short_circuit_current * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE))
+  end
+
+  alias :isc :short_circuit_current
+
+  # Imp
+  def maximum_power_point_current
+    @pv_module.maximum_power_point_current_0 * (@pv_module.performance_coefficients[0] * @effective_solar_irradiance + @pv_module.performance_coefficients[1] * @effective_solar_irradiance.power(2)) * 
+    (1 + @pv_module.normalized_temperature_coefficient_for_maximum_power_point_current * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE))
+  end
+
+  alias :imp :maximum_power_point_current
+
+  # Voc
+  def open_circuit_voltage
+    @pv_module.open_circuit_voltage_0 + @pv_module.number_of_cells_in_series * thermal_voltage_per_cell * BigMath.log(@effective_solar_irradiance, 4) +
+    temperature_coefficient_for_open_circuit_voltage * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE)
+  end
+
+  alias :voc :open_circuit_voltage
+
+  # Vmp
+  def maximum_power_point_voltage
+    @pv_module.maximum_power_point_voltage_0 + @pv_module.performance_coefficients[2] * @pv_module.number_of_cells_in_series * thermal_voltage_per_cell * BigMath.log(@effective_solar_irradiance, 4) + 
+    @pv_module.performance_coefficients[3] * @pv_module.number_of_cells_in_series * (thermal_voltage_per_cell * BigMath.log(@effective_solar_irradiance, 4)).power(2) + 
+    temperature_coefficient_for_maximum_power_point_voltage * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE)
+  end
+
+  alias :vmp :maximum_power_point_voltage
+
+  # Ix
+  def fourth_point_current
+    @pv_module.fourth_point_current_0 * (@pv_module.performance_coefficients[4] * @effective_solar_irradiance + @pv_module.performance_coefficients[5] * @effective_solar_irradiance.power(2)) * 
+    (1 + @pv_module.normalized_temperature_coefficient_for_short_circuit_current * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE))
+  end
+
+  alias :ix :fourth_point_current
+
+  # Ixx
+  def fifth_point_current
+    @pv_module.fifth_point_current_0 * (@pv_module.performance_coefficients[6] * @effective_solar_irradiance + @pv_module.performance_coefficients[7] * @effective_solar_irradiance.power(2)) * 
+    (1 + @pv_module.normalized_temperature_coefficient_for_maximum_power_point_current * (@estimated_cell_temperature - REFERECE_CELL_TEMPERATURE))
+  end
+
+  alias :ixx :fifth_point_current
+
+  private
+    # Ee
+    def empirical_effective_solar_irradiance(absolute_air_mass, angle_of_incidence, beam_irradiance, ground_irradiance, sky_diffuse_irradiance)
+      diffuse_irradiance = ground_irradiance + sky_diffuse_irradiance
+
+      empirical_solar_spectral_influence_on_short_circuit_current(absolute_air_mass) * 
+      ((beam_irradiance * empirical_opitical_influence_on_short_circuit_current(angle_of_incidence) + @pv_module.diffuse_irradiance_factor * diffuse_irradiance)/REFERECE_SOLAR_IRRADIANCE)
+    end
+
+    # f1(AMa)
+    def empirical_solar_spectral_influence_on_short_circuit_current(absolute_air_mass)
+      @pv_module.air_mass_variation_polynomial_coefficients[4] +
+      @pv_module.air_mass_variation_polynomial_coefficients[3] * absolute_air_mass + 
+      @pv_module.air_mass_variation_polynomial_coefficients[2] * absolute_air_mass.power(2) +
+      @pv_module.air_mass_variation_polynomial_coefficients[1] * absolute_air_mass.power(3) +
+      @pv_module.air_mass_variation_polynomial_coefficients[0] * absolute_air_mass.power(4)
+    end
+
+    # f2(AOI)
+    def empirical_opitical_influence_on_short_circuit_current(angle_of_incidence)
+      @pv_module.angle_of_incidence_polynominal_coefficients[5] + 
+      @pv_module.angle_of_incidence_polynominal_coefficients[4] * angle_of_incidence + 
+      @pv_module.angle_of_incidence_polynominal_coefficients[3] * angle_of_incidence.power(2) + 
+      @pv_module.angle_of_incidence_polynominal_coefficients[2] * angle_of_incidence.power(3) + 
+      @pv_module.angle_of_incidence_polynominal_coefficients[1] * angle_of_incidence.power(4) + 
+      @pv_module.angle_of_incidence_polynominal_coefficients[0] * angle_of_incidence.power(5)
+    end
+
+    # delta(Tc)
+    def thermal_voltage_per_cell
+      @pv_module.diode_factor * BOLTZMANN_CONSTANT * (@estimated_cell_temperature + BigDecimal('273.15')) / ELEMENTARY_CHARGE
+    end
+
+    # BetaVoc
+    def temperature_coefficient_for_open_circuit_voltage
+      @pv_module.temperature_coefficient_for_open_circuit_voltage + @pv_module.irradiance_dependent_temperature_coefficient_for_open_circuit_voltage * (1-@effective_solar_irradiance)
+    end
+
+    # BetaVmp
+    def temperature_coefficient_for_maximum_power_point_voltage
+      @pv_module.temperature_coefficient_for_maximum_power_point_voltage + @pv_module.irradiance_dependent_temperature_coefficient_for_maximum_power_point_voltage * (1-@effective_solar_irradiance)
+    end
+
+end

data/lib/PVLIB_Ruby/pv_temperature.rb

@@ -0,0 +1,26 @@
+require 'bigdecimal'
+
+# This class calculates cell and module temperatures using Sandia PV Array Performance Model, which is based on:
+# King, D. et al, 2004, "Sandia Photovoltaic Array Performance Model", SAND Report 3535, Sandia National Laboratories, Albuquerque, NM
+class PvTemperature
+  include CalculationHelper
+
+  def initialize(pv_module, solar_irradiance_incident_on_module_surface, reference_solar_irradiance, wind_speed, air_temperature)
+    @pv_module = pv_module
+    @solar_irradiance_incident_on_module_surface = solar_irradiance_incident_on_module_surface
+    @reference_solar_irradiance = reference_solar_irradiance
+    @wind_speed = wind_speed
+    @air_temperature = air_temperature
+  end
+
+  def cell_temperature
+    module_temperature + 
+    (@solar_irradiance_incident_on_module_surface / @reference_solar_irradiance) * @pv_module.cell_module_back_surface_temperature_difference
+  end
+
+  def module_temperature
+    back_surface_module_temperature = @solar_irradiance_incident_on_module_surface * bigdecimal_exp(@pv_module.temperature_upper_limit_coefficient + @pv_module.wind_speed_decrease_rate_coefficient * @wind_speed) + 
+                                      @air_temperature
+  end
+
+end

data/lib/PVLIB_Ruby/solar_ephemeris.rb

@@ -0,0 +1,145 @@
+require 'bigdecimal'
+
+class SolarEphemeris
+  include CalculationHelper
+
+  def initialize(time, location, pressure: BigDecimal('101325'), temperature: BigDecimal('12'))
+    @time = time
+    @location = location
+    @pressure = pressure
+    @temperature = temperature
+
+    utc_offset_in_hours = BigDecimal(@time.utc_offset.to_s) / BigDecimal('3600')
+    dec_hours = BigDecimal(@time.hour.to_s) + BigDecimal(@time.min.to_s) / BigDecimal('60') + BigDecimal(@time.sec.to_s) / BigDecimal('3600')
+    @universal_date = BigDecimal(@time.yday.to_s) + ((dec_hours - utc_offset_in_hours) / BigDecimal('24')).floor
+    @universal_hour = (dec_hours - utc_offset_in_hours).modulo(BigDecimal('24')) 
+  end
+
+  def sun_azimuth
+    calculated_sun_azimuth = radian_to_degree(bigdecimal_atan2(BigDecimal('-1') * bigdecimal_sin(degree_to_radian(hour_angle)), 
+                                              bigdecimal_cos(degree_to_radian(@location.latitude)) * bigdecimal_tan(degree_to_radian(declination)) - bigdecimal_sin(degree_to_radian(@location.latitude)) * bigdecimal_cos(degree_to_radian(hour_angle))))
+    sun_azimuth = calculated_sun_azimuth < BigDecimal('0') ? calculated_sun_azimuth + BigDecimal('360') : calculated_sun_azimuth
+    sun_azimuth
+  end
+
+  def sun_elevation
+    # Caching because it's used several times for the calculation of apparent_sun_elevation
+    @sun_elevation ||= radian_to_degree(bigdecimal_asin(bigdecimal_cos(degree_to_radian(@location.latitude)) * bigdecimal_cos(degree_to_radian(declination)) * bigdecimal_cos(degree_to_radian(hour_angle)) + 
+                                                        bigdecimal_sin(degree_to_radian(@location.latitude)) * bigdecimal_sin(degree_to_radian(declination))))
+  end
+
+  def apparent_sun_elevation
+    sun_elevation + reflaction
+  end
+
+  def solar_time
+    # BigDecimal#sign returns 2 for positive value, -2 for negative value, 1 for 0 and -1 for -0. 
+    # Hence, decided to use the following statement: 
+    hour_angle_sign = hour_angle < BigDecimal('0') ? BigDecimal('-1') : BigDecimal('1')
+
+    adjusted_hour_angle = hour_angle.abs > BigDecimal('180') ? hour_angle - (BigDecimal('360') * hour_angle_sign) : hour_angle
+
+    (BigDecimal('180') + adjusted_hour_angle) / BigDecimal('15')
+  end
+
+  private
+    def hour_angle
+      loc_ast - rt_ascen
+    end
+
+    def declination
+      radian_to_degree(bigdecimal_asin(bigdecimal_sin(degree_to_radian(obliquity)) * bigdecimal_sin(degree_to_radian(ec_lon))))
+    end
+
+    # For hour_angle calculation
+    # LocAST in PVLIB_MatLab. Couldn't find what it stands for. 
+    # Also for other variables names, couldn't find what they stand for. 
+    def loc_ast
+      t = e_zero / 36525
+
+      gmst0 = BigDecimal('6') / BigDecimal('24') + BigDecimal('38') / BigDecimal('1440') + 
+              (BigDecimal('45.836') + BigDecimal('8640184.542') * t + BigDecimal('0.0929') * t.power(2) ) / BigDecimal('86400')
+      gmst0 = BigDecimal('360') * (gmst0 - gmst0.floor)
+      gmsti = (gmst0 + BigDecimal('360') * (BigDecimal('1.0027379093') * @universal_hour / BigDecimal('24'))).modulo(BigDecimal('360')) 
+
+      (BigDecimal('360') + gmsti + @location.longitude).modulo(BigDecimal('360')) 
+    end
+
+    # For hour_angle calculation
+    # RtAscen in PVLIB_MatLab. Couldn't find what it stands for. 
+    # Also for other variables names, couldn't find what they stand for. 
+    def rt_ascen
+      radian_to_degree(bigdecimal_atan2(bigdecimal_cos(degree_to_radian(obliquity)) * bigdecimal_sin(degree_to_radian(ec_lon)), bigdecimal_cos(degree_to_radian(ec_lon))))
+    end
+
+    def obliquity
+      BigDecimal('23.452294') - BigDecimal('0.0130125') * t1 - BigDecimal('0.00000164') * t1.power(2) + BigDecimal('0.000000503') * t1.power(3)
+    end
+
+    def ec_lon
+      (ml_perigee + true_anom).modulo(BigDecimal('360')) - abber
+    end
+
+    def ml_perigee
+      BigDecimal('281.22083') + BigDecimal('0.0000470684') * epoch_date + BigDecimal('0.000453') * t1.power(2) + BigDecimal('0.000003') * t1.power(3)
+    end
+
+    def true_anom
+      # Cannot find a proper name for these values:
+      temp1 = ((BigDecimal('1') + eccen) / (BigDecimal('1') - eccen)).power(0.5) * bigdecimal_tan(degree_to_radian(eccen_anom / BigDecimal('2')))
+      temp2 = radian_to_degree(bigdecimal_atan2(temp1, BigDecimal('1')))
+      BigDecimal('2') * temp2.modulo(BigDecimal('360'))
+    end
+
+    def abber
+      BigDecimal('20') / BigDecimal('3600')
+    end
+
+    def eccen
+      BigDecimal('0.01675104') - BigDecimal('0.0000418') * t1 - BigDecimal('0.000000126') * t1.power(2)
+    end
+
+    def eccen_anom
+      mean_anom = (BigDecimal('358.47583') + BigDecimal('0.985600267') * epoch_date - BigDecimal('0.00015') * t1.power(2) - BigDecimal('0.000003') * t1.power(3)).modulo(BigDecimal('360'))
+
+      eccen_anom = mean_anom
+      e = BigDecimal('0')
+      while (mean_anom - e).abs > BigDecimal('0.0001')
+        e = eccen_anom
+        eccen_anom = mean_anom + radian_to_degree(eccen) * bigdecimal_sin(degree_to_radian(e))
+      end
+      eccen_anom
+    end
+
+    def t1
+      epoch_date / BigDecimal('36525')
+    end
+
+    def epoch_date
+      e_zero + @universal_hour / BigDecimal('24')
+    end 
+
+    def e_zero
+      year_starting_from_1900 = BigDecimal(@time.year.to_s) - BigDecimal('1900')
+      year_begin = BigDecimal('365') * year_starting_from_1900 + ((year_starting_from_1900 - BigDecimal('1')) / BigDecimal('4')).floor - BigDecimal('0.5')
+      year_begin + @universal_date
+    end
+
+    # For apparent_sun_elevation calculation
+    def reflaction
+      tangent_of_elevation = bigdecimal_tan(degree_to_radian(sun_elevation))
+
+      raw_reflaction = if sun_elevation > BigDecimal('5') && sun_elevation <= BigDecimal('85')
+                         BigDecimal('58.1') / tangent_of_elevation - BigDecimal('0.07') / (tangent_of_elevation.power(3)) + BigDecimal('.000086') / (tangent_of_elevation.power(5))
+                       elsif sun_elevation > BigDecimal('-0.575') && sun_elevation <= BigDecimal('5')
+                         sun_elevation * (BigDecimal('-518.2') + sun_elevation * (BigDecimal('103.4') + sun_elevation * (BigDecimal('-12.79') + sun_elevation * BigDecimal('0.711')))) + BigDecimal('1735')
+                       elsif sun_elevation > BigDecimal('-1') && sun_elevation <= BigDecimal('-0.575')
+                         BigDecimal('-20.774') / tangent_of_elevation
+                       else
+                         BigDecimal('0')
+                       end
+ 
+      raw_reflaction * (BigDecimal('283') / (BigDecimal('273') + @temperature)) * @pressure / BigDecimal('101325') / BigDecimal('3600')
+    end
+
+end

data/lib/PVLIB_Ruby/version.rb

@@ -0,0 +1,3 @@
+module PvlibRuby
+  VERSION = "0.1.0"
+end

data/lib/PVLIB_Ruby.rb

@@ -0,0 +1,17 @@
+require 'pry'
+require 'pvlib_ruby/version'
+require 'pvlib_ruby/models/active_csv'
+require 'pvlib_ruby/models/pv_module'
+require 'pvlib_ruby/models/inverter'
+require 'pvlib_ruby/pv_performance_characterization'
+require 'pvlib_ruby/dc_to_ac_conversion'
+require 'pvlib_ruby/calculation_helper'
+require 'pvlib_ruby/air_mass'
+require 'pvlib_ruby/pv_temperature'
+require 'pvlib_ruby/plain_of_array_irradiance'
+require 'pvlib_ruby/models/location'
+require 'pvlib_ruby/solar_ephemeris'
+
+module PvlibRuby
+  # Your code goes here...
+end

metadata

@@ -0,0 +1,130 @@
+--- !ruby/object:Gem::Specification
+name: PVLIB_Ruby
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- Tadatoshi Takahashi
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2015-06-17 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: activemodel
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.9'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.9'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: PVLIB_MatLab translated to Ruby code in Object-Oriented way.
+email:
+- tadatoshi@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".ruby-gemset"
+- ".ruby-version"
+- ".travis.yml"
+- CODE_OF_CONDUCT.md
+- Gemfile
+- LICENSE.txt
+- PVLIB_Ruby.gemspec
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- examples/data/sandia_inverter_example.csv
+- examples/data/sandia_module_example.csv
+- examples/pv_power_generation_example.rb
+- lib/PVLIB_Ruby.rb
+- lib/PVLIB_Ruby/air_mass.rb
+- lib/PVLIB_Ruby/calculation_helper.rb
+- lib/PVLIB_Ruby/dc_to_ac_conversion.rb
+- lib/PVLIB_Ruby/models/active_csv.rb
+- lib/PVLIB_Ruby/models/inverter.rb
+- lib/PVLIB_Ruby/models/location.rb
+- lib/PVLIB_Ruby/models/pv_module.rb
+- lib/PVLIB_Ruby/plain_of_array_irradiance.rb
+- lib/PVLIB_Ruby/pv_performance_characterization.rb
+- lib/PVLIB_Ruby/pv_temperature.rb
+- lib/PVLIB_Ruby/solar_ephemeris.rb
+- lib/PVLIB_Ruby/version.rb
+homepage: https://github.com/tadatoshi/PVLIB_Ruby
+licenses:
+- MIT
+metadata:
+  allowed_push_host: https://rubygems.org
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.4.7
+signing_key: 
+specification_version: 4
+summary: PVLIB written in Ruby.
+test_files: []