mk_time 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 061938b212f49cb24231b13987da223ddc6abf41
+  data.tar.gz: 957e47cc7988230237b9ab31192e26f7e7d78e8c
+SHA512:
+  metadata.gz: 10b62a77ebb9991e75d82811e50d75e74115164fe9a0cca570da9a504df475efd9008a0856f8b9dc52db80510f83c6fd41c23ffa33b690c1d67d676681f7c301
+  data.tar.gz: 2e97a1543253eb803a0e393e1720e3b294678ba6b71a6a66966e41fd000794df937ef80aea194bff5b83ebc9b0f72809ebba7c47c81fe70e3a3fe48c907fb4a3

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_time.gemspec
+gemspec
+
+group :development do
+  gem "guard"
+  gem "guard-rspec", "~> 4.7.0"
+end
+

data/Guardfile

@@ -0,0 +1,42 @@
+# 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)
+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,67 @@
+# MkTime
+
+## Introduction
+
+This is a gem library which converts from UTC time to another time system.
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'mk_time'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mk_time
+
+## Usage
+
+### Instantiation
+
+    t = MkTime.new("20160726")
+    t = MkTime.new("20160726123456")
+    e = MkTime.new
+
+* You can set TDB 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.
+
+### Conversion
+
+    puts "     UTC: #{t.utc}"
+    puts "     JST: #{t.jst}"
+    puts "      JD: #{t.jd} day"
+    puts "       T: #{t.t} century"
+    puts " UTC-TAI: #{t.utc_tai} sec"
+    puts "LEAP_SEC: #{t.leap_sec} sec"  # = tuc_tai
+    puts "    DUT1: #{t.dut1} sec"
+    puts "  deltaT: #{t.dt} sec"
+    puts "     TAI: #{t.tai.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "     UT1: #{t.ut1.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "      TT: #{t.tt .instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "     TCG: #{t.tcg.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "     TCB: #{t.tcb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+    puts "     TDB: #{t.tdb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+
+* To check times under 1 second, I use `instance_eval` method.
+
+## 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_time.
+
+
+## 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_time"
+
+# 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_time

@@ -0,0 +1,20 @@
+#!/usr/bin/env ruby
+
+require "mk_time"
+
+t = MkTime.new(ARGV[0])
+puts "     UTC: #{t.utc}"
+puts "     JST: #{t.jst}"
+puts "      JD: #{t.jd} day"
+puts "       T: #{t.t} century"
+puts " UTC-TAI: #{t.utc_tai} sec"
+puts "LEAP_SEC: #{t.leap_sec} sec"  # = tuc_tai
+puts "    DUT1: #{t.dut1} sec"
+puts "  deltaT: #{t.dt} sec"
+puts "     TAI: #{t.tai.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "     UT1: #{t.ut1.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "      TT: #{t.tt .instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "     TCG: #{t.tcg.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "     TCB: #{t.tcb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+puts "     TDB: #{t.tdb.instance_eval { '%s.%03d' % [strftime('%Y-%m-%d %H:%M:%S'), (usec / 1000.0).round] }}"
+

data/lib/mk_time.rb

@@ -0,0 +1,15 @@
+require "date"
+require "mk_time/version"
+require "mk_time/consts"
+require "mk_time/argument"
+require "mk_time/compute"
+require "mk_time/calc"
+
+module MkTime
+  def self.new(arg)
+    arg ||= Time.now.strftime("%Y%m%d%H%M%S")
+    utc = MkTime::Argument.new(arg).get_utc
+    return unless utc
+    return MkTime::Calc.new(utc)
+  end
+end

data/lib/mk_time/argument.rb

@@ -0,0 +1,24 @@
+module MkTime
+  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_time/calc.rb

