equationoftime 4.1.2 → 4.1.4

data/examples/test_celes.rb

@@ -1,51 +0,0 @@
-#require_relative 'celes_core'
-require 'celes'
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-
-eot = Eot.new
-
-# Celes.methods.sort.each {|m| p m }
-
-dt = DateTime.now.to_time.utc.to_datetime
-jd = Celes.cal2jd(dt.year, dt.month, dt.day)
-
-ajd = jd[0] + jd[1] + dt.day_fraction
-
-p eot.ajd = ajd
-
-p eot.delta_equinox()[0]
-p Celes.nut06a(ajd,0)[0]
-
-p eot.delta_equinox()[1]
-p Celes.nut06a(ajd,0)[1]
-
-p eot.delta_equinox()[2]
-p Celes.falp03(eot.ta)
-
-p eot.delta_equinox()[3]
-p Celes.faom03(eot.ta)
-
-require 'benchmark'
-
-n = 5_000
-Benchmark.bm do |x|
-  #~ x.report("celes")  { n.times { Celes.cal2jd(dt.year, dt.month, dt.day) } }
-  #~ x.report("ruby")  { n.times { DateTime.new(dt.year, dt.month, dt.day).jd } }
-  x.report("celes")  { n.times { Celes.nut06a(ajd,0)[0] }}
-  x.report("celes")  { n.times { Celes.nut06a(ajd,0)[1] }}
-  x.report("celes")  { n.times { Celes.falp03(eot.ta) }}
-  x.report("celes")  { n.times { Celes.faom03(eot.ta) }}
-
-  x.report("ruby")  { n.times { eot.delta_equinox()[0] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[1] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[2] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[3] }}
-end
-

data/examples/test_ceot.rb

@@ -1,55 +0,0 @@
-require 'benchmark'
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-
-require 'eot'
-
-eot = Eot.new
-eot.ajd = DateTime.now.to_time.utc.to_datetime.ajd
-eot.ajd_to_datetime(eot.ajd)
-
-p eot.eccentricity_Earth()
-p eot.eoe(eot.ta)
-
-n = 500_000
-Benchmark.bm do |x|
-  x.report("reot")  { n.times {eot.eccentricity_Earth()} }
-  x.report("ceot")  { n.times {eot.eoe(eot.ta)} }  
-end
-puts
-
-p eot.center() * Eot::R2D 
-p eot.eqc(eot.ma, eot.ta) * Eot::R2D
-
-Benchmark.bm do |x|
-  x.report("reot") {n.times {eot.center()}}
-  x.report("ceot") {n.times {eot.eqc(eot.ma, eot.ta)}}
-end
-puts 
-
-p eot.tl_Sun()
-p eot.tl(eot.ma, eot.ta)
-
-n = 500_000
-Benchmark.bm do |x|
-  x.report("reot")  { n.times {eot.tl_Sun()} }
-  x.report("ceot")  { n.times {eot.tl(eot.ma, eot.ta)} }  
-end
-puts
-
-p eot.al_Sun()
-p eot.al(eot.ma, eot.ta, eot.omega())
-
-n = 500_000
-Benchmark.bm do |x|
-  x.report("reot")  { n.times {eot.al_Sun()} }
-  x.report("ceot")  { n.times {eot.al(eot.ma, eot.ta, eot.omega())} }  
-end
-puts
-

data/examples/test_poly_eval.rb

@@ -1,32 +0,0 @@
-# test_poly_eval.rb
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-
-# 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,29 +0,0 @@
-# eot_methods_list.rb
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-$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,53 +0,0 @@
-# a list of sunrise and sunset times for a year of dates.
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-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,26 +0,0 @@
-# usage_example.rb
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-
-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,222 +0,0 @@
-# use_angles.rb
-
-begin
-  require 'eot'
-rescue LoadError
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  require 'eot' 
-end
-
-
-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

data/lib/eot/displays.rb

