equationoftime 4.1.1 → 4.1.2

data/equationoftime.gemspec

@@ -1,7 +1,8 @@
 # coding: utf-8
-lib = File.expand_path('../lib', __FILE__)
+
+lib = File.expand_path('./lib', __FILE__)
 $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-require 'eot/version'
+require 'eot'
 
 Gem::Specification.new do |spec|
      # Metadata
@@ -21,20 +22,27 @@ Gem::Specification.new do |spec|
      spec.files         = Dir.glob("ext/**/*.{c,h,rb, so}") + Dir.glob("**") +
                           Dir.glob("lib/**/*.rb") + Dir.glob(".*")
      #spec.files         = `git ls-files -z`.split($/)#split("\x0")
-     spec.extensions    << "ext/ceot/extconf.rb"
+     spec.extensions    << "ext/eot/extconf.rb"
      spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
      spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
      spec.require_paths = ["lib"]
      spec.platform      = Gem::Platform::RUBY
+# note: the last good build was with these development dependencies below
      spec.add_development_dependency "bundler", "~> 1.5"
-     spec.add_development_dependency "rake"
-     spec.add_development_dependency "rake-compiler"
-     spec.add_development_dependency "minitest"	 
-     spec.add_development_dependency "rspec"
-     spec.add_development_dependency "yard"
+     #spec.add_development_dependency "rake",, "~> 10.3.2"
+     #spec. add_development_dependency 'rake-compiler', '~> 0.9.3'
+     #spec.add_development_dependency "hoe", "~> 3.12.0"
+     #spec.add_development_dependency "guard", "~> 2.6.1"
+     #spec.add_development_dependency "guard-minitest", "~> 2.3.2"
+     #spec.add_development_dependency "minitest", "~> 5.4.1"
+     #spec.add_development_dependency "ZenTest", "~> 4.10.1" 
+     #spec.add_development_dependency "rspec", "~> 3.1.0"
+     #spec.add_development_dependency "yard", "~> 0.8.7.4"
+     #spec.add_development_dependency "rdoc", "~> 4.1.2"
+
      spec.add_runtime_dependency "celes"
      #spec.add_runtime_dependency "astro-algo"
-     spec.add_runtime_dependency "multi_xml"
+     #spec.add_runtime_dependency "multi_xml"
      spec.add_runtime_dependency "rest-client"
      #spec.add_runtime_dependency "safe_yaml"	 
      

data/examples/analemma_data_generator.rb

@@ -0,0 +1,58 @@
+# analemma_data_generator.rb
+
+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'
+require 'safe_yaml'
+require 'time'
+
+eot    = Eot.new
+
+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  = start_time.to_datetime.jd
+finish_jd = finish_time.to_datetime.jd
+
+fstr = "%b %d"
+
+@data, @group, @group_id = ["2014", eot.addr, eot.latitude, eot.longitude], {}, 1
+
+(start_jd..finish_jd).each do |jd|     
+  date    = Date.jd(jd + 0.5).strftime(fstr)
+  
+  eot.ajd = jd  
+  delta_t = eot.show_minutes(eot.time_eot())
+  trans   = eot.display_time_string(12 - eot.time_eot()/60.0)
+  e1      = eot.show_minutes(eot.time_delta_orbit())
+  e2      = eot.show_minutes(eot.time_delta_oblique())
+  decline = eot.degrees_to_s(eot.dec_Sun())
+  @group  = { "id"     => "#{@group_id}", 
+              "date"   => date,
+              "julian" => "#{jd}",
+                "rise" => eot.sunrise_dt().to_json, 
+                 "eot" => delta_t,                   
+             "transit" => eot.local_noon_dt().to_json,
+         "declination" => "#{decline}",
+                 "set" => eot.sunset_dt().to_json             
+             }  
+  @data << @group
+  @group_id  += 1
+end
+
+file_path = "analemma_data.yml"
+File::open( file_path, "w" ) do |f|
+  YAML.dump( @data, f )
+end
+
+puts "File analemma_data.yml processed"

data/examples/check_date_type.rb

