lunation 0.1.4 → 0.1.6
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- checksums.yaml +4 -4
- data/.rubocop.yml +18 -0
- data/CHANGELOG.md +12 -0
- data/config/periodic_terms_earth_nutation.yml +63 -567
- data/config/periodic_terms_moon_latitude.yml +60 -300
- data/config/periodic_terms_moon_longitude_distance.yml +60 -360
- data/lib/lunation/calculation/angle.rb +9 -6
- data/lib/lunation/calculation/earth_position_vsop87.rb +27 -14
- data/lib/lunation/calculation/moon_illuminated_fraction.rb +3 -2
- data/lib/lunation/calculation/moon_position.rb +47 -43
- data/lib/lunation/calculation/nutation_and_obliquity.rb +19 -18
- data/lib/lunation/calculation/sun_position.rb +6 -3
- data/lib/lunation/calculation/timekeeping.rb +14 -13
- data/lib/lunation/calculation.rb +5 -1
- data/lib/lunation/version.rb +1 -1
- metadata +3 -3
|
@@ -4,12 +4,14 @@ module Lunation
|
|
|
4
4
|
class Calculation
|
|
5
5
|
module MoonPosition
|
|
6
6
|
# rubocop:disable Layout/LineLength
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
7
|
+
MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_longitude_distance.yml", __dir__).freeze
|
|
8
|
+
MOON_LATITUDE_PERIOD_TERMS_PATH = File.expand_path("../../../config/periodic_terms_moon_latitude.yml", __dir__).freeze
|
|
9
|
+
MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS = YAML.load_file(MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS_PATH).freeze
|
|
10
|
+
MOON_LATITUDE_PERIOD_TERMS = YAML.load_file(MOON_LATITUDE_PERIOD_TERMS_PATH).freeze
|
|
11
|
+
MOON_MEAN_LONGITUDE_CONSTANTS = [218.3164477, 481_267.88123421, -0.0015786, 1.fdiv(538_841), -1.fdiv(65_194_000)].freeze
|
|
12
|
+
MOON_MEAN_ELONGATION_CONSTANTS = [297.8501921, 445_267.1114034, -0.0018819, 1.fdiv(545_868), -1.fdiv(113_065_000)].freeze
|
|
13
|
+
MOON_MEAN_ANOMALY_CONSTANTS = [134.9633964, 477_198.8675055, 0.0087414, 1.fdiv(69_699), -1.fdiv(14_712_000)].freeze
|
|
14
|
+
MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.2720950, 483_202.0175233, -0.0036539, -1.fdiv(3_526_000), 1.fdiv(863_310_000)].freeze
|
|
13
15
|
CORRECTION_ECCENTRICITY_OF_EARTH_CONSTANTS = [1, -0.002516, -0.0000074].freeze
|
|
14
16
|
FIXED_DISTANCE_BETWEEN_EARTH_AND_MOON = 385_000.56
|
|
15
17
|
RADIUS_OF_EARTH = 6378.14
|
|
@@ -75,22 +77,22 @@ module Lunation
|
|
|
75
77
|
# UNIT: decimal degrees
|
|
76
78
|
def moon_heliocentric_longitude
|
|
77
79
|
@moon_heliocentric_longitude ||= begin
|
|
78
|
-
result =
|
|
80
|
+
result = MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS.inject(0.0) do |acc, elem|
|
|
79
81
|
sine_argument = Angle.from_decimal_degrees(
|
|
80
|
-
elem[
|
|
81
|
-
elem[
|
|
82
|
-
elem[
|
|
83
|
-
elem[
|
|
82
|
+
elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
|
|
83
|
+
elem[1] * sun_mean_anomaly2.decimal_degrees +
|
|
84
|
+
elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
|
|
85
|
+
elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
|
|
84
86
|
)
|
|
85
87
|
|
|
86
|
-
if elem[
|
|
88
|
+
if elem[4].nil?
|
|
87
89
|
next acc
|
|
88
|
-
elsif [1, -1].include?(elem[
|
|
89
|
-
acc + elem[
|
|
90
|
-
elsif [-2, 2].include?(elem[
|
|
91
|
-
acc + elem[
|
|
90
|
+
elsif [1, -1].include?(elem[1])
|
|
91
|
+
acc + elem[4] * correction_eccentricity_of_earth * sine_argument.sin
|
|
92
|
+
elsif [-2, 2].include?(elem[1])
|
|
93
|
+
acc + elem[4] * correction_eccentricity_of_earth**2 * sine_argument.sin
|
|
92
94
|
else
|
|
93
|
-
acc + elem[
|
|
95
|
+
acc + elem[4] * sine_argument.sin
|
|
94
96
|
end
|
|
95
97
|
end + 3958 * correction_venus.sin +
|
|
96
98
|
1962 * (moon_mean_longitude - moon_argument_of_latitude_high_precision).sin +
|
|
@@ -103,22 +105,22 @@ module Lunation
|
|
|
103
105
|
# UNIT: decimal degrees
|
|
104
106
|
def moon_heliocentric_latitude
|
|
105
107
|
@moon_heliocentric_latitude ||= begin
|
|
106
|
-
result =
|
|
108
|
+
result = MOON_LATITUDE_PERIOD_TERMS.inject(0.0) do |acc, elem|
|
|
107
109
|
sine_argument = Angle.from_decimal_degrees(
|
|
108
|
-
elem[
|
|
109
|
-
elem[
|
|
110
|
-
elem[
|
|
111
|
-
elem[
|
|
110
|
+
elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
|
|
111
|
+
elem[1] * sun_mean_anomaly2.decimal_degrees +
|
|
112
|
+
elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
|
|
113
|
+
elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
|
|
112
114
|
)
|
|
113
115
|
|
|
114
|
-
if elem[
|
|
116
|
+
if elem[4].nil?
|
|
115
117
|
next acc
|
|
116
|
-
elsif [1, -1].include?(elem[
|
|
117
|
-
acc + elem[
|
|
118
|
-
elsif [-2, 2].include?(elem[
|
|
119
|
-
acc + elem[
|
|
118
|
+
elsif [1, -1].include?(elem[1])
|
|
119
|
+
acc + elem[4] * correction_eccentricity_of_earth * sine_argument.sin
|
|
120
|
+
elsif [-2, 2].include?(elem[1])
|
|
121
|
+
acc + elem[4] * correction_eccentricity_of_earth**2 * sine_argument.sin
|
|
120
122
|
else
|
|
121
|
-
acc + elem[
|
|
123
|
+
acc + elem[4] * sine_argument.sin
|
|
122
124
|
end
|
|
123
125
|
end - 2235 * moon_mean_longitude.sin +
|
|
124
126
|
382 * correction_latitude.sin +
|
|
@@ -134,22 +136,22 @@ module Lunation
|
|
|
134
136
|
# UNIT: 1000 km
|
|
135
137
|
def moon_heliocentric_distance
|
|
136
138
|
@moon_heliocentric_distance ||= begin
|
|
137
|
-
result =
|
|
139
|
+
result = MOON_LONGITUDE_AND_DISTANCE_PERIOD_TERMS.inject(0.0) do |acc, elem|
|
|
138
140
|
cosine_argument = Angle.from_decimal_degrees(
|
|
139
|
-
elem[
|
|
140
|
-
elem[
|
|
141
|
-
elem[
|
|
142
|
-
elem[
|
|
141
|
+
elem[0] * moon_mean_elongation_from_sun_high_precision.decimal_degrees +
|
|
142
|
+
elem[1] * sun_mean_anomaly2.decimal_degrees +
|
|
143
|
+
elem[2] * moon_mean_anomaly_high_precision.decimal_degrees +
|
|
144
|
+
elem[3] * moon_argument_of_latitude_high_precision.decimal_degrees
|
|
143
145
|
)
|
|
144
146
|
|
|
145
|
-
if elem[
|
|
147
|
+
if elem[5].nil?
|
|
146
148
|
next acc
|
|
147
|
-
elsif [1, -1].include?(elem[
|
|
148
|
-
acc + elem[
|
|
149
|
-
elsif [-2, 2].include?(elem[
|
|
150
|
-
acc + elem[
|
|
149
|
+
elsif [1, -1].include?(elem[1])
|
|
150
|
+
acc + elem[5] * correction_eccentricity_of_earth * cosine_argument.cos
|
|
151
|
+
elsif [-2, 2].include?(elem[1])
|
|
152
|
+
acc + elem[5] * correction_eccentricity_of_earth**2 * cosine_argument.cos
|
|
151
153
|
else
|
|
152
|
-
acc + elem[
|
|
154
|
+
acc + elem[5] * cosine_argument.cos
|
|
153
155
|
end
|
|
154
156
|
end
|
|
155
157
|
result.round
|
|
@@ -160,7 +162,8 @@ module Lunation
|
|
|
160
162
|
# UNIT: Angle
|
|
161
163
|
def moon_ecliptic_longitude
|
|
162
164
|
@moon_ecliptic_longitude ||= Angle.from_decimal_degrees(
|
|
163
|
-
moon_mean_longitude.decimal_degrees +
|
|
165
|
+
moon_mean_longitude.decimal_degrees +
|
|
166
|
+
moon_heliocentric_longitude.fdiv(1_000_000)
|
|
164
167
|
)
|
|
165
168
|
end
|
|
166
169
|
|
|
@@ -168,7 +171,7 @@ module Lunation
|
|
|
168
171
|
# UNIT: Angle
|
|
169
172
|
def moon_ecliptic_latitude
|
|
170
173
|
@moon_ecliptic_latitude ||= Angle.from_decimal_degrees(
|
|
171
|
-
moon_heliocentric_latitude
|
|
174
|
+
moon_heliocentric_latitude.fdiv(1_000_000),
|
|
172
175
|
normalize: false
|
|
173
176
|
)
|
|
174
177
|
end
|
|
@@ -178,7 +181,7 @@ module Lunation
|
|
|
178
181
|
def distance_between_earth_and_moon
|
|
179
182
|
@distance_between_earth_and_moon ||= begin
|
|
180
183
|
result = FIXED_DISTANCE_BETWEEN_EARTH_AND_MOON +
|
|
181
|
-
(
|
|
184
|
+
moon_heliocentric_distance.fdiv(1_000)
|
|
182
185
|
result.round(1)
|
|
183
186
|
end
|
|
184
187
|
end
|
|
@@ -194,7 +197,8 @@ module Lunation
|
|
|
194
197
|
# (apparent λ) Moon apparent longitude (A.A. p. 343)
|
|
195
198
|
# UNIT: Angle
|
|
196
199
|
def moon_apparent_ecliptic_longitude
|
|
197
|
-
@moon_apparent_ecliptic_longitude ||=
|
|
200
|
+
@moon_apparent_ecliptic_longitude ||=
|
|
201
|
+
moon_ecliptic_longitude + nutation_in_longitude
|
|
198
202
|
end
|
|
199
203
|
|
|
200
204
|
# (α) geocentric (apparent) right ascension of the moon (13.3 A.A. p. 93)
|
|
@@ -2,13 +2,14 @@ module Lunation
|
|
|
2
2
|
class Calculation
|
|
3
3
|
module NutationAndObliquity
|
|
4
4
|
# rubocop:disable Layout/LineLength
|
|
5
|
-
|
|
5
|
+
NUTATION_IN_OBLIQUITY_PERIODIC_TERMS_PATH = File.expand_path("../../../config/periodic_terms_earth_nutation.yml", __dir__)
|
|
6
|
+
NUTATION_IN_OBLIQUITY_PERIODIC_TERMS = YAML.load_file(NUTATION_IN_OBLIQUITY_PERIODIC_TERMS_PATH).freeze
|
|
6
7
|
ECLIPTIC_MEAN_OBLIQUITY_CONSTANTS = [21.448, -4680.93, -1.55, 1_999.25, -51.38, -249.67, -39.05, 7.12, 27.87, 5.79, 2.45].freeze
|
|
7
|
-
MOON_MEAN_ELONGATION_FROM_SUN_CONSTANTS = [297.85036, 445_267.111480, -0.0019142, 1
|
|
8
|
-
SUN_MEAN_ANOMALY_CONSTANTS = [357.52772, 35_999.050340, -0.0001603, -1
|
|
9
|
-
MOON_MEAN_ANOMALY_CONSTANTS = [134.96298, 477_198.867398, 0.0086972, 1
|
|
10
|
-
MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.27191, 483_202.017538, -0.0036825, 1
|
|
11
|
-
LONGITUDE_OF_THE_ASCENDING_NODE_CONSTANTS = [125.04452, -1934.136261, 0.0020708, 1
|
|
8
|
+
MOON_MEAN_ELONGATION_FROM_SUN_CONSTANTS = [297.85036, 445_267.111480, -0.0019142, 1.fdiv(189_474)].freeze
|
|
9
|
+
SUN_MEAN_ANOMALY_CONSTANTS = [357.52772, 35_999.050340, -0.0001603, -1.fdiv(300_000)].freeze
|
|
10
|
+
MOON_MEAN_ANOMALY_CONSTANTS = [134.96298, 477_198.867398, 0.0086972, 1.fdiv(56_250)].freeze
|
|
11
|
+
MOON_ARGUMENT_OF_LATITUDE_CONSTANTS = [93.27191, 483_202.017538, -0.0036825, 1.fdiv(327_270)].freeze
|
|
12
|
+
LONGITUDE_OF_THE_ASCENDING_NODE_CONSTANTS = [125.04452, -1934.136261, 0.0020708, 1.fdiv(450_000)].freeze
|
|
12
13
|
# rubocop:enable Layout/LineLength
|
|
13
14
|
|
|
14
15
|
# (D) Mean elongation of the moon from the sun (A.A. p. 144)
|
|
@@ -63,13 +64,13 @@ module Lunation
|
|
|
63
64
|
@nutation_in_longitude ||= begin
|
|
64
65
|
result = NUTATION_IN_OBLIQUITY_PERIODIC_TERMS.inject(0.0) do |acc, elem|
|
|
65
66
|
argument = Angle.from_decimal_degrees(
|
|
66
|
-
elem[
|
|
67
|
-
elem[
|
|
68
|
-
elem[
|
|
69
|
-
elem[
|
|
70
|
-
elem[
|
|
67
|
+
elem[0] * moon_mean_elongation_from_sun.decimal_degrees +
|
|
68
|
+
elem[1] * sun_mean_anomaly.decimal_degrees +
|
|
69
|
+
elem[2] * moon_mean_anomaly.decimal_degrees +
|
|
70
|
+
elem[3] * moon_argument_of_latitude.decimal_degrees +
|
|
71
|
+
elem[4] * longitude_of_ascending_node.decimal_degrees
|
|
71
72
|
)
|
|
72
|
-
coefficient = elem[
|
|
73
|
+
coefficient = elem[5] + elem[6] * time
|
|
73
74
|
acc + coefficient * argument.sin
|
|
74
75
|
end / 10_000.0
|
|
75
76
|
Angle.from_decimal_arcseconds(result)
|
|
@@ -82,13 +83,13 @@ module Lunation
|
|
|
82
83
|
@nutation_in_obliquity ||= begin
|
|
83
84
|
result = NUTATION_IN_OBLIQUITY_PERIODIC_TERMS.inject(0.0) do |acc, elem|
|
|
84
85
|
argument = Angle.from_decimal_degrees(
|
|
85
|
-
elem[
|
|
86
|
-
elem[
|
|
87
|
-
elem[
|
|
88
|
-
elem[
|
|
89
|
-
elem[
|
|
86
|
+
elem[0] * moon_mean_elongation_from_sun.decimal_degrees +
|
|
87
|
+
elem[1] * sun_mean_anomaly.decimal_degrees +
|
|
88
|
+
elem[2] * moon_mean_anomaly.decimal_degrees +
|
|
89
|
+
elem[3] * moon_argument_of_latitude.decimal_degrees +
|
|
90
|
+
elem[4] * longitude_of_ascending_node.decimal_degrees
|
|
90
91
|
)
|
|
91
|
-
coefficient = elem[
|
|
92
|
+
coefficient = elem[7] + elem[8] * time
|
|
92
93
|
acc + coefficient * argument.cos
|
|
93
94
|
end / 10_000.0
|
|
94
95
|
Angle.from_decimal_arcseconds(result)
|
|
@@ -58,7 +58,10 @@ module Lunation
|
|
|
58
58
|
# UNIT: Astronomical Units (AU)
|
|
59
59
|
def distance_between_earth_and_sun_in_astronomical_units
|
|
60
60
|
@distance_between_earth_and_sun_in_astronomical_units ||= begin
|
|
61
|
-
result = 1.000001018 *
|
|
61
|
+
result = 1.000001018 *
|
|
62
|
+
(1 - earth_orbit_eccentricity**2) / (
|
|
63
|
+
(1 + earth_orbit_eccentricity * sun_anomaly.cos)
|
|
64
|
+
)
|
|
62
65
|
result.round(7)
|
|
63
66
|
end
|
|
64
67
|
end
|
|
@@ -70,8 +73,8 @@ module Lunation
|
|
|
70
73
|
(distance_between_earth_and_sun_in_astronomical_units * 149_597_870).floor
|
|
71
74
|
end
|
|
72
75
|
|
|
73
|
-
# (Ω) Longitude of the ascending node of the Moon's mean orbit on the ecliptic
|
|
74
|
-
# A.A. p. 164
|
|
76
|
+
# (Ω) Longitude of the ascending node of the Moon's mean orbit on the ecliptic
|
|
77
|
+
# (low precision) A.A. p. 164
|
|
75
78
|
# UNIT: Angle
|
|
76
79
|
def longitude_of_ascending_node_low_precision
|
|
77
80
|
@longitude_of_ascending_node_low_precision ||=
|
|
@@ -7,14 +7,14 @@ module Lunation
|
|
|
7
7
|
# rubocop:disable Layout/LineLength
|
|
8
8
|
DELTA_T_500BC_500AD = [10_583.6, -1_014.41, 33.78311, -5.952053, -0.1798452, 0.022174192, 0.0090316521].freeze
|
|
9
9
|
DELTA_T_500AD_1600AD = [1_574.2, -556.01, 71.23472, 0.319781, -0.8503463, -0.005050998, 0.0083572073].freeze
|
|
10
|
-
DELTA_T_1600AD_1700AD = [120, -0.9808, -0.01532, 1
|
|
11
|
-
DELTA_T_1700AD_1800AD = [8.83, 0.1603, -0.0059285, 0.00013336, -1
|
|
10
|
+
DELTA_T_1600AD_1700AD = [120, -0.9808, -0.01532, 1.fdiv(7_129)].freeze
|
|
11
|
+
DELTA_T_1700AD_1800AD = [8.83, 0.1603, -0.0059285, 0.00013336, -1.fdiv(1_174_000)].freeze
|
|
12
12
|
DELTA_T_1800AD_1860AD = [13.72, -0.332447, 0.0068612, 0.0041116, -0.00037436, 0.0000121272, -0.0000001699, 0.000000000875].freeze
|
|
13
|
-
DELTA_T_1860AD_1900AD = [7.62, 0.5737, -0.251754, 0.01680668, -0.0004473624, 1
|
|
13
|
+
DELTA_T_1860AD_1900AD = [7.62, 0.5737, -0.251754, 0.01680668, -0.0004473624, 1.fdiv(233_174)].freeze
|
|
14
14
|
DELTA_T_1900AD_1920AD = [-2.79, 1.494119, -0.0598939, 0.0061966, -0.000197].freeze
|
|
15
15
|
DELTA_T_1920AD_1941AD = [21.20, 0.84493, -0.076100, 0.0020936].freeze
|
|
16
|
-
DELTA_T_1941AD_1961AD = [29.07, 0.407, -1
|
|
17
|
-
DELTA_T_1961AD_1986AD = [45.45, 1.067, -1
|
|
16
|
+
DELTA_T_1941AD_1961AD = [29.07, 0.407, -1.fdiv(233), 1.fdiv(2_547)].freeze
|
|
17
|
+
DELTA_T_1961AD_1986AD = [45.45, 1.067, -1.fdiv(260), -1.fdiv(718)].freeze
|
|
18
18
|
DELTA_T_1986AD_2005AD = [63.86, 0.3345, -0.060374, 0.0017275, 0.000651814, 0.00002373599].freeze
|
|
19
19
|
DELTA_T_2005AD_2050AD = [62.92, 0.32217, 0.005589].freeze
|
|
20
20
|
# rubocop:enable Layout/LineLength
|
|
@@ -23,7 +23,7 @@ module Lunation
|
|
|
23
23
|
# UNIT: centuries from the Epoch J2000.0
|
|
24
24
|
def time
|
|
25
25
|
@time ||=
|
|
26
|
-
(
|
|
26
|
+
(julian_ephemeris_day - EPOCH_J2000_JDE).fdiv(JULIAN_CENTURY_IN_DAYS).round(12)
|
|
27
27
|
end
|
|
28
28
|
|
|
29
29
|
# ΔT Difference between Dynamical Time (TD) and Universal Time (UT, more commonly
|
|
@@ -36,12 +36,12 @@ module Lunation
|
|
|
36
36
|
@delta_t ||= begin
|
|
37
37
|
result = case datetime.year
|
|
38
38
|
when -1_999...-500
|
|
39
|
-
elapsed_years = (@decimal_year - 1_820)
|
|
39
|
+
elapsed_years = (@decimal_year - 1_820).fdiv(100)
|
|
40
40
|
-20 + 32 * elapsed_years**2
|
|
41
41
|
when -500...500
|
|
42
|
-
Horner.compute(@decimal_year
|
|
42
|
+
Horner.compute(@decimal_year.fdiv(100), DELTA_T_500BC_500AD)
|
|
43
43
|
when 500...1_600
|
|
44
|
-
Horner.compute((@decimal_year - 1_000.0)
|
|
44
|
+
Horner.compute((@decimal_year - 1_000.0).fdiv(100), DELTA_T_500AD_1600AD)
|
|
45
45
|
when 1_600...1_700
|
|
46
46
|
Horner.compute(@decimal_year - 1_600.0, DELTA_T_1600AD_1700AD)
|
|
47
47
|
when 1_700...1_800
|
|
@@ -63,7 +63,8 @@ module Lunation
|
|
|
63
63
|
when 2_005...2_050
|
|
64
64
|
Horner.compute(@decimal_year - 2_000, DELTA_T_2005AD_2050AD)
|
|
65
65
|
when 2_050...2_150
|
|
66
|
-
-20 + 32 * ((@decimal_year - 1_820) / 100)**2 - 0.5628 *
|
|
66
|
+
-20 + 32 * ((@decimal_year - 1_820) / 100)**2 - 0.5628 *
|
|
67
|
+
(2_150 - @decimal_year)
|
|
67
68
|
when 2_150..3_000
|
|
68
69
|
-20 + 32 * @decimal_year**2
|
|
69
70
|
end
|
|
@@ -74,7 +75,7 @@ module Lunation
|
|
|
74
75
|
# (TD) Dynamical time (A.A. p. 77)
|
|
75
76
|
# UNIT: ISO 8601 Date and time with offset
|
|
76
77
|
def dynamical_time
|
|
77
|
-
@dynamical_time ||= datetime + delta_t.round.
|
|
78
|
+
@dynamical_time ||= datetime + delta_t.round.fdiv(SECONDS_PER_DAY)
|
|
78
79
|
end
|
|
79
80
|
|
|
80
81
|
# (JDE) Julian ephemeris day (A.A. p. 59)
|
|
@@ -86,13 +87,13 @@ module Lunation
|
|
|
86
87
|
# (t) time, measured in Julian millennia from the epoch J2000.0 (32.1, A.A. p. 218)
|
|
87
88
|
# UNIT: millennia from the Epoch J2000.0
|
|
88
89
|
def time_millennia
|
|
89
|
-
@time_millennia ||= (
|
|
90
|
+
@time_millennia ||= time.fdiv(10).round(12)
|
|
90
91
|
end
|
|
91
92
|
|
|
92
93
|
# (U) Time measured in units of 10_000 Julian years from J2000.0 (A.A. p. 147)
|
|
93
94
|
# UNIT: 10_000 years (myriads) from the Epoch J2000.0
|
|
94
95
|
def time_myriads
|
|
95
|
-
@time_myriads ||= time
|
|
96
|
+
@time_myriads ||= time.fdiv(100)
|
|
96
97
|
end
|
|
97
98
|
end
|
|
98
99
|
end
|
data/lib/lunation/calculation.rb
CHANGED
|
@@ -25,11 +25,14 @@ module Lunation
|
|
|
25
25
|
|
|
26
26
|
def to_h
|
|
27
27
|
{
|
|
28
|
+
# rubocop:disable Layout/LineLength
|
|
28
29
|
corrected_obliquity_of_ecliptic: corrected_obliquity_of_ecliptic,
|
|
29
30
|
correction_eccentricity_of_earth: correction_eccentricity_of_earth,
|
|
30
31
|
correction_jupiter: correction_jupiter,
|
|
31
32
|
correction_latitude: correction_latitude,
|
|
32
33
|
correction_venus: correction_venus,
|
|
34
|
+
datetime: datetime,
|
|
35
|
+
decimal_year: decimal_year,
|
|
33
36
|
delta_t: delta_t,
|
|
34
37
|
distance_between_earth_and_moon: distance_between_earth_and_moon,
|
|
35
38
|
distance_between_earth_and_sun_in_astronomical_units: distance_between_earth_and_sun_in_astronomical_units,
|
|
@@ -78,11 +81,12 @@ module Lunation
|
|
|
78
81
|
time_millennia: time_millennia,
|
|
79
82
|
time_myriads: time_myriads,
|
|
80
83
|
time: time,
|
|
84
|
+
# rubocop:enable Layout/LineLength
|
|
81
85
|
}
|
|
82
86
|
end
|
|
83
87
|
|
|
84
88
|
def to_s
|
|
85
|
-
|
|
89
|
+
<<~HEREDOC
|
|
86
90
|
DATE AND TIME (UT): #{datetime}
|
|
87
91
|
----------------------------------------------------------------------------------
|
|
88
92
|
DECIMAL YEAR: #{decimal_year}
|
data/lib/lunation/version.rb
CHANGED
metadata
CHANGED
|
@@ -1,14 +1,14 @@
|
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
|
2
2
|
name: lunation
|
|
3
3
|
version: !ruby/object:Gem::Version
|
|
4
|
-
version: 0.1.
|
|
4
|
+
version: 0.1.6
|
|
5
5
|
platform: ruby
|
|
6
6
|
authors:
|
|
7
7
|
- Ivo Kalverboer
|
|
8
8
|
autorequire:
|
|
9
9
|
bindir: exe
|
|
10
10
|
cert_chain: []
|
|
11
|
-
date: 2024-12-
|
|
11
|
+
date: 2024-12-31 00:00:00.000000000 Z
|
|
12
12
|
dependencies: []
|
|
13
13
|
description: Lunation offers a Ruby implementation of Meeus's Astronomical Algorithms.
|
|
14
14
|
email:
|
|
@@ -59,7 +59,7 @@ licenses:
|
|
|
59
59
|
metadata:
|
|
60
60
|
homepage_uri: https://www.github.com/valerius/lunation
|
|
61
61
|
source_code_uri: https://www.github.com/valerius/lunation
|
|
62
|
-
changelog_uri: https://
|
|
62
|
+
changelog_uri: https://github.com/valerius/lunation/blob/main/CHANGELOG.md
|
|
63
63
|
post_install_message:
|
|
64
64
|
rdoc_options: []
|
|
65
65
|
require_paths:
|