equationoftime 4.1.1 → 4.1.2

data/lib/eot/init.rb

@@ -8,11 +8,11 @@ require_relative 'times'
 class Eot
       
 
-  # From init.rb:<br>
+  # From init.rb:  
   # address used for GeoLatLng.addr
   attr_accessor :addr
   
-  # From init.rb:<br>
+  # From init.rb:  
   # Astronomical Julian Day Number is an instance of DateTime class.
   # When new Equation of Time class is initialized 
   # @ajd = DateTime.now.to_time.utc.to_datetime.jd.to_f
@@ -29,53 +29,53 @@ class Eot
     @ma = Celes.falp03(@ta)
   end
 
-  # From init.rb:<br>
+  # From init.rb:  
   # Nutation Data is an instance of Array class.
   # @data = nutation_table5_3a.yaml 
   # YAML File loaded when new Eot class is initialized.
   #attr_reader :data  
 
-  # From init.rb:<br>
+  # From init.rb:  
   # @date is an instance of DateTime class.  
   # When new Eot class is initialized @date = ajd_to_datetime(@ajd)  
   attr_accessor :date
 
-  # From init.rb:<br>
+  # From init.rb:  
   # Julian Day Number is an instance of DateTime class.
   # When new Eot class is initialized @jd = @ajd
   attr_accessor :jd
 
-  # From init.rb:<br>
+  # From init.rb:  
   # Latitude input is an instance of Float class.
   # When new Eot class is initialized @latitude = 0.0
   # May use GeoLatLng class to set it also but please comment that out  
   # if no internet connection is present via proxies or firewalls :D  
   attr_accessor :latitude
   
-  # From init.rb:<br>
+  # From init.rb:  
   # Longitude input is an instance of Float class.
   # When new Eot class is initialized @longitude = 0.0
   # May use GeoLatLng class to set it also but please comment that out 
   # if no internet connection is present via proxies or firewalls :D   
   attr_accessor :longitude
 
-  # From init.rb:<br>
+  # From init.rb:  
   # Mean Anomaly gets called a lot so class attribute saves it.
   # @ma = Celes.falp03(@ta) 
   attr_accessor :ma
   
 
-  # From init.rb:<br>
+  # From init.rb:  
   # JCT gets called a lot so class attribute it.
   # Setting @ajd will set this and @ma
   # @ta = (( @ajd - DJ00 ) / DJC).to_f  
   attr_accessor :ta
       
-  # From init.rb:<br>
+  # From init.rb:  
   # Initialize to set attributes 
   # You may use GeoLatLng to set up @latitude and @longitude but you need to have
   # internet so if not please comment it out for now.    
-  def initialize(addr=nil)
+  def initialize()
 
     # file_path     = File.expand_path( File.dirname( __FILE__ ) + "/nutation_table5_3a.yaml" )
     # @data         = YAML::load( File.open( file_path, 'r'), :safe => true  ).freeze
@@ -88,11 +88,11 @@ class Eot
     @ma = Celes.falp03(@ta)                                      
     
     # comment out below if you do not have internet connection
-    geo = GeoLatLng.new()   
-    @addr = geo.addr           
+    geo = GeoLatLng.new()      
     geo.get_coordinates_from_address
-    @latitude.nil?  ? @latitude = geo.lat : @latitude
-    @longitude.nil? ? @longitude = geo.lng : @longitude
+    @addr = geo.addr
+    @latitude = geo.lat
+    @longitude = geo.lng
   end
   
 end
@@ -106,7 +106,7 @@ if __FILE__ == $PROGRAM_NAME
   require 'eot'
   eot = Eot.new
 
-  p eot.ajd
+  # p eot.ajd
   p eot.date
   p eot.jd
 
@@ -116,8 +116,9 @@ if __FILE__ == $PROGRAM_NAME
   p eot.latitude
   p eot.longitude
   
-  spec = File.expand_path('../../../tests/minitest', __FILE__)
-  $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
-  require 'init_spec'
+  #~ spec = File.expand_path('../../../tests/minitest', __FILE__)
+  #~ $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
+  #~ require 'init_spec'
+  system 'rake test'
 
 end

data/lib/eot/nutation.rb

@@ -5,7 +5,7 @@
 
 class Eot    
  
-  # From nutation.rb<br>
+  # From nutation.rb:   
   # Returns array with [ delta_eps, delta_psi, ma_sun, omega]
   # celes gem is used here now for performance. It is a Ruby wrapper for
   # see http://www.iausofa.org/

