when_exe 0.3.4 → 0.3.5

data/lib/when_exe/region/coptic.rb

@@ -0,0 +1,88 @@
+# -*- coding: utf-8 -*-
+=begin
+  Copyright (C) 2011-2014 Takashi SUGA
+
+  You may use and/or modify this file according to the license described in the LICENSE.txt file included in this archive.
+=end
+
+module When
+
+  class BasicTypes::M17n
+
+    CopticTerms = [self, [
+      "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
+      "locale:[=en:, ja=ja:, alias]",
+      "names:[CopticTerms=]",
+      "[Coptic=en:Coptic_calendar,       コプト暦    ]",
+      "[Ethiopian=en:Ethiopian_calendar, エチオピア暦=en:Ethiopian_calendar]",
+
+      [self,
+        "names:[EgyptianMonth=, 月=ja:%%<月_(暦)>]",
+        "[tut=,      トート=      ]",
+        "[baba=,     バーバ=      ]",
+        "[hatur=,    ハートール=  ]",
+        "[kiyahak=,  キヤハーク=  ]",
+        "[tuba=,     トーバ=      ]",
+        "[amshir=,   アムシール=  ]",
+        "[baramhat=, バラムハート=]",
+        "[barmuda=,  バルムーダ=  ]",
+        "[bashans=,  バシャンス=  ]",
+        "[ba'una=,   バウーナ=    ]",
+        "[abib=,     アビーブ=    ]",
+        "[misra=,    ミスラー=    ]",
+        "[epagomen=, エパゴメネ=  ]"
+      ],
+
+      [self,
+        "names:[EthiopianMonth=, 月=ja:%%<月_(暦)>]",
+        "[Mäskäräm=, マスカラム=  ]",
+        "[Ṭəqəmt=,   テケルト=    ]",
+        "[Ḫədar=,    ヘダル=      ]",
+        "[Taḫśaś=,   ターサス=    ]",
+        "[Ṭərr=,     テル=        ]",
+        "[Yäkatit=,  イェカティト=]",
+        "[Mägabit=,  メガビト=    ]",
+        "[Miyazya=,  ミアジア=    ]",
+        "[Gənbot=,   ゲエンポト=  ]",
+        "[Säne=,     セネ=        ]",
+        "[Ḥamle=,    ハムレ=      ]",
+        "[Nähase=,   ネハッセ=    ]",
+        "[Ṗagʷəmen=, パゴウメン=  ]"
+      ]
+    ]]
+  end
+
+  module CalendarTypes
+
+    _egyptian_month_indices  = [
+      When::Coordinates::Index.new({:unit =>13, :trunk=>When::Parts::Resource._instance('_m:CopticTerms::EgyptianMonth::*')}),
+      When::Coordinates::Index.new
+    ]
+
+    _ethiopian_month_indices = [
+      When::Coordinates::Index.new({:unit =>13, :trunk=>When::Parts::Resource._instance('_m:CopticTerms::EthiopianMonth::*')}),
+      When::Coordinates::Index.new
+    ]
+
+    #
+    # Coptic Calendar in Egypt and Ethiopia
+    #
+    Coptic =  [{'Epoch'=>{'284Y'=>{'origin_of_MSC' =>   1,
+                                   'label'         => Parts::Resource._instance('_m:CopticTerms::Coptic'),
+                                   'indices'       => _egyptian_month_indices},
+                            '8Y'=>{'origin_of_MSC' => 277,
+                                   'label'         => Parts::Resource._instance('_m:CopticTerms::Ethiopian'),
+                                   'indices'       => _ethiopian_month_indices}}}, CyclicTableBased, {
+      'label'         => Parts::Resource._instance('_m:CopticTerms::Coptic'),
+      'origin_of_LSC' => 1825030,
+      'origin_of_MSC' =>       1,
+      'epoch_in_CE'   =>     285,
+      'indices'       => _egyptian_month_indices,
+      'rule_table' => {
+        'T' => {'Rule'  =>[365,365,366,365]},
+        365 => {'Length'=>[30]*12+[5]},
+        366 => {'Length'=>[30]*12+[6]}
+      }
+    }]
+  end
+end

