mk_greenwich 0.1.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 577c15044a40fd368c0dc1922bdcfad86ce2adbc
+  data.tar.gz: ab54e375e1c5a4a0aaafebdbb96bc0fbc317057c
+SHA512:
+  metadata.gz: 630c3cfe10dc9758910c0a625adcdb7042023479fc617931aae646df771172e4977fa9d30bad5ec51f74299afca6dab5ef175e6336be4c3323d2ae8fd47c101d
+  data.tar.gz: 5148bcaff4f812c041e3419752b9522de2e262a41defc3429e04b4e7b24be1649377f0e8f691f019f4ea741e74d3235380769bbc915f21457466b6511bcd7525

data/.gitignore

@@ -0,0 +1,9 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/.rspec

@@ -0,0 +1,2 @@
+--format documentation
+--color

data/.travis.yml

@@ -0,0 +1,5 @@
+sudo: false
+language: ruby
+rvm:
+  - 2.3.1
+before_install: gem install bundler -v 1.12.5

data/Gemfile

@@ -0,0 +1,10 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in mk_greenwich.gemspec
+gemspec
+
+group :development do
+  gem "guard"
+  gem "guard-rspec", "~> 4.7.0"
+end
+

data/Guardfile

@@ -0,0 +1,70 @@
+# A sample Guardfile
+# More info at https://github.com/guard/guard#readme
+
+## Uncomment and set this to only include directories you want to watch
+# directories %w(app lib config test spec features) \
+#  .select{|d| Dir.exists?(d) ? d : UI.warning("Directory #{d} does not exist")}
+
+## Note: if you are using the `directories` clause above and you are not
+## watching the project directory ('.'), then you will want to move
+## the Guardfile to a watched dir and symlink it back, e.g.
+#
+#  $ mkdir config
+#  $ mv Guardfile config/
+#  $ ln -s config/Guardfile .
+#
+# and, you'll have to watch "config/Guardfile" instead of "Guardfile"
+
+# Note: The cmd option is now required due to the increasing number of ways
+#       rspec may be run, below are examples of the most common uses.
+#  * bundler: 'bundle exec rspec'
+#  * bundler binstubs: 'bin/rspec'
+#  * spring: 'bin/rspec' (This will use spring if running and you have
+#                          installed the spring binstubs per the docs)
+#  * zeus: 'zeus rspec' (requires the server to be started separately)
+#  * 'just' rspec: 'rspec'
+
+guard :rspec, cmd: "bundle exec rspec" do
+  require "guard/rspec/dsl"
+  dsl = Guard::RSpec::Dsl.new(self)
+
+  # Feel free to open issues for suggestions and improvements
+
+  # RSpec files
+  rspec = dsl.rspec
+  watch(rspec.spec_helper) { rspec.spec_dir }
+  watch(rspec.spec_support) { rspec.spec_dir }
+  watch(rspec.spec_files)
+
+  # Ruby files
+  ruby = dsl.ruby
+  dsl.watch_spec_files_for(ruby.lib_files)
+
+  ## Rails files
+  #rails = dsl.rails(view_extensions: %w(erb haml slim))
+  #dsl.watch_spec_files_for(rails.app_files)
+  #dsl.watch_spec_files_for(rails.views)
+
+  #watch(rails.controllers) do |m|
+  #  [
+  #    rspec.spec.call("routing/#{m[1]}_routing"),
+  #    rspec.spec.call("controllers/#{m[1]}_controller"),
+  #    rspec.spec.call("acceptance/#{m[1]}")
+  #  ]
+  #end
+
+  ## Rails config changes
+  #watch(rails.spec_helper)     { rspec.spec_dir }
+  #watch(rails.routes)          { "#{rspec.spec_dir}/routing" }
+  #watch(rails.app_controller)  { "#{rspec.spec_dir}/controllers" }
+
+  ## Capybara features specs
+  #watch(rails.view_dirs)     { |m| rspec.spec.call("features/#{m[1]}") }
+  #watch(rails.layouts)       { |m| rspec.spec.call("features/#{m[1]}") }
+
+  ## Turnip features and steps
+  #watch(%r{^spec/acceptance/(.+)\.feature$})
+  #watch(%r{^spec/acceptance/steps/(.+)_steps\.rb$}) do |m|
+  #  Dir[File.join("**/#{m[1]}.feature")][0] || "spec/acceptance"
+  #end
+end

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2016 Masaru Koizumi
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,64 @@
+# MkGreenwich
+
+## Installation
+
+This is the gem library which calculates Greenwich time.
+
+```ruby
+gem 'mk_greenwich'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mk_greenwich
+
+## Usage
+
+### Instantiation
+
+    require 'mk_greenwich'
+    
+    g = MkGreenwich.new("20160906")
+    g = MkGreenwich.new("20160906123456")
+    g = MkGreenwich.new
+
+* You can set UTC formatted "YYYYMMDD" or "YYYYMMDDHHMMSS" as an argument.
+* If you don't set an argument, this class considers the system time to have been set as an argument.
+
+### Calculation
+
+    puts " UTC: #{g.utc.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "  TT: #{g.tt.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts " UT1: #{g.ut1.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts " TDB: #{g.tdb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts " ERA: #{g.era}"
+    puts "  EO: #{g.eo}"
+    puts "GAST: #{g.gast} rad"
+    puts "    = #{g.gast_deg} °"
+    puts "    = #{g.gast_hms}"
+    puts "GMST: #{g.gmst} rad"
+    puts "    = #{g.gmst_deg} °"
+    puts "    = #{g.gmst_hms}"
+    puts "  EE: #{g.ee} rad"
+    puts "    = #{g.ee_deg} °"
+    puts "    = #{g.ee_hms}"
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/komasaru/mk_greenwich.
+
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](http://opensource.org/licenses/MIT).
+

data/Rakefile

@@ -0,0 +1,6 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "mk_greenwich"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/exe/mk_greenwich

@@ -0,0 +1,21 @@
+#!/usr/bin/env ruby
+
+require "mk_greenwich"
+
+g = MkGreenwich.new(ARGV[0])  # argument should be UT
+puts " UTC: #{g.utc.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "  TT: #{g.tt.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts " UT1: #{g.ut1.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts " TDB: #{g.tdb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts " ERA: #{g.era}"
+puts "  EO: #{g.eo}"
+puts "GAST: #{g.gast} rad"
+puts "    = #{g.gast_deg} °"
+puts "    = #{g.gast_hms}"
+puts "GMST: #{g.gmst} rad"
+puts "    = #{g.gmst_deg} °"
+puts "    = #{g.gmst_hms}"
+puts "  EE: #{g.ee} rad"
+puts "    = #{g.ee_deg} °"
+puts "    = #{g.ee_hms}"
+

data/lib/mk_greenwich/argument.rb

@@ -0,0 +1,24 @@
+module MkGreenwich
+  class Argument
+    def initialize(arg)
+      @arg = arg
+    end
+
+    #=========================================================================
+    # 引数取得
+    #
+    # * コマンドライン引数を取得して日時の妥当性チェックを行う
+    # * コマンドライン引数無指定なら、現在日時とする。
+    #
+    # @return: utc  (Time Object)
+    #=========================================================================
+    def get_utc
+      (puts Const::MSG_ERR_1; return) unless @arg =~ /^\d{8}$|^\d{14}$/
+      year, month, day = @arg[ 0, 4].to_i, @arg[ 4, 2].to_i, @arg[ 6, 2].to_i
+      hour, min,   sec = @arg[ 8, 2].to_i, @arg[10, 2].to_i, @arg[12, 2].to_i
+      (puts Const::MSG_ERR_2; return) unless Date.valid_date?(year, month, day)
+      (puts Const::MSG_ERR_2; return) if hour > 23 || min > 59 || sec > 59
+      return Time.new(year, month, day, hour, min, sec, "+00:00")
+    end
+  end
+end

data/lib/mk_greenwich/cip_cio.rb

@@ -0,0 +1,92 @@
+require "mk_greenwich/fundamental_argument"
+
+module MkGreenwich
+  class CipCio
+    include MkGreenwich::FundamentalArgument
+
+    #-------------------------------------------------------------------------
+    # Initialization
+    #
+    # @param:  t  (Julian Century)
+    #-------------------------------------------------------------------------
+    def initialize(t)
+      @t = t
+    end
+
+    #-------------------------------------------------------------------------
+    # Extract from the bias-precession-nutation matrix the X,Y coordinates
+    # of the Celestial Intermediate Pole.
+    #
+    # @param:  r       (Rotation Matrix)
+    # @return: [x, y]  (x, y cordinates of CIP)
+    #-------------------------------------------------------------------------
+    def bpn2xy(r)
+      return [r[2][0], r[2][1]]
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # The CIO locator s, positioning the Celestial Intermediate Origin on
+    # the equator of the Celestial Intermediate Pole, given the CIP's X,Y
+    # coordinates.  Compatible with IAU 2006/2000A precession-nutation.
+    #
+    # @param:  x  (x coordinate of CIP)
+    # @param:  y  (y coordinate of CIP)
+    # @return: s  (CIO locator (Unit: rad))
+    #-------------------------------------------------------------------------
+    def calc_s_06(x, y)
+      # Fundamental Arguments (from IERS Conventions 2003)
+      fa = [
+        # Mean anomaly of the Moon.(Ref: iauFal03(t))
+        calc_l_iers2003(@t),
+        # Mean anomaly of the Sun.(Ref: iauFalp03(t))
+        calc_p_iers2003(@t),
+        # Mean longitude of the Moon minus that of the ascending node.(Ref: iauFaf03(t))
+        calc_f_iers2003(@t),
+        # Mean elongation of the Moon from the Sun.(Ref: iauFad03(t))
+        calc_d_iers2003(@t),
+        # Mean longitude of the ascending node of the Moon.(Ref: iauFaom03(t))
+        calc_om_iers2003(@t),
+        # Mean longitude of Venus.(Ref: iauFave03(t))
+        calc_ve_iers2003(@t),
+        # Mean longitude of Earth.(Ref: iauFae03(t))
+        calc_ea_iers2003(@t),
+        # General precession in longitude.(Ref: iauFapa03(t))
+        calc_pa_iers2003(@t)
+      ]
+      # Evaluate s.
+      w_0, w_1, w_2, w_3, w_4, w_5 = Const::SP
+      (Const::S_0.size - 1).downto(0) do |i|
+        a = 0.0
+        0.upto(7) { |j| a += Const::S_0[i][0][j] * fa[j] }
+        w_0 += Const::S_0[i][1] * Math.sin(a) + Const::S_0[i][2] * Math.cos(a)
+      end
+      (Const::S_1.size - 1).downto(0) do |i|
+        a = 0.0
+        0.upto(7) { |j| a += Const::S_1[i][0][j] * fa[j] }
+        w_1 += Const::S_1[i][1] * Math.sin(a) + Const::S_1[i][2] * Math.cos(a)
+      end
+      (Const::S_2.size - 1).downto(0) do |i|
+        a = 0.0
+        0.upto(7) { |j| a += Const::S_2[i][0][j] * fa[j] }
+        w_2 += Const::S_2[i][1] * Math.sin(a) + Const::S_2[i][2] * Math.cos(a)
+      end
+      (Const::S_3.size - 1).downto(0) do |i|
+        a = 0.0
+        0.upto(7) { |j| a += Const::S_3[i][0][j] * fa[j] }
+        w_3 += Const::S_3[i][1] * Math.sin(a) + Const::S_3[i][2] * Math.cos(a)
+      end
+      (Const::S_4.size - 1).downto(0) do |i|
+        a = 0.0
+        0.upto(7) { |j| a += Const::S_4[i][0][j] * fa[j] }
+        w_4 += Const::S_4[i][1] * Math.sin(a) + Const::S_4[i][2] * Math.cos(a)
+      end
+      return (w_0 + (w_1 + (w_2 + (w_3 + (w_4 + w_5 \
+             * @t) * @t) * @t) * @t) * @t) * Const::AS2R - x * y / 2.0
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_greenwich/const.rb

@@ -0,0 +1,97 @@
+module MkGreenwich
+  module Const
+    MSG_ERR_1 = "[ERROR] Format: YYYYMMDD or YYYYMMDDHHMMSS"
+    MSG_ERR_2 = "[ERROR] Invalid date-time!"
+    J2000   = 2451545.0  # Reference epoch (J2000.0), Julian Date
+    JC      = 36525.0    # Days per Julian century
+    # ==== for CipCio ====
+    # Polynomial coefficients
+    SP = [94.00e-6, 3808.65e-6, -122.68e-6, -72574.11e-6, 27.98e-6, 15.62e-6]
+    # Terms of order t^0
+    S_0 = [
+      [[0,  0,  0,  0,  1,  0,  0,  0], -2640.73e-6,  0.39e-6],
+      [[0,  0,  0,  0,  2,  0,  0,  0],   -63.53e-6,  0.02e-6],
+      [[0,  0,  2, -2,  3,  0,  0,  0],   -11.75e-6, -0.01e-6],
+      [[0,  0,  2, -2,  1,  0,  0,  0],   -11.21e-6, -0.01e-6],
+      [[0,  0,  2, -2,  2,  0,  0,  0],     4.57e-6,  0.00e-6],
+      [[0,  0,  2,  0,  3,  0,  0,  0],    -2.02e-6,  0.00e-6],
+      [[0,  0,  2,  0,  1,  0,  0,  0],    -1.98e-6,  0.00e-6],
+      [[0,  0,  0,  0,  3,  0,  0,  0],     1.72e-6,  0.00e-6],
+      [[0,  1,  0,  0,  1,  0,  0,  0],     1.41e-6,  0.01e-6],
+      [[0,  1,  0,  0, -1,  0,  0,  0],     1.26e-6,  0.01e-6],
+      [[1,  0,  0,  0, -1,  0,  0,  0],     0.63e-6,  0.00e-6],
+      [[1,  0,  0,  0,  1,  0,  0,  0],     0.63e-6,  0.00e-6],
+      [[0,  1,  2, -2,  3,  0,  0,  0],    -0.46e-6,  0.00e-6],
+      [[0,  1,  2, -2,  1,  0,  0,  0],    -0.45e-6,  0.00e-6],
+      [[0,  0,  4, -4,  4,  0,  0,  0],    -0.36e-6,  0.00e-6],
+      [[0,  0,  1, -1,  1, -8, 12,  0],     0.24e-6,  0.12e-6],
+      [[0,  0,  2,  0,  0,  0,  0,  0],    -0.32e-6,  0.00e-6],
+      [[0,  0,  2,  0,  2,  0,  0,  0],    -0.28e-6,  0.00e-6],
+      [[1,  0,  2,  0,  3,  0,  0,  0],    -0.27e-6,  0.00e-6],
+      [[1,  0,  2,  0,  1,  0,  0,  0],    -0.26e-6,  0.00e-6],
+      [[0,  0,  2, -2,  0,  0,  0,  0],     0.21e-6,  0.00e-6],
+      [[0,  1, -2,  2, -3,  0,  0,  0],    -0.19e-6,  0.00e-6],
+      [[0,  1, -2,  2, -1,  0,  0,  0],    -0.18e-6,  0.00e-6],
+      [[0,  0,  0,  0,  0,  8,-13, -1],     0.10e-6, -0.05e-6],
+      [[0,  0,  0,  2,  0,  0,  0,  0],    -0.15e-6,  0.00e-6],
+      [[2,  0, -2,  0, -1,  0,  0,  0],     0.14e-6,  0.00e-6],
+      [[0,  1,  2, -2,  2,  0,  0,  0],     0.14e-6,  0.00e-6],
+      [[1,  0,  0, -2,  1,  0,  0,  0],    -0.14e-6,  0.00e-6],
+      [[1,  0,  0, -2, -1,  0,  0,  0],    -0.14e-6,  0.00e-6],
+      [[0,  0,  4, -2,  4,  0,  0,  0],    -0.13e-6,  0.00e-6],
+      [[0,  0,  2, -2,  4,  0,  0,  0],     0.11e-6,  0.00e-6],
+      [[1,  0, -2,  0, -3,  0,  0,  0],    -0.11e-6,  0.00e-6],
+      [[1,  0, -2,  0, -1,  0,  0,  0],    -0.11e-6,  0.00e-6]
+    ]
+    # Terms of order t^1
+    S_1 = [
+      [[0,  0,  0,  0,  2,  0,  0,  0], -0.07e-6,  3.57e-6],
+      [[0,  0,  0,  0,  1,  0,  0,  0],  1.73e-6, -0.03e-6],
+      [[0,  0,  2, -2,  3,  0,  0,  0],  0.00e-6,  0.48e-6]
+    ]
+    # Terms of order t^2
+    S_2 = [
+      [[0,  0,  0,  0,  1,  0,  0,  0], 743.52e-6, -0.17e-6],
+      [[0,  0,  2, -2,  2,  0,  0,  0],  56.91e-6,  0.06e-6],
+      [[0,  0,  2,  0,  2,  0,  0,  0],   9.84e-6, -0.01e-6],
+      [[0,  0,  0,  0,  2,  0,  0,  0],  -8.85e-6,  0.01e-6],
+      [[0,  1,  0,  0,  0,  0,  0,  0],  -6.38e-6, -0.05e-6],
+      [[1,  0,  0,  0,  0,  0,  0,  0],  -3.07e-6,  0.00e-6],
+      [[0,  1,  2, -2,  2,  0,  0,  0],   2.23e-6,  0.00e-6],
+      [[0,  0,  2,  0,  1,  0,  0,  0],   1.67e-6,  0.00e-6],
+      [[1,  0,  2,  0,  2,  0,  0,  0],   1.30e-6,  0.00e-6],
+      [[0,  1, -2,  2, -2,  0,  0,  0],   0.93e-6,  0.00e-6],
+      [[1,  0,  0, -2,  0,  0,  0,  0],   0.68e-6,  0.00e-6],
+      [[0,  0,  2, -2,  1,  0,  0,  0],  -0.55e-6,  0.00e-6],
+      [[1,  0, -2,  0, -2,  0,  0,  0],   0.53e-6,  0.00e-6],
+      [[0,  0,  0,  2,  0,  0,  0,  0],  -0.27e-6,  0.00e-6],
+      [[1,  0,  0,  0,  1,  0,  0,  0],  -0.27e-6,  0.00e-6],
+      [[1,  0, -2, -2, -2,  0,  0,  0],  -0.26e-6,  0.00e-6],
+      [[1,  0,  0,  0, -1,  0,  0,  0],  -0.25e-6,  0.00e-6],
+      [[1,  0,  2,  0,  1,  0,  0,  0],   0.22e-6,  0.00e-6],
+      [[2,  0,  0, -2,  0,  0,  0,  0],  -0.21e-6,  0.00e-6],
+      [[2,  0, -2,  0, -1,  0,  0,  0],   0.20e-6,  0.00e-6],
+      [[0,  0,  2,  2,  2,  0,  0,  0],   0.17e-6,  0.00e-6],
+      [[2,  0,  2,  0,  2,  0,  0,  0],   0.13e-6,  0.00e-6],
+      [[2,  0,  0,  0,  0,  0,  0,  0],  -0.13e-6,  0.00e-6],
+      [[1,  0,  2, -2,  2,  0,  0,  0],  -0.12e-6,  0.00e-6],
+      [[0,  0,  2,  0,  0,  0,  0,  0],  -0.11e-6,  0.00e-6]
+    ]
+    # Terms of order t^3
+    S_3 = [
+      [[0,  0,  0,  0,  1,  0,  0,  0],  0.30e-6, -23.42e-6],
+      [[0,  0,  2, -2,  2,  0,  0,  0], -0.03e-6,  -1.46e-6],
+      [[0,  0,  2,  0,  2,  0,  0,  0], -0.01e-6,  -0.25e-6],
+      [[0,  0,  0,  0,  2,  0,  0,  0],  0.00e-6,   0.23e-6]
+    ]
+    # Terms of order t^4
+    S_4 = [
+      [[0,  0,  0,  0,  1,  0,  0,  0], -0.26e-6, -0.01e-6]
+    ]
+    AS2R    = 4.848136811095359935899141e-6  # Arcseconds to radians
+    PI      = 3.141592653589793238462643     # PI
+    PI2     = 6.283185307179586476925287     # 2 * PI
+    PI_180  = 0.017453292519943295           # PI / 180
+    TURNAS  = 1296000.0                      # Arcseconds in a full circle
+  end
+end

data/lib/mk_greenwich/era_eors.rb

@@ -0,0 +1,116 @@
+module MkGreenwich
+  #===========================================================================
+  # Class for
+  #   ERA(Earth rotation angle (IAU 2000 model), 地球回転角),
+  #   EORS(Equation of the origins, 原点差)
+  #   GMST(Greenwich mean sidereal time, グリニッジ平均恒星時)
+  #   GAST(Greenwich apparent sidereal time, グリニッジ視恒星時)
+  #   EE(Equation of Equinoxes, 分点均差)
+  #===========================================================================
+  class EraEors
+    #-------------------------------------------------------------------------
+    # Initialization
+    #
+    # @param:  jd  (Julian Day)
+    #-------------------------------------------------------------------------
+    def initialize(jd)
+      @jd, @t = jd, jd - Const::J2000  # JD, JD2000.0
+    end
+
+    #-------------------------------------------------------------------------
+    # Earth rotation angle (IAU 2000 model).
+    #
+    # @param:  <none>
+    # @return: ERA  (Earth rotation angle (Unit: rad, Range: 0-2pi), 地球回転角)
+    #-------------------------------------------------------------------------
+    def calc_era(jd, t)
+      # Fractional part of T (days).
+      f = jd % 1.0
+      # Earth rotation angle at this UT1.
+      return norm_angle((f + 0.7790572732640 + 0.00273781191135448 * t) * Const::PI2)
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Equation of the origins, given the classical NPB matrix and the
+    # quantity s.
+    #
+    # @param:  r   (Rotation matrix)
+    # @param:  s   (CIO locator)
+    # @return: EO  (Equation of the origin (Unit: rad), 原点差)
+    #-------------------------------------------------------------------------
+    def calc_eo(r_mtx, s)
+      x = r_mtx[2][0]
+      ax = x / (1.0 + r_mtx[2][2])
+      xs = 1.0 - ax * x
+      ys = -ax * r_mtx[2][1]
+      zs = -x
+      p = r_mtx[0][0] * xs + r_mtx[0][1] * ys + r_mtx[0][2] * zs
+      q = r_mtx[1][0] * xs + r_mtx[1][1] * ys + r_mtx[1][2] * zs
+      return (p != 0 || q != 0) ? s - Math.atan2(q, p) : s
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Greenwich apparent sidereal time
+    #
+    # @param:  era   (Earth rotation angle)
+    # @param:  eo    (Equation of the origin)
+    # @return: GAST  (Greenwich apparent sidereal time (Unit: rad), グリニッジ視恒星時)
+    #-------------------------------------------------------------------------
+    def calc_gast(era, eo)
+      return norm_angle(era - eo)
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Greenwich mean sidereal time, IAU 2006.
+    #
+    # @param:  gast  (Greenwich apparent sidereal time, グリニッジ視恒星時)
+    # @param:  t     (Julian Century)
+    # @return: GMST  (Greenwich mean sidereal time (Unit: rad), グリニッジ平均恒星時)
+    #-------------------------------------------------------------------------
+    def calc_gmst(gast, t)
+      return norm_angle(gast +
+        (   0.014506    + \
+        (4612.156534    + \
+        (   1.3915817   + \
+        (  -0.00000044  + \
+        (  -0.000029956 + \
+        (  -0.0000000368) \
+        * t) * t) * t) * t) * t) * Const::AS2R)
+    rescue => e
+      raise
+    end
+    #-------------------------------------------------------------------------
+    # Equation of Equinoxes
+    #
+    # @param:  gast  (Greenwich apparent sidereal time, グリニッジ視恒星時)
+    # @param:  gmst  (Greenwich mean sidereal time, グリニッジ平均恒星時)
+    # @return: EE    (Equation of Equinoxes (Unit: rad), 分点均差)
+    #-------------------------------------------------------------------------
+    def calc_ee(gast, gmst)
+      return gast - gmst
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Normalize angle into the range 0 <= a < 2pi.
+    #
+    # @param:  angle  (Before normalized)
+    # @return: angle  (Normalized angle)
+    #-------------------------------------------------------------------------
+    def norm_angle(angle)
+      while angle < 0;          angle += Const::PI2; end
+      while angle > Const::PI2; angle -= Const::PI2; end
+      return angle
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_greenwich/fundamental_argument.rb

@@ -0,0 +1,130 @@
+module MkGreenwich
+  #===========================================================================
+  # Module for Fundamental arguments
+  #===========================================================================
+  module FundamentalArgument
+    module_function
+
+    #-------------------------------------------------------------------------
+    # Mean anomaly of the Moon (IERS 2003)
+    #
+    # @param:  t  (Julian Centry)
+    # @return: l  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_l_iers2003(t)
+      return ((    485868.249036  + \
+              (1717915923.2178    + \
+              (        31.8792    + \
+              (         0.051635  + \
+              (        -0.00024470) \
+              * t) * t) * t) * t) % Const::TURNAS) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Mean anomaly of the Sun (IERS 2003)
+    #
+    # @param:  t  (Julian Centry)
+    # @return: p  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_p_iers2003(t)
+      return ((   1287104.793048  + \
+              ( 129596581.0481    + \
+              (       - 0.5532    + \
+              (         0.000136  + \
+              (       - 0.00001149) \
+              * t) * t) * t) * t) % Const::TURNAS) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Mean longitude of the Moon minus that of the ascending node (IERS 2003)
+    #
+    # @param:  t  (Julian Centry)
+    # @return: f  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_f_iers2003(t)
+      return ((     335779.526232 + \
+              (1739527262.8478    + \
+              (       -12.7512    + \
+              (        -0.001037  + \
+              (         0.00000417) \
+              * t) * t) * t) * t) % Const::TURNAS) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # mean elongation of the Moon from the Sun (IERS 2003)
+    #
+    # @param:  t  (Julian Centry)
+    # @return: d  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_d_iers2003(t)
+      return ((   1072260.703692  + \
+              (1602961601.2090    + \
+              (       - 6.3706    + \
+              (         0.006593  + \
+              (       - 0.00003169) \
+              * t) * t) * t) * t) % Const::TURNAS ) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Mean longitude of the ascending node of the Moon (IERS 2003)
+    #
+    # @param:  t   (Julian Centry)
+    # @return: om  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_om_iers2003(t)
+     return ((    450160.398036  + \
+             (  -6962890.5431    + \
+             (         7.4722    + \
+             (         0.007702  + \
+             (        -0.00005939) \
+             * t) * t) * t) * t) % Const::TURNAS) * Const::AS2R
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Venus longitudes (IERS 2003)
+    #
+    # @param:  t   (Julian Centry)
+    # @return: ve  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_ve_iers2003(t)
+      return (3.176146697 + 1021.3285546211 * t) % Const::PI2
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # Earth longitudes (IERS 2003)
+    #
+    # @param:  t   (Julian Centry)
+    # @return: ea  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_ea_iers2003(t)
+      return (1.753470314 + 628.3075849991 * t) % Const::PI2
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # General accumulated precession in longitude (IERS 2003)
+    #
+    # @param:  t   (Julian Centry)
+    # @return: pa  (Unit: rad)
+    #-------------------------------------------------------------------------
+    def calc_pa_iers2003(t)
+      return (0.024381750 + 0.00000538691 * t) * t
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_greenwich/greenwich.rb

@@ -0,0 +1,213 @@
+module MkGreenwich
+  class Greenwich
+    attr_reader :utc, :tt, :ut1, :tdb
+
+    def initialize(utc)
+      @utc = utc                                              # 協定世界時
+      t_utc = MkTime.new(@utc.strftime("%Y%m%d%H%M%S"))
+      @tt  = t_utc.tt                                         # 地球時(for UTC)
+      @ut1 = t_utc.ut1                                        # 世界時1(for TT)
+      @tdb = t_utc.tdb                                        # 太陽系力学時(for TT)
+      @jd  = t_utc.jd                                         # ユリウス日(for TT)
+      @jc = calc_jc(@jd)                                      # ユリウス世紀数(for TT)
+      @jd_ut1 = MkTime.new(@ut1.strftime("%Y%m%d%H%M%S")).jd  # ユリウス日(for UT1)
+      @t = @jd_ut1 - Const::J2000
+      bpn = EphBpn.new(@tdb.strftime("%Y%m%d%H%M%S"))
+      @r_mtx = prod_mtx(bpn.r_bias_prec, bpn.r_nut)
+      cc = CipCio.new(@jc)
+      x, y = cc.bpn2xy(@r_mtx)
+      @s = cc.calc_s_06(x, y)
+      @e = EraEors.new(@jd_ut1)
+    end
+
+    #=========================================================================
+    # Earth rotation angle
+    #
+    # * IAU 2000 model
+    #
+    # @param:  <none>
+    # @return: @era
+    #=========================================================================
+    def era
+      @era = @e.calc_era(@jd, @t)
+      return @era
+    end
+
+    #=========================================================================
+    # Equation of the origins, given the classical NPB matrix and the
+    # quantity s.
+    #
+    # @param:  <none>
+    # @return: @eo
+    #=========================================================================
+    def eo
+      @eo = @e.calc_eo(@r_mtx, @s)
+      return @eo
+    end
+
+    #=========================================================================
+    # Greenwich apparent sidereal time (Unit: rad)
+    #
+    # @param:  <none>
+    # @return: @gast
+    #=========================================================================
+    def gast
+      @era  = @e.calc_era(@jd, @t)   unless @era
+      @eo   = @e.calc_eo(@r_mtx, @s) unless @eo
+      @gast = @e.calc_gast(@era, @eo)
+      return @gast
+    end
+
+    #=========================================================================
+    # Greenwich apparent sidereal time (Unit: deg)
+    #
+    # @param:  <none>
+    # @return: @gast_deg
+    #=========================================================================
+    def gast_deg
+      @gast = gast unless @gast
+      @gast_deg = @gast / Const::PI_180
+      return @gast_deg
+    end
+
+    #=========================================================================
+    # Greenwich apparent sidereal time (Unit: HMS)
+    #
+    # @param:  <none>
+    # @return: @gast_hms
+    #=========================================================================
+    def gast_hms
+      @gast_deg = gast_deg unless @gast_deg
+      @gast_hms = deg2hms(@gast_deg)
+      return @gast_hms
+    end
+
+    #=========================================================================
+    # Greenwich mean sidereal time, IAU 2006. (Unit: rad)
+    #
+    # @param:  <none>
+    # @return: @gmst
+    #=========================================================================
+    def gmst
+      @era  = @e.calc_era(@jd, @t)   unless @era
+      @gmst = @e.calc_gmst(@era, @jc)
+      return @gmst
+    end
+
+    #=========================================================================
+    # Greenwich mean sidereal time, IAU 2006. (Unit: deg)
+    #
+    # @param:  <none>
+    # @return: @gmst_geg
+    #=========================================================================
+    def gmst_deg
+      @gmst = gmst unless @gast
+      @gmst_deg = @gmst / Const::PI_180
+      return @gmst_deg
+    end
+
+    #=========================================================================
+    # Greenwich mean sidereal time, IAU 2006. (Unit: HMS)
+    #
+    # @param:  <none>
+    # @return: @gmst_hms
+    #=========================================================================
+    def gmst_hms
+      @gmst_deg = gmst_deg unless @gmst_deg
+      @gmst_hms = deg2hms(@gmst_deg)
+      return @gmst_hms
+    end
+
+    #=========================================================================
+    # Equation of Equinoxes (Unit: rad)
+    #
+    # @param:  <none>
+    # @return: @ee
+    #=========================================================================
+    def ee
+      @gast = gast unless @gast
+      @gmst = gmst unless @gmst
+      @ee = @e.calc_ee(@gast, @gmst)
+      return @ee
+    end
+
+    #=========================================================================
+    # Equation of Equinoxes (Unit: deg)
+    #
+    # @param:  <none>
+    # @return: @ee_deg
+    #=========================================================================
+    def ee_deg
+      @ee = ee unless @ee
+      @ee_deg = @ee / Const::PI_180
+      return @ee_deg
+    end
+
+    #=========================================================================
+    # Equation of Equinoxes (Unit: HMS)
+    #
+    # @param:  <none>
+    # @return: @ee_hms
+    #=========================================================================
+    def ee_hms
+      @ee_deg = ee_deg unless @ee_deg
+      @ee_hms = deg2hms(@ee_deg)
+      return @ee_hms
+    end
+
+    private
+
+    #-------------------------------------------------------------------------
+    # ユリウス世紀数の計算
+    #
+    # @param:  jd  (Julian Day)
+    # @return: jc  (Julian Centry)
+    #-------------------------------------------------------------------------
+    def calc_jc(jd)
+      return (jd - Const::J2000) / Const::JC
+    rescue => e
+      raise
+    end
+
+    #-------------------------------------------------------------------------
+    # 行列の積 (3x3)
+    #
+    # @param:  r_1  (3x3-matrix)
+    # @param:  r_2  (3x3-matrix)
+    # @return: r    (3x3-matrix)
+    #-------------------------------------------------------------------------
+    def prod_mtx(r_1, r_2)
+      r = Array.new(3).map { |a| Array.new(3, 0.0) }
+
+      begin
+        0.upto(2) do |i|
+          0.upto(2) do |j|
+            r[i][j] = (0..2).inject(0.0) { |s, k| s + r_1[i][k] * r_2[k][j] }
+          end
+        end
+        return r
+      rescue => e
+        raise
+      end
+    end
+
+    def deg2hms(deg)
+      sign = ""
+
+      begin
+        h = (deg / 15.0).truncate
+        _m = (deg - h * 15.0) * 4.0
+        m = _m.truncate
+        s = (_m - m) * 60.0
+        if s < 0
+          s *= -1
+          sign = "-"
+        end
+        return sprintf("%s%d h %02d m %06.3f s", sign, h, m, s)
+      rescue => e
+        raise
+      end
+    end
+  end
+end
+

data/lib/mk_greenwich/version.rb

@@ -0,0 +1,3 @@
+module MkGreenwich
+  VERSION = "0.1.1"
+end

data/lib/mk_greenwich.rb

@@ -0,0 +1,18 @@
+require "date"
+require "eph_bpn"
+require "mk_time"
+require "mk_greenwich/version"
+require "mk_greenwich/argument"
+require "mk_greenwich/const"
+require "mk_greenwich/cip_cio"
+require "mk_greenwich/era_eors"
+require "mk_greenwich/greenwich"
+
+module MkGreenwich
+  def self.new(arg)
+    arg ||= Time.now.strftime("%Y%m%d%H%M%S")
+    ut = MkGreenwich::Argument.new(arg).get_utc
+    return unless ut
+    return MkGreenwich::Greenwich.new(ut)
+  end
+end

data/mk_greenwich.gemspec

@@ -0,0 +1,35 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mk_greenwich/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "mk_greenwich"
+  spec.version       = MkGreenwich::VERSION
+  spec.authors       = ["komasaru"]
+  spec.email         = ["masaru@mk-mode.com"]
+
+  spec.summary       = %q{Greenwich time calculation tool.}
+  spec.description   = %q{MkGreenwich is a Greenwich time calculation tool.}
+  spec.homepage      = "https://github.com/komasaru/mk_greenwich"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org. To allow pushes either set the 'allowed_push_host'
+  # to allow pushing to a single host or delete this section to allow pushing to any host.
+  #if spec.respond_to?(:metadata)
+  #  spec.metadata['allowed_push_host'] = "TODO: Set to 'http://mygemserver.com'"
+  #else
+  #  raise "RubyGems 2.0 or newer is required to protect against public gem pushes."
+  #end
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_dependency "eph_bpn"
+  spec.add_dependency "mk_time"
+  spec.add_development_dependency "bundler", "~> 1.12"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+end

metadata

@@ -0,0 +1,135 @@
+--- !ruby/object:Gem::Specification
+name: mk_greenwich
+version: !ruby/object:Gem::Version
+  version: 0.1.1
+platform: ruby
+authors:
+- komasaru
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2016-09-06 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: eph_bpn
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: mk_time
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.12'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.12'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+description: MkGreenwich is a Greenwich time calculation tool.
+email:
+- masaru@mk-mode.com
+executables:
+- mk_greenwich
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- Guardfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- exe/mk_greenwich
+- lib/mk_greenwich.rb
+- lib/mk_greenwich/argument.rb
+- lib/mk_greenwich/cip_cio.rb
+- lib/mk_greenwich/const.rb
+- lib/mk_greenwich/era_eors.rb
+- lib/mk_greenwich/fundamental_argument.rb
+- lib/mk_greenwich/greenwich.rb
+- lib/mk_greenwich/version.rb
+- mk_greenwich.gemspec
+homepage: https://github.com/komasaru/mk_greenwich
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.6.6
+signing_key: 
+specification_version: 4
+summary: Greenwich time calculation tool.
+test_files: []