equationoftime 4.0.0 → 4.1.1

data/examples/ptime.rb

@@ -1,162 +0,0 @@
-require 'date'
-require 'bigdecimal'
-include Math
-
-
-def degrees_to_radians(degrees)
-    bd(degrees * PI / 180.0)
-end
-
-def radians_to_degrees(radians)
-    bd(radians * 180.0 / PI)
-end
-
-# Truncate large angles
-def truncate(x)
-    bd(360.0 * ( x / 360.0 - Integer( x / 360.0)))
-end
-
-def bd(var)
-    BigDecimal(var.to_s).round(16)# * 1.0
-end
-
-def calc_time_julian_century(t)
-# Julian Day Number j(2000) subtracted
-  (t - 2451545.0)/36525.0
-# Time in fractional centurey
-end
-
-def calc_eccentricity_earth_orbit(t)
-    0.016708617 - t * (0.000042037 + 0.0000001235 * t)#unitless
-    #~ puts "EccentricityEarthOrbit = #{eccentricity}"
-end
-
-def mean_anomaly(t)
-  (357.527723333 + truncate(35640.0 * t) \
-         + t * (359.05034 - t * (0.0001602777778 \
-             + t / 300000.0))) % 360.0
-#~ t2 = t **2
-#~ t3 = t **3
-#~ t4 = t **4
-#~ @mean_anomaly = 1287104.79305 + 129596581.0481 * t - 0.5532 * t2 + 0.000136 * t3 - 0.00001149 * t4
-#~ @mean_anomaly = 357.527723333 + 359.05034 * t + trunc(35640.0 * t) - 0.0001602777778 * t2 - t3/300000.0
-     #~ @mean_anomaly = 357.52910 + t * 35999.05034 - t2 * 0.0001559 - t3 * 0.00000048
-     #~ puts "GeomMeanAnomalySun = #{mean_anomaly%360}"
-     #~ return mean_anomaly % 360 # in degrees
-end
-
-def calc_geom_mean_long_sun(t)
-    #~ //Mean longitude of the sun
-    t2 = t * t
-    t3 = t2 * t
-    t4 = t3 * t
-    t5 = t4 * t
-    truncate(280.4664567 + t * 36000.76982779 \
-        + t2 * 0.0003032028 + t3/49931.0 \
-        - t4/15299.0 - t5/1988000.0)
-end
-
-def sun_eq_of_center(t)
-    #t2 = t **2
-    mean_anomaly = mean_anomaly(t)
-    #~ puts "GeomMeanAnomalySun = #{m}"
-    marad = degrees_to_radians(mean_anomaly)
-    sinmarad = sin(marad)
-    eccentricity = calc_eccentricity_earth_orbit(t)
-    sin2marad = sin(2.0 * marad)
-    sin3marad = sin(3.0 * marad)
-    sinmarad * (1.914602 - t * (0.004817 + 0.000014 * t))\
-         + sin2marad * (0.019993 - 0.000101 * t)\
-         + sin3marad * 0.000289
-end
- 
-def true_anomaly(t)
-    mean_anomaly = mean_anomaly(t)
-    eqc = sun_eq_of_center(t)
-    bd(mean_anomaly + eqc) * 1.0
-end
-
-def calc_sun_true_long(t)
-    mean_lon = calc_geom_mean_long_sun(t)
-    eqc = sun_eq_of_center(t)
-    mean_lon + eqc - 0.01397 * (Time.now.year - 2000)
-end
-
-def mean_obliquity_of_ecliptic(t)
-    bd(23.4392911111 - t * (0.0130041666667 + t * (1.63888888889e-07 - (5.03611111111e-07 * t)))) * 1.0
-end
-
-def omega(t)
-    bd(125.04455501 - truncate(1800.0 * t)   \
-		    - t * (134.136261972     \
-		    - t * (0.00207561111111  \
-		    + t * (2.13944444444e-06 \
-		    - t * 1.64972222222e-08  \
-    )))) * 1.0
-    bd(125.04455501 - truncate(1800.0 * t)              \
-		    - t * 134.136261972                 \
-		    + t * t * 0.00207561111111          \
-		    + t * t * t * 2.13944444444e-06     \
-		    - t * t * t * t * 1.64972222222e-08 \
-       ) * 1.0
-end
-  
-def obliquity_correction(t)
-    eps0 = mean_obliquity_of_ecliptic(t)
-    omega = omega(t)
-    delta_eps = bd(0.00256 * cos(degrees_to_radians(omega)))
-    bd(eps0 + delta_eps)
-end
-
-def sun_right_ascension(t)
-    lambda = calc_sun_true_long(t)
-    epsilon = obliquity_correction(t)
-    y0 = bd(sin(degrees_to_radians(lambda)))
-    x0 = bd(cos(degrees_to_radians(lambda)))
-    y0 = bd(y0 * cos(degrees_to_radians(epsilon)))
-    bd(radians_to_degrees(atan2(-y0, -x0)) + 180.0) * 1.0
-end
-
-def mean_long_aries(t)
-    d = bd(t * 36525.0)
-    bd(truncate(280.46061666 + d * 360.98564736629\
-                             + t * (t * 0.000387933\
-			     - t * (t / 38710000.0)))) * 1.0
-end
-
-pdate = Date.new(2014, 1, 3)
-pd = pdate.ajd * 1.0
-#puts pd
-daysecs = 2 * 84600.0
-
-for seconds in 1..daysecs
-  #puts seconds / 84600.0
-  jdt = calc_time_julian_century(pd + seconds / 84600.0)
-  mas = bd(mean_anomaly(jdt)).round(3) * 1.0
-  #puts mas
-  tas = bd(true_anomaly(jdt)).round(3) * 1.0
-  #puts tas
-  if (mas + tas) == 720
-    d = pd + seconds / 84600.0
-    #puts d
-    pt = jdt % 1.0
-    break
-  #~ else
-    #~ d = 0
-    #~ pt = 0
-  end
-end
-
-#~ puts d
-#~ puts pt
-puts "2014 perihelion occurred at #{d}"
-tls = calc_sun_true_long(pt)
-puts "lambda perihelion = #{tls}"
-sra = sun_right_ascension(pt)
-puts "Right ascesion = #{sra}"
-ma = mean_long_aries(pt)
-puts "GHA Aries = #{ma}"
-
-
-
-

