equationoftime 3.0.0 → 4.0.0

data/examples/rbeot.c

@@ -0,0 +1,24 @@
+#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/test_eot.rb

@@ -0,0 +1,17 @@
+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

@@ -3,6 +3,7 @@
 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
@@ -18,21 +19,9 @@ $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
     # end
 # end
 
-require 'eot'
-eot = Eot.new
+my_array   = [1.0, -0.5, 3.0]
+my_array.inject(0.0) {|sum, n| p sum * 2.0 + n }
 
-@current = DateTime.now.to_time.utc
-# comparing the difference for two methods of day_fraction
-df1 = @current.to_datetime.day_fraction.to_f
-df2 = @current.hour / 24.0 + @current.min / 1_440.0 + @current.sec / 86_400.0 + @current.usec / 86_400_000_000.0
-
-ajd = @current.to_datetime.ajd.to_f
-
-# setting up to test poly_eval in each angle method using it.
-tjc = eot.time_julian_century(2456734.00110875)
-#t = (Date.today.jd - 2451545.0)/365250.0
-
-eot.eccentricity_Earth( tjc )
-eot.gml_Sun( tjc )
-p eot.ml_Aries( tjc )
-eot.mo_Earth( tjc )
+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

@@ -11,13 +11,10 @@ eot = Eot.new
 # puts eot.public_methods(false).sort
 # puts eot.nil?
 loop do
-p date = DateTime.now.to_time.utc.to_datetime
-#~ dt = Astro.delta_T(date) / 86400.0
-ajd = date.ajd
-#~ ajd1 = ajd + dt
-p date2 = DateTime.jd(ajd + 0.5)
-sleep 0.7
-puts
+  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)

data/examples/times_year.rb

@@ -0,0 +1,50 @@
+# 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

@@ -6,8 +6,18 @@ $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
 require 'eot'
 
 eot = Eot.new()
-p eot.date
-p eot.longitude
-p eot.latitude
-p eot.sunrise_dt
-p eot.sunset_dt
+# 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

@@ -6,150 +6,214 @@ $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
 require 'eot'
 
 
-file = <<EOS
+file = <<DISPLAY
 
-This will explain how the equation of time is derived by showing
-each method that comprises the formula.
+This will explain how the equation of time is derived by showing each method that comprises the formula.
 
 eot = Eot.new
-#{eot = Eot.new}
 
-Set the ajd instance to the deisired date ex : eot.ajd = Date.parse("yyyy, mm, dd").jd
+Instance of Eot class eot = #{eot = Eot.new}
 
-eot.ajd = Date.today.jd
-#{eot.ajd = Date.today.jd}
-The ajd has been set to #{eot.ajd}
 
+Set the class attribute ajd for the deisired date ex : eot.ajd = Date.parse("yyyy, mm, dd").jd.to_f
 
-The Astronomical Julian Day Number plus the fractional day time now will be used to calculate current angles.
+eot.ajd = Date.today.jd.to_f
 
-The current time expressed as an AJDN : 
-#{eot.ajd}
+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.
 
-mu is the solar Mean anomaly. It is an angle with respect to perihelion of Earth orbit around Sun.
-It is calculated as if the Earths orbit were circular.
+The Astronomical Julian Day Number is a half day less than Julian Day Number with the fractional day time included when we DateTime class.
 
-mu = eot.ma
-#{mu = eot.ma} deg.
+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.
 
-nu is the True anomaly. It 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.
 
-eqc = eot.center()
-#{eqc = eot.center()} deg.
+eot.center() = #{eot.center()} degrees.
 
-nu = eot.ta_Sun()
-#{nu = eot.ta_Sun()} deg. 
+nu = equation of center + M or mu or mean anomaly
 
-We find the delta or difference of these first two angles.
-Notice it is counter sign of equation of center. It just may have taken longer to calculate.
+The delta is M - nu.
 
-delta_eccentric = mu - nu
-#{delta_eccentric = mu - nu} deg. 
 
-The next delta is calculated as follows.
+The method inside the gem is called delta_orbit()
 
-lambda = mean longitude + equation of center
+eot.delta_orbit = #{eot.delta_orbit} degrees
 
-mean_longitude = eot.gml_Sun()
-#{mean_longitude = eot.gml_Sun()} deg.
+We could just change the sign of eot.center().  
 
-eqc = eot.center(eot.time_julian_century())
-#{eqc = eot.center(eot.time_julian_century())} deg
+In fact let's just do that now. a1 = - eot.center()
+  
+- eot.center() = #{- eot.center()}
 
-true longitude =  mean longitude + equation of center
-#{mean_longitude + eqc}
 
-Lambda is the \"ecliptic coordinate system\" angle of the sun.
+The next delta is calculated as follows.
 