@@ -0,0 +1,25 @@
+module MkTime
+  class Calc
+    attr_reader :utc, :jst, :jd, :t, :utc_tai, :leap_sec, :dut1,
+                :ut1, :tai, :dt, :tt, :tcg, :tcb, :tdb
+    include Compute
+    alias leap_sec utc_tai
+
+    def initialize(arg)
+      @utc      = arg
+      @jst      = utc2jst(@utc)
+      @jd       = gc2jd(@utc)
+      @t        = jd2t(@jd)
+      @utc_tai  = get_utc_tai(@utc)
+      @dut1     = get_dut1(@utc)
+      @ut1      = utc2ut1(@utc)
+      @tai      = utc2tai(@utc)
+      @dt       = calc_dt(@utc)
+      @tt       = tai2tt(@tai)
+      @tcg      = tt2tcg(@tt)
+      @tcb      = tt2tcb(@tt)
+      @tdb      = tcb2tdb(@tcb)
+    end
+  end
+end
+

data/lib/mk_time/compute.rb

@@ -0,0 +1,524 @@
+module MkTime
+  module Compute
+    module_function
+
+    #=========================================================================
+    # UTC(協定世界時) -> JST(日本標準時)
+    #
+    # * JST = UTC + 09:00
+    #
+    # @param:  utc  (Time Object)
+    # @return: jst  (Time Object)
+    #=========================================================================
+    def utc2jst(utc)
+      jst = utc + Const::JST_OFFSET * 60 * 60
+      return Time.new(
+        jst.year, jst.month, jst.day,
+        jst.hour, jst.min, jst.sec,
+        sprintf("+%02d:00", Const::JST_OFFSET)
+      )
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # 年月日(グレゴリオ暦) -> JD(ユリウス日)
+    #
+    # * フリーゲルの公式を使用する
+    #   JD = int(365.25 * year)
+    #      + int(year / 400)
+    #      - int(year / 100)
+    #      + int(30.59 (month - 2))
+    #      + day
+    #      + 1721088
+    # * 上記の int(x) は厳密には、 x を超えない最大の整数
+    # * 「ユリウス日」でなく「準ユリウス日」を求めるなら、
+    #   `+ 1721088` を `- 678912` とする。
+    #
+    # @param:  t   (Time Object)
+    # @return: jd  (ユリウス日)
+    #=========================================================================
+    def gc2jd(t)
+      year, month, day = t.year, t.month, t.day
+      hour, min, sec   = t.hour, t.min, t.sec
+
+      begin
+        # 1月,2月は前年の13月,14月とする
+        if month < 3
+          year  -= 1
+          month += 12
+        end
+        # 日付(整数)部分計算
+        jd  = (365.25 * year).floor \
+            + (year / 400.0).floor \
+            - (year / 100.0).floor \
+            + (30.59 * (month - 2)).floor \
+            + day \
+            + 1721088.5
+        # 時間(小数)部分計算
+        t = (sec / 3600.0 + min / 60.0 + hour) / 24.0
+        return jd + t
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # JD(ユリウス日) -> T(ユリウス世紀数)
+    #
+    # * t = (JD - 2451545) / 36525.0
+    #
+    # @param:  jd  (ユリウス日)
+    # @return: t   (ユリウス世紀数)
+    #=========================================================================
+    def jd2t(jd)
+      return (jd - Const::J2000) / (Const::JY * 100)
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # UTC(協定世界時) - TAI(国際原子時) (= うるう秒総和) の取得
+    #
+    # * Ref: http://jjy.nict.go.jp/QandA/data/leapsec.html
+    #
+    # @param:  utc      (Time Object)
+    # @return: utc_tai  (sec)
+    #=========================================================================
+    def get_utc_tai(utc)
+      utc_tai = 0
+      target = utc.strftime("%Y%m%d")
+
+      begin
+        Const::LEAP_SECS.reverse.each do |date, sec|
+          if date <= target
+            utc_tai = sec
+            break
+          end
+        end
+        return utc_tai
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # DUT1 (= UT1(世界時1) - UTC(協定世界時)) の取得
+    #
+    # * Ref: http://jjy.nict.go.jp/QandA/data/dut1.html
+    #
+    # @param:  utc   (Time Object)
+    # @return: dut1  (sec)
+    #=========================================================================
+    def get_dut1(utc)
+      dut1 = 0
+      target = utc.strftime("%Y%m%d")
+
+      begin
+        Const::DUT1S.reverse.each do |date, sec|
+          if date <= target
+            dut1 = sec
+            break
+          end
+        end
+        return dut1
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # UTC(協定世界時) -> UT1(世界時1)
+    #
+    # * UT1 = UTC + DUT1
+    #
+    # @param:  utc  (Time Object)
+    # @return: ut1  (Time Object)
+    #=========================================================================
+    def utc2ut1(utc)
+      return utc + Rational(@dut1)
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # UTC(協定世界時) -> TAI(国際原子時)
+    #
+    # * TAI = UTC - UTC_TAI
+    #
+    # @param:  utc  (Time Object)
+    # @return: tai  (Time Object)
+    #=========================================================================
+    def utc2tai(utc)
+      return utc - @utc_tai
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # ΔT の計算
+    #
+    # * うるう秒調整が明確な場合は、うるう秒総和を使用した計算
+    # * そうでない場合は、NASA の計算式による計算
+    #
+    # @param:  utc  (Time Object)
+    # @return: dt   (sec)
+    #=========================================================================
+    def calc_dt(utc)
+      return Rational(Const::TT_TAI, 1000) - @utc_tai - @dut1 unless @utc_tai == 0
+      y = utc.year + (utc.month - 0.5) / 12
+      case
+      when                     utc.year <  -500; dt = calc_dt_before_m500(y)
+      when -500 <= utc.year && utc.year <   500; dt = calc_dt_before_500(y)
+      when  500 <= utc.year && utc.year <  1600; dt = calc_dt_before_1600(y)
+      when 1600 <= utc.year && utc.year <  1700; dt = calc_dt_before_1700(y)
+      when 1700 <= utc.year && utc.year <  1800; dt = calc_dt_before_1800(y)
+      when 1800 <= utc.year && utc.year <  1860; dt = calc_dt_before_1860(y)
+      when 1860 <= utc.year && utc.year <  1900; dt = calc_dt_before_1900(y)
+      when 1900 <= utc.year && utc.year <  1920; dt = calc_dt_before_1920(y)
+      when 1920 <= utc.year && utc.year <  1941; dt = calc_dt_before_1941(y)
+      when 1941 <= utc.year && utc.year <  1961; dt = calc_dt_before_1961(y)
+      when 1961 <= utc.year && utc.year <  1986; dt = calc_dt_before_1986(y)
+      when 1986 <= utc.year && utc.year <  2005; dt = calc_dt_before_2005(y)
+      when 2005 <= utc.year && utc.year <  2050; dt = calc_dt_before_2050(y)
+      when 2050 <= utc.year && utc.year <= 2150; dt = calc_dt_until_2150(y)
+      when 2150 <  utc.year                    ; dt = calc_dt_after_2150(y)
+      end
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (year < -500)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_m500(y)
+      t = (y - 1820) / 100.0
+      dt = -20 + 32 * t ** 2
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (-500 <= year && year < 500)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_500(y)
+      t = y / 100.0
+      dt  = 10583.6
+           (-1014.41        + \
+           (   33.78311     + \
+           (   -5.952053    + \
+           (   -0.1798452   + \
+           (    0.022174192 + \
+           (    0.0090316521) \
+           * t) * t) * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (500 <= year && year < 1600)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1600(y)
+      t = (y - 1000) / 100.0
+      dt  = 1574.2         + \
+           (-556.01        + \
+           (  71.23472     + \
+           (   0.319781    + \
+           (  -0.8503463   + \
+           (  -0.005050998 + \
+           (   0.0083572073) \
+           * t) * t) * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1600 <= year && year < 1700)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1700(y)
+      t = y - 1600
+      dt  = 120           + \
+           ( -0.9808      + \
+           ( -0.01532     + \
+           (  1.0 / 7129.0) \
+           * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1700 <= year && year < 1800)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1800(y)
+      t = y - 1700
+      dt  =  8.83           + \
+           ( 0.1603         + \
+           (-0.0059285      + \
+           ( 0.00013336     + \
+           (-1.0 / 1174000.0) \
+           * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1800 <= year && year < 1860)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1860(y)
+      t = y - 1800
+      dt  = 13.72          + \
+           (-0.332447      + \
+           ( 0.0068612     + \
+           ( 0.0041116     + \
+           (-0.00037436    + \
+           ( 0.0000121272  + \
+           (-0.0000001699  + \
+           ( 0.000000000875) \
+           * t) * t) * t) * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1860 <= year && year < 1900)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1900(y)
+      t = y - 1860
+      dt  =  7.62          + \
+           ( 0.5737        + \
+           (-0.251754      + \
+           ( 0.01680668    + \
+           (-0.0004473624  + \
+           ( 1.0 / 233174.0) \
+           * t) * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1900 <= year && year < 1920)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1920(y)
+      t = y - 1900
+      dt  = -2.79      + \
+           ( 1.494119  + \
+           (-0.0598939 + \
+           ( 0.0061966 + \
+           (-0.000197  ) \
+           * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1920 <= year && year < 1941)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1941(y)
+      t = y - 1920
+      dt  = 21.20     + \
+           ( 0.84493  + \
+           (-0.076100 + \
+           ( 0.0020936) \
+           * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1941 <= year && year < 1961)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1961(y)
+      t = y - 1950
+      dt  = 29.07      + \
+           ( 0.407     + \
+           (-1 / 233.0 + \
+           ( 1 / 2547.0) \
+           * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1961 <= year && year < 1986)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_1986(y)
+      t = y - 1975
+      dt = 45.45      + \
+          ( 1.067     + \
+          (-1 / 260.0 + \
+          (-1 / 718.0)  \
+          * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (1986 <= year && year < 2005)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_2005(y)
+      t = y - 2000
+      dt  = 63.86         + \
+           ( 0.3345       + \
+           (-0.060374     + \
+           ( 0.0017275    + \
+           ( 0.000651814  + \
+           ( 0.00002373599) \
+           * t) * t) * t) * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (2005 <= year && year < 2050)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_before_2050(y)
+      t = y - 2000
+      dt  = 62.92    + \
+           ( 0.32217 + \
+           ( 0.005589) \
+           * t) * t
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (2050 <= year && year <= 2150)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_until_2150(y)
+      dt  = -20 \
+          + 32 * ((y - 1820) / 100.0) ** 2
+          - 0.5628 * (2150 - y)
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # deltaT (2150 < year)
+    #
+    # @param:  y   (Year coefficient for deltaT calculation)
+    # @return: dt  (sec)
+    #=========================================================================
+    def calc_dt_after_2150(y)
+      t = (y - 1820) / 100.0
+      dt  = -20 + 32 * t ** 2
+      return dt
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # TT(地球時)
+    #
+    # * TT = TAI + TT_TAI
+    #      = UT1 + ΔT
+    #=========================================================================
+    def tai2tt(tai)
+      return tai + Rational(Const::TT_TAI, 1000)
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # TCG(地球重心座標時)
+    #
+    # * TCG = TT + L_G * (JD - T_0) * 86,400
+    #   （JD: ユリウス日,
+    #     L_G = 6.969290134 * 10^(-10), T_0 = 2,443,144.5003725）
+    #
+    # @param:  tt   (Time Object)
+    # @return: tcg  (Time Object)
+    #=========================================================================
+    def tt2tcg(tt)
+      return tt + Const::L_G * (@jd - Const::T_0) * 86400
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # TCB(太陽系重心座標時)
+    #
+    # * TCB = TT + L_B * (JD - T_0) * 86400
+    #
+    # @param:  tt   (Time Object)
+    # @return: tcb  (Time Object)
+    #=========================================================================
+    def tt2tcb(tt)
+      return tt + Const::L_B * (@jd - Const::T_0) * 86400
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # TCB(太陽系重心座標時) -> TDB(太陽系力学時)
+    #
+    # * TDB = TCB - L_B * (JD_TCB - T_0) * 86400 + TDB_0
+    #
+    # @param:  tcb  (Time Object)
+    # @return: tdb  (Time Object)
+    #=========================================================================
+    def tcb2tdb(tcb)
+      jd_tcb = gc2jd(tcb)
+      return tcb - Const::L_B * (jd_tcb - Const::T_0) * 86400 + Const::TDB_0
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_time/consts.rb

@@ -0,0 +1,179 @@
+module MkTime
+  module Const
+    MSG_ERR_1  = "[ERROR] Format: YYYYMMDD or YYYYMMDDHHMMSS"
+    MSG_ERR_2  = "[ERROR] Invalid date-time!"
+    JST_OFFSET = 9                # JST offset from UTC
+    J2000      = 2451545          # Julian Day of 2000-01-01 12:00:00
+    JY         = 365.25           # 1 Julian Year
+    TT_TAI     = 32184            # (TT - TAI) * 1000 (計算時の誤差抑制のため 1000 倍している)
+    L_G        = 6.969290134e-10  # for TCG
+    L_B        = 1.550519768e-8   # for TCG, TDB
+    T_0        = 2443144.5003725  # for TCG, TDB, TCB
+    TDB_0      = -6.55e-5         # for TDB
+    LEAP_SECS  = [
+                   ["19720101", -10],
+                   ["19720701", -11],
+                   ["19730101", -12],
+                   ["19740101", -13],
+                   ["19750101", -14],
+                   ["19760101", -15],
+                   ["19770101", -16],
+                   ["19780101", -17],
+                   ["19790101", -18],
+                   ["19800101", -19],
+                   ["19810701", -20],
+                   ["19820701", -21],
+                   ["19830701", -22],
+                   ["19850701", -23],
+                   ["19880101", -24],
+                   ["19900101", -25],
+                   ["19910101", -26],
+                   ["19920701", -27],
+                   ["19930701", -28],
+                   ["19940701", -29],
+                   ["19960101", -30],
+                   ["19970701", -31],
+                   ["19990101", -32],
+                   ["20060101", -33],
+                   ["20090101", -34],
+                   ["20120701", -35],
+                   ["20150701", -36],
+                   ["20170101", -37],
+                   ["20190101",   0]  # (<= Provisional end point)
+                 ].freeze  # Leap Second's adjustment
+    DUT1S      = [
+                   ["19880317",  0.2],
+                   ["19880512",  0.1],
+                   ["19880825",  0.0],
+                   ["19881110", -0.1],
+                   ["19890119", -0.2],
+                   ["19890406", -0.3],
+                   ["19890608", -0.4],
+                   ["19890921", -0.5],
+                   ["19891116", -0.6],
+                   ["19900101",  0.3],
+                   ["19900301",  0.2],
+                   ["19900412",  0.1],
+                   ["19900510",  0.0],
+                   ["19900726", -0.1],
+                   ["19900920", -0.2],
+                   ["19901101", -0.3],
+                   ["19910101",  0.6],
+                   ["19910207",  0.5],
+                   ["19910321",  0.4],
+                   ["19910425",  0.3],
+                   ["19910620",  0.2],
+                   ["19910822",  0.1],
+                   ["19911017",  0.0],
+                   ["19911121", -0.1],
+                   ["19920123", -0.2],
+                   ["19920227", -0.3],
+                   ["19920402", -0.4],
+                   ["19920507", -0.5],
+                   ["19920701",  0.4],
+                   ["19920903",  0.3],
+                   ["19921022",  0.2],
+                   ["19921126",  0.1],
+                   ["19930114",  0.0],
+                   ["19930218", -0.1],
+                   ["19930401", -0.2],
+                   ["19930506", -0.3],
+                   ["19930701",  0.6],
+                   ["19930819",  0.5],
+                   ["19931007",  0.4],
+                   ["19931118",  0.3],
+                   ["19931230",  0.2],
+                   ["19940210",  0.1],
+                   ["19940317",  0.0],
+                   ["19940421", -0.1],
+                   ["19940609", -0.2],
+                   ["19940701",  0.8],
+                   ["19940811",  0.7],
+                   ["19941006",  0.6],
+                   ["19941117",  0.5],
+                   ["19941222",  0.4],
+                   ["19950223",  0.3],
+                   ["19950316",  0.2],
+                   ["19950413",  0.1],
+                   ["19950525",  0.0],
+                   ["19950713", -0.1],
+                   ["19950907", -0.2],
+                   ["19951026", -0.3],
+                   ["19951130", -0.4],
+                   ["19960101",  0.5],
+                   ["19960222",  0.4],
+                   ["19960411",  0.3],
+                   ["19960516",  0.2],
+                   ["19960808",  0.1],
+                   ["19961003",  0.0],
+                   ["19961205", -0.1],
+                   ["19970206", -0.2],
+                   ["19970320", -0.3],
+                   ["19970508", -0.4],
+                   ["19970626", -0.5],
+                   ["19970701",  0.5],
+                   ["19970918",  0.4],
+                   ["19971030",  0.3],
+                   ["19971218",  0.2],
+                   ["19980219",  0.1],
+                   ["19980326",  0.0],
+                   ["19980507", -0.1],
+                   ["19980813", -0.2],
+                   ["19981126", -0.3],
+                   ["19990101",  0.7],
+                   ["19990304",  0.6],
+                   ["19990527",  0.5],
+                   ["19991014",  0.4],
+                   ["20000106",  0.3],
+                   ["20000413",  0.2],
+                   ["20001019",  0.1],
+                   ["20010301",  0.0],
+                   ["20011004", -0.1],
+                   ["20020214", -0.2],
+                   ["20021024", -0.3],
+                   ["20030403", -0.4],
+                   ["20040429", -0.5],
+                   ["20050317", -0.6],
+                   ["20060101",  0.3],
+                   ["20060427",  0.2],
+                   ["20060928",  0.1],
+                   ["20061222",  0.0],
+                   ["20070315", -0.1],
+                   ["20070614", -0.2],
+                   ["20071129", -0.3],
+                   ["20080313", -0.4],
+                   ["20080807", -0.5],
+                   ["20081120", -0.6],
+                   ["20090101",  0.4],
+                   ["20090312",  0.3],
+                   ["20090611",  0.2],
+                   ["20091112",  0.1],
+                   ["20100311",  0.0],
+                   ["20100603", -0.1],
+                   ["20110106", -0.2],
+                   ["20110512", -0.3],
+                   ["20111104", -0.4],
+                   ["20120209", -0.5],
+                   ["20120510", -0.6],
+                   ["20120701",  0.4],
+                   ["20121025",  0.3],
+                   ["20130131",  0.2],
+                   ["20130411",  0.1],
+                   ["20130822",  0.0],
+                   ["20131121", -0.1],
+                   ["20140220", -0.2],
+                   ["20140508", -0.3],
+                   ["20140925", -0.4],
+                   ["20141225", -0.5],
+                   ["20150319", -0.6],
+                   ["20150528", -0.7],
+                   ["20150701",  0.3],
+                   ["20150917",  0.2],
+                   ["20151126",  0.1],
+                   ["20160131",  0.0],
+                   ["20160324", -0.1],
+                   ["20160519", -0.2],
+                   ["20160819",  0.0]  # (<= Provisional end point)
+                 ].freeze  # DUT1 adjustment
+  end
+end

data/lib/mk_time/version.rb

@@ -0,0 +1,3 @@
+module MkTime
+  VERSION = "0.1.0"
+end

data/mk_time.gemspec

@@ -0,0 +1,33 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mk_time/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "mk_time"
+  spec.version       = MkTime::VERSION
+  spec.authors       = ["komasaru"]
+  spec.email         = ["masaru@mk-mode.com"]
+
+  spec.summary       = %q{Time system conversion library.}
+  spec.description   = %q{MkTime is a time system conversion library.}
+  spec.homepage      = "https://github.com/komasaru/mk_time"
+  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_development_dependency "bundler", "~> 1.12"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+end

metadata

@@ -0,0 +1,105 @@
+--- !ruby/object:Gem::Specification
+name: mk_time
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- komasaru
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2016-07-26 00:00:00.000000000 Z
+dependencies:
+- !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: MkTime is a time system conversion library.
+email:
+- masaru@mk-mode.com
+executables:
+- mk_time
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- Guardfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- exe/mk_time
+- lib/mk_time.rb
+- lib/mk_time/argument.rb
+- lib/mk_time/calc.rb
+- lib/mk_time/compute.rb
+- lib/mk_time/consts.rb
+- lib/mk_time/version.rb
+- mk_time.gemspec
+homepage: https://github.com/komasaru/mk_time
+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: Time system conversion library.
+test_files: []