data/examples/rbeot.c

@@ -1,24 +0,0 @@
-#include "ruby.h"
-
-#ifndef DBL2NUM
-# define DBL2NUM(dbl) rb_float_new(dbl)
-#endif
-
-VALUE cEot;
-ID id_status;
-
-static VALUE
-ma_Sun(VALUE self, VALUE vt) {
-
-  rb_ivar_set(self, id_status, INT2FIX(0));
-  return DBL2NUM(maSun(NUM2DBL(vt)));
-}
-
-
-void
-Init_rbeot(void) {
-
-    cEot = rb_define_class("EOT", rb_cObject);
-    id_status = rb_intern("@status");
-    rb_define_method(cEot, "ma", ma_Sun, 1);
-}

data/examples/suntimes.rb

@@ -1,28 +0,0 @@
-# suntimes.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-# The following gem is not in rubygems.org yet. Please use the gem on github. Thanks!
-require 'eot'
-
-eot = Eot.new
-# Change these for your geo location
-@latitude = 41.9474  # + in Northern Hemishere,  - in Southern Hemisphere
-@longitude = -88.74467  # + East of zulu time zone, - West of zulu time zone
-@zone = -5 # Time offset for zone.
-
-
-eot.longitude= @longitude
-eot.latitude= @latitude
-puts
-@date = Date.today
-# get the times for a few more days.
-for i in 0..10
-
-  eot.ajd = (@date.jd + i) 
-
-  puts "Sunrise #{(eot.sunrise_dt().to_time).strftime("%c")}"
-  puts "Sunset #{(eot.sunset_dt().to_time).strftime("%c")}"
-  puts
-end

data/examples/suntimes_test.rb

