equationoftime 3.0.0 → 4.0.0

data/examples/celes_parts/zpv.c

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

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

@@ -11,11 +11,9 @@ output =<<EOO
 
 eot = Eot.new()
 
-eot.date = #{eot.date} = #{eot.date.class}
+eot.ajd = #{eot.ajd}
 
-eot.ajd = eot.date.ajd.to_f = #{eot.ajd = eot.date.ajd.to_f}
-
-eot.eot = #{eot.eot()} # getting default today
+eot.eot = #{eot.eot()} # this is radians for eot today.
 
 eot.date = #{eot.date = "2014-10-14"} = #{eot.date.class}
 

data/examples/compare_geoc_long_ra.rb

@@ -7,25 +7,33 @@ require 'eot'
 
 eot = Eot.new()
 
-start  = "2013-3-20"
-finish = "2013-3-20"
+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 min in 0..59
-  h      = 10 / 24.0
-  fd     = h + min / 1440.0
-  timejc = eot.time_julian_century(ajd + fd) 
-  ml     = eot.geometric_mean_longitude(timejc)
-  tl     = eot.true_longitude(timejc)
-  ra     = eot.right_ascension(timejc)
-  puts "time\t\t\tmean longitude\t\tright ascension"
-	puts "10:#{min}\t\t\t#{ml}\t\t#{ra}\t\t#{}"
+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,40 @@
+# 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

@@ -2,12 +2,37 @@
 
 require 'bigdecimal'
 
-factor = 1 / (1 - 1 / 365.25)
-bd_factor = BigDecimal.new("#{1 / (1 - 1 / 365.25)}") * 1.0
+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
+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

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

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

data/examples/eot_methods_list.rb

@@ -7,10 +7,37 @@ 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.public_methods(false).sort
-# puts eot.nil?
+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}
 
-date = DateTime.now.to_time.utc.to_datetime
 puts date
 

data/examples/eot_suntimes.rb

@@ -1,116 +1,117 @@
 #!/usr/bin/env ruby
-# eot_suntimes.rb
+# eot_suntimes.rb leftover from building the parts in the gem so just left it here.
 #
 lib = File.expand_path('../../lib', __FILE__)
 $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
 
-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
+#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
-eot.longitude = loc.lng
-eot.ajd = Date.today.jd
-eot.latitude = loc.lat
-
-puts "Using this files methods"
-puts sunrise_jd()
-
-# 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
+#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.date}
-#{eot.date.jd}
+
+#{eot.latitude}
+#{eot.longitude}
 
 #{eot.sunrise_jd()}
-#{DateTime.jd(eot.sunrise_jd() + 0.5)}
-#{DateTime.jd(eot.sunrise_jd() + 0.5).to_time}
+#{eot.ajd_to_datetime(eot.sunrise_jd())}
+#{eot.ajd_to_datetime(eot.sunrise_jd()).to_time}
 
 #{eot.sunset_jd()}
-#{DateTime.jd(eot.sunset_jd() + 0.5)}
-#{DateTime.jd(eot.sunset_jd() + 0.5).to_time}
+#{eot.ajd_to_datetime(eot.sunset_jd())}
+#{eot.ajd_to_datetime(eot.sunset_jd()).to_time}
 
 #{Date.jd(Date.today.jd)}
 #{Date.today.jd}
@@ -135,6 +136,9 @@ Using eot gem methods
 
 #{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