@@ -0,0 +1,60 @@
+# check_date_type.rb
+
+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()
+
+output =<<EOO
+
+eot = Eot.new()
+
+eot.ajd = #{eot.ajd}
+
+eot.eot = #{eot.eot()} # this is radians for eot today.
+
+eot.date = #{eot.date = "2014-10-14"} = #{eot.date.class}
+
+eot.ajd = Date.parse(eot.date).jd.to_f = #{eot.ajd = Date.parse(eot.date).jd.to_f}
+
+eot.eot = #{eot.eot} # getting default today
+
+eot.equation_of_time(eot.time_julian_century(eot.ajd)) = #{eot.eot()}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()} with lat and lng 0.
+
+eot.longitude = #{eot.longitude = 0}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+
+eot.latitude = #{eot.latitude = 0}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+
+eot.longitude = #{eot.longitude = 180.0}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+
+eot.longitude = #{eot.longitude = -180.0}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+
+eot.longitude = #{eot.longitude = 0}
+eot.latitude = #{eot.latitude = 90}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+
+eot.longitude = #{eot.longitude = 0}
+eot.latitude = #{eot.latitude = -90}
+
+eot.sunrise_dt() = #{eot.sunrise_dt()}
+EOO
+
+puts output

data/examples/compare_geoc_long_ra.rb

@@ -0,0 +1,44 @@
+# compare_geoc_long_ra.rb
+
+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()
+
+start  = "2014-6-21"
+finish = "2014-6-21"
+
+start_date  = Date.parse(start)
+finish_date = Date.parse(finish)
+
+puts "start julian date  = #{start_date.jd} date  = #{start}"
+puts "finish julian date = #{finish_date.jd} date = #{finish}"
+puts
+
+ajd = start_date.ajd * 1.0
+
+for psec in 630..631
+  h       = 10 / 24.0
+  fd      = h + 45 / 1440.0 + 31 / 86400.0 + 691 / 86400.0 / 1e3 + 244 / 86400.0 / 1e6 + 989 / 86400.0 / 1e9 + psec / 86400.0 / 1e12 
+  eot.ajd = ajd + fd
+  ml      = eot.geometric_mean_longitude()
+  tl      = eot.true_longitude()
+  da      = ml - tl
+  al      = eot.apparent_longitude()
+  ra      = eot.right_ascension()
+  db      = al - ra
+#  puts "time\t\t\tmean longitude\t\t\ttrue_longitude\t\t\tdifference"
+#  puts "10:45:31.691_#{usec}\t\t\t#{ml}\t\t#{tl}\t\t#{da}"
+#  puts
+  puts "time\t\t\tapparent longitude\t\tright ascension\t\t\tdifference"
+	puts "10:45:31.691_244_989_#{psec}\t\t\t#{al}\t\t#{ra}\t\t#{db}"
+	puts
+end
+#~ end

data/examples/data_table_for_astro_dog.rb

@@ -0,0 +1,45 @@
+# class that goes with equationoftime.herokuapp.com/analemma
+class AnalemmaDataTable
+
+  attr_accessor  :data, :finish, :start, :table, :html
+
+  def initialize 
+    @start        = Time.utc( 2014, "jan", 1, 12, 0, 0 ).to_s
+    @start_jd     = Date.parse(@start).jd.to_s
+    @finish       = Time.utc( 2014, "dec", 31, 12, 0, 0 ).to_s
+    @finish_jd    = Date.parse(@finish).jd.to_s
+    @span         = Date.parse( @finish ).jd - Date.parse( @start ).jd
+    @file_path    = File.expand_path( File.dirname( __FILE__ ) + "/analemma_data.yml" )
+    @data         = YAML::load( File.open( @file_path, 'r'), :safe => true ).freeze   
+    @table        = ""
+    @html         = ""
+
+    (0..@span).each do |i|
+      p @data[i]
+#      @jd          = @data[i].fetch ":jd"
+#      @date        = @data[i].fetch ":date"
+#      @delta_1     = @data[i].fetch ":delta_1"
+#      @delta_2     = @data[i].fetch ":delta_2"
+#      @delta_t     = @data[i].fetch ":delta_t"
+#      @declination = @data[i].fetch ":declination"
+#      @transit     = @data[i].fetch ":transit"
+    end      
+  end
+end
+
+if __FILE__ == $PROGRAM_NAME
+  
+  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 'eot'
+  
+  adt = AnalemmaDataTable.new
+
+end

data/examples/earth_rotation.rb

