equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/zpv.c

@@ -1,79 +0,0 @@
-#include "sofam.h"
-
-void iauZpv(double pv[2][3])
-/*
-**  - - - - - - -
-**   i a u Z p v
-**  - - - - - - -
-**
-**  Zero a pv-vector.
-**
-**  Status:  vector/matrix support function.
-**
-**  Returned:
-**     pv       double[2][3]      pv-vector
-**
-**  Called:
-**     iauZp        zero p-vector
-**
-**  This revision:  2008 May 11
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   iauZp(pv[0]);
-   iauZp(pv[1]);
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/zr.c

@@ -1,83 +0,0 @@
-#include "sofam.h"
-
-void iauZr(double r[3][3])
-/*
-**  - - - - - -
-**   i a u Z r
-**  - - - - - -
-**
-**  Initialize an r-matrix to the null matrix.
-**
-**  Status:  vector/matrix support function.
-**
-**  Returned:
-**     r        double[3][3]    r-matrix
-**
-**  This revision:  2012 April 3
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   r[0][0] = 0.0;
-   r[0][1] = 0.0;
-   r[0][2] = 0.0;
-   r[1][0] = 0.0;
-   r[1][1] = 0.0;
-   r[1][2] = 0.0;
-   r[2][0] = 0.0;
-   r[2][1] = 0.0;
-   r[2][2] = 0.0;
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/check_date_type.rb

@@ -1,55 +0,0 @@
-# check_date_type.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-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

@@ -1,39 +0,0 @@
-# compare_geoc_long_ra.rb
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-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

@@ -1,40 +0,0 @@
-# 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
-  
-  lib = File.expand_path('../../lib', __FILE__)
-  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-  
-  require 'safe_yaml'
-  require 'eot'
-  
-  adt = AnalemmaDataTable.new
-
-end

data/examples/earth_rotation.rb

@@ -1,38 +0,0 @@
-# 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
-
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-require 'eot'
-
-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.c

@@ -1,15 +0,0 @@
-#include "eot.h"
-double maSun(double t)
-{
-	   double a;
-
-
-	/* Mean anomaly of the Sun (IERS Conventions 2003). */
-	   a = fmod(         1287104.793048 +
-	             t * (      129596581.0481   +
-	             t * (                    - 0.5532   +
-	             t * (                    0.000136 +
-	             t * (       - 0.00001149 ) ) ) ), 1296000.0 ) * 4.848136811095359935899141e-6;
-
-	   return a;
-}

data/examples/eot.h

@@ -1,2 +0,0 @@
-#include "math.h"
-double maSun(double t);

data/examples/eot_methods_list.rb

@@ -1,43 +0,0 @@
-# eot_methods_list.rb
-#
-lib = File.expand_path('../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-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

@@ -1,57 +0,0 @@
-
-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()
-