@@ -1,47 +0,0 @@
-# suntimes_test.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require 'eot'
-
-eot = Eot.new
-
-# Change these for your geographic location Ex: of mine.
-# Note the two coordinate parameters must be passed as floats.
-# All parameters may be passed at any time.
-# eot.date = @date
-# eot.longitude = @longitude
-# eot.latitude = @latitude
-@latitude = 0.0 # sign in Northern Hemishere is none,  - (minus)for Southern Hemisphere
-@longitude = 0.0  # East of UTC is none, - (minus)West of UTC
-@date = "2000-01-01"
-eot.jd = Date.parse(@date).jd
-
-print "Enter your longitude in decimal format Ex: -125.099 \n"
-
-lng = gets
-
-lng.nil? ? eot.longitude = @longitude : eot.longitude = lng.chomp.to_f		
-
-puts "Your longitude is #{eot.longitude}"
-
-print "Enter your latitude in decimal format Ex: 39.3339 \n"
-
-lat = gets
-
-lat.nil? ? eot.latitude = @latitude : eot.latitude = lat.chomp.to_f
-
-puts "Your latitude is #{eot.latitude}"
-
-print "Enter the date you want to compute sun times for Ex: yyyy-mm-dd \n"
-
-date = gets
-
-date.nil? ? eot.date = @date : eot.date = date.chomp.to_f
-
-puts "The date is #{eot.date}"
-
-puts "Sunrise #{eot.sunrise_dt.to_time}"
-puts "Sunset #{eot.sunset_dt.to_time}"
-

data/examples/test_eot.rb

@@ -1,17 +0,0 @@
-require_relative 'rbeot'
-require_relative '../lib/eot'
-
-rbeot = EOT.new
-eot = Eot.new
-
-p eot.ta
-p rbeot.ma(eot.ta)
-#~ p eot.ma_Sun
-
-require 'benchmark'
-
-n = 5_000
-Benchmark.bm do |x|
-  x.report("eot")  { n.times { rbeot.ma(eot.ta) } }
-  #~ x.report("ruby")  { n.times {eot.ma_Sun() } }
-end

data/examples/test_poly_eval.rb

@@ -1,27 +0,0 @@
-# test_poly_eval.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-# From astro-algo gem 
-# class Array    
-    # Evaluate polynomial using Horner's method.
-    # Array consists of coefficients of a polynomial, the
-    # coefficient of the highest order term first, the
-    # constant coefficient last.
-    # Returns evaluation of polynomial at +x+.
-    #
-    # Example: evaluate the polynomial x**2 - 0.5*x + 3.0 where x = 2.0
-    #    [1.0, -0.5, 3.0].poly_eval(2.0)   # => 6.0
-    # def poly_eval(x)
-        # self.inject(0.0) {|p, a|		
-        # p*x + a}
-    # end
-# end
-
-my_array   = [1.0, -0.5, 3.0]
-my_array.inject(0.0) {|sum, n| p sum * 2.0 + n }
-
-ecc = [-0.0000001235, -0.000042037, 0.016708617]
-ta  = 0.14
-ecc.inject(0.0) {|sum, n| p sum * ta + n }  

data/examples/time_scales.rb

@@ -1,26 +0,0 @@
-# eot_methods_list.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require  'eot'
-$DEBUG and set_trace_func proc { |event, file, line, id, binding, classname|
-       printf "%8s %s:%-2d %10s %8s\n", event, file, line, id, classname
-    }
-eot = Eot.new
-# puts eot.public_methods(false).sort
-# puts eot.nil?
-loop do
-  puts DateTime.now.to_time.utc.nsec.inspect
-  puts DateTime.jd(DateTime.now.to_time.utc.to_datetime.ajd + 0.5).to_time.nsec.inspect
-  sleep 0.7
-  puts
-end
-
-# puts Astro.solar_longitude(date)
-# tjc = eot.time_julian_century( ajd)
-# puts eot.tl_Sun(tjc)
-
-# tjc = eot.time_julian_century( ajd + dt )
-# puts eot.tl_Sun(tjc)
-

data/examples/times_year.rb

@@ -1,50 +0,0 @@
-# a list of sunrise and sunset times for a year of dates.
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require 'eot'
-require 'safe_yaml'
-require 'time'
-
-eot    = Eot.new
-
-# set your latitude and longitude first.
-eot.latitude  = 41.9474
-eot.longitude = -88.74467
-
-# make some start and finish dates.
-start  = "2014-1-1"
-finish = "2014-12-31"
-
-start_time  = Time.utc( 2014, "jan", 1, 12, 0, 0 )
-finish_time = Time.utc( 2014, "dec", 31, 12, 0, 0 )
-
-start_jd  = Date.parse(start).jd
-finish_jd = Date.parse(finish).jd
-
-fstr = "%b %d"
-
-@data, @group, @group_id = [], {}, 1
-
-(start_jd..finish_jd).each do |jd|     
-#  date    = DateTime.jd(jd + 0.5).to_date
-  eot.ajd = jd
-  rise    = eot.sunrise_dt().to_time.to_json
-  trans   = eot.local_noon_dt().to_time.to_json
-  set     = eot.sunset_dt().to_time.to_json
-  @group  = {  "id" => @group_id, 
-             "rise" => rise,
-             "noon" => trans,  
-             "sset" => set
-           }  
-  @data << @group
-  @group_id  += 1
-end
-
-file_path = "rise_set_data.yml"
-File::open( file_path, "w" ) do |f|
-  YAML.dump( @data, f )
-end
-
-puts "File rise_set_data.yml processed"             