@@ -0,0 +1,42 @@
+# earth_rotation.rb
+
+require 'bigdecimal'
+
+factor = 1 / (1 - 1 / 366.0)
+bd_factor = BigDecimal.new("#{1 / (1 - 1 / 366.0)}") * 1.0
+puts "360 degrees x #{factor} = #{360 * factor} degrees"
+puts "360 degrees x #{bd_factor} = #{360 * bd_factor} degrees"
+
+sidereal_minutes = 4 / factor
+bd_sidereal_minutes = 4 / bd_factor
+puts sidereal_minutes
+puts sidereal_minutes
+
+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()
+# From angles.rb:<br>
+# angle factor for daily sidereal time (experimental)  
+def factor eot
+  eot.ajd  =        Date.parse("2000-01-01").jd    
+  tlaa  =           eot.tl_Aries()
+  eot.ajd  =        eot.ajd + 1
+  tlab  =           eot.tl_Aries()
+  dif   =           (tlab - tlaa) * Eot::R2D 
+  f1    =           dif / 360.0 + 1 
+  1 / f1
+end
+
+p factor(eot) * 4
+p factor(eot) * 4 * 360.0
+p factor(eot) * 4 * 360.0 / 60.0
+p factor(eot) * 4 * 360.0 / 60.0 * 366
+shd = factor(eot) * 4 * 360.0 / 60.0
+sdnl = shd * 366 * 4
+p sdnl / (4 * 365 + 1)

data/examples/eot_methods_list.rb

@@ -0,0 +1,48 @@
+# 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
+
+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
+    }
+    
+date = DateTime.now.to_time.utc.to_datetime
+puts date
+ 
+eot = Eot.new
+puts eot.nil?
+puts eot.public_methods(false).count
+list = eot.public_methods(false).sort
+sym_list = []
+list.each {|i|sym_list << i.to_sym}
+#p sym_list
+my_list =[:al_Sun, :angle_delta_oblique, :angle_delta_orbit, :angle_delta_psi, :angle_equation_of_time, 
+          :apparent_longitude, :apparent_longitude_string, :center, :cosine_al_Sun, :cosine_apparent_longitude, 
+          :cosine_tl_Sun, :cosine_to_Earth, :cosine_true_longitude, :cosine_true_obliquity,  :dec_Sun, 
+          :declination, :declination_string, :degrees_to_s, :delta_epsilon, :delta_equinox, :delta_oblique, :delta_orbit,
+          :delta_psi, :delta_t_ecliptic, :delta_t_elliptic, :display_equation_of_time, :display_time_string,
+          :eccentricity_Earth, :eccentricity_earth_orbit, :ecliptic_longitude, :eot, :eot_jd, :eq_of_equinox, 
+          :equation_of_center, :geometric_mean_longitude, :gml_Sun, :ha_Sun, :horizon_angle,
+          :jd_to_date_string, :julian_period_day_fraction_to_time, :latitude, :latitude=, :local_noon_dt, 
+          :ma_Sun, :mean_anomaly, :mean_anomaly_string, :mean_local_noon_dt, :mean_longitude_aries, 
+          :mean_obliquity, :mean_obliquity_of_ecliptic, :ml_Aries, :mo_Earth, :mod_360, :now,
+          :obliquity_correction, :omega, :ra_Sun, :right_ascension, :right_ascension_string, :show_minutes, 
+          :show_now, :sine_al_Sun, :sine_apparent_longitude, :sine_tl_Sun, :sine_to_Earth, 
+          :sine_true_longitude, :string_al_Sun, :string_day_fraction_to_time, :string_dec_Sun, 
+          :string_delta_oblique, :string_delta_orbit, :string_eot, :string_eqc, :string_jd_to_date, 
+          :string_ma_Sun, :string_ra_Sun, :string_ta_Sun, :string_time, :string_tl_Sun, :string_to_Earth, 
+          :sunrise_dt, :sunrise_jd, :sunset_dt, :sunset_jd, :ta_Sun, :time_delta_oblique, :time_delta_orbit, 
+          :time_eot, :time_julian_centurey, :time_julian_century, :tl_Aries, :tl_Sun, :to_Earth, 
+          :true_anomaly, :true_anomaly_string, :true_longitude, :true_longitude_aries, 
+          :true_longitude_string, :true_obliquity, :true_obliquity_string, :truncate]
+my_list.collect {|e| p e}
+
+puts date
+

data/examples/eot_plot.r

