equationoftime 3.0.0

data/lib/eot/displays.rb

@@ -0,0 +1,213 @@
+# displays.rb
+
+class Eot
+
+  # From displays.rb<br>
+  # String formatter for d:m:s display 
+  def degrees_to_s( degrees = 0.0 )
+    degrees.nil? ? degrees = 0.0 : degrees
+    degrees < 0 ? sign_string      = "-" : sign_string = "+"    
+                  absolute_degrees = degrees.abs
+          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 = bd( 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<br>
+  # 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<br>
+  # String for time now
+  def show_now(now = now(Time.now.utc))
+    show_minutes(now)
+  end  
+  
+  # From displays.rb<br>
+  # String format of apparent longitude 
+  def string_al_Sun( ta = A2000 )
+    degrees_to_s( al_Sun( ta ) )
+  end
+  alias_method :apparent_longitude_string, :string_al_Sun
+
+  # From displays.rb<br>
+  # 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<br>
+  # String format of declination 
+  def string_dec_Sun( ta = A2000 )
+    degrees_to_s( dec_Sun( ta ) )
+  end
+  alias_method :declination_string, :string_dec_Sun
+  
+  # From displays.rb<br>
+  # String format for delta oblique
+  def string_delta_oblique( ta = A2000 )
+    show_minutes(delta_oblique( ta ))
+  end
+  
+  # From displays.rb<br>
+  # String format for delta orbit
+  def string_delta_orbit( ta = A2000 )
+    show_minutes(delta_orbit( ta ))
+  end
+  
+  # From displays.rb<br>
+  # String format for centre
+  def string_eqc( ta = A2000 )
+    degrees_to_s( center( ta ))
+  end
+  
+  # From displays.rb<br>
+  # 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<br>
+  # String format conversion of jd to date
+  def string_jd_to_date( jd = J2000 )
+    jd = check_jd_zero( jd )
+    Date.jd( jd ).to_s
+  end 
+  alias_method :jd_to_date_string, :string_jd_to_date 
+
+  # From displays.rb<br>
+  # 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<br>
+  # String format of right ascension
+  def string_ra_Sun( ta = A2000 )
+    degrees_to_s( ra_Sun( ta ) )
+  end
+  alias_method :right_ascension_string, :string_ra_Sun     
+
+  # From displays.rb<br>
+  # String format of true anomaly  
+  def string_ta_Sun( ta = A2000 )
+    degrees_to_s( ta_Sun( ta ) )
+  end
+  alias_method :true_anomaly_string, :string_ta_Sun
+
+  # From displays.rb<br>
+  # 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 = bd( 60.0 * ( decimal - hours )) * 1.0
+      minutes = Integer( mindecimal )
+      seconds = bd( 60.0 * ( mindecimal - minutes )) * 1.0    
+      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<br>
+  # String format of true longitude 
+  def string_tl_Sun( ta = A2000 )
+    degrees_to_s( tl_Sun( ta ) )
+  end
+  alias_method :true_longitude_string, :string_tl_Sun 
+
+  # From displays.rb<br>
+  # String format of true obliquity 
+  def string_to_Earth( ta = A2000 )
+    degrees_to_s( to_Earth( ta ) )
+  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

data/lib/eot/geo_lat_lng_smt.rb

@@ -0,0 +1,80 @@
+# geo_lat_lng_smt.rb
+
+require 'multi_xml'
+require 'rest-client'
+
+  # class for location lookup 
+  # in geo_lat_lng_smt.rb 
+class GeoLatLng
+
+  # Base address for Google maps api
+  attr_reader :base
+
+  # Default US set to PARCS
+  attr_reader :default_us
+  
+  # Default International set to GMT Museum
+  attr_reader :default_int
+  
+  # Address entered
+  attr_accessor :addr
+  
+  # Latitude returned
+  attr_accessor :lat
+  
+  # Longitude returned
+  attr_accessor :lng
+
+  # Base Google URL
+  attr_reader :base
+
+  # Instance variables
+  def initialize
+       
+    @base           = "http://maps.googleapis.com/maps/api/geocode/xml?sensor=false&address="
+    @default_us     = "3333 Coyote Hill Road, Palo Alto, CA, 94304, USA"#do you copy? :D
+    @default_int    = "Blackheath Ave, London SE10 8XJ, UK"
+    @lat            = 0.0
+    @lng            = 0.0
+    MultiXml.parser = :rexml#:libxml#:ox # :nokogiri
+  
+  end   
+
+  # set address  
+  def addr=(addr = @default_int)
+    @addr = addr
+  end  
+
+  # coordinates lookup   
+  def get_coordinates_from_address       
+    res        = RestClient.get(
+                                URI.encode(
+                                           "#{ @base }#{ @addr }"
+                                          )
+                               )    
+    parsed_res = MultiXml.parse( res )
+    result     = parsed_res[ "GeocodeResponse" ][ "result" ]
+    status     = parsed_res[ "GeocodeResponse" ][ "status" ] 
+    if status != "OK"      
+      @default      
+    else      
+      if result.count != 4          
+        @default          
+      else
+        @lat   = parsed_res[ "GeocodeResponse" ][ "result" ][ "geometry" ][ "location" ][ "lat" ]
+        @lng   = parsed_res[ "GeocodeResponse" ][ "result" ][ "geometry" ][ "location" ][ "lng" ]
+      end
+    end
+    
+  end
+
+end
+
+if __FILE__ == $PROGRAM_NAME
+
+  spec = File.expand_path('../../../tests/minitest', __FILE__)
+  $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
+  require 'geo_spec'
+  
+end
+

data/lib/eot/init.rb

@@ -0,0 +1,93 @@
+# init.rb
+
+class Eot    
+
+  # From init.rb:<br>
+  # Astronomical Julian Day Number is an instance of Date class.
+  # When new Equation of Time class is initialized @ajd = jd = today DateTime.jd
+  # Used for getting the equation of time now if it is set for that. There is 
+  # always a slight delay in the computation though.
+  attr_accessor :ajd
+
+  # From init.rb:<br>
+  # 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>
+  # @date is an instance of Date class.  
+  # When new Eot class is initialized @date = today  
+  attr_accessor :date
+
+  # From init.rb:<br>
+  # Julian Day Number is an instance of Date class.
+  # When new Eot class is initialized @jd = jd today  
+  attr_accessor :jd
+
+  # From init.rb:<br>
+  # 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 is commented out or will fail if no 
+  # internet connection is present.  
+  attr_accessor :latitude
+  
+  # From init.rb:<br>
+  # 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 is commented out or will fail if no 
+  # internet connection is present.  
+  attr_accessor :longitude
+
+  # From init.rb:<br>
+  # Mean Anomaly gets called a lot so attribute accessor saves time 
+  attr_accessor :ma
+
+  # From init.rb:<br>
+  # address used for GeoLatLng.addr when used.(commented out) 
+  attr_accessor :addr
+
+      
+  # From init.rb:<br>
+  # Initialize loads nutation data with safe_yaml and is frozen, atrributes are set.
+  # You may use GeoLatLng to set up @latitude and @longitude but you need to have
+  # internet so that is commented out for now.    
+  def initialize(addr=nil)
+
+    file_path     = File.expand_path( File.dirname( __FILE__ ) + "/nutation_table5_3a.yaml" )
+    @data         = YAML::load( File.open( file_path, 'r'), :safe => true  ).freeze
+ 
+    @ajd.nil? ? @ajd = DateTime.now.to_time.utc.to_datetime.jd.to_f : @ajd = self.ajd      
+    @date.nil? ? @date = DateTime.now.to_time.utc.to_date : @date = self.date
+    @jd.nil? ? @jd  = DateTime.now.to_time.utc.to_datetime.jd.to_f : @jd = self.jd
+    @latitude.nil? ? @latitude = 0.0 : @latitude = self.latitude
+    @longitude.nil? ? @longitude = 0.0 : @longitude = self.longitude    
+    @ma.nil? ? @ma = ma_Sun() : @ma
+
+    #geo = GeoLatLng.new
+    @addr = addr
+    #addr.nil? ? geo.addr=(geo.default_int) : geo.addr=addr
+# uncomment below if you have internet connection
+    #geo.get_coordinates_from_address
+    #@latitude = geo.lat.to_f
+    #@longitude = geo.lng.to_f
+  end
+  
+end
+
+
+# we can run some tests from inside this file.
+if __FILE__ == $PROGRAM_NAME
+
+  lib = File.expand_path('../../../lib', __FILE__)
+  puts "Loading gem from #{lib}/eot.rb"
+  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+  require 'eot'
+  e = Eot.new
+  p e.date
+
+  spec = File.expand_path('../../../tests/minitest', __FILE__)
+  $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
+  require 'init_spec'
+
+end

data/lib/eot/nutation.rb

@@ -0,0 +1,70 @@
+# nutation.rb
+
+require 'mutex_m'
+
+class Eot    
+ 
+  # From nutation.rb<br>
+  # nutation data terms are used here.
+  # arc seconds are used for these terms
+  # see http://aa.usno.navy.mil/publications/docs/Circular_179.pdf page 46 (5.19) 
+  def delta_equinox( ta = A2000 )
+    ta = check_jct_zero( ta )    
+    
+    ma_moon = [-0.00024470, 0.051635, 31.8792, 1717915923.2178, 485868.249036].inject(0.0) {|p, a| p * ta[0] + a}		
+
+    ma_sun  = [-0.00001149, 0.000136, -0.5532, 129596581.0481, 1287104.79305].inject(0.0) {|p, a| p * ta[0] + a}	
+                    
+    md_moon = [0.00000417, -0.001037, -12.7512, 1739527262.8478, 335779.526232].inject(0.0) {|p, a| p * ta[0] + a}
+            
+    me_moon = [-0.00003169, 0.006593, -6.3706, 1602961601.2090, 1072260.70369].inject(0.0) {|p, a| p * ta[0] + a} 
+           
+    omega   = [-0.00005939, 0.007702, 7.4722, -6962890.5431, 450160.398036].inject(0.0) {|p, a| p * ta[0] + a}            
+    
+    # declare and clear these two variables for the sigma loop
+    delta_psi, delta_eps = 0, 0
+
+    lines = data.size - 1
+    for i in 0..lines
+      fma_sun    = data[i][0].to_i
+      fma_moon   = data[i][1].to_i  	
+      fmd_moon   = data[i][2].to_i
+      fme_moon   = data[i][3].to_i  
+      fomega     = data[i][4].to_i
+
+      sine       = sind(fma_moon * ma_moon +
+                        fma_sun  * ma_sun  +
+                        fmd_moon * md_moon +
+                        fme_moon * me_moon +
+                        fomega   * omega)
+          
+      cosine     = cosd(fma_moon * ma_moon +
+                        fma_sun  * ma_sun  +
+                        fmd_moon * md_moon +
+                        fme_moon * me_moon +
+                        fomega   * omega)
+          
+      delta_psi += (data[i][6].to_f                  + 
+                    data[i][7].to_f  * ta[0]) * sine +
+                    data[i][10].to_f * cosine
+          
+      delta_eps += (data[i][8].to_f                     + 
+                    data[i][9].to_f   * ta[0]) * cosine +				 
+                    data[i][12].to_f  * sine						
+                         
+    end
+
+    delta_eps  = to_deg( delta_eps ) / 1000.0
+    delta_psi  = to_deg( delta_psi ) / 1000.0
+
+    [ delta_eps, delta_psi, ma_sun, omega, ma_moon, md_moon, me_moon]
+  end
+
+end
+if __FILE__ == $PROGRAM_NAME
+
+  spec = File.expand_path('../../../tests/minitest', __FILE__)
+  $LOAD_PATH.unshift(spec) unless $LOAD_PATH.include?(spec)
+  require 'nutation_spec'
+
+end