solar_geometry_calculation 0.1.4

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 900bb37f7b3f5cdbe8dc83bd45585cc66ff45f02
+  data.tar.gz: d1704fa9695996338eb4d0d01f95f8deb47dbcba
+SHA512:
+  metadata.gz: 1dc72c21d597a19454cbc2a83197deda1eb428775dfc456b078bd2d138e8f206239fa218298fc0c2e78d6a4b9fe5d41fabb9ea45599411d681e690e40d8dd963
+  data.tar.gz: 9a2faaef30525a75e90696d0c7b9675a75cec9c5ad411912d646f8967a3d5d0500409e0ddafcdcb2f4bebaaec968659fb77bc1ea58c98015141c97fe45720d64

data/.gitignore

@@ -0,0 +1,4 @@
+*.gem
+.bundle
+Gemfile.lock
+pkg/*

data/.rspec

@@ -0,0 +1 @@
+--color

data/.ruby-gemset

@@ -0,0 +1 @@
+solar_geometry_calculation

data/.ruby-version

@@ -0,0 +1 @@
+2.1.2

data/Gemfile

@@ -0,0 +1,4 @@
+source "http://rubygems.org"
+
+# Specify your gem's dependencies in solar_geometry_calculation.gemspec
+gemspec

data/LICENSE

@@ -0,0 +1,20 @@
+Copyright (c) 2010-2014 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/README.md

@@ -0,0 +1,30 @@
+# Solar Geometry Calculation
+
+Calculate solar geometry, such as height or azimuth of the Sun for the purpose of obtaining information for solar panel installation. 
+
+
+## Installation
+
+Add to your Gemfile and run the `bundle` command to install it.
+
+```ruby
+gem "solar_geometry_calculation"
+```
+
+**Requires Ruby 2.0.0 or later.**
+
+
+## Usage
+
+Instantiate SolarPositionCalculation class.
+
+```ruby
+solar_position_calculation = SolarPositionCalculation.new(latitude: "45.5")
+```
+
+
+## Development
+
+Questions or problems? Please post them on the [issue tracker](https://github.com/tadatoshi/solar_geometry_calculation/issues). You can contribute changes by forking the project and submitting a pull request. You can ensure the tests passing by running `bundle` and `rake`.
+
+This gem is created by Tadatoshi Takahashi and is under the MIT License.

data/Rakefile

@@ -0,0 +1,8 @@
+require 'bundler'
+Bundler::GemHelper.install_tasks
+
+require 'rspec/core/rake_task'
+
+RSpec::Core::RakeTask.new(:spec)
+
+task default: :spec

data/lib/solar_geometry_calculation.rb

@@ -0,0 +1,9 @@
+require 'bigdecimal'
+
+module SolarGeometryCalculation
+
+  autoload :SolarPositionCalculation, 'solar_geometry_calculation/solar_position_calculation'
+  autoload :SolarTime, 'solar_geometry_calculation/solar_time'
+  autoload :CalculationHelper, 'solar_geometry_calculation/calculation_helper'
+
+end

data/lib/solar_geometry_calculation/calculation_helper.rb

@@ -0,0 +1,45 @@
+module SolarGeometryCalculation
+
+  module CalculationHelper
+  
+    def day_of_the_year(year, month, day, hour = 0, minute = 0)
+      if hour == 0 && minute == 0
+        BigDecimal((Date.new(year.to_i, month.to_i, day.to_i).yday).to_s)
+      else
+        day_of_the_year_one_day_before = BigDecimal((Date.new(year.to_i, month.to_i, day.to_i).yday).to_s) - 1
+        decimal_hour_minute = BigDecimal(hour.to_s) + BigDecimal(minute.to_s) / 60
+        (day_of_the_year_one_day_before * 24 + decimal_hour_minute + 12) / 24
+      end
+    end 
+  
+    def degree_to_radian(degree)
+      degree * BigDecimal(Math::PI.to_s) / 180
+    end
+
+    def radian_to_degree(radian)
+      radian * 180 / BigDecimal(Math::PI.to_s)
+    end  
+    
+    def big_decimal_sin(angle_in_radian)
+      BigDecimal(Math.sin(angle_in_radian).to_s)
+    end
+
+    def big_decimal_cos(angle_in_radian)
+      BigDecimal(Math.cos(angle_in_radian).to_s)
+    end 
+
+    def big_decimal_tan(angle_in_radian)
+      BigDecimal(Math.tan(angle_in_radian).to_s)
+    end   
+
+    def big_decimal_asin(angle_in_radian)
+      BigDecimal(Math.asin(angle_in_radian).to_s)
+    end
+
+    def big_decimal_acos(angle_in_radian)
+      BigDecimal(Math.acos(angle_in_radian).to_s)
+    end    
+    
+  end
+
+end

data/lib/solar_geometry_calculation/solar_position_calculation.rb

@@ -0,0 +1,116 @@
+require 'active_support/core_ext/object/blank'
+
+module SolarGeometryCalculation
+
+  class SolarPositionCalculation
+  	include CalculationHelper
+
+    attr_accessor :id, :latitude, :meridian, :longitude, :year, :month, :day, :hour, :minute, :surface_inclination, :surface_azimuth, :solar_time
+    attr_writer :timezone_identifier	
+
+    def initialize(attributes = {})
+      attributes.each { |attribute_name, attribute_value| self.send("#{attribute_name}=".to_sym, attribute_value) }
+      get_solar_time
+    end
+
+
+    def declination(with_hour_and_minutes: true)
+      day_of_the_year = with_hour_and_minutes ? day_of_the_year(@solar_time.year, @solar_time.month, @solar_time.day, @solar_time.hour, @solar_time.minute) : day_of_the_year(@solar_time.year, @solar_time.month, @solar_time.day)
+      radian = degree_to_radian(BigDecimal("360") / BigDecimal("365") * (day_of_the_year + 284))
+      (23.45 * big_decimal_sin(radian)).round(2)
+    end
+
+    def hour_angle
+      if @solar_time.blank?
+        hour_for_calculation = BigDecimal(@hour.to_s)
+        minute_for_calculation = BigDecimal(@minute.to_s)
+      else
+        hour_for_calculation = @solar_time.hour
+        minute_for_calculation = @solar_time.minute
+      end
+   
+      hour_with_minutes_in_digit = hour_for_calculation + minute_for_calculation.div(60, 2)
+    
+      ((hour_with_minutes_in_digit - 12) * 15).round(2)
+    end
+
+    def solar_elevation
+      latitude_in_radian = degree_to_radian(BigDecimal(self.latitude.to_s))
+      declination_in_radian = degree_to_radian(self.declination)
+      hour_angle_in_radian = degree_to_radian(self.hour_angle)
+    
+      sine_solar_elevation = big_decimal_sin(latitude_in_radian) * big_decimal_sin(declination_in_radian) + 
+                             big_decimal_cos(latitude_in_radian) * big_decimal_cos(declination_in_radian) * big_decimal_cos(hour_angle_in_radian)
+      radian_to_degree(big_decimal_asin(sine_solar_elevation)).round(2)
+    end  
+
+    def solar_azimuth
+      declination_in_radian = degree_to_radian(self.declination)
+      hour_angle_in_radian = degree_to_radian(self.hour_angle)
+      solar_elevation_in_radian = degree_to_radian(self.solar_elevation)
+    
+      sine_solar_azimuth = big_decimal_cos(declination_in_radian) * big_decimal_sin(hour_angle_in_radian) / big_decimal_cos(solar_elevation_in_radian)
+      radian_to_degree(big_decimal_asin(sine_solar_azimuth)).round(2)
+    end  
+
+    def sunrise
+      sunrise_time_with_munutes_in_digit = BigDecimal((12 - (sunrise_sunset_base_angle / 15)).to_s)    
+      sunrise_sunset_time(sunrise_time_with_munutes_in_digit)
+    end
+  
+    def sunset
+      sunset_time_with_munutes_in_digit = BigDecimal((12 + (sunrise_sunset_base_angle / 15)).to_s)    
+      sunrise_sunset_time(sunset_time_with_munutes_in_digit)
+    end  
+
+    def angle_of_incidence
+      latitude_in_radian = degree_to_radian(BigDecimal(self.latitude.to_s))
+      declination_in_radian = degree_to_radian(self.declination)
+      hour_angle_in_radian = degree_to_radian(self.hour_angle)
+      surface_inclination_in_radian = degree_to_radian(BigDecimal(@surface_inclination.to_s))
+      surface_azimuth_in_radian = degree_to_radian(BigDecimal(@surface_azimuth.to_s))
+
+      first_value = big_decimal_sin(declination_in_radian) * big_decimal_sin(latitude_in_radian) * big_decimal_cos(surface_inclination_in_radian)
+      second_value = big_decimal_sin(declination_in_radian) * big_decimal_cos(latitude_in_radian) * big_decimal_sin(surface_inclination_in_radian) * big_decimal_cos(surface_azimuth_in_radian)
+      third_value = big_decimal_cos(declination_in_radian) * big_decimal_cos(latitude_in_radian) * big_decimal_cos(surface_inclination_in_radian) * big_decimal_cos(hour_angle_in_radian)
+      fourth_value = big_decimal_cos(declination_in_radian) * big_decimal_sin(latitude_in_radian) * big_decimal_sin(surface_inclination_in_radian) * big_decimal_cos(surface_azimuth_in_radian) * big_decimal_cos(hour_angle_in_radian)
+      fifth_value = big_decimal_cos(declination_in_radian) * big_decimal_sin(surface_inclination_in_radian) * big_decimal_sin(surface_azimuth_in_radian) * big_decimal_sin(hour_angle_in_radian)
+      cosine_of_angle_of_incidence = first_value - second_value + third_value + fourth_value + fifth_value
+    
+      radian_to_degree(big_decimal_acos(cosine_of_angle_of_incidence)).round(2)
+    end          
+
+    private 
+      def get_solar_time
+        if @year.present? && @month.present? && @day.present?
+          if @address.present? && @place.present?
+            @solar_time = SolarTime.create_by_standard_time(BigDecimal(@year.to_s), BigDecimal(@month.to_s), BigDecimal(@day.to_s), BigDecimal(@hour.to_s), BigDecimal(@minute.to_s), @place.meridian, @place.longitude, @place.timezone_identifier) 
+          elsif @meridian.present? && @longitude.present? && @timezone_identifier.present?
+            @solar_time = SolarTime.create_by_standard_time(BigDecimal(@year.to_s), BigDecimal(@month.to_s), BigDecimal(@day.to_s), BigDecimal(@hour.to_s), BigDecimal(@minute.to_s), @meridian, @longitude, @timezone_identifier)           
+          end
+        end
+      end 
+
+      def sunrise_sunset_base_angle
+        latitude_in_radian = degree_to_radian(BigDecimal(self.latitude.to_s))
+        declination_in_radian = degree_to_radian(self.declination(with_hour_and_minutes: false))
+
+        sunrise_sunset_base_angle_in_radian = big_decimal_acos(-(big_decimal_tan(latitude_in_radian) * big_decimal_tan(declination_in_radian)))
+        radian_to_degree(sunrise_sunset_base_angle_in_radian)
+      end  
+
+      def sunrise_sunset_time(sunrise_or_sunset_time_with_munute_in_digit)
+        sunrise_or_sunset_hour, sunrise_or_sunset_minute_in_10_base = sunrise_or_sunset_time_with_munute_in_digit.divmod(1)   
+        sunrise_or_sunset_solar_time = SolarTime.create_by_solar_time_values(@solar_time.year, @solar_time.month, @solar_time.day, sunrise_or_sunset_hour, (sunrise_or_sunset_minute_in_10_base * 60).round, @meridian, @longitude, @solar_time.timezone_identifier)
+        "#{sunrise_or_sunset_solar_time.standard_time_hour.to_i}:#{self.class.format_minute(sunrise_or_sunset_solar_time.standard_time_minute.to_i)}"      
+      end  
+
+    class << self
+      def format_minute(minute)
+        minute.to_s.length == 1 ? "0#{minute.to_s}" : minute.to_s
+      end
+    end         
+
+  end
+
+end

data/lib/solar_geometry_calculation/solar_time.rb

@@ -0,0 +1,206 @@
+require 'tzinfo'
+require 'active_support/core_ext/time/calculations'
+require 'active_support/core_ext/date/calculations'
+
+module SolarGeometryCalculation
+
+  class SolarTime
+  
+    attr_accessor :year, :month, :day, :hour, :minute, :meridian, :longitude, :timezone_identifier
+    attr_writer :standard_time
+  
+    # Methods for standard_time provided because the method to get minute for Time is min and that naming is not consistent with minute method for this class: 
+    def standard_time_year
+      @standard_time.year
+    end
+  
+    def standard_time_month
+      @standard_time.month
+    end
+  
+    def standard_time_day
+      @standard_time.day
+    end
+  
+    def standard_time_hour
+      @standard_time.hour
+    end
+  
+    def standard_time_minute
+      @standard_time.min
+    end
+  
+    class << self
+    
+      # Moved from app/models/factories/solar_time_factory.rb because "rails server" in development environment cannot find this module. 
+      # -- At least after I tried rack-mini-profiler. Even after I commented it out in Gemfile, the problem still persists.
+      def create_by_standard_time(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+        solar_time = SolarTime.new
+
+        # Uses Time because DateTime doesn't take daylight savings time even with the one from Rails:
+        solar_time.standard_time = Time.new(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute)
+
+        solar_time.meridian = meridian
+        solar_time.longitude = longitude
+        solar_time.timezone_identifier = timezone_identifier 
+
+        standard_time_year_without_offset, standard_time_month_without_offset, standard_time_day_without_offset, standard_time_hour_without_offset, standard_time_minute_without_offset = 
+          SolarTime.standard_time_without_daylight_savings_time_offset(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)    
+
+        solar_time.year, solar_time.month, solar_time.day, solar_time.hour, solar_time.minute = 
+          SolarTime.standard_time_without_daylight_savings_time_offset_to_solar_time(standard_time_year_without_offset, standard_time_month_without_offset, standard_time_day_without_offset, standard_time_hour_without_offset, standard_time_minute_without_offset, meridian, longitude, timezone_identifier)
+
+        solar_time
+      end
+    
+      # Moved from app/models/factories/solar_time_factory.rb because "rails server" in development environment cannot find this module. 
+      # -- At least after I tried rack-mini-profiler. Even after I commented it out in Gemfile, the problem still persists.
+      def create_by_solar_time_values(solar_time_year, solar_time_month, solar_time_day, solar_time_hour, solar_time_minute, meridian, longitude, timezone_identifier)
+        solar_time = SolarTime.new
+        solar_time.year = solar_time_year
+        solar_time.month = solar_time_month
+        solar_time.day = solar_time_day
+        solar_time.hour = solar_time_hour
+        solar_time.minute = solar_time_minute
+        solar_time.meridian = meridian
+        solar_time.longitude = longitude    
+        solar_time.timezone_identifier = timezone_identifier    
+
+        standard_time_year_without_offset, standard_time_month_without_offset, standard_time_day_without_offset, standard_time_hour_without_offset, standard_time_minute_without_offset = 
+          SolarTime.solar_time_to_standard_time_without_daylight_savings_time_offset(solar_time.year, solar_time.month, solar_time.day, solar_time.hour, solar_time.minute, meridian, longitude, timezone_identifier)
+
+        standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute = 
+          SolarTime.standard_time_with_daylight_savings_time_offset(standard_time_year_without_offset, standard_time_month_without_offset, standard_time_day_without_offset, standard_time_hour_without_offset, standard_time_minute_without_offset, meridian, longitude, timezone_identifier)
+
+        # Uses Time because DateTime doesn't take daylight savings time even with the one from Rails:
+        solar_time.standard_time = Time.new(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute)
+
+        solar_time
+      end    
+    
+      def standard_time_without_daylight_savings_time_offset(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+        NegativeDaylightSavingsTimeConverter.new.convert(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+      end
+    
+      def standard_time_with_daylight_savings_time_offset(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+        PositiveDaylightSavingsTimeConverter.new.convert(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+      end    
+
+      def standard_time_without_daylight_savings_time_offset_to_solar_time(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+        ToSolarTimeConverter.new.convert(standard_time_year, standard_time_month, standard_time_day, standard_time_hour, standard_time_minute, meridian, longitude, timezone_identifier)
+      end   
+    
+      def solar_time_to_standard_time_without_daylight_savings_time_offset(solar_time_year, solar_time_month, solar_time_day, solar_time_hour, solar_time_minute, meridian, longitude, timezone_identifier)
+        FromSolarTimeConverter.new.convert(solar_time_year, solar_time_month, solar_time_day, solar_time_hour, solar_time_minute, meridian, longitude, timezone_identifier)
+      end    
+
+    end
+  
+  end
+
+  module TimeConverter
+  
+    def convert(year, month, day, hour, minute, meridian, longitude, timezone_identifier)
+
+      # Creating Time object as UTC since Time.local gets Time object in machine's timezone. 
+      time = Time.utc(year, month, day, hour, minute)
+
+      time = convert_by_offset(time, meridian, longitude, timezone_identifier)
+
+      [BigDecimal(time.year.to_s), BigDecimal(time.month.to_s), BigDecimal(time.day.to_s), BigDecimal(time.hour.to_s), BigDecimal(time.min.to_s)]
+
+    end
+  
+    private
+      def convert_by_offset(time, meridian, longitude, timezone_identifier)
+        raise "must be implemented by a class"
+      end
+  
+      def target_time(time, offset_in_seconds)
+        raise "must be implemented by a class"
+      end
+    
+  end  
+
+  module DaylightSavingsTimeConverter 
+    include TimeConverter
+  
+    private
+      def convert_by_offset(time, meridian, longitude, timezone_identifier)
+        timezone = TZInfo::Timezone.get(timezone_identifier)
+        timezone_period = timezone.period_for_local(time, true)
+
+        if timezone_period.dst?
+          time = target_time(time, timezone_period.std_offset)
+        end
+      
+        time
+      end    
+  
+      def target_time(time, offset_in_seconds)
+        raise "must be implemented by a class"
+      end
+  
+  end
+
+  module SolarTimeConverter 
+    include TimeConverter
+    include SolarGeometryCalculation::CalculationHelper
+  
+    private
+      def convert_by_offset(time, meridian, longitude, timezone_identifier)
+        minutes_difference = minutes_difference(time.year, time.month, time.day, meridian, longitude)
+        time = target_time(time, minutes_difference * 60)
+      end    
+
+      def target_time(time, offset_in_seconds)
+        raise "must be implemented by a class"
+      end
+  
+      def minutes_difference(year, month, day, meridian, longitude)
+        -(4 * (BigDecimal(meridian.to_s) - BigDecimal(longitude.to_s))) + equation_of_time(year, month, day).round
+      end    
+
+      def equation_of_time(year, month, day)
+        b_in_radian = degree_to_radian((day_of_the_year(year, month, day) - 1) * BigDecimal("360") / BigDecimal("365"))
+  
+        BigDecimal("229.2") * (BigDecimal("0.000075") + BigDecimal("0.001868") * big_decimal_cos(b_in_radian) - BigDecimal("0.032077") * big_decimal_sin(b_in_radian) - BigDecimal("0.014615") * big_decimal_cos(2 * b_in_radian) - BigDecimal("0.04089") * big_decimal_sin(2 * b_in_radian))                  
+      end
+  
+  end
+
+  module OffsetAdvanceConverter
+    private
+      def target_time(time, offset_in_seconds)
+        time.advance(seconds: offset_in_seconds)
+      end    
+  end
+
+  module OffsetSubtractionConverter
+    private
+      def target_time(time, offset_in_seconds)
+        time.ago(offset_in_seconds)
+      end    
+  end  
+
+  class NegativeDaylightSavingsTimeConverter
+    include DaylightSavingsTimeConverter
+    include OffsetSubtractionConverter
+  end
+
+  class PositiveDaylightSavingsTimeConverter
+    include DaylightSavingsTimeConverter
+    include OffsetAdvanceConverter 
+  end  
+
+  class ToSolarTimeConverter
+    include SolarTimeConverter
+    include OffsetAdvanceConverter
+  end
+
+  class FromSolarTimeConverter
+    include SolarTimeConverter
+    include OffsetSubtractionConverter
+  end  
+
+end

data/lib/solar_geometry_calculation/version.rb

@@ -0,0 +1,3 @@
+module SolarGeometryCalculation
+  VERSION = "0.1.4"
+end

data/solar_geometry_calculation.gemspec

@@ -0,0 +1,26 @@
+# -*- encoding: utf-8 -*-
+$:.push File.expand_path("../lib", __FILE__)
+require "solar_geometry_calculation/version"
+
+Gem::Specification.new do |s|
+  s.name        = "solar_geometry_calculation"
+  s.version     = SolarGeometryCalculation::VERSION
+  s.platform    = Gem::Platform::RUBY
+  s.authors     = ["Tadatoshi Takahashi"]
+  s.email       = ["tadatoshi@gmail.com"]
+  s.homepage    = "https://github.com/tadatoshi/solar_geometry_calculation"
+  s.summary     = %q{Performs calculation for solar geometry}
+  s.description = %q{Mathematical part of obtaining solar geometry information}
+
+  s.rubyforge_project = "solar_geometry_calculation"
+
+  s.files         = `git ls-files`.split("\n")
+  s.test_files    = `git ls-files -- {test,spec,features}/*`.split("\n")
+  s.executables   = `git ls-files -- bin/*`.split("\n").map{ |f| File.basename(f) }
+  s.require_paths = ["lib"]
+
+  s.add_dependency 'tzinfo'
+  s.add_dependency 'activesupport', ">= 4.1.1"
+
+  s.add_development_dependency "rspec"
+end

data/spec/solar_geometry_calculation/calculation_helper_spec.rb

@@ -0,0 +1,40 @@
+require 'spec_helper'
+
+describe SolarGeometryCalculation::CalculationHelper do
+  
+  class CalculationHelperExtender; extend SolarGeometryCalculation::CalculationHelper; end
+  
+  context "Day of the year" do
+  
+    it "should calculate the day of the year" do
+    
+      CalculationHelperExtender.day_of_the_year(2012, 2, 19).should == 50
+    
+    end
+    
+    it "should calculate the day of the year taking into account of the time in the day" do
+      
+      CalculationHelperExtender.day_of_the_year(2012, 2, 19, 6, 50).should be_within(0.01).of(49.78)
+      CalculationHelperExtender.day_of_the_year(2012, 2, 19, 17, 10).should be_within(0.01).of(50.22)
+      
+    end
+  
+  end
+  
+  context "Radian degree" do
+    
+    it "should convert a degree to Radian because Math.sin expects Radian" do
+      
+      CalculationHelperExtender.degree_to_radian(1).should be_within(0.0001).of(0.0175)
+      
+    end
+    
+    it "should convert radian to degree" do
+      
+      CalculationHelperExtender.radian_to_degree(1).should be_within(0.0001).of(57.2958)
+      
+    end    
+    
+  end  
+  
+end

data/spec/solar_geometry_calculation/solar_position_calculation_spec.rb

@@ -0,0 +1,176 @@
+require 'spec_helper'
+
+describe SolarGeometryCalculation::SolarPositionCalculation do
+
+  context "declination" do
+    
+    it "should calculate declination" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.declination.should == -0.19
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.declination.should == -20.14
+      
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_3.declination.should == 19.83        
+      
+    end
+    
+    context "(solar energy and applications course - homework 2)" do
+      
+      it "should calculate declination with latitude 45 degrees 28 minutes" do
+
+        solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+        solar_position_calculation_1.declination.should be_within(0.1).of(-11.95)
+                
+        # At sunrise 6:50 (6.83h):
+        solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 6, :minute => 50, :surface_inclination => "0", :surface_azimuth => "0")
+        solar_position_calculation_2.declination.should be_within(0.1).of(-12.02)
+
+        # At sunset 17:10 (17.17h):
+        solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 2, :day => 19, :hour => 17, :minute => 10, :surface_inclination => "0", :surface_azimuth => "0")
+        solar_position_calculation_3.declination.should be_within(0.1).of(-11.87)     
+
+      end      
+      
+    end
+    
+  end
+
+  context "hour angle" do
+    
+    it "should calculate hour angle" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.hour_angle.should be_within(0.1).of(-42.03)
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.hour_angle.should be_within(0.15).of(-1.40)
+
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 13, :minute => 15, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_3.hour_angle.should be_within(0.1).of(3.37)           
+      
+    end
+    
+  end
+
+  context "height of the sun" do
+    
+    it "should calculate solar elevation" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.solar_elevation.should be_within(0.1).of(43.09) 
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.solar_elevation.should == 24.35
+      
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.5", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_3.solar_elevation.should be_within(0.1).of(61.37)             
+      
+    end
+    
+  end
+  
+  context "Solar azimuth" do
+    
+    it "should calculate solar azimuth" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.solar_azimuth.should be_within(0.1).of(-51.58)
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.solar_azimuth.should be_within(0.1).of(-31.03)
+      
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_3.solar_azimuth.should be_within(0.1).of(-70.83)   
+
+    end
+    
+  end  
+  
+  context "Sunrise" do
+    
+    it "should calculate sunrise" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.sunrise.should == "6:49"
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.sunrise.should == "7:34"
+      
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_3.sunrise.should == "5:36"      
+      
+    end
+
+    it "should calculate the same sunrise time throughout the same day" do
+
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 3, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.sunrise.should == "6:49"
+
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.46", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 21, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.sunrise.should == "6:49" 
+
+      solar_position_calculation_1.sunrise.should == solar_position_calculation_2.sunrise    
+
+    end
+    
+    context "(solar energy and applications course - homework 2)" do
+      
+      it "should calculate sunrise with latitude 45 degrees 28 minutes" do
+
+        solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+        solar_position_calculation_1.sunrise.should == "6:59"        
+
+      end      
+      
+    end        
+    
+  end
+
+  context "Sunset" do
+
+    it "should calculate the same sunset time throughout the same day" do
+
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 2, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.sunset.should == "17:19"
+
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 22, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_2.sunset.should == "17:19" 
+
+      solar_position_calculation_1.sunset.should == solar_position_calculation_2.sunset    
+
+    end
+    
+    context "(solar energy and applications course - homework 2)" do
+      
+      it "should calculate sunset with latitude 45 degrees 28 minutes" do
+
+        solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2012, :month => 2, :day => 19, :hour => 12, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+        solar_position_calculation_1.sunset.should == "17:19"        
+
+      end      
+      
+    end
+    
+  end  
+  
+  context "Angle of incidence" do
+    
+    it "should calculate angle of incidence" do
+      
+      solar_position_calculation_1 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 9, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "0", :surface_azimuth => "0")
+      solar_position_calculation_1.angle_of_incidence.should be_within(0.1).of(58.73)
+      
+      solar_position_calculation_2 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 1, :day => 21, :hour => 10, :minute => 0, :surface_inclination => "70", :surface_azimuth => "-45")
+      solar_position_calculation_2.angle_of_incidence.should be_within(0.1).of(13.28)
+      
+      solar_position_calculation_3 = SolarGeometryCalculation::SolarPositionCalculation.new(:latitude => "45.47", :meridian => "-75", :longitude => "-73.75", :timezone_identifier => "America/Montreal", :year => 2010, :month => 7, :day => 24, :hour => 10, :minute => 0, :surface_inclination => "70", :surface_azimuth => "-45")
+      solar_position_calculation_3.angle_of_incidence.should == 32.78     
+      
+    end
+    
+  end
+
+end

data/spec/solar_geometry_calculation/solar_time_spec.rb

@@ -0,0 +1,118 @@
+require 'spec_helper'
+
+describe SolarGeometryCalculation::SolarTime do
+  
+  context "Factory" do
+    
+    it "should create an instance based on the standard time" do
+      
+      solar_time = SolarGeometryCalculation::SolarTime.create_by_standard_time(2010, 9, 21, 12, 0, -75, -73.75, "America/Los_Angeles")
+      
+      solar_time.standard_time_year.should == 2010
+      solar_time.standard_time_month.should == 9
+      solar_time.standard_time_day.should == 21
+      solar_time.standard_time_hour.should == 12
+      solar_time.standard_time_minute.should == 0
+      
+      solar_time.year.should == 2010
+      solar_time.month.should == 9
+      solar_time.day.should == 21
+      solar_time.hour.should == 11
+      solar_time.minute.should == 12 
+      
+    end
+    
+    it "should create an instance based on the solar time values" do
+      
+      solar_time = SolarGeometryCalculation::SolarTime.create_by_solar_time_values(2010, 9, 21, 11, 12, -75, -73.75, "America/Los_Angeles")
+      
+      solar_time.year.should == 2010
+      solar_time.month.should == 9
+      solar_time.day.should == 21
+      solar_time.hour.should == 11
+      solar_time.minute.should == 12 
+      
+      solar_time.standard_time_year.should == 2010
+      solar_time.standard_time_month.should == 9
+      solar_time.standard_time_day.should == 21
+      solar_time.standard_time_hour.should == 12
+      solar_time.standard_time_minute.should == 0 
+      
+    end    
+  
+  end  
+
+  context "helper" do
+  
+    it "should get the time without daylight saving time offset from the time with daylight saving time offset" do
+    
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset, 2012, 07, 13, 10, 12, -75, -78.63, "America/Los_Angeles").should == [2012, 7, 13, 9, 12]
+      
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset, 2012, 8, 1, 0, 30, -75, -78.63, "America/Los_Angeles").should == [2012, 7, 31, 23, 30]      
+    
+    end
+    
+    it "should get the time with daylight saving time offset from the time without daylight saving time offset" do
+    
+      SolarGeometryCalculation::SolarTime.send(:standard_time_with_daylight_savings_time_offset, 2012, 07, 13, 9, 12, -75, -78.63, "America/Los_Angeles").should == [2012, 7, 13, 10, 12]
+      
+      SolarGeometryCalculation::SolarTime.send(:standard_time_with_daylight_savings_time_offset, 2012, 7, 31, 23, 30, -75, -78.63, "America/Los_Angeles").should == [2012, 8, 1, 0, 30]      
+    
+    end    
+
+    it "should calculate solar time" do
+
+      # Data for Toronto:
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 2, 19, 7, 14, -75, -78.63, "America/Los_Angeles").should == [2010, 2, 19, 6, 45]        
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 9, 21, 10, 0, -75, -73.55, "America/Los_Angeles").should == [2010, 9, 21, 10, 12]
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 9, 21, 10, 0, BigDecimal("-75"), BigDecimal("-73.55"), "America/Los_Angeles").should == [2010, 9, 21, 10, 12]
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 9, 21, 10, 0, -75, -78.30, "America/Los_Angeles").should == [2010, 9, 21, 9, 53]
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 8, 1, 0, 0, -75, -78.30, "America/Los_Angeles").should == [2010, 7, 31, 23, 39]    
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 7, 31, 23, 59, -75, -78.30, "America/Los_Angeles").should == [2010, 7, 31, 23, 38]    
+
+      SolarGeometryCalculation::SolarTime.send(:standard_time_without_daylight_savings_time_offset_to_solar_time, 2010, 9, 21, 12, 0, -75, -73.75, "America/Los_Angeles").should == [2010, 9, 21, 12, 12]    
+
+    end    
+    
+  end
+
+  context "(solar energy and applications course - homework 2)" do
+    
+    it "should calculate average solar hour" do
+
+      solar_time_1 = SolarGeometryCalculation::SolarTime.create_by_standard_time(2012, 02, 19, 12, 30, -75, -73.75, "America/Los_Angeles")
+      solar_time_1.year.should == 2012
+      solar_time_1.month.should == 2
+      solar_time_1.day.should == 19
+      solar_time_1.hour.should == 12
+      solar_time_1.minute.should be_within(1).of(21)
+      
+      solar_time_2 = SolarGeometryCalculation::SolarTime.create_by_standard_time(2012, 02, 19, 12, 30, BigDecimal("-75"), BigDecimal("-73.75"), "America/Los_Angeles")
+      solar_time_1.year.should == 2012
+      solar_time_1.month.should == 2
+      solar_time_1.day.should == 19
+      solar_time_1.hour.should == 12
+      solar_time_1.minute.should be_within(1).of(21)     
+
+    end    
+    
+  end
+  
+  context "converter" do
+    
+    class TempSolarTimeConverter; include SolarGeometryCalculation::SolarTimeConverter; end
+
+    it "should calculate the adjustment by the Equation of time" do
+    
+      TempSolarTimeConverter.new.send(:equation_of_time, 2012, 02, 19).should be_within(0.01).of(-14.11)
+  
+    end
+  
+  end
+
+end

data/spec/spec_helper.rb

@@ -0,0 +1,2 @@
+require 'solar_geometry_calculation'
+

metadata

@@ -0,0 +1,107 @@
+--- !ruby/object:Gem::Specification
+name: solar_geometry_calculation
+version: !ruby/object:Gem::Version
+  version: 0.1.4
+platform: ruby
+authors:
+- Tadatoshi Takahashi
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-08-18 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: tzinfo
+  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: activesupport
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 4.1.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 4.1.1
+- !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: Mathematical part of obtaining solar geometry information
+email:
+- tadatoshi@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".ruby-gemset"
+- ".ruby-version"
+- Gemfile
+- LICENSE
+- README.md
+- Rakefile
+- lib/solar_geometry_calculation.rb
+- lib/solar_geometry_calculation/calculation_helper.rb
+- lib/solar_geometry_calculation/solar_position_calculation.rb
+- lib/solar_geometry_calculation/solar_time.rb
+- lib/solar_geometry_calculation/version.rb
+- solar_geometry_calculation.gemspec
+- spec/solar_geometry_calculation/calculation_helper_spec.rb
+- spec/solar_geometry_calculation/solar_position_calculation_spec.rb
+- spec/solar_geometry_calculation/solar_time_spec.rb
+- spec/spec_helper.rb
+homepage: https://github.com/tadatoshi/solar_geometry_calculation
+licenses: []
+metadata: {}
+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: solar_geometry_calculation
+rubygems_version: 2.2.2
+signing_key: 
+specification_version: 4
+summary: Performs calculation for solar geometry
+test_files:
+- spec/solar_geometry_calculation/calculation_helper_spec.rb
+- spec/solar_geometry_calculation/solar_position_calculation_spec.rb
+- spec/solar_geometry_calculation/solar_time_spec.rb
+- spec/spec_helper.rb