equationoftime 3.0.0 → 4.0.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (251) hide show
  1. checksums.yaml +4 -4
  2. data/.buildpath +5 -0
  3. data/.project +17 -0
  4. data/.rvmrc +1 -0
  5. data/.settings/org.eclipse.ltk.core.refactoring.prefs +2 -0
  6. data/Makefile +238 -0
  7. data/Rakefile +40 -3
  8. data/analemma_data.yml +2925 -0
  9. data/equationoftime.gemspec +6 -3
  10. data/examples/analemma_data_generator.rb +12 -12
  11. data/examples/celes_parts/.RUBYARCHDIR.time +0 -0
  12. data/examples/celes_parts/.cproject +56 -0
  13. data/examples/celes_parts/.project +78 -0
  14. data/examples/celes_parts/a2af.c +120 -0
  15. data/examples/celes_parts/a2tf.c +116 -0
  16. data/examples/celes_parts/af2a.c +107 -0
  17. data/examples/celes_parts/anp.c +82 -0
  18. data/examples/celes_parts/anpm.c +82 -0
  19. data/examples/celes_parts/bi00.c +116 -0
  20. data/examples/celes_parts/bp00.c +173 -0
  21. data/examples/celes_parts/bp06.c +136 -0
  22. data/examples/celes_parts/bpn2xy.c +100 -0
  23. data/examples/celes_parts/c2i00a.c +139 -0
  24. data/examples/celes_parts/c2i00b.c +139 -0
  25. data/examples/celes_parts/c2i06a.c +136 -0
  26. data/examples/celes_parts/c2ibpn.c +142 -0
  27. data/examples/celes_parts/c2ixy.c +131 -0
  28. data/examples/celes_parts/c2ixys.c +123 -0
  29. data/examples/celes_parts/c2s.c +96 -0
  30. data/examples/celes_parts/c2t00a.c +154 -0
  31. data/examples/celes_parts/c2t00b.c +150 -0
  32. data/examples/celes_parts/c2t06a.c +152 -0
  33. data/examples/celes_parts/c2tcio.c +122 -0
  34. data/examples/celes_parts/c2teqx.c +122 -0
  35. data/examples/celes_parts/c2tpe.c +167 -0
  36. data/examples/celes_parts/c2txy.c +159 -0
  37. data/examples/celes_parts/cal2jd.c +139 -0
  38. data/examples/celes_parts/celes_core.c +2522 -0
  39. data/examples/celes_parts/cp.c +80 -0
  40. data/examples/celes_parts/cpv.c +82 -0
  41. data/examples/celes_parts/cr.c +83 -0
  42. data/examples/celes_parts/d2dtf.c +206 -0
  43. data/examples/celes_parts/d2tf.c +160 -0
  44. data/examples/celes_parts/dat.c +289 -0
  45. data/examples/celes_parts/depend +25 -0
  46. data/examples/celes_parts/dtdb.c +1213 -0
  47. data/examples/celes_parts/dtf2d.c +196 -0
  48. data/examples/celes_parts/ee00.c +128 -0
  49. data/examples/celes_parts/ee00a.c +135 -0
  50. data/examples/celes_parts/ee00b.c +141 -0
  51. data/examples/celes_parts/ee06a.c +122 -0
  52. data/examples/celes_parts/eect00.c +282 -0
  53. data/examples/celes_parts/eform.c +147 -0
  54. data/examples/celes_parts/eo06a.c +131 -0
  55. data/examples/celes_parts/eors.c +108 -0
  56. data/examples/celes_parts/epb.c +96 -0
  57. data/examples/celes_parts/epb2jd.c +91 -0
  58. data/examples/celes_parts/epj.c +93 -0
  59. data/examples/celes_parts/epj2jd.c +91 -0
  60. data/examples/celes_parts/epv00.c +2589 -0
  61. data/examples/celes_parts/eqeq94.c +131 -0
  62. data/examples/celes_parts/era00.c +136 -0
  63. data/examples/celes_parts/extconf.rb +4 -0
  64. data/examples/celes_parts/fad03.c +103 -0
  65. data/examples/celes_parts/fae03.c +102 -0
  66. data/examples/celes_parts/faf03.c +106 -0
  67. data/examples/celes_parts/faju03.c +102 -0
  68. data/examples/celes_parts/fal03.c +103 -0
  69. data/examples/celes_parts/falp03.c +103 -0
  70. data/examples/celes_parts/fama03.c +102 -0
  71. data/examples/celes_parts/fame03.c +102 -0
  72. data/examples/celes_parts/fane03.c +99 -0
  73. data/examples/celes_parts/faom03.c +104 -0
  74. data/examples/celes_parts/fapa03.c +103 -0
  75. data/examples/celes_parts/fasa03.c +102 -0
  76. data/examples/celes_parts/faur03.c +99 -0
  77. data/examples/celes_parts/fave03.c +102 -0
  78. data/examples/celes_parts/fk52h.c +143 -0
  79. data/examples/celes_parts/fk5hip.c +126 -0
  80. data/examples/celes_parts/fk5hz.c +160 -0
  81. data/examples/celes_parts/fw2m.c +134 -0
  82. data/examples/celes_parts/fw2xy.c +120 -0
  83. data/examples/celes_parts/gc2gd.c +134 -0
  84. data/examples/celes_parts/gc2gde.c +200 -0
  85. data/examples/celes_parts/gd2gc.c +134 -0
  86. data/examples/celes_parts/gd2gce.c +138 -0
  87. data/examples/celes_parts/gmst00.c +145 -0
  88. data/examples/celes_parts/gmst06.c +136 -0
  89. data/examples/celes_parts/gmst82.c +151 -0
  90. data/examples/celes_parts/gst00a.c +138 -0
  91. data/examples/celes_parts/gst00b.c +146 -0
  92. data/examples/celes_parts/gst06.c +140 -0
  93. data/examples/celes_parts/gst06a.c +131 -0
  94. data/examples/celes_parts/gst94.c +131 -0
  95. data/examples/celes_parts/h2fk5.c +148 -0
  96. data/examples/celes_parts/hfk5z.c +175 -0
  97. data/examples/celes_parts/ir.c +83 -0
  98. data/examples/celes_parts/jd2cal.c +155 -0
  99. data/examples/celes_parts/jdcalf.c +161 -0
  100. data/examples/celes_parts/num00a.c +121 -0
  101. data/examples/celes_parts/num00b.c +121 -0
  102. data/examples/celes_parts/num06a.c +125 -0
  103. data/examples/celes_parts/numat.c +109 -0
  104. data/examples/celes_parts/nut00a.c +2047 -0
  105. data/examples/celes_parts/nut00b.c +372 -0
  106. data/examples/celes_parts/nut06a.c +153 -0
  107. data/examples/celes_parts/nut80.c +325 -0
  108. data/examples/celes_parts/nutm80.c +117 -0
  109. data/examples/celes_parts/obl06.c +118 -0
  110. data/examples/celes_parts/obl80.c +118 -0
  111. data/examples/celes_parts/p06e.c +321 -0
  112. data/examples/celes_parts/p2pv.c +83 -0
  113. data/examples/celes_parts/p2s.c +91 -0
  114. data/examples/celes_parts/pap.c +139 -0
  115. data/examples/celes_parts/pas.c +96 -0
  116. data/examples/celes_parts/pb06.c +144 -0
  117. data/examples/celes_parts/pdp.c +84 -0
  118. data/examples/celes_parts/pfw06.c +165 -0
  119. data/examples/celes_parts/plan94.c +514 -0
  120. data/examples/celes_parts/pm.c +83 -0
  121. data/examples/celes_parts/pmat00.c +118 -0
  122. data/examples/celes_parts/pmat06.c +122 -0
  123. data/examples/celes_parts/pmat76.c +141 -0
  124. data/examples/celes_parts/pmp.c +85 -0
  125. data/examples/celes_parts/pn.c +109 -0
  126. data/examples/celes_parts/pn00.c +177 -0
  127. data/examples/celes_parts/pn00a.c +162 -0
  128. data/examples/celes_parts/pn00b.c +162 -0
  129. data/examples/celes_parts/pn06.c +187 -0
  130. data/examples/celes_parts/pn06a.c +152 -0
  131. data/examples/celes_parts/pnm00a.c +121 -0
  132. data/examples/celes_parts/pnm00b.c +121 -0
  133. data/examples/celes_parts/pnm06a.c +124 -0
  134. data/examples/celes_parts/pnm80.c +126 -0
  135. data/examples/celes_parts/pom00.c +115 -0
  136. data/examples/celes_parts/ppp.c +85 -0
  137. data/examples/celes_parts/ppsp.c +94 -0
  138. data/examples/celes_parts/pr00.c +142 -0
  139. data/examples/celes_parts/prec76.c +148 -0
  140. data/examples/celes_parts/pv2p.c +81 -0
  141. data/examples/celes_parts/pv2s.c +144 -0
  142. data/examples/celes_parts/pvdpv.c +102 -0
  143. data/examples/celes_parts/pvm.c +86 -0
  144. data/examples/celes_parts/pvmpv.c +87 -0
  145. data/examples/celes_parts/pvppv.c +87 -0
  146. data/examples/celes_parts/pvstar.c +207 -0
  147. data/examples/celes_parts/pvu.c +93 -0
  148. data/examples/celes_parts/pvup.c +88 -0
  149. data/examples/celes_parts/pvxpv.c +107 -0
  150. data/examples/celes_parts/pxp.c +94 -0
  151. data/examples/celes_parts/rm2v.c +111 -0
  152. data/examples/celes_parts/rv2m.c +118 -0
  153. data/examples/celes_parts/rx.c +110 -0
  154. data/examples/celes_parts/rxp.c +99 -0
  155. data/examples/celes_parts/rxpv.c +86 -0
  156. data/examples/celes_parts/rxr.c +99 -0
  157. data/examples/celes_parts/ry.c +110 -0
  158. data/examples/celes_parts/rz.c +110 -0
  159. data/examples/celes_parts/s00.c +371 -0
  160. data/examples/celes_parts/s00a.c +143 -0
  161. data/examples/celes_parts/s00b.c +143 -0
  162. data/examples/celes_parts/s06.c +368 -0
  163. data/examples/celes_parts/s06a.c +145 -0
  164. data/examples/celes_parts/s2c.c +85 -0
  165. data/examples/celes_parts/s2p.c +88 -0
  166. data/examples/celes_parts/s2pv.c +103 -0
  167. data/examples/celes_parts/s2xpv.c +87 -0
  168. data/examples/celes_parts/sepp.c +105 -0
  169. data/examples/celes_parts/seps.c +93 -0
  170. data/examples/celes_parts/sofa.h +379 -0
  171. data/examples/celes_parts/sofam.h +155 -0
  172. data/examples/celes_parts/sp00.c +118 -0
  173. data/examples/celes_parts/starpm.c +205 -0
  174. data/examples/celes_parts/starpv.c +264 -0
  175. data/examples/celes_parts/sxp.c +84 -0
  176. data/examples/celes_parts/sxpv.c +85 -0
  177. data/examples/celes_parts/taitt.c +110 -0
  178. data/examples/celes_parts/taiut1.c +112 -0
  179. data/examples/celes_parts/taiutc.c +182 -0
  180. data/examples/celes_parts/tcbtdb.c +132 -0
  181. data/examples/celes_parts/tcgtt.c +109 -0
  182. data/examples/celes_parts/tdbtcb.c +137 -0
  183. data/examples/celes_parts/tdbtt.c +122 -0
  184. data/examples/celes_parts/test_celes.rb +48 -0
  185. data/examples/celes_parts/tf2a.c +107 -0
  186. data/examples/celes_parts/tf2d.c +107 -0
  187. data/examples/celes_parts/tr.c +93 -0
  188. data/examples/celes_parts/trxp.c +93 -0
  189. data/examples/celes_parts/trxpv.c +93 -0
  190. data/examples/celes_parts/tttai.c +110 -0
  191. data/examples/celes_parts/tttcg.c +112 -0
  192. data/examples/celes_parts/tttdb.c +121 -0
  193. data/examples/celes_parts/ttut1.c +110 -0
  194. data/examples/celes_parts/ut1tai.c +111 -0
  195. data/examples/celes_parts/ut1tt.c +110 -0
  196. data/examples/celes_parts/ut1utc.c +193 -0
  197. data/examples/celes_parts/utctai.c +163 -0
  198. data/examples/celes_parts/utcut1.c +151 -0
  199. data/examples/celes_parts/xy06.c +2758 -0
  200. data/examples/celes_parts/xys00a.c +133 -0
  201. data/examples/celes_parts/xys00b.c +133 -0
  202. data/examples/celes_parts/xys06a.c +133 -0
  203. data/examples/celes_parts/zp.c +77 -0
  204. data/examples/celes_parts/zpv.c +79 -0
  205. data/examples/celes_parts/zr.c +83 -0
  206. data/examples/check_date_type.rb +2 -4
  207. data/examples/compare_geoc_long_ra.rb +19 -11
  208. data/examples/data_table_for_astro_dog.rb +40 -0
  209. data/examples/earth_rotation.rb +28 -3
  210. data/examples/eot.c +15 -0
  211. data/examples/eot.h +2 -0
  212. data/examples/eot_methods_list.rb +30 -3
  213. data/examples/eot_suntimes.rb +99 -95
  214. data/examples/extconf.rb +2 -0
  215. data/examples/{read_nutation_data.rb → gmst_gast_non_sofa.rb} +103 -98
  216. data/examples/nutation_table5_3a.yaml +1658 -1658
  217. data/examples/rbeot.c +24 -0
  218. data/examples/test_eot.rb +17 -0
  219. data/examples/test_poly_eval.rb +6 -17
  220. data/examples/time_scales.rb +4 -7
  221. data/examples/times_year.rb +50 -0
  222. data/examples/usage_example.rb +15 -5
  223. data/examples/use_angles.rb +151 -87
  224. data/lib/eot/angles.rb +130 -177
  225. data/lib/eot/constants.rb +41 -44
  226. data/lib/eot/displays.rb +40 -37
  227. data/lib/eot/geo_lat_lng_smt.rb +7 -6
  228. data/lib/eot/init.rb +62 -39
  229. data/lib/eot/nutation.rb +53 -45
  230. data/lib/eot/times.rb +32 -41
  231. data/lib/eot/utilities.rb +11 -90
  232. data/lib/eot/version.rb +3 -2
  233. data/rise_set_data.yml +1461 -0
  234. data/run_tests_eclipse.rb +1 -0
  235. data/tests/minitest/aliased_angles_spec.rb +119 -165
  236. data/tests/minitest/aliased_displays_spec.rb +41 -39
  237. data/tests/minitest/aliased_utilities_spec.rb +11 -22
  238. data/tests/minitest/angles_spec.rb +134 -181
  239. data/tests/minitest/constants_spec.rb +4 -10
  240. data/tests/minitest/delta_epsilon_spec.rb +13 -16
  241. data/tests/minitest/displays_spec.rb +39 -37
  242. data/tests/minitest/geo_spec.rb +15 -11
  243. data/tests/minitest/init_spec.rb +30 -16
  244. data/tests/minitest/nutation_spec.rb +19 -13
  245. data/tests/minitest/times_spec.rb +77 -75
  246. data/tests/minitest/utilities_spec.rb +10 -94
  247. data/wiki.md +7 -9
  248. metadata +218 -26
  249. data/examples/data_table.rb +0 -26
  250. data/lib/eot/nutation_table5_3a.yaml +0 -9532
  251. data/tests/minitest/aliased_times_spec.rb +0 -36