@@ -0,0 +1,57 @@
+
+R version 3.0.0 (2013-04-03) -- "Masked Marvel"
+Copyright (C) 2013 The R Foundation for Statistical Computing
+Platform: i386-w64-mingw32/i386 (32-bit)
+
+R is free software and comes with ABSOLUTELY NO WARRANTY.
+You are welcome to redistribute it under certain conditions.
+Type 'license()' or 'licence()' for distribution details.
+
+R is a collaborative project with many contributors.
+Type 'contributors()' for more information and
+'citation()' on how to cite R or R packages in publications.
+
+Type 'demo()' for some demos, 'help()' for on-line help, or
+'help.start()' for an HTML browser interface to help.
+Type 'q()' to quit R.
+
+#formulas, equations and notation from http://homeweb2.unifr.ch/hungerbu/pub/sonnenuhr/sundial.ps
+# and by the paper of C. Blatter (http://dz-srv1.sub.uni-goettingen.de/sub/digbib/loader?did=D241775)
+ 
+epsilon=23.45*2*pi/360
+alpha  =78.5 *2*pi/360
+kappa  =0.016722
+
+mu<-function(t,epsilon,alpha,kappa) {
+zaehl= -sin(delta(t,kappa))+tan(epsilon/2)^2*sin(2*(t-alpha)+delta(t,kappa))
+nenn =  cos(delta(t,kappa))+tan(epsilon/2)^2*cos(2*(t-alpha)+delta(t,kappa))
+return( atan(zaehl/nenn) )
+}
+
+# approximation by Blatter
+psi<-function(t) {
+   return(t+delta(t))
+ }
+ 
+ delta<-function(t,kappa) {
+   return( 2*sin(t)*kappa+5/4*sin(2*t)*kappa^2 )
+ }
+ 
+ png(filename="Equation of time.jpg", width=1024, height=768, pointsize=12)
+ par(bg="whitesmoke")
+ time=seq(-0.6,2*pi+0.4,length=365+(0.6+0.4)/(2*pi)*365)
+ da<-paste("1.",as.character(seq(1,12,by=2)),".07",sep="")
+ da<-as.Date(c(da,"1.1.08"),"%d.%m.%y")
+ #umrechungsfaktor: 24h=2*pi
+ plot(time,mu(time,epsilon,alpha,kappa)*24*60/(2*pi),type="l",col="red",xlab="Tag",ylab="Time in Minutes",main="Equation of time: true solar time  -  mean solar time",cex.main=1.5,lwd=3,xaxt="n",xaxs="i",xlim=range(time))
+ lines(time,mu(time,epsilon=0,alpha,kappa)*24*60/(2*pi),lwd=2,col="navy",lty="dotdash")
+ lines(time,mu(time,epsilon,alpha,kappa=0)*24*60/(2*pi),lwd=2,col="purple",lty="dashed")
+ #grid
+ abline(h=c(-15,-10,-5,5,10,15),lty="dotted",col="grey")
+ abline(h=0,lty="solid",col="grey33")
+ tage=cumsum(c(0,31,28,31,30,31,30,31,31,30,31,30,31))/365-3/365
+ abline(v=tage*2*pi,lty="dotted",col="grey")
+ axis(1,at=seq(0,2*pi,length=7)-3*2*pi/365,labels=as.character(da,"1. %b"))
+ legend(x=5.5,y=-10.9,legend=c("Equation of time","Ecliptic","Earth Orbit"),col=c("red","navy","purple"),lwd=c(3,2,2),lty=c("solid","dotdash","dashed"))
+ dev.off()
+

data/examples/eot_suntimes.rb