data/lib/eot/times.rb

@@ -2,35 +2,35 @@
 
 class Eot    
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Pass in an AJD number
   # Returns a DateTime object
   def ajd_to_datetime(ajd)
     DateTime.jd(ajd + 0.5)
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns EOT as an AJD Julian number
   def eot_jd
     time_eot() / DAY_MINUTES 
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd and @longitude attributes
   # Returns DateTime object of local noon or solar transit
   def local_noon_dt()
     ajd_to_datetime(@ajd - @longitude / 360.0 - eot_jd())
   end
   
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd and @longitude attributes
   # Returns DateTime object of local mean noon or solar transit
   def mean_local_noon_dt()
     ajd_to_datetime(@ajd - @longitude / 360.0)
   end
   
-  # From times.rb:<br>
+  # From times.rb:  
   # sets @ajd to DateTime.now
   # Returns EOT (equation of time) now in decimal minutes form
   def now()
@@ -39,56 +39,56 @@ class Eot
     time_eot()
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns a DateTime object of local sunrise
   def sunrise_dt()                    
     ajd_to_datetime(sunrise_jd())     
   end
 
-  # From times.rb:<br> 
+  # From times.rb:   
   # Uses @ajd attribute
   # Returns Sunrise as a Julian Day Number
   def sunrise_jd()    
     local_noon_dt().ajd - ha_Sun() * R2D / 360.0
   end 	
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns a DateTime object of local sunset
   def sunset_dt()                                                   
     ajd_to_datetime(sunset_jd())      
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns Sunset as a Julian Day Number
   def sunset_jd()
     local_noon_dt().ajd + ha_Sun() * R2D / 360.0 
   end
   
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns Oblique component of EOT in decimal minutes time  
   def time_delta_oblique()         
     (tl_Sun() - ra_Sun()) * R2D * SM        
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns Orbit component of EOT in decimal minutes time 
   def time_delta_orbit()      
     (@ma - ta_Sun()) * R2D * SM
   end  
   
-  # From times.rb:<br>
+  # From times.rb:  
   # Uses @ajd attribute
   # Returns EOT as a float for decimal minutes time
   def time_eot()
     eot() * R2D * SM	  
   end
 
-  # From times.rb:<br>
+  # From times.rb:  
   # All calculations with ( ta )  were based on this.
   # Julian Century Time is a fractional century
   # Julian Day Number DJ00 is subtracted from the JDN or AJDN and then divided
@@ -96,15 +96,15 @@ class Eot
   # Deprecated 
   def time_julian_century()
     t1 = ( @ajd - DJ00 ) / DJC
-    t2 = t1 * t1
-    t3 = t1 * t2
-    t4 = t2 * t2
-    t5 = t2 * t3
-    t6 = t3 * t3
-    t7 = t3 * t4
-    t8 = t4 * t4
-    t9 = t4 * t5
-    t10 = t5 * t5
+    # t2 = t1 * t1
+    # t3 = t1 * t2
+    # t4 = t2 * t2
+    # t5 = t2 * t3
+    # t6 = t3 * t3
+    # t7 = t3 * t4
+    # t8 = t4 * t4
+    # t9 = t4 * t5
+    # t10 = t5 * t5
 #    @ta = [ t1, t2, t3, t4, t5, t6, t7, t8, t9, t10 ] 
     @ta = t1     
   end

data/lib/eot/utilities.rb

@@ -1,35 +1,35 @@
 # utilities.rb
 
 class Eot
-  # From utilities.rb:<br>
+  # From utilities.rb:  
   # A check for default J2000
   # sets default when arg is nil
   def check_jd_nil( jd = DJ00 )
     jd.nil? ? jd = DJ00 : jd
   end
 
-  # From utilities.rb:<br>
+  # From utilities.rb:  
   # A check for default J2000
   # sets default when arg is zero
   def check_jd_zero( jd = DJ00 )
     jd == 0 ? jd = DJ00 : jd = check_jd_nil( jd )
   end
 
-  # From utilities.rb:<br>
+  # From utilities.rb:  
   # A check for default DT2000
   # sets default when arg is nil
   def check_t_nil( dt = DT2000 )
     dt.nil? ? dt = DT2000 : dt
   end
 
-  # From utilities.rb:<br>
+  # From utilities.rb:  
   # A check for default DT2000
   # sets default when arg is zero
   def check_t_zero( dt = DT2000 )
     dt == 0 ? dt = DT2000 : dt = check_t_nil( dt )
   end
 
