sun 0.1.0 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: dc8811fc890b22bef67635691e8c41411c57d2d1
-  data.tar.gz: 50804066785010d6a19a051631b0ee0582e88fcc
+  metadata.gz: 62b638fc2309bc08dff55796cee15932fd4d791d
+  data.tar.gz: c6c556fafabec3ea2d959d1c2c13c80299d15836
 SHA512:
-  metadata.gz: 57dbb11a8ad58ba70cfff77d5378aab61f931a5fc13bbe36ea156021a76ddde0e2ea0cfbac064ea9eb5d8a7527558d812f45b2b47177cdb232c282a640a22030
-  data.tar.gz: 7db0301164c792ebf20d4fd02aabea6199cf1564c47edc7fc344dac5882381bca1f3174aa092822bbd38263b4934d71130b3a0fa02c745e0ee5d7fda5a8e9448
+  metadata.gz: 0b8411da7c981e71efb2584e05278c7f001bb3ca069feb5281dcb7af8782b58d5f7d2ff83f7156ed14a2eb29c0412f8e91e789fbe3d6aaa9283a3637f54539ae
+  data.tar.gz: 7747c1543b1769044b0430c7ed6257b88f4d903eb688e04aed82bc5521c3b016e2bdc3e6e942cf6ac77832b1d17473160273b5ce3fe187572d152883bf9a8b2c

data/lib/sun.rb

@@ -1,5 +1,121 @@
 require 'sun/version'
 
 module Sun
-  # Your code goes here...
+  SOLAR_ZENITH = 90.833
+
+  # Noon on January 1st, 2000
+  JULIAN_CONSTANT = 2451545.to_r
+
+  # Sun times
+  def self.solar_noon(time, latitude, longitude)
+    minutes = solar_noon_minutes(time, latitude, longitude)
+    date_at_time(time.to_date, minutes)
+  end
+
+  def self.sunrise(time, latitude, longitude)
+    minutes = sunrise_minutes(time, latitude, longitude)
+    date_at_time(time.to_date, minutes)
+  end
+
+  def self.sunset(time, latitude, longitude)
+    minutes = sunset_minutes(time, latitude, longitude)
+    date_at_time(time.to_date, minutes)
+  end
+
+  # Helpers
+  def self.degrees(radians)
+    180.0 * radians / Math::PI
+  end
+
+  def self.radians(degrees)
+    Math::PI * degrees / 180.0
+  end
+
+  def self.date_to_unix_time(date)
+    Time.utc(date.year, date.month, date.day).to_i
+  end
+
+  def self.date_at_time(date, minutes)
+    date_to_unix_time(date) + minutes * 60
+  end
+
+  def self.julian_days(time)
+    time.utc.to_datetime.ajd
+  end
+
+  def self.julian_century(time)
+    (julian_days(time) - JULIAN_CONSTANT) / 36525
+  end
+
+  def self.mean_obliquity_of_ecliptic(julian_century)
+    23 + (26 + ((21.448 - julian_century * (46.815 + julian_century * (0.00059 - julian_century * 0.001813)))) / 60) / 60
+  end
+
+  def self.oblique_correction(julian_century)
+    mean_obliquity_of_ecliptic(julian_century) + 0.00256 * Math.cos(radians(125.04 - 1934.136 * julian_century))
+  end
+
+  def self.geometric_mean_anomoly(julian_century)
+    357.52911 + julian_century * (35999.05029 - 0.0001537 * julian_century)
+  end
+
+  # MOD(280.46646+G2*(36000.76983 + G2*0.0003032),360)
+  def self.geometric_mean_longitude(julian_century)
+    280.46646 + julian_century * (36000.76983 + julian_century * 0.0003032) % 360
+  end
+
+  def self.y(oblique_correction)
+    Math.tan(radians(Rational(oblique_correction, 2))) * Math.tan(radians(Rational(oblique_correction, 2)))
+  end
+
+  def self.eccentricity_of_earth_orbit(julian_century)
+    0.016708634 - julian_century * (0.000042037 + 0.0000001267 * julian_century)
+  end
+
+  def self.equation_of_center(julian_century)
+    geometric_mean_anomoly = geometric_mean_anomoly(julian_century)
+    Math.sin(radians(geometric_mean_anomoly)) * (1.914602 - julian_century * (0.004817 + 0.000014 * julian_century)) + Math.sin(radians(2 * geometric_mean_anomoly)) * (0.019993 - 0.000101 * julian_century) + Math.sin(radians(3 * geometric_mean_anomoly)) * 0.00028
+    # =SIN(RADIANS(J2))*(1.914602-G2*(0.004817+0.000014*G2))+SIN(RADIANS(2*J2))*(0.019993-0.000101*G2)+SIN(RADIANS(3*J2))*0.00028
+  end
+
+  def self.true_longitude(julian_century)
+    geometric_mean_longitude(julian_century) + equation_of_center(julian_century)
+  end
+
+  def self.apparent_longitude(julian_century)
+    true_longitude(julian_century) - 0.00569 - 0.00478 * Math.sin(radians(125.04 - 1934.136 * julian_century))
+  end
+
+  def self.declination(oblique_correction, julian_century)
+    degrees(Math.asin(Math.sin(radians(oblique_correction)) * Math.sin(radians(apparent_longitude(julian_century)))))
+  end
+
+  def self.equation_of_time(time, longitude)
+    julian_century = julian_century(time)
+    oblique_correction = oblique_correction(julian_century)
+    geometric_mean_anomoly = geometric_mean_anomoly(julian_century)
+    geometric_mean_longitude = geometric_mean_longitude(julian_century)
+    eccentricity_of_earth_orbit = eccentricity_of_earth_orbit(julian_century)
+    y = y(oblique_correction)
+    4 * degrees(y * Math.sin(2 * radians(geometric_mean_longitude)) - 2 * eccentricity_of_earth_orbit * Math.sin(radians(geometric_mean_anomoly)) + 4 * eccentricity_of_earth_orbit * y * Math.sin(radians(geometric_mean_anomoly)) * Math.cos(2 * radians(geometric_mean_longitude)) - 0.5 * y * y * Math.sin(4 * radians(geometric_mean_longitude)) - 1.25 * eccentricity_of_earth_orbit * eccentricity_of_earth_orbit * Math.sin(2 * radians(geometric_mean_anomoly)))
+  end
+
+  def self.hour_angle(time, latitude)
+    julian_century = julian_century(time)
+    oblique_correction = oblique_correction(julian_century)
+    declination = declination(oblique_correction, julian_century)
+    degrees(Math.acos(Math.cos(radians(90.833)) / (Math.cos(radians(latitude)) * Math.cos(radians(declination))) - Math.tan(radians(latitude)) * Math.tan(radians(declination))))
+  end
+
+  def self.solar_noon_minutes(time, latitude, longitude)
+    720 - (4 * longitude) - equation_of_time(time, longitude)
+  end
+
+  def self.sunrise_minutes(time, latitude, longitude)
+    solar_noon_minutes(time, latitude, longitude) - 4 * hour_angle(time, latitude)
+  end
+
+  def self.sunset_minutes(time, latitude, longitude)
+    solar_noon_minutes(time, latitude, longitude) + 4 * hour_angle(time, latitude)
+  end
 end