@@ -0,0 +1,149 @@
+#!/usr/bin/env ruby
+# eot_suntimes.rb leftover from building the parts in the gem so just left it here.
+#
+begin
+  require 'eot'
+rescue LoadError
+  lib = File.expand_path('../../lib', __FILE__)
+  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+  require 'eot' 
+end
+
+#def sunrise_jd(arg)
+#
+#	arg.date.nil? ? result = "date parameter not set" : result = "date parameter OK"
+#	
+#	arg.longitude.nil? ? result = "longitude parameter not set" : result = "longitude parameter OK"
+#	
+#	arg.latitude.nil? ? result = "latitude parameter not set" : result = "latitude parameter OK"	
+#	
+#	utc_noon_jd = arg.date.jd
+#	
+#	lng_jd = arg.longitude / 360.0
+#	
+#	mean_utc_noon_jd = utc_noon_jd - lng_jd
+#	
+#	tjc_mean_local_noon_array = arg.time_julian_century(mean_utc_noon_jd)
+#	
+#	eot = arg.equation_of_time(tjc_mean_local_noon_array) 
+#	
+#	eot_jd = eot / 1440.0
+#	
+#	true_utc_noon_jd = mean_utc_noon_jd - eot_jd
+#	
+#	tjc_true_local_noon_array = arg.time_julian_century(true_utc_noon_jd)
+#	
+#	mean_anomaly = arg.ma_Sun(tjc_true_local_noon_array)
+#	
+#	lha = arg.ha_Sun(tjc_true_local_noon_array, mean_anomaly)
+#	
+#	lha_jd = lha / 360.0
+#	
+#	# rise_jd
+#	true_utc_noon_jd - lha_jd	
+#	
+#end
+#
+#def sunset_jd(arg)
+#
+#  arg.date.nil? ? result = "date parameter not set" : result = "date parameter OK"
+#	
+#  arg.longitude.nil? ? result = "longitude parameter not set" : result = "longitude parameter OK"
+#	
+#  arg.latitude.nil? ? result = "latitude parameter not set" : result = "latitude parameter OK"  
+#  
+#  utc_noon_jd = arg.date.jd
+#  
+#  lng_jd = arg.longitude / 360.0
+#  
+#  mean_utc_noon_jd = utc_noon_jd - lng_jd
+#  
+#  tjc_mean_local_noon = arg.time_julian_century(mean_utc_noon_jd)
+#  
+#  eot_jd = arg.equation_of_time(tjc_mean_local_noon) / 1440.0
+#  
+#  true_utc_noon_jd = mean_utc_noon_jd - eot_jd
+#  
+#  tjc_true_local_noon = arg.time_julian_century(true_utc_noon_jd)
+#  
+#  mean_anomaly = arg.ma_Sun(tjc_true_local_noon)
+#  
+#  lha_jd = arg.ha_Sun(tjc_true_local_noon, mean_anomaly) / 360.0
+#  
+#  # sunset jd
+#  true_utc_noon_jd  + lha_jd
+#  
+#end
+
+
+require 'eot'
+
+addr = "8000 South Michigan Ave., Chicago, IL"
+#loc = GeoLatLng.new
+#loc.addr = addr
+#loc.get_coordinates_from_address
+#puts loc.lat, loc.lng
+
+eot = Eot.new(addr)
+#eot.longitude = loc.lng
+#eot.ajd = Date.today.jd
+#eot.latitude = loc.lat
+
+#puts "Using this files methods"
+#puts sunrise_jd(Date.today)
+
+## Note: DateTime.jd() renders time from midnight not noon
+##       so 12 hours need to be added for correct time.
+#puts DateTime.jd(sunrise_jd() + 0.5)
+#puts DateTime.jd(sunrise_jd() + 0.5).to_time.utc
+#
+#puts
+#puts sunset_jd()
+#
+#puts DateTime.jd(sunset_jd() + 0.5)
+#puts DateTime.jd(sunset_jd() + 0.5).to_time.utc
+
+OUT =<<EOS
+
+Using eot gem methods
+
+#{eot.latitude}
+#{eot.longitude}
+
+#{eot.sunrise_jd()}
+#{eot.ajd_to_datetime(eot.sunrise_jd())}
+#{eot.ajd_to_datetime(eot.sunrise_jd()).to_time}
+
+#{eot.sunset_jd()}
+#{eot.ajd_to_datetime(eot.sunset_jd())}
+#{eot.ajd_to_datetime(eot.sunset_jd()).to_time}
+
+#{Date.jd(Date.today.jd)}
+#{Date.today.jd}
+
+#{DateTime.jd(Date.today.jd)}
+#{DateTime.jd(Date.today.jd + 0.5)}
+
+#{DateTime.jd(Date.today.jd).to_time.utc}
+#{DateTime.jd(Date.today.jd + 0.5).to_time.utc}
+
+#{Time.now.to_datetime}
+#{Time.now.to_datetime.jd}
+
+#{Time.now.utc.to_datetime}
+#{Time.now.utc.to_date.jd}
+
+#{Date.jd(Time.now.to_date.jd)}
+#{Date.jd(Time.now.utc.to_date.jd)}
+
+#{DateTime.jd(Time.now.to_datetime.jd)}
+#{DateTime.jd(Time.now.utc.to_datetime.jd)}
+
+#{DateTime.jd(Time.now.to_datetime.jd).to_time}
+#{DateTime.jd(Time.now.utc.to_datetime.jd).to_time}
+
+#{DateTime.jd(Time.now.to_datetime.jd + 0.5).to_time}
+#{DateTime.jd(Time.now.utc.to_datetime.jd + 0.5).to_time}
+EOS
+
+puts OUT