@@ -0,0 +1,112 @@
1
+ #include "sofam.h"
2
+
3
+ int iauTaiut1(double tai1, double tai2, double dta,
4
+ double *ut11, double *ut12)
5
+ /*
6
+ ** - - - - - - - - - -
7
+ ** i a u T a i u t 1
8
+ ** - - - - - - - - - -
9
+ **
10
+ ** Time scale transformation: International Atomic Time, TAI, to
11
+ ** Universal Time, UT1.
12
+ **
13
+ ** Status: canonical.
14
+ **
15
+ ** Given:
16
+ ** tai1,tai2 double TAI as a 2-part Julian Date
17
+ ** dta double UT1-TAI in seconds
18
+ **
19
+ ** Returned:
20
+ ** ut11,ut12 double UT1 as a 2-part Julian Date
21
+ **
22
+ ** Returned (function value):
23
+ ** int status: 0 = OK
24
+ **
25
+ ** Notes:
26
+ **
27
+ ** 1) tai1+tai2 is Julian Date, apportioned in any convenient way
28
+ ** between the two arguments, for example where tai1 is the Julian
29
+ ** Day Number and tai2 is the fraction of a day. The returned
30
+ ** UT11,UT12 follow suit.
31
+ **
32
+ ** 2) The argument dta, i.e. UT1-TAI, is an observed quantity, and is
33
+ ** available from IERS tabulations.
34
+ **
35
+ ** Reference:
36
+ **
37
+ ** Explanatory Supplement to the Astronomical Almanac,
38
+ ** P. Kenneth Seidelmann (ed), University Science Books (1992)
39
+ **
40
+ ** This revision: 2011 May 14
41
+ **
42
+ ** Original version 2012-03-01
43
+ **
44
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
45
+ **
46
+ */
47
+ {
48
+ double dtad;
49
+
50
+
51
+ /* Result, safeguarding precision. */
52
+ dtad = dta / DAYSEC;
53
+ if ( tai1 > tai2 ) {
54
+ *ut11 = tai1;
55
+ *ut12 = tai2 + dtad;
56
+ } else {
57
+ *ut11 = tai1 + dtad;
58
+ *ut12 = tai2;
59
+ }
60
+
61
+ /* Status (always OK). */
62
+ return 0;
63
+
64
+ /*----------------------------------------------------------------------
65
+ **
66
+ ** Celes is a wrapper of the SOFA Library for Ruby.
67
+ **
68
+ ** This file is redistributed and relicensed in accordance with
69
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
70
+ **
71
+ ** The original library is available from IAU Standards of
72
+ ** Fundamental Astronomy (http://www.iausofa.org/).
73
+ **
74
+ **
75
+ **
76
+ **
77
+ **
78
+ ** Copyright (C) 2013, Naoki Arita
79
+ ** All rights reserved.
80
+ **
81
+ ** Redistribution and use in source and binary forms, with or without
82
+ ** modification, are permitted provided that the following conditions
83
+ ** are met:
84
+ **
85
+ ** 1 Redistributions of source code must retain the above copyright
86
+ ** notice, this list of conditions and the following disclaimer.
87
+ **
88
+ ** 2 Redistributions in binary form must reproduce the above copyright
89
+ ** notice, this list of conditions and the following disclaimer in
90
+ ** the documentation and/or other materials provided with the
91
+ ** distribution.
92
+ **
93
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
94
+ ** the International Astronomical Union nor the names of its
95
+ ** contributors may be used to endorse or promote products derived
96
+ ** from this software without specific prior written permission.
97
+ **
98
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
99
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
100
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
101
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
102
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
103
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
104
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
105
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
106
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
107
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
108
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
109
+ ** POSSIBILITY OF SUCH DAMAGE.
110
+ **
111
+ **--------------------------------------------------------------------*/
112
+ }
@@ -0,0 +1,182 @@
1
+ #include "sofam.h"
2
+
3
+ int iauTaiutc(double tai1, double tai2, double *utc1, double *utc2)
4
+ /*
5
+ ** - - - - - - - - - -
6
+ ** i a u T a i u t c
7
+ ** - - - - - - - - - -
8
+ **
9
+ ** Time scale transformation: International Atomic Time, TAI, to
10
+ ** Coordinated Universal Time, UTC.
11
+ **
12
+ ** Status: canonical.
13
+ **
14
+ ** Given:
15
+ ** tai1,tai2 double TAI as a 2-part Julian Date (Note 1)
16
+ **
17
+ ** Returned:
18
+ ** utc1,utc2 double UTC as a 2-part quasi Julian Date (Notes 1-3)
19
+ **
20
+ ** Returned (function value):
21
+ ** int status: +1 = dubious year (Note 4)
22
+ ** 0 = OK
23
+ ** -1 = unacceptable date
24
+ **
25
+ ** Notes:
26
+ **
27
+ ** 1) tai1+tai2 is Julian Date, apportioned in any convenient way
28
+ ** between the two arguments, for example where tai1 is the Julian
29
+ ** Day Number and tai2 is the fraction of a day. The returned utc1
30
+ ** and utc2 form an analogous pair, except that a special convention
31
+ ** is used, to deal with the problem of leap seconds - see the next
32
+ ** note.
33
+ **
34
+ ** 2) JD cannot unambiguously represent UTC during a leap second unless
35
+ ** special measures are taken. The convention in the present
36
+ ** function is that the JD day represents UTC days whether the
37
+ ** length is 86399, 86400 or 86401 SI seconds.
38
+ **
39
+ ** 3) The function iauD2dtf can be used to transform the UTC quasi-JD
40
+ ** into calendar date and clock time, including UTC leap second
41
+ ** handling.
42
+ **
43
+ ** 4) The warning status "dubious year" flags UTCs that predate the
44
+ ** introduction of the time scale and that are too far in the future
45
+ ** to be trusted. See iauDat for further details.
46
+ **
47
+ ** Called:
48
+ ** iauJd2cal JD to Gregorian calendar
49
+ ** iauDat delta(AT) = TAI-UTC
50
+ ** iauCal2jd Gregorian calendar to JD
51
+ **
52
+ ** References:
53
+ **
54
+ ** McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
55
+ ** IERS Technical Note No. 32, BKG (2004)
56
+ **
57
+ ** Explanatory Supplement to the Astronomical Almanac,
58
+ ** P. Kenneth Seidelmann (ed), University Science Books (1992)
59
+ **
60
+ ** This revision: 2011 May 14
61
+ **
62
+ ** Original version 2012-03-01
63
+ **
64
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
65
+ */
66
+ {
67
+ int big1;
68
+ int i, iy, im, id, js;
69
+ double a1, a2, d1, dats1, d2, fd, ddats, dats2, datd, as1, as2, da;
70
+
71
+
72
+ /* Put the two parts of the TAI into big-first order. */
73
+ big1 = ( tai1 >= tai2 );
74
+ if ( big1 ) {
75
+ a1 = tai1;
76
+ a2 = tai2;
77
+ } else {
78
+ a1 = tai2;
79
+ a2 = tai1;
80
+ }
81
+
82
+ /* See if the TAI can possibly be in a leap-second day. */
83
+ d1 = a1;
84
+ dats1 = 0.0;
85
+ for ( i = -1; i <= 3; i++ ) {
86
+ d2 = a2 + (double) i;
87
+ if ( iauJd2cal(d1, d2, &iy, &im, &id, &fd) ) return -1;
88
+ js = iauDat(iy, im, id, 0.0, &dats2);
89
+ if ( js < 0 ) return -1;
90
+ if ( i == -1 ) dats1 = dats2;
91
+ ddats = dats2 - dats1;
92
+ datd = dats1 / DAYSEC;
93
+ if ( fabs(ddats) >= 0.5 ) {
94
+
95
+ /* Yes. Get TAI for the start of the UTC day that */
96
+ /* ends in a leap. */
97
+ if ( iauCal2jd(iy, im, id, &d1, &d2) ) return -1;
98
+ as1 = d1;
99
+ as2 = d2 - 1.0 + datd;
100
+
101
+ /* Is the TAI after this point? */
102
+ da = a1 - as1;
103
+ da = da + ( a2 - as2 );
104
+ if ( da > 0 ) {
105
+
106
+ /* Yes: fraction of the current UTC day that has elapsed. */
107
+ fd = da * DAYSEC / ( DAYSEC + ddats );
108
+
109
+ /* Ramp TAI-UTC to bring about SOFA's JD(UTC) convention. */
110
+ datd += ddats * ( fd <= 1.0 ? fd : 1.0 ) / DAYSEC;
111
+ }
112
+
113
+ /* Done. */
114
+ break;
115
+ }
116
+ dats1 = dats2;
117
+ }
118
+
119
+ /* Subtract the (possibly adjusted) TAI-UTC from TAI to give UTC. */
120
+ a2 -= datd;
121
+
122
+ /* Return the UTC result, preserving the TAI order. */
123
+ if ( big1 ) {
124
+ *utc1 = a1;
125
+ *utc2 = a2;
126
+ } else {
127
+ *utc1 = a2;
128
+ *utc2 = a1;
129
+ }
130
+
131
+ /* Status. */
132
+ return js;
133
+
134
+ /*----------------------------------------------------------------------
135
+ **
136
+ ** Celes is a wrapper of the SOFA Library for Ruby.
137
+ **
138
+ ** This file is redistributed and relicensed in accordance with
139
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
140
+ **
141
+ ** The original library is available from IAU Standards of
142
+ ** Fundamental Astronomy (http://www.iausofa.org/).
143
+ **
144
+ **
145
+ **
146
+ **
147
+ **
148
+ ** Copyright (C) 2013, Naoki Arita
149
+ ** All rights reserved.
150
+ **
151
+ ** Redistribution and use in source and binary forms, with or without
152
+ ** modification, are permitted provided that the following conditions
153
+ ** are met:
154
+ **
155
+ ** 1 Redistributions of source code must retain the above copyright
156
+ ** notice, this list of conditions and the following disclaimer.
157
+ **
158
+ ** 2 Redistributions in binary form must reproduce the above copyright
159
+ ** notice, this list of conditions and the following disclaimer in
160
+ ** the documentation and/or other materials provided with the
161
+ ** distribution.
162
+ **
163
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
164
+ ** the International Astronomical Union nor the names of its
165
+ ** contributors may be used to endorse or promote products derived
166
+ ** from this software without specific prior written permission.
167
+ **
168
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
169
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
170
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
171
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
172
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
173
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
174
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
175
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
176
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
177
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
178
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
179
+ ** POSSIBILITY OF SUCH DAMAGE.
180
+ **
181
+ **--------------------------------------------------------------------*/
182
+ }
@@ -0,0 +1,132 @@
1
+ #include "sofam.h"
2
+
3
+ int iauTcbtdb(double tcb1, double tcb2, double *tdb1, double *tdb2)
4
+ /*
5
+ ** - - - - - - - - - -
6
+ ** i a u T c b t d b
7
+ ** - - - - - - - - - -
8
+ **
9
+ ** Time scale transformation: Barycentric Coordinate Time, TCB, to
10
+ ** Barycentric Dynamical Time, TDB.
11
+ **
12
+ ** Status: canonical.
13
+ **
14
+ ** Given:
15
+ ** tcb1,tcb2 double TCB as a 2-part Julian Date
16
+ **
17
+ ** Returned:
18
+ ** tdb1,tdb2 double TDB as a 2-part Julian Date
19
+ **
20
+ ** Returned (function value):
21
+ ** int status: 0 = OK
22
+ **
23
+ ** Notes:
24
+ **
25
+ ** 1) tcb1+tcb2 is Julian Date, apportioned in any convenient way
26
+ ** between the two arguments, for example where tcb1 is the Julian
27
+ ** Day Number and tcb2 is the fraction of a day. The returned
28
+ ** tdb1,tdb2 follow suit.
29
+ **
30
+ ** 2) The 2006 IAU General Assembly introduced a conventional linear
31
+ ** transformation between TDB and TCB. This transformation
32
+ ** compensates for the drift between TCB and terrestrial time TT,
33
+ ** and keeps TDB approximately centered on TT. Because the
34
+ ** relationship between TT and TCB depends on the adopted solar
35
+ ** system ephemeris, the degree of alignment between TDB and TT over
36
+ ** long intervals will vary according to which ephemeris is used.
37
+ ** Former definitions of TDB attempted to avoid this problem by
38
+ ** stipulating that TDB and TT should differ only by periodic
39
+ ** effects. This is a good description of the nature of the
40
+ ** relationship but eluded precise mathematical formulation. The
41
+ ** conventional linear relationship adopted in 2006 sidestepped
42
+ ** these difficulties whilst delivering a TDB that in practice was
43
+ ** consistent with values before that date.
44
+ **
45
+ ** 3) TDB is essentially the same as Teph, the time argument for the
46
+ ** JPL solar system ephemerides.
47
+ **
48
+ ** Reference:
49
+ **
50
+ ** IAU 2006 Resolution B3
51
+ **
52
+ ** This revision: 2011 May 14
53
+ **
54
+ ** Original version 2012-03-01
55
+ **
56
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
57
+ */
58
+ {
59
+
60
+ /* 1977 Jan 1 00:00:32.184 TT, as two-part JD */
61
+ static const double t77td = DJM0 + DJM77;
62
+ static const double t77tf = TTMTAI/DAYSEC;
63
+
64
+ /* TDB (days) at TAI 1977 Jan 1.0 */
65
+ static const double tdb0 = TDB0/DAYSEC;
66
+
67
+ double d;
68
+
69
+
70
+ /* Result, safeguarding precision. */
71
+ if ( tcb1 > tcb2 ) {
72
+ d = tcb1 - t77td;
73
+ *tdb1 = tcb1;
74
+ *tdb2 = tcb2 + tdb0 - ( d + ( tcb2 - t77tf ) ) * ELB;
75
+ } else {
76
+ d = tcb2 - t77td;
77
+ *tdb1 = tcb1 + tdb0 - ( d + ( tcb1 - t77tf ) ) * ELB;
78
+ *tdb2 = tcb2;
79
+ }
80
+
81
+ /* Status (always OK). */
82
+ return 0;
83
+
84
+ /*----------------------------------------------------------------------
85
+ **
86
+ ** Celes is a wrapper of the SOFA Library for Ruby.
87
+ **
88
+ ** This file is redistributed and relicensed in accordance with
89
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
90
+ **
91
+ ** The original library is available from IAU Standards of
92
+ ** Fundamental Astronomy (http://www.iausofa.org/).
93
+ **
94
+ **
95
+ **
96
+ **
97
+ **
98
+ ** Copyright (C) 2013, Naoki Arita
99
+ ** All rights reserved.
100
+ **
101
+ ** Redistribution and use in source and binary forms, with or without
102
+ ** modification, are permitted provided that the following conditions
103
+ ** are met:
104
+ **
105
+ ** 1 Redistributions of source code must retain the above copyright
106
+ ** notice, this list of conditions and the following disclaimer.
107
+ **
108
+ ** 2 Redistributions in binary form must reproduce the above copyright
109
+ ** notice, this list of conditions and the following disclaimer in
110
+ ** the documentation and/or other materials provided with the
111
+ ** distribution.
112
+ **
113
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
114
+ ** the International Astronomical Union nor the names of its
115
+ ** contributors may be used to endorse or promote products derived
116
+ ** from this software without specific prior written permission.
117
+ **
118
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
119
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
120
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
121
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
122
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
123
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
124
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
125
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
126
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
127
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
128
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
129
+ ** POSSIBILITY OF SUCH DAMAGE.
130
+ **
131
+ **--------------------------------------------------------------------*/
132
+ }
@@ -0,0 +1,109 @@
1
+ #include "sofam.h"
2
+
3
+ int iauTcgtt(double tcg1, double tcg2, double *tt1, double *tt2)
4
+ /*
5
+ ** - - - - - - - - -
6
+ ** i a u T c g t t
7
+ ** - - - - - - - - -
8
+ **
9
+ ** Time scale transformation: Geocentric Coordinate Time, TCG, to
10
+ ** Terrestrial Time, TT.
11
+ **
12
+ ** Status: canonical.
13
+ **
14
+ ** Given:
15
+ ** tcg1,tcg2 double TCG as a 2-part Julian Date
16
+ **
17
+ ** Returned:
18
+ ** tt1,tt2 double TT as a 2-part Julian Date
19
+ **
20
+ ** Returned (function value):
21
+ ** int status: 0 = OK
22
+ **
23
+ ** Note:
24
+ **
25
+ ** tcg1+tcg2 is Julian Date, apportioned in any convenient way
26
+ ** between the two arguments, for example where tcg1 is the Julian
27
+ ** Day Number and tcg22 is the fraction of a day. The returned
28
+ ** tt1,tt2 follow suit.
29
+ **
30
+ ** References:
31
+ **
32
+ ** McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),.
33
+ ** IERS Technical Note No. 32, BKG (2004)
34
+ **
35
+ ** IAU 2000 Resolution B1.9
36
+ **
37
+ ** This revision: 2010 May 14
38
+ **
39
+ ** Original version 2012-03-01
40
+ **
41
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
42
+ */
43
+ {
44
+
45
+ /* 1977 Jan 1 00:00:32.184 TT, as MJD */
46
+ static const double t77t = DJM77 + TTMTAI/DAYSEC;
47
+
48
+
49
+ /* Result, safeguarding precision. */
50
+ if ( tcg1 > tcg2 ) {
51
+ *tt1 = tcg1;
52
+ *tt2 = tcg2 - ( ( tcg1 - DJM0 ) + ( tcg2 - t77t ) ) * ELG;
53
+ } else {
54
+ *tt1 = tcg1 - ( ( tcg2 - DJM0 ) + ( tcg1 - t77t ) ) * ELG;
55
+ *tt2 = tcg2;
56
+ }
57
+
58
+ /* OK status. */
59
+ return 0;
60
+
61
+ /*----------------------------------------------------------------------
62
+ **
63
+ ** Celes is a wrapper of the SOFA Library for Ruby.
64
+ **
65
+ ** This file is redistributed and relicensed in accordance with
66
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
67
+ **
68
+ ** The original library is available from IAU Standards of
69
+ ** Fundamental Astronomy (http://www.iausofa.org/).
70
+ **
71
+ **
72
+ **
73
+ **
74
+ **
75
+ ** Copyright (C) 2013, Naoki Arita
76
+ ** All rights reserved.
77
+ **
78
+ ** Redistribution and use in source and binary forms, with or without
79
+ ** modification, are permitted provided that the following conditions
80
+ ** are met:
81
+ **
82
+ ** 1 Redistributions of source code must retain the above copyright
83
+ ** notice, this list of conditions and the following disclaimer.
84
+ **
85
+ ** 2 Redistributions in binary form must reproduce the above copyright
86
+ ** notice, this list of conditions and the following disclaimer in
87
+ ** the documentation and/or other materials provided with the
88
+ ** distribution.
89
+ **
90
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
91
+ ** the International Astronomical Union nor the names of its
92
+ ** contributors may be used to endorse or promote products derived
93
+ ** from this software without specific prior written permission.
94
+ **
95
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
96
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
97
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
98
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
99
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
100
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
101
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
102
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
103
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
104
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
105
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
106
+ ** POSSIBILITY OF SUCH DAMAGE.
107
+ **
108
+ **--------------------------------------------------------------------*/
109
+ }