data/lib/when_exe/region/ephemeric_notes.rb

@@ -8,346 +8,361 @@
 #
 # Ephemeris を用いる暦注
 #
-module When::CalendarTypes
-
-  class CalendarNote
-    #
-    # 太陽と月の位置によるイベント
-    #
-    class LuniSolarPositions < self
-
-      # 座標の分子
-      #
-      # @return [Numeric]
-      #
-      attr_reader :num
-
-      # 座標の分母
-      #
-      # @return [Numeric]
-      #
-      attr_reader :den
-
-      # 計算アルゴリズム
-      #
-      # @return [When::Ephemeris::Formula]
-      #
-      attr_reader :formula
-
-      # enumerator の周期
-      #
-      # @return [Numeric]
-      #
-      attr_reader :delta
-
-      # 没滅計算用の補正
-      #
-      # @return [Numeric]
-      #
-      attr_reader :margin
-
-      # イベントの日時
-      #
-      # @param [When::TM::TemporalPosition] date イベントを探す基準とする日時
-      # @param [Array<Numeric>] parameter 座標の分子と分母( num, den)
-      #
-      #   num 座標の分子 (デフォルト @num)
-      #
-      #   den 座標の分母 (デフォルト @den)
-      #
-      # @param [String] parameter  座標の分子と分母("#{ num }/#{ den }" の形式)
-      # @param [Integer] precision 取得したい時間位置の分解能(デフォルト date の分解能)
-      #
-      # @return [When::TM::CalDate] date またはその直後のイベントの日時
-      #
-      def term(date, parameter=nil, precision=date.precision)
-        precision = nil if precision == When::SYSTEM
-        num, den  = parameter.kind_of?(String) ? parameter.split(/\//, 2) : parameter
-        num = (num || @num).to_f
-        den = (den || @den).to_f
-        date      = date.floor(precision) if precision
-        options   = date._attr
-        quot, mod = (@formula.time_to_cn(date)*30.0).divmod(den)
-        cycle     = quot * den + num
-        cycle    += den if mod > (num % den)
-        time      = When::TM::JulianDate._d_to_t(@formula.cn_to_time(cycle/30.0))
-        time      = date.time_standard.from_dynamical_time(time) if @formula.is_dynamical
-        date      = date.frame.jul_trans(When::TM::JulianDate.universal_time(time), options)
-        precision ? date.floor(precision) : date
-      end
+class When::CalendarNote
+  #
+  # 太陽と月の位置によるイベント
+  #
+  class LuniSolarPositions < self
 
-      # 日付に対応する座標
-      #
-      # @param [When::TM::TemporalPosition] date 日付
-      # @param [Numeric] delta 周期の補正(土用の時刻の補正に使用,デフォルト 0)
-      #
-      # @return [Array<Integer>] Array< Integer, 0 or 1 or 2 >
-      #
-      #   [Integer]     対応する座標
-      #
-      #   [0 or 1 or 2] 座標の進み(0 なら 没, 2 なら滅)
-      #
-      def position(date, delta=0)
-        date   = date.floor
-        p0, p1 = [date, date.succ].map {|d| (30.0 * @formula.time_to_cn(d) - @margin + delta).floor}
-        [p1 % @den, p1 - p0]
-      end
+    # 座標の分子
+    #
+    # @return [Numeric]
+    #
+    attr_reader :num
 
-      #
-      # イベントの標準的な間隔を返す
-      #
-      # @param [String] parameter 座標の分子と分母("#{ num }/#{ den }" の形式)
-      #
-      # @return [When::TM::IntervalLength]
-      def term_delta(parameter=nil)
-        return @delta unless parameter
-        num, den = parameter.split(/\//, 2)
-        When::TM::IntervalLength.new([(den || @den).to_f-1,1].max, 'day')
-      end
-    end
+    # 座標の分母
+    #
+    # @return [Numeric]
+    #
+    attr_reader :den
 
+    # 計算アルゴリズム
     #
-    # 二十四節気
+    # @return [When::Ephemeris::Formula]
     #
-    class SolarTerms < LuniSolarPositions
+    attr_reader :formula
 
-      private
+    # enumerator の周期
+    #
+    # @return [Numeric]
+    #
+    attr_reader :delta
 
-      # オブジェクトの正規化
-      #   num     - 太陽黄経/度の分子 (デフォルト   0 - 春分)
-      #   den     - 太陽黄経/度の分母 (デフォルト 360 - 1年)
-      #   formula - 計算アルゴリズム(デフォルト '_ep:Formula?formula=12S')
-      #   delta   - enumerator の周期 (デフォルト (den/360)年)
-      #   margin  - 没滅計算用の補正  (デフォルト 1E-8)
-      def _normalize(args=[], options={})
-        num, den, formula, delta, margin = args
-        @num     = (num || @num  ||   0).to_f
-        @den     = (den || @den  || 360).to_f
-        @formula = When.Resource(formula || @formula ||'Formula?formula=12S', '_ep:')
-        @delta   = When.Duration(delta   || @delta   || When::TM::IntervalLength.new(@den/360, 'year'))
-        @margin  = (margin || @margin    || 1E-8).to_f
-        super
-      end
-    end
+    # 没滅計算用の補正
+    #
+    # @return [Numeric]
+    #
+    attr_reader :margin
 
+    # イベントの日時
+    #
+    # @param [When::TM::TemporalPosition] date イベントを探す基準とする日時
+    # @param [Array<Numeric>] parameter 座標の分子と分母( num, den)
+    #
+    #   num 座標の分子 (デフォルト @num)
     #
-    # 月の位相
+    #   den 座標の分母 (デフォルト @den)
     #
-    class LunarPhases < LuniSolarPositions
+    # @param [String] parameter  座標の分子と分母("#{ num }/#{ den }" の形式)
+    # @param [Integer] precision 取得したい時間位置の分解能(デフォルト date の分解能)
+    #
+    # @return [When::TM::CalDate] date またはその直後のイベントの日時
+    #
+    def event(date, parameter=nil, precision=date.precision)
+      num, den  = parameter.kind_of?(String) ? parameter.split(/\//, 2) : parameter
+      num = (num || @num).to_f
+      den = (den || @den).to_f
+      date      = date.floor(precision) if precision < date.precision
+      options   = date._attr
+      is_date_and_time = options.key?(:clock) || precision > When::DAY
+      options[:precision] = precision
+      options[:clock]   ||= date.frame.time_basis || When::TM::Clock.local_time
+      quot, mod = (@formula.time_to_cn(date)*30.0).divmod(den)
+      cycle     = quot * den + num
+      cycle    += den if mod > (num % den)
+      time      = When::TM::JulianDate._d_to_t(@formula.cn_to_time(cycle/30.0))
+      time      = date.time_standard.from_dynamical_time(time) if @formula.is_dynamical
+      event = date.frame.jul_trans(When::TM::JulianDate.universal_time(time), options)
+      is_date_and_time ? event : event.to_cal_date
+    end
 
-      private
+    # 日付に対応する座標
+    #
+    # @param [When::TM::TemporalPosition] date 日付
+    # @param [Numeric] delta 周期の補正(土用の時刻の補正に使用,デフォルト 0)
+    #
+    # @return [Array<Integer>] Array< Integer, 0 or 1 or 2 >
+    #
+    #   [Integer]     対応する座標
+    #
+    #   [0 or 1 or 2] 座標の進み(0 なら 没, 2 なら滅)
+    #
+    def position(date, delta=0)
+      date   = date.floor
+      p0, p1 = [date, date.succ].map {|d| (30.0 * @formula.time_to_cn(d) - @margin + delta).floor}
+      [p1 % @den, p1 - p0]
+    end
 
-      # オブジェクトの正規化
-      #   num     - 月の位相/12度の分子 (デフォルト  0 - 朔)
-      #   den     - 月の位相/12度の分母 (デフォルト 30 - 1月)
-      #   formula - 計算アルゴリズム(デフォルト '_ep:Formula?formula=1L')
-      #   delta   - enumerator の周期 (デフォルト (den/30)月)
-      #   margin  - 没滅計算用の補正  (デフォルト 1E-8)
-      def _normalize(args=[], options={})
-        num, den, formula, delta, margin = args
-        @num     = (num || @num  ||  0).to_f
-        @den     = (den || @den  || 30).to_f
-        @formula = When.Resource(formula || @formula ||'Formula?formula=1L', '_ep:')
-        @delta   = When.Duration(delta   || @delta   || When::TM::IntervalLength.new(@den/30, 'month'))
-        @margin  = (margin || @margin    || 1E-8).to_f
-        super
-      end
+    #
+    # イベントの標準的な間隔を返す
+    #
+    # @param [String] parameter 座標の分子と分母("#{ num }/#{ den }" の形式)
+    #
+    # @return [When::TM::IntervalLength]
+    def event_delta(parameter=nil)
+      return @delta unless parameter
+      num, den = parameter.split(/\//, 2)
+      When::TM::IntervalLength.new([(den || @den).to_f,1].max*0.9, 'day')
     end
+  end
+
+  #
+  # 二十四節気
+  #
+  class SolarTerms < LuniSolarPositions
 
+    # 二十四節気のための event の別名
     #
-    # 天体暦の暦注
+    # @return [When::TM::CalDate] date またはその直後のイベントの日時
+    alias :term :event
+
+    private
+
+    # 二十四節気のための event_delta の別名
+    alias :term_delta :event_delta
+
+    # オブジェクトの正規化
+    #   num     - 太陽黄経/度の分子 (デフォルト   0 - 春分)
+    #   den     - 太陽黄経/度の分母 (デフォルト 360 - 1年)
+    #   formula - 計算アルゴリズム(デフォルト '_ep:Formula?formula=12S')
+    #   delta   - enumerator の周期 (デフォルト (den/360)年)
+    #   margin  - 没滅計算用の補正  (デフォルト 1E-8)
+    def _normalize(args=[], options={})
+      num, den, formula, delta, margin = args
+      @num     = (num || @num  ||   0).to_f
+      @den     = (den || @den  || 360).to_f
+      @formula = When.Resource(formula || @formula ||'Formula?formula=12S', '_ep:')
+      @delta   = When.Duration(delta   || @delta   || When::TM::IntervalLength.new(@den/360, 'year'))
+      @margin  = (margin || @margin    || 1E-8).to_f
+      super
+    end
+  end
+
+  #
+  # 月の位相
+  #
+  class LunarPhases < LuniSolarPositions
+
+    # 月の位相のための event の別名
     #
-    class EphemericNote < self
+    # @return [When::TM::CalDate] date またはその直後のイベントの日時
+    alias :phase :event
+
+    private
+
+    # 月の位相のための event_delta の別名
+    alias :phase_delta :event_delta
+
+    # オブジェクトの正規化
+    #   num     - 月の位相/12度の分子 (デフォルト  0 - 朔)
+    #   den     - 月の位相/12度の分母 (デフォルト 30 - 1月)
+    #   formula - 計算アルゴリズム(デフォルト '_ep:Formula?formula=1L')
+    #   delta   - enumerator の周期 (デフォルト (den/30)月)
+    #   margin  - 没滅計算用の補正  (デフォルト 1E-8)
+    def _normalize(args=[], options={})
+      num, den, formula, delta, margin = args
+      @num     = (num || @num  ||  0).to_f
+      @den     = (den || @den  || 30).to_f
+      @formula = When.Resource(formula || @formula ||'Formula?formula=1L', '_ep:')
+      @delta   = When.Duration(delta   || @delta   || When::TM::IntervalLength.new(@den/30, 'month'))
+      @margin  = (margin || @margin    || 1E-8).to_f
+      super
+    end
+  end
 
-      NoteObjects = [When::BasicTypes::M17n, [
-        "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
-        "locale:[=en:, ja=ja:, alias=ja:]",
-        "names:[Ephemeric]",
+  #
+  # 天体暦の暦注
+  #
+  class EphemericNote < self
 
-        # 年の暦注 ----------------------------
-        [When::BasicTypes::M17n,
-          "names:[year]"
-        ],
+    NoteObjects = [When::BasicTypes::M17n, [
+      "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
+      "locale:[=en:, ja=ja:, alias=ja:]",
+      "names:[Ephemeric]",
+
+      # 年の暦注 ----------------------------
+      [When::BasicTypes::M17n,
+        "names:[year]"
+      ],
 
-        # 月の暦注 ----------------------------
+      # 月の暦注 ----------------------------
+      [When::BasicTypes::M17n,
+        "names:[month]",
         [When::BasicTypes::M17n,
-          "names:[month]",
+          "names:[Month]"
+        ]
+      ],
+
+      # 日の暦注 ----------------------------
+      [When::BasicTypes::M17n,
+        "names:[day]",
+          "[Sunrise,   日の出          ]", # 日の出
+          "[Sunset,    日の入り        ]", # 日の入り
+          [When::Coordinates::Residue,
+            "label:[Moon_Age=, 正午月齢=ja:%%<月齢>]",
+            "divisor:60",
+            "format:[%s(%4.1f)]"
+          ],
+          "[Moonrise,  月の出          ]", # 月の出
+          "[Moonset=,  月の入り        ]", # 月の入り
           [When::BasicTypes::M17n,
-            "names:[Month]"
+            "names:[Tide, 潮汐]",            # 満潮干潮日時
+            "[High_Tide=en:Tide, 満潮=ja:%%<潮汐>]",
+            "[Low_Tide=en:Tide,  干潮=ja:%%<潮汐>]"
           ]
-        ],
+      ]
+    ]]
 
-        # 日の暦注 ----------------------------
-        [When::BasicTypes::M17n,
-          "names:[day]",
-            "[Sunrise,   日の出          ]", # 日の出
-            "[Sunset,    日の入り        ]", # 日の入り
-            [When::Coordinates::Residue,
-              "label:[Moon_Age=, 正午月齢=ja:%%<月齢>]",
-              "divisor:60",
-              "format:[%s(%4.1f)]"
-            ],
-            "[Moonrise,  月の出          ]", # 月の出
-            "[Moonset=,  月の入り        ]", # 月の入り
-            [When::BasicTypes::M17n,
-              "names:[Tide, 潮汐]",            # 満潮干潮日時
-              "[High_Tide=en:Tide, 満潮=ja:%%<潮汐>]",
-              "[Low_Tide=en:Tide,  干潮=ja:%%<潮汐>]"
-            ]
-        ]
-      ]]
-
-      #
-      # 日の出
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options dummy
-      #
-      # @return [When::TM::TemporalPosition] 日の出の時刻をイベント時刻とする
-      #
-      def sunrise(date, options={})
-        event = formula(date.location.iri).sunrise(date)
-        event.events = [@root['Sunrise']]
-        event
-      rescue
-        nil
-      end
+    #
+    # 日の出
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options dummy
+    #
+    # @return [When::TM::TemporalPosition] 日の出の時刻をイベント時刻とする
+    #
+    def sunrise(date, options={})
+      event = formula(date.location.iri).sunrise(date)
+      event.events = [@root['Sunrise']]
+      event
+    rescue
+      nil
+    end
 
-      #
-      # 日の入り
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options dummy
-      #
-      # @return [When::TM::TemporalPosition] 日の入りの時刻をイベント時刻とする
-      #
-      def sunset(date, options={})
-        event = formula(date.location.iri).sunset(date)
-        event.events = [@root['Sunset']]
-        event
-      rescue
-        nil
-      end
+    #
+    # 日の入り
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options dummy
+    #
+    # @return [When::TM::TemporalPosition] 日の入りの時刻をイベント時刻とする
+    #
+    def sunset(date, options={})
+      event = formula(date.location.iri).sunset(date)
+      event.events = [@root['Sunset']]
+      event
+    rescue
+      nil
+    end
 
-      #
-      # 正午月齢
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options dummy
-      #
-      # @return [Numeric] 正午における朔からの経過日数
-      #
-      def moon_age(date, options={})
-        @phase ||= When.CalendarNote('LunarPhases')
-        noon = date.floor(When::DAY,When::SYSTEM) + 0.5
-        @root['Moon_Age'][noon.to_f - @phase.term(noon, [-30.0,30.0]).to_f]
-      end
+    #
+    # 正午月齢
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options dummy
+    #
+    # @return [Numeric] 正午における朔からの経過日数
+    #
+    def moon_age(date, options={})
+      @phase ||= When.CalendarNote('LunarPhases')
+      noon = date.floor(When::DAY,When::SYSTEM) + 0.5
+      @root['Moon_Age'][noon.to_f - @phase.phase(noon, [-30.0,30.0]).to_f]
+    end
 
-      #
-      # 月の出
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options dummy
-      #
-      # @return [When::TM::TemporalPosition] 月の出の時刻をイベント時刻とする
-      #
-      def moonrise(date, options={})
-        event = formula(date.location.iri).moonrise(date)
-        event.events = [@root['Moonrise']]
-        event
-      rescue
-        nil
-      end
+    #
+    # 月の出
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options dummy
+    #
+    # @return [When::TM::TemporalPosition] 月の出の時刻をイベント時刻とする
+    #
+    def moonrise(date, options={})
+      event = formula(date.location.iri).moonrise(date)
+      event.events = [@root['Moonrise']]
+      event
+    rescue
+      nil
+    end
 
-      #
-      # 月の入り
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options dummy
-      #
-      # @return [When::TM::TemporalPosition] 月の入りの時刻をイベント時刻とする
-      #
-      def moonset(date, options={})
-        event = formula(date.location.iri).moonset(date)
-        event.events = [@root['Moonset']]
-        event
-      rescue
-        nil
-      end
+    #
+    # 月の入り
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options dummy
+    #
+    # @return [When::TM::TemporalPosition] 月の入りの時刻をイベント時刻とする
+    #
+    def moonset(date, options={})
+      event = formula(date.location.iri).moonset(date)
+      event.events = [@root['Moonset']]
+      event
+    rescue
+      nil
+    end
 
-      # 干潮・満潮の日時
-      #
-      # @param [When::TM::TemporalPosition] date
-      # @param [Hash] options
-      # @option options [String] tide 潮汐計算方式 'Horizontal' - 地平高度基準, 'Equatorial' - 子午線通過基準(デフォルト)
-      #
-      # @return [Array<Array<Integer, When::TM::TemporalPosotion>>] 干潮・満潮の日時の Array
-      #
-      #   [Integer] +1:満潮, -1:干潮
-      #
-      #   [When::TM::TemporalPosotion] 干潮・満潮の日時
-      #
-      #
-      def tide(date, options={})
-        return nil unless @interval
-        @target  ||= When.Resource('_ep:Moon')
-        events     = @root['Tide']
-        form       = formula(date.location.iri)
-        type       = options[:tide] =~ /^horizon/i ? nil : 0
-
-        now        = +date
-        high_tides = []
-        5.times do |i|
-          high_tide = form.day_event(now + i - 2, type, @target) + @interval
-          high_tides << high_tide if high_tides.size == 0 || high_tide > high_tides[-1] + 0.5
-        end
-        tides = []
-
-        (high_tides.size-1).times do |i|
-          tides << [0, high_tides[i]]
-          tides << [1, 0.75*high_tides[i] + 0.25*high_tides[i+1]]
-          tides << [0, 0.50*high_tides[i] + 0.50*high_tides[i+1]]
-          tides << [1, 0.25*high_tides[i] + 0.75*high_tides[i+1]]
-        end
-        tides << [0, high_tides[-1]]
-
-        today = +date.floor(When::DAY)...+date.ceil(When::DAY)
-        seed  = date._attr
-        seed[:clock] ||= When::TM::Clock.local_time
-        tides.select {|x| today.include?(x[1])}.map {|x|
-          d = form._to_seed_type(x[1], seed)
-          d.events = [events[x[0]]]
-          d
-        }
-      rescue
-        nil
+    # 干潮・満潮の日時
+    #
+    # @param [When::TM::TemporalPosition] date
+    # @param [Hash] options 以下の通り
+    # @option options [String] tide 潮汐計算方式 'Horizontal' - 地平高度基準, 'Equatorial' - 子午線通過基準(デフォルト)
+    #
+    # @return [Array<Array<Integer, When::TM::TemporalPosotion>>] 干潮・満潮の日時の Array
+    #
+    #   [Integer] +1:満潮, -1:干潮
+    #
+    #   [When::TM::TemporalPosotion] 干潮・満潮の日時
+    #
+    #
+    def tide(date, options={})
+      return nil unless @interval
+      @target  ||= When.Resource('_ep:Moon')
+      events     = @root['Tide']
+      form       = formula(date.location.iri)
+      type       = options[:tide] =~ /^horizon/i ? nil : 0
+
+      now        = +date
+      high_tides = []
+      5.times do |i|
+        high_tide = form.day_event(now + i - 2, type, @target) + @interval
+        high_tides << high_tide if high_tides.size == 0 || high_tide > high_tides[-1] + 0.5
       end
+      tides = []
 
-      private
-
-      # オブジェクトの正規化
-      #   long     - 計算に用いる経度 / 度
-      #   lat      - 計算に用いる緯度 / 度
-      #   alt      - 計算に用いる高度 / m
-      #   interval - 高潮間隔(月の子午線通過から満潮までの時間) / 時間
-      def _normalize(args=[], options={})
-        if @location
-          @location = When.Resource(@location)
-          @formula  = When.Resource(@formula || "Formula?location=(#{@location.iri})", '_ep:')
-        else
-          @formula  = {}
-        end
-        @interval  = @interval.sub('@','.').to_f / 24 if @interval
-        @root      = When.CalendarNote('EphemericNote/NoteObjects::day')
-        @prime   ||= [%w(Month), %w(Sunrise Sunset Moon_Age)]
-        super
+      (high_tides.size-1).times do |i|
+        tides << [0, high_tides[i]]
+        tides << [1, 0.75*high_tides[i] + 0.25*high_tides[i+1]]
+        tides << [0, 0.50*high_tides[i] + 0.50*high_tides[i+1]]
+        tides << [1, 0.25*high_tides[i] + 0.75*high_tides[i+1]]
       end
+      tides << [0, high_tides[-1]]
+
+      today = +date.floor(When::DAY)...+date.ceil(When::DAY)
+      seed  = date._attr
+      seed[:clock] ||= When::TM::Clock.local_time
+      tides.select {|x| today.include?(x[1])}.map {|x|
+        d = form._to_seed_type(x[1], seed)
+        d.events = [events[x[0]]]
+        d
+      }
+    rescue
+      nil
+    end
 
-      # 計算に用いる Ephemeris
-      def formula(location)
-        return @formula unless @formula.kind_of?(Hash)
-        @formula[location] ||= When.Resource("Formula?location=(#{location})", '_ep:')
+    private
+
+    # オブジェクトの正規化
+    #   long     - 計算に用いる経度 / 度
+    #   lat      - 計算に用いる緯度 / 度
+    #   alt      - 計算に用いる高度 / m
+    #   interval - 高潮間隔(月の子午線通過から満潮までの時間) / 時間
+    def _normalize(args=[], options={})
+      if @location
+        @location = When.Resource(@location)
+        @formula  = When.Resource(@formula || "Formula?location=(#{@location.iri})", '_ep:')
+      else
+        @formula  = {}
       end
+      @interval  = @interval.sub('@','.').to_f / 24 if @interval
+      @root      = When.CalendarNote('EphemericNote/NoteObjects::day')
+      @prime   ||= [%w(Month), %w(Sunrise Sunset Moon_Age)]
+      super
+    end
+
+    # 計算に用いる Ephemeris
+    def formula(location)
+      return @formula unless @formula.kind_of?(Hash)
+      @formula[location] ||= When.Resource("Formula?location=(#{location})", '_ep:')
     end
   end
 end