@@ -1,216 +0,0 @@
-# displays.rb
-
-class Eot
-
-  # From displays.rb  
-  # String formatter for d:m:s display 
-  def degrees_to_s( radians = 0.0 )
-    radians.nil? ? radians = 0.0 : radians
-    radians < 0 ? sign_string      = "-" : sign_string = "+"    
-                  absolute_degrees = radians.abs * R2D
-          absolute_degrees_integer = Integer( absolute_degrees )
-          absolute_decimal_minutes = 60.0 * 
-                                     (
-                                      absolute_degrees - 
-                                      absolute_degrees_integer
-                                     )
-                                     
-          absolute_minutes_integer = Integer( absolute_decimal_minutes )
-          
-          absolute_decimal_seconds = 60.0 * 
-                                     (
-                                      absolute_decimal_minutes - 
-                                      absolute_minutes_integer
-                                      )
-                                      
-          absolute_seconds_integer = Integer( absolute_decimal_seconds )
-          
-    absolute_milli_seconds_integer = Integer(1000.0 *
-	                                           (  
-                                              absolute_decimal_seconds - 
-                                              absolute_seconds_integer
-                                             )
-                                            )
-                            sign_string +
-      "%03d" % absolute_degrees_integer +
-                                    ":" + 
-      "%02d" % absolute_minutes_integer +
-                                    ":" + 
-      "%02d" % absolute_seconds_integer +
-                                    "." + 
-      "%3.3d" % absolute_milli_seconds_integer
-  end  
-  
-  # From displays.rb  
-  # String formatter for + and - time 
-  def show_minutes(min = 0.0)
-    min.nil? ? min = 0.0 : min
-    time = Time.utc(1, 1, 1, 0, 0, 0, 0.0)
-    time = time + (min.abs * 60.0)
-    if min < 0.0
-      sign = "-"
-    else
-      sign = "+"
-    end
-    time.strftime("#{sign}%M:%S.%3N")
-  end
-
-  # From displays.rb  
-  # String for time now
-  def show_now(now = now(Time.now.utc))
-    show_minutes(now)
-  end  
-  
-  # From displays.rb  
-  # String format of apparent longitude 
-  def string_al_Sun()
-    degrees_to_s( al_Sun() )
-  end
-  alias_method :apparent_longitude_string, :string_al_Sun
-
-  # From displays.rb  
-  # String formatter for fraction of Julian day number
-  def string_day_fraction_to_time( jpd_time = 0.0 )
-    jpd_time.nil? ? jpd_time = 0.0 : jpd_time
-    fraction = jpd_time + 0.5 - Integer( jpd_time )
-    hours = Integer( fraction * DAY_HOURS )
-    minutes = Integer(( fraction - hours / DAY_HOURS ) * DAY_MINUTES )
-    seconds = Integer(( fraction - hours / 24.0 - minutes / DAY_MINUTES ) * DAY_SECONDS )
-    "%02d" % hours   +
-    ":"              +
-    "%02d" % minutes +
-    ":"              +
-    "%02d" % seconds
-  end 
-  alias_method :julian_period_day_fraction_to_time, :string_day_fraction_to_time 
-
-  # From displays.rb  
-  # String format of declination 
-  def string_dec_Sun()
-    degrees_to_s( dec_Sun() )
-  end
-  alias_method :declination_string, :string_dec_Sun
-  
-  # From displays.rb  
-  # String format for delta oblique
-  def string_delta_oblique()
-    show_minutes(delta_oblique())
-  end
-  
-  # From displays.rb  
-  # String format for delta orbit
-  def string_delta_orbit()
-    show_minutes(delta_orbit())
-  end
-  
-  # From displays.rb  
-  # String format for centre
-  def string_eqc()
-    degrees_to_s( center())
-  end
-  
-  # From displays.rb  
-  # Equation of time output for minutes and seconds
-  def string_eot()    
-    eot = time_eot()  
-    min_eot = eot
-    if min_eot < 0.0
-      sign = "-"
-    else
-      sign = "+"
-    end
-    eot = min_eot.abs
-    minutes = Integer( eot )
-    seconds = ( eot - minutes ) * 60.0
-    decimal_seconds = ( seconds - Integer( seconds )) * 100.0
-    min = "%02d" % minutes
-    sec = "%02d" %  seconds
-    dec_sec = "%01d" % decimal_seconds
-    sign << min << "m, " << sec << "." << dec_sec << "s"
-  end 
-  alias_method :display_equation_of_time, :string_eot 
-  
-  # From displays.rb  
-  # String format conversion of jd to date
-  def string_jd_to_date( jd = DJ00 )
-    jd = check_jd_zero( jd )
-    Date.jd( jd ).to_s
-  end 
-  alias_method :jd_to_date_string, :string_jd_to_date 
-
-  # From displays.rb  
-  # String format of mean anomaly    
-  def string_ma_Sun()
-    degrees_to_s( @ma )
-  end
-  alias_method :mean_anomaly_string, :string_ma_Sun 
-
-  # From displays.rb  
-  # String format of right ascension
-  def string_ra_Sun()
-    degrees_to_s( ra_Sun() )
-  end
-  alias_method :right_ascension_string, :string_ra_Sun     
-
-  # From displays.rb  
-  # String format of true anomaly  
-  def string_ta_Sun( )
-    degrees_to_s( ta_Sun() )
-  end
-  alias_method :true_anomaly_string, :string_ta_Sun
-
-  # From displays.rb  
-  # String formatter for h:m:s display 
-  def string_time( dt = DT2000 )
-    dt = check_t_zero( dt )
-    
-    if dt.class == DateTime      
-      hours   = dt.hour
-      minutes = dt.min
-      seconds = dt.sec
-      intsecs = Integer( seconds )
-      decsecs = Integer(( seconds - intsecs ).round( 3 ) * 1000.0 )
-    else
-      decimal    = dt % DAY_HOURS
-      hours      = Integer( decimal )
-      mindecimal = 60.0 * ( decimal - hours )
-      minutes    = Integer( mindecimal )
-      seconds    = 60.0 * ( mindecimal - minutes )    
-      intsecs    = Integer( seconds )
-      decsecs    = Integer(( seconds - intsecs ).round( 3 ) * 1000.0 )
-    end
-    
-    "%02d" % hours   +
-                 ":" + 
-    "%02d" % minutes + 
-                 ":" + 
-    "%02d" % intsecs +
-                 "." +
-    "%3.3d" % decsecs
-  end
-  alias_method :display_time_string, :string_time    
-
-  # From displays.rb  
-  # String format of true longitude 
-  def string_tl_Sun()
-    degrees_to_s( tl_Sun() )
-  end
-  alias_method :true_longitude_string, :string_tl_Sun 
-
-  # From displays.rb  
-  # String format of true obliquity 
-  def string_to_Earth()
-    degrees_to_s( to_Earth() )
-  end
-  alias_method :true_obliquity_string, :string_to_Earth 
-
-end
-
-if __FILE__ == $PROGRAM_NAME
-
-  spec = File.expand_path('../../../tests/minitest', __FILE__)
-  $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
-  require 'displays_spec'
-  require 'aliased_displays_spec'
-
-end