-  # From utilities.rb:<br>
+  # From utilities.rb:  
   # Keeps large angles in range of 360.0
   # aliased by truncate
   def mod_360( x = 0.0 )

data/lib/eot/version.rb

@@ -1,7 +1,3 @@
-
-## version.rb
-
-class Eot
-  # From version.rb<br>
-  VERSION = "4.1.1"
+class Eot
+  VERSION = "4.1.2"
 end

data/test/aliased_angles_spec.rb

@@ -0,0 +1,239 @@
+# aliased_angles_spec.rb
+#
+
+require_relative 'spec_config'
+
+lib = File.expand_path('../../../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'eot'
+
+Eot_aliased_angles = Eot.new()
+
+describe 'tests ajd of 2456885.0' do
+  
+  before(:each) do  
+    Eot_aliased_angles.ajd  = 2456885.0
+    ajd = Eot_aliased_angles.ajd
+    # somtimes need date to check values somewhere else
+    Eot_aliased_angles.date = Eot_aliased_angles.ajd_to_datetime(ajd)
+  end 
+
+  it 'expected    2456885.0 for Eot_aliased_angles.ajd'do
+    assert_equal( 2456885.0, Eot_aliased_angles.ajd )    
+  end 
+  
+  it 'expected    "2014-08-15T12:00:00+00:00" for Eot_aliased_angles.date'.to_s do
+    assert_equal( "2014-08-15T12:00:00+00:00", Eot_aliased_angles.date.to_s )    
+  end 
+
+  it 'expected    3.8508003966038915 for Eot_aliased_angles.ma'do
+    assert_equal( 3.8508003966038915, Eot_aliased_angles.ma )
+  end
+
+  it 'expected    2.4887103398436143 returned by Eot_aliased_angles.apparent_longitude()? ' do
+    assert_equal( 2.4887103398436143, Eot_aliased_angles.apparent_longitude() )    
+  end 
+  
+  it 'expected    -0.7943361570447028 returned by Eot_aliased_angles.cosine_apparent_longitude()? ' do
+    assert_equal( -0.7943361570447028, Eot_aliased_angles.cosine_apparent_longitude() ) 
+  end  
+  
+  it 'expected    -0.7943772759574919 returned by Eot_aliased_angles.cosine_true_longitude()? ' do
+    assert_equal( -0.7943772759574919, Eot_aliased_angles.cosine_true_longitude() ) 
+  end
+  
+  it 'expected    0.9175115346811911 returned by Eot_aliased_angles.cosine_true_obliquity()? ' do
+    assert_equal( 0.9175115346811911, Eot_aliased_angles.cosine_true_obliquity() )
+  end
+  
+  it 'expected     0.24401410218543554 returned by Eot_aliased_angles.declination()? ' do
+     assert_equal( 0.24401410218543554, Eot_aliased_angles.declination() ) 
+  end  
+  
+  it 'expected    -0.04103082558803539 returned by Eot_aliased_angles.delta_t_ecliptic()? ' do
+    assert_equal( -0.04103082558803539, Eot_aliased_angles.delta_t_ecliptic() )
+  end
+  
+  it 'expected    0.021413249720702243 returned by Eot_aliased_angles.delta_t_elliptic()? ' do
+    assert_equal( 0.021413249720702243, Eot_aliased_angles.delta_t_elliptic() ) 
+  end   
+  
+  it 'expected    0.016702468499021204 returned by Eot_aliased_angles.eccentricity_earth_orbit()? ' do
+    assert_equal( 0.016702468499021204, Eot_aliased_angles.eccentricity_earth_orbit() )
+  end
+  
+  it 'expected    -0.021413249720702462 returned by Eot_aliased_angles.equation_of_center()? ' do
+    assert_equal( -0.021413249720702462, Eot_aliased_angles.equation_of_center() ) 	
+  end
+  
+  it 'expected    2.5101912804141424 returned by Eot_aliased_angles.geometric_mean_longitude()? ' do
+    assert_equal( 2.5101912804141424, Eot_aliased_angles.geometric_mean_longitude() )
+  end
+
+  it 'expected    1.9143229567859146 returned by Eot_aliased_angles.horizon_angle()? ' do
+    assert_equal( 1.9143229567859146, Eot_aliased_angles.horizon_angle() ) 
+  end
+
+  it 'expected    3.8508003966038915 returned by Eot_aliased_angles.mean_anomaly()? ' do
+    assert_equal( 3.8508003966038915, Eot_aliased_angles.mean_anomaly() ) 
+  end  
+  
+  it 'expected     2.510089864980358 returned by Eot_aliased_angles.mean_longitude_aries()? ' do
+    assert_equal(  2.510089864980358, Eot_aliased_angles.mean_longitude_aries() )
+  end
+  
+  it 'expected    0.40905940254265843 returned by Eot_aliased_angles.mean_obliquity()? ' do
+    assert_equal( 0.40905940254265843, Eot_aliased_angles.mean_obliquity() )
+  end
+  
+  it 'expected    0.40905940254265843 returned by Eot_aliased_angles.mean_obliquity_of_ecliptic()? ' do
+    assert_equal( 0.40905940254265843, Eot_aliased_angles.mean_obliquity_of_ecliptic() )
+  end
+  
+  it 'expected    0.40901870461547685 returned by Eot_aliased_angles.obliquity_correction()? ' do
+    assert_equal( 0.40901870461547685, Eot_aliased_angles.obliquity_correction() )
+  end
+ 
+  it 'expected    2.5297411654316497 returned by Eot_aliased_angles.right_ascension()? ' do
+    assert_equal( 2.5297411654316497, Eot_aliased_angles.right_ascension() )
+  end
+  
+  it 'expected    0.6074784519729512 returned by Eot_aliased_angles.sine_apparent_longitude()? ' do
+    assert_equal( 0.6074784519729512, Eot_aliased_angles.sine_apparent_longitude() )
+  end  
+
+  it 'expected    0.6074246812917259 returned by Eot_aliased_angles.sine_true_longitude()? ' do
+    assert_equal( 0.6074246812917259, Eot_aliased_angles.sine_true_longitude() )
+  end
+
+  it 'expected    3.8293871468831893 returned by Eot_aliased_angles.true_anomaly()? ' do
+    assert_equal( 3.8293871468831893, Eot_aliased_angles.true_anomaly() ) 
+  end  
+
+  it 'expected    2.48877803069344 returned by Eot_aliased_angles.true_longitude()? ' do
+    assert_equal( 2.48877803069344, Eot_aliased_angles.true_longitude() )
+  end
+
+  it 'expected     2.5101242776531474 returned by Eot_aliased_angles.true_longitude_aries()? ' do
+    assert_equal(  2.5101242776531474, Eot_aliased_angles.true_longitude_aries() )
+  end    
+
+  it 'expected    0.40901870461547685 returned by Eot_aliased_angles.true_obliquity()? ' do
+    assert_equal( 0.40901870461547685, Eot_aliased_angles.true_obliquity() )
+  end    
+end
+
+describe 'tests ajd of 2455055.5' do
+
+  before(:each) do
+    Eot_aliased_angles.ajd                   = 2455055.0     
+    ajd = Eot_aliased_angles.ajd
+    # check date for this ajd when needed.
+    Eot_aliased_angles.date = Eot_aliased_angles.ajd_to_datetime(ajd)    
+  end
+  
+  it 'expected    455055.0, returned by Eot_aliased_angles.' do
+    assert_equal( 2455055.0, Eot_aliased_angles.ajd )
+  end
+  
+  it 'expected    "2009-08-11T12:00:00+00:00", returned by Eot_aliased_angles.date.to_s' do
+    assert_equal( "2009-08-11T12:00:00+00:00", Eot_aliased_angles.date.to_s )
+  end
+  
+  it 'expected    3.7871218188949207, returned by Eot_aliased_angles.ma' do
+    assert_equal( 3.7871218188949207, Eot_aliased_angles.ma ) 
+  end 
+
+  it 'expected    3.7871218188949207 Eot_aliased_angles.ma  returned by Eot_angles.mean_anomaly()' do
+    assert_equal( 3.7871218188949207, Eot_aliased_angles.mean_anomaly() )
+  end   
+  
+  it 'expected    2.4252140645725033 returned by Eot_aliased_angles.apparent_longitude()? ' do
+    assert_equal( 2.4252140645725033, Eot_aliased_angles.apparent_longitude() )
+  end
+  
+  it 'expected    -0.7541886969975007 returned by Eot_aliased_angles.cosine_apparent_longitude()? ' do
+    assert_equal( -0.7541886969975007, Eot_aliased_angles.cosine_apparent_longitude() )
+  end  
+  
+  it 'expected    -0.7542060769936684 returned by Eot_aliased_angles.cosine_true_longitude()? ' do
+    assert_equal( -0.7542060769936684, Eot_aliased_angles.cosine_true_longitude() )
+  end
+  
+  it 'expected    0.9174818088112336 returned by Eot_aliased_angles.cosine_true_obliquity()? ' do
+    assert_equal( 0.9174818088112336, Eot_aliased_angles.cosine_true_obliquity() )
+  end
+  
+  it 'expected     0.26426912722944046 returned by Eot_aliased_angles.declination()? ' do
+     assert_equal( 0.26426912722944046, Eot_aliased_angles.declination() )
+  end
+  
+  it 'expected    -0.04234904897476355 returned by Eot_aliased_angles.delta_t_ecliptic()? ' do
+    assert_equal( -0.04234904897476355, Eot_aliased_angles.delta_t_ecliptic() )    	
+  end
+  
+  it 'expected    0.019768413456709766 returned by Eot_aliased_angles.delta_t_elliptic()? ' do
+    assert_equal( 0.019768413456709766, Eot_aliased_angles.delta_t_elliptic() )
+  end 
+  
+  it 'expected     0.016704576164208475 returned by Eot_aliased_angles.eccentricity_earth_orbit()? ' do
+     assert_equal( 0.016704576164208475, Eot_aliased_angles.eccentricity_earth_orbit() )
+  end
+  
+  it 'expected    -0.019768413456709915 returned by Eot_aliased_angles.equation_of_center()? ' do
+    assert_equal( -0.019768413456709915, Eot_aliased_angles.equation_of_center() )
+  end
+
+  it 'expected    2.445008945789877 returned by Eot_aliased_angles.geometric_mean_longitude()? ' do
+    assert_equal( 2.445008945789877, Eot_aliased_angles.geometric_mean_longitude() )
+  end
+
+  it 'expected    1.9434600427973594 returned by Eot_aliased_angles.horizon_angle()? ' do
+    assert_equal( 1.9434600427973594, Eot_aliased_angles.horizon_angle() ) 
+  end
+
+  it 'expected    3.7871218188949207 returned by Eot_aliased_angles.mean_anomaly()? ' do
+    assert_equal( 3.7871218188949207, Eot_aliased_angles.mean_anomaly() ) 
+  end  
+  
+  it 'expected    2.444907382260759 returned by Eot_aliased_angles.mean_longitude_aries()? ' do
+    assert_equal( 2.444907382260759, Eot_aliased_angles.mean_longitude_aries() )
+  end
+  
+  it 'expected    0.4090707793981491 returned by Eot_aliased_angles.mean_obliquity()? ' do
+    assert_equal( 0.4090707793981491, Eot_aliased_angles.mean_obliquity() )
+  end
+  
+  it 'expected    0.4090934409048494 returned by Eot_aliased_angles.obliquity_correction()? ' do
+    assert_equal( 0.4090934409048494, Eot_aliased_angles.obliquity_correction() )
+  end  
+  
+  it 'expected    2.467563113547267 returned by Eot_aliased_angles.right_ascension()? ' do
+    assert_equal( 2.467563113547267, Eot_aliased_angles.right_ascension() ) 
+  end
+  
+  it 'expected    0.6566577566139093 returned by Eot_aliased_angles.sine_apparent_longitude()? ' do
+    assert_equal( 0.6566577566139093, Eot_aliased_angles.sine_apparent_longitude() ) 
+  end  
+
+  it 'expected    0.6566377946979757 returned by Eot_aliased_angles.sine_true_longitude()? ' do
+    assert_equal( 0.6566377946979757, Eot_aliased_angles.sine_true_longitude() )
+  end
+
+  it 'expected    3.767353405438211 returned by Eot_aliased_angles.true_anomaly()? ' do
+    assert_equal( 3.767353405438211, Eot_aliased_angles.true_anomaly() )
+  end
+  
+  it 'expected    2.4252405323331674 returned by Eot_aliased_angles.true_longitude()? ' do
+    assert_equal( 2.4252405323331674, Eot_aliased_angles.true_longitude() )
+  end 
+  
+  it 'expected    2.4449774607872907 returned by Eot_aliased_angles.true_longitude_aries()? ' do
+    assert_equal( 2.4449774607872907, Eot_aliased_angles.true_longitude_aries() ) 
+  end
+
+  it 'expected    0.4090934409048494 returned by Eot_aliased_angles.true_obliquity()? ' do
+    assert_equal( 0.4090934409048494, Eot_aliased_angles.true_obliquity() )	
+  end      
+  
+end