data/examples/usage_example.rb

@@ -1,23 +0,0 @@
-# usage_example.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require 'eot'
-
-eot = Eot.new()
-# find what the defaults are set to
-show = <<RAN
-ajd       = #{eot.ajd}
-date      = #{eot.date}
-jd        = #{eot.jd}
-longitude = #{eot.longitude}
-latitude  = #{eot.latitude}
-sunrise   = #{eot.sunrise_dt}
-trnasit   = #{eot.local_noon_dt}
-min eot   = #{eot.time_eot}
-sunset    = #{eot.sunset_dt}
-gasa      = #{eot.tl_Aries}
-RAN
-
-puts show

data/examples/use_angles.rb

@@ -1,219 +0,0 @@
-# use_angles.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require 'eot'
-
-
-file = <<DISPLAY
-
-This will explain how the equation of time is derived by showing each method that comprises the formula.
-
-eot = Eot.new
-
-Instance of Eot class eot = #{eot = Eot.new}
-
-
-Set the class attribute ajd for the deisired date ex : eot.ajd = Date.parse("yyyy, mm, dd").jd.to_f
-
-eot.ajd = Date.today.jd.to_f
-
-eot.ajd = #{eot.ajd = Date.today.jd.to_f}
-
-The class attribute ajd has been set to a nice whole number even though it is a Float class now.
-
-The Astronomical Julian Day Number is a half day less than Julian Day Number with the fractional day time included when we DateTime class.
-
-DateTime.now.day_fraction.to_f = #{DateTime.now.day_fraction.to_f}
-
-DateTime class has the method day_fraction and Date class does not.
-
-We used Date class which doesn't give us the UTC time of date. It yields your time zone date.
-
-Let's just prove this.
-
-Date.today.to_time.utc = #{Date.today.to_time.utc}
-
-It looks like that was your midnight in UTC time.
-
-Maybe we could get away with getting the ajd from that.
-
-Date.today.to_time.utc.to_date.ajd.to_f = #{Date.today.to_time.utc.to_date.ajd.to_f}
-
-No. That's a half day less than the Julian Day Number.
-
-So let's start over setting the adj correctly for UTC using the DateTime class.
-
-eot.ajd = DateTime.now.to_time.utc.to_datetime.ajd.to_f
-
-eot.ajd = #{eot.ajd = DateTime.now.to_time.utc.to_datetime.ajd.to_f}
-
-
-Because this gem is going to calculate the sunrise and sunset times it's okay to set the day for noon and is encouraged.
-
-So go ahead and set the ajd using the Date class and jd method as was first done.
-
-The class init.rb sets everything to defaults and one of them is today's datetime.jd 
-
-But for right now we won't do that so we get current angles of the present time UTC.
-
-
-
-eot.ma is the solar Mean anomaly mu(M) and is called quite often hence the attribute.
-
-eot.ma was set using eot.ma_Sun() method.
-
-It is an angle with respect to perihelion of Earth orbit around Sun in an assumed circular orbit.
-
-M is calculated within the nutation.rb file under delta_equinox() method and returned as the third element of an array.
-
-eot.ma = #{eot.ma}
-
-eot.ma_Sun() = #{eot.ma_Sun()}
-
-
-angle one = eot.ma_Sun() - eot.ta_Sun()
-
-a1 = #{eot.ma_Sun() - eot.ta_Sun()} degrees
-
-a1 is a delta angle difference for mean and true anomalies in degrees.
-
-True anomaly (nu) is the actual angle with respect to perihelion in an eliptical orbit.
-
-It was calculated by using the equation of center formula which you may access also.
-
-eot.center() = #{eot.center()} degrees.
-
-nu = equation of center + M or mu or mean anomaly
-
-The delta is M - nu.
-
-
-The method inside the gem is called delta_orbit()
-
-eot.delta_orbit = #{eot.delta_orbit} degrees
-
-We could just change the sign of eot.center().  
-
-In fact let's just do that now. a1 = - eot.center()
-  
-- eot.center() = #{- eot.center()}
-
-
-The next delta is calculated as follows.
-
-lambda or apparent longitude of the Sun minus alpha or right ascension of the Sun.
-
-a2 = lambda - alpha
-
-lambda = mean longitude + equation of center
-
-eot.gml_Sun() + eot.center() = #{eot.gml_Sun() + eot.center()} or
-
-eot.tl_Sun = #{eot.tl_Sun}
-
-Lambda is the \"ecliptic coordinate system\" angle of the sun.
-
-Alpha is the Right Ascension of the sun in the celestial coordinate system.
-
-alpha = eot.ra_Sun()
-
-alpha = #{eot.ra_Sun()}
-
-lambda - alpha = #{eot.tl_Sun - eot.ra_Sun()}
-
-The method inside the gem is called delta_oblique.
-
-eot.delta_oblique() = #{eot.delta_oblique()}
-
-
-The sum of these two delta angles is the equation of time angle in degrees not time.
-
-eot.eot() = #{eot.eot()}
-
-That was for the time UTC now which was easy.
-
-
-Now using your longitude to compute your mean solar transit from longitude / 15.0
-
-my longitude is -88.74467 deg. so I'll set the longitude attribute in Eot class.
-
-eot.longitude = #{eot.longitude = -88.74467}
-
-
-So my mean solar transit UTC is 12.0 - my_longitude / 15.0
-
-my_mean_noon = #{12.0 - -88.74467 / 15.0} hr.
-
-Be sure to set everything for UTC noon before adding it to your longitude because we were doing angles for time now.
-  
-eot.ajd = #{eot.ajd = Date.today.jd.to_f}
-  
-
-
-Using the Eot class I may also compute it with eot.mean_local_noon_dt() method if I have the ajd set right.
-
-eot.mean_local_noon_dt() = #{eot.mean_local_noon_dt()}
-
-I'm showing what time mean noon for my_longitude normally is without the equation of time considered.
-
-This time is true only four days a year. Lets get the right one for today to be sure.
-
-
-
-We can add longitude plus equation of time as angles. 
-
-eot.longitude + eot.eot() = #{eot.longitude} + #{eot.eot()}
-
-What's the time of solar transit here now?
-
-12.0 - (eot.longitude + eot.eot()) / 15.0 = #{12.0 - (eot.longitude + eot.eot()) / 15.0}
-  
-That looks about right.
-
-
-
-Let's make a time out of that by converting it to a fractional day.
-
-(12.0 - (eot.longitude + eot.eot()) / 15.0) / 24.0 = #{(12.0 - (eot.longitude + eot.eot()) / 15.0) / 24.0}
-  
-We'll add that to the ajd then.
-
-ajd = Date.today.ajd.to_f
-
-ajd = #{Date.today.ajd.to_f}
-
-#{Date.today.ajd.to_f} + #{(12.0 - (eot.longitude + eot.eot()) / 15.0) / 24.0} = #{Date.today.ajd.to_f + (12.0 - (eot.longitude + eot.eot()) / 15.0) / 24.0}
-
-
-There is a method in Eot class to make that a DateTime object.
-
-eot.ajd_to_datetime(ajd + transit)
-
-#{eot.ajd_to_datetime(Date.today.ajd.to_f + (12.0 - (eot.longitude + eot.eot()) / 15.0) / 24.0)}
-
-
-Be sure to set your latitude attribute as well when you want to compute sunrise and sunset times for the days of interest.
-
-eot.latitude = my_latitude
-
-
-Your longitude and the GHA of the Sun will match at True Solar Transit time.
-
-Check it here http://douglasallen.github.io/planets/ just put in the correct time and date then hit calculate.
-
-Does GHA Sun = your longitude?
-
-It's off a little bit because that site uses a different formula that is commonly used in older programs.
-
-I get about 3 seconds difference in the EOT is all.
-DISPLAY
-
-
-
-
-
-
-
-puts file