data/lib/sun/version.rb

@@ -1,3 +1,3 @@
 module Sun
-  VERSION = "0.1.0"
+  VERSION = "0.1.1"
 end

data/sun.gemspec

@@ -26,7 +26,10 @@ Gem::Specification.new do |spec|
   spec.executables = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
   spec.require_paths = ['lib']
 
-  spec.add_development_dependency 'bundler', '~> 1.10'
-  spec.add_development_dependency 'rake', '~> 10.0'
+  spec.add_development_dependency 'bundler'
+  spec.add_development_dependency 'byebug'
   spec.add_development_dependency 'minitest'
+  spec.add_development_dependency 'minitest-reporters'
+  spec.add_development_dependency 'rake'
+  spec.add_development_dependency 'tzinfo'
 end

metadata

@@ -1,43 +1,43 @@
 --- !ruby/object:Gem::Specification
 name: sun
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.1
 platform: ruby
 authors:
 - Cory Kaufman-Schofield
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2015-08-26 00:00:00.000000000 Z
+date: 2015-09-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ~>
+    - - '>='
       - !ruby/object:Gem::Version
-        version: '1.10'
+        version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ~>
+    - - '>='
       - !ruby/object:Gem::Version
-        version: '1.10'
+        version: '0'
 - !ruby/object:Gem::Dependency
-  name: rake
+  name: byebug
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ~>
+    - - '>='
       - !ruby/object:Gem::Version
-        version: '10.0'
+        version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ~>
+    - - '>='
       - !ruby/object:Gem::Version
-        version: '10.0'
+        version: '0'
 - !ruby/object:Gem::Dependency
   name: minitest
   requirement: !ruby/object:Gem::Requirement
@@ -52,6 +52,48 @@ dependencies:
     - - '>='
       - !ruby/object:Gem::Version
         version: '0'
+- !ruby/object:Gem::Dependency
+  name: minitest-reporters
+  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'
+- !ruby/object:Gem::Dependency
+  name: rake
+  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'
+- !ruby/object:Gem::Dependency
+  name: tzinfo
+  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: 
 email:
 - cory@corykaufman.com