-lambda_1 = eot.tl_Sun(eot.time_julian_century())
-#{lambda_1 = eot.tl_Sun(eot.time_julian_century())} deg
+lambda or apparent longitude of the Sun minus alpha or right ascension of the Sun.
 
-This may also be used and is a very small difference to lambda 1
+a2 = lambda - alpha
+
+lambda = mean longitude + equation of center
 
-lambda_2 = eot.al_Sun(eot.time_julian_century())
-#{lambda_2 = eot.al_Sun(eot.time_julian_century())} deg
+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(eot.time_julian_century())
-#{alpha = eot.ra_Sun(eot.time_julian_century())} deg
+alpha = eot.ra_Sun()
+
+alpha = #{eot.ra_Sun()}
+
+lambda - alpha = #{eot.tl_Sun - eot.ra_Sun()}
 
-Delta ecliptic is the difference in these two angles.
+The method inside the gem is called delta_oblique.
 
-delta_ecliptic = lambda_1 - alpha
-#{delta_ecliptic = lambda_1 - alpha} deg.
+eot.delta_oblique() = #{eot.delta_oblique()}
 
-The method inside the gem is called angle_delta_oblique.
 
-The sum of these two delta angles is
+The sum of these two delta angles is the equation of time angle in degrees not time.
 
-delta_total = delta_eccentric + delta_ecliptic
-#{delta_total = delta_eccentric + delta_ecliptic} deg.
+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 #{my_longitude = -88.74467}
-#{my_longitude} deg. note: -deg. because I'm west of UTC or 0.0 deg.
 
-So my mean solar transit UTC is
+my longitude is -88.74467 deg. so I'll set the longitude attribute in Eot class.
 
-mean_noon = 12 - my_longitude / 15.0
-#{mean_noon = 12 - my_longitude / 15.0} hr.
+eot.longitude = #{eot.longitude = -88.74467}
 
-using the equation of time library method get_time(Date.today.to_datetime + mean_noon) 
-I'm showing what time mean noon for my_longitude normally is without
-the equation of time considered and this is true only four days a year.
 
-Be sure to set your longitude and latitude as attributes.
-eot.longitude = my_longitude
-eot.latitude = my_latitude
-#{eot.longitude = my_longitude}
+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.mean_local_noon_dt()
-#{eot.mean_local_noon_dt()}
+eot.longitude + eot.eot() = #{eot.longitude} + #{eot.eot()}
 
-Now we need to calculate the delta total for our mean solar transit.
-We will use the mean noon time to get equation of time.
-We will be using mean_noon as a time object.
-First create a time object for the date at 12:00 and subtract 
-(my_longitude / 15.0) * 3600.0 to get the time object.
-This is what get_time did without the longitude because it was provided before.
+What's the time of solar transit here now?
 
-mean_noon = Time.utc(DateTime.now.to_time.utc.to_datetime.year, DateTime.now.to_time.utc.to_datetime.month, DateTime.now.to_time.utc.to_datetime.day, 12, 0, 0) - (my_longitude / 15.0) * 3600.0}
-#{mean_noon = Time.utc(DateTime.now.to_time.utc.to_datetime.year, DateTime.now.to_time.utc.to_datetime.month, DateTime.now.to_time.utc.to_datetime.day, 12, 0, 0) - (my_longitude / 15.0) * 3600.0}
+12.0 - (eot.longitude + eot.eot()) / 15.0 = #{12.0 - (eot.longitude + eot.eot()) / 15.0}
+  
+That looks about right.
 
-delta_total = eot.eot()
-#{delta_total = eot.eot()} min
 
-Let us convert that to an angle.
 
-delta_deg = delta_total / 4.0 * 360.0 / 360.98564736629  Note: not quite 4 minutes per degree.
-#{delta_deg = delta_total / 4.0 * 360.0 / 360.98564736629} deg.
+Let's make a time out of that by converting it to a fractional day.
 
-Now we can add that from our longitude or mean noon angle
+(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.
 
-noon_angle = my_longitude + delta_deg
-#{noon_angle = my_longitude + delta_deg}
+ajd = Date.today.ajd.to_f
 
-Your longitude and the GHA of the Sun will match atTrue Solar Transit time.
+ajd = #{Date.today.ajd.to_f}
 
-So lets convert that angle to a time now.
+#{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}
 
-12 - noon_angle / 15.0
-#{12 - noon_angle / 15.0}
 
-my_transit = eot.ajd + (12 - noon_angle / 15.0) / 24.0
-#{my_transit = eot.ajd + (12 - noon_angle / 15.0) / 24.0}
+There is a method in Eot class to make that a DateTime object.
 
-DateTime.jd(my_transit)
-#{DateTime.jd(my_transit)}
+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.
 
-check it here http://douglasallen.github.io/planets/ just put in the correct time and date then hit calculate.
 Does GHA Sun = your longitude?
-EOS
+
+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