| line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
|
1
|
|
|
|
|
|
|
=head1 NAME |
|
2
|
|
|
|
|
|
|
|
|
3
|
|
|
|
|
|
|
Time::TCB - Barycentric Coordinate Time |
|
4
|
|
|
|
|
|
|
|
|
5
|
|
|
|
|
|
|
=head1 SYNOPSIS |
|
6
|
|
|
|
|
|
|
|
|
7
|
|
|
|
|
|
|
use Time::TCB qw(tcb_instant_to_mjd tcb_mjd_to_instant); |
|
8
|
|
|
|
|
|
|
|
|
9
|
|
|
|
|
|
|
$mjd = tcb_instant_to_mjd($instant); |
|
10
|
|
|
|
|
|
|
$instant = tcb_mjd_to_instant($mjd); |
|
11
|
|
|
|
|
|
|
|
|
12
|
|
|
|
|
|
|
=head1 DESCRIPTION |
|
13
|
|
|
|
|
|
|
|
|
14
|
|
|
|
|
|
|
Barycentric Coordinate Time (TCB) is a coordinate time scale representing |
|
15
|
|
|
|
|
|
|
time in the Sol system. Specifically, it is the proper time experienced |
|
16
|
|
|
|
|
|
|
by a distant clock comoving with the barycentre of the Sol system. |
|
17
|
|
|
|
|
|
|
|
|
18
|
|
|
|
|
|
|
This module represents instants on the TCB time scale as a scalar number |
|
19
|
|
|
|
|
|
|
of SI seconds since an epoch. This is an appropriate form for all manner |
|
20
|
|
|
|
|
|
|
of calculations. TCB is defined with a well-known point at TAI instant |
|
21
|
|
|
|
|
|
|
1977-01-01T00:00:00.0 at the Terran geocentre. This point is assigned the |
|
22
|
|
|
|
|
|
|
scalar value -460_080_000, putting the epoch at approximately the date |
|
23
|
|
|
|
|
|
|
at which the resolution defining TCB was adopted by the International |
|
24
|
|
|
|
|
|
|
Astronomical Union. This epoch is deliberately very different from |
|
25
|
|
|
|
|
|
|
those used for Geocentric Coordinate Time (TCG) in L and for |
|
26
|
|
|
|
|
|
|
Terrestrial Time (TT) in L, to avoid confusion between them. |
|
27
|
|
|
|
|
|
|
|
|
28
|
|
|
|
|
|
|
There is also a conventional way to represent TCB instants using day-based |
|
29
|
|
|
|
|
|
|
notations associated with planetary rotation `time' scales. The `day' |
|
30
|
|
|
|
|
|
|
of TCB is a nominal period of exactly 86400 SI seconds, which is slightly |
|
31
|
|
|
|
|
|
|
shorter than an actual Terran day. The well-known point at TAI instant |
|
32
|
|
|
|
|
|
|
1977-01-01T00:00:00.0 is assigned the label 1977-01-01T00:00:32.184 |
|
33
|
|
|
|
|
|
|
(MJD 43144.0003725). Because TCB is not connected to Terran rotation, |
|
34
|
|
|
|
|
|
|
and so has no inherent concept of a day, it is somewhat misleading to |
|
35
|
|
|
|
|
|
|
use such day-based notations. Conversion between this notation and |
|
36
|
|
|
|
|
|
|
the linear count of seconds is supported by this module. The day-based |
|
37
|
|
|
|
|
|
|
notations for TT, TCG, and TCB instants yield very similar values for |
|
38
|
|
|
|
|
|
|
corresponding instants, so care must be taken to avoid confusion. |
|
39
|
|
|
|
|
|
|
|
|
40
|
|
|
|
|
|
|
=cut |
|
41
|
|
|
|
|
|
|
|
|
42
|
|
|
|
|
|
|
package Time::TCB; |
|
43
|
|
|
|
|
|
|
|
|
44
|
1
|
|
|
1
|
|
22756
|
{ use 5.006; } |
|
|
1
|
|
|
|
|
3
|
|
|
|
1
|
|
|
|
|
45
|
|
|
45
|
1
|
|
|
1
|
|
5
|
use warnings; |
|
|
1
|
|
|
|
|
2
|
|
|
|
1
|
|
|
|
|
28
|
|
|
46
|
1
|
|
|
1
|
|
11
|
use strict; |
|
|
1
|
|
|
|
|
11
|
|
|
|
1
|
|
|
|
|
34
|
|
|
47
|
|
|
|
|
|
|
|
|
48
|
1
|
|
|
1
|
|
1179
|
use Math::BigRat 0.13; |
|
|
1
|
|
|
|
|
88445
|
|
|
|
1
|
|
|
|
|
6
|
|
|
49
|
|
|
|
|
|
|
|
|
50
|
|
|
|
|
|
|
our $VERSION = "0.002"; |
|
51
|
|
|
|
|
|
|
|
|
52
|
1
|
|
|
1
|
|
2244
|
use parent "Exporter"; |
|
|
1
|
|
|
|
|
287
|
|
|
|
1
|
|
|
|
|
7
|
|
|
53
|
|
|
|
|
|
|
our @EXPORT_OK = qw(tcb_instant_to_mjd tcb_mjd_to_instant); |
|
54
|
|
|
|
|
|
|
|
|
55
|
|
|
|
|
|
|
=head1 FUNCTIONS |
|
56
|
|
|
|
|
|
|
|
|
57
|
|
|
|
|
|
|
=over |
|
58
|
|
|
|
|
|
|
|
|
59
|
|
|
|
|
|
|
=item tcb_instant_to_mjd(INSTANT) |
|
60
|
|
|
|
|
|
|
|
|
61
|
|
|
|
|
|
|
Converts from a count of seconds to a Modified Julian Date in the manner |
|
62
|
|
|
|
|
|
|
conventional for TCB. The MJD can be further converted to other forms of |
|
63
|
|
|
|
|
|
|
day-based date using other modules. The input must be a C |
|
64
|
|
|
|
|
|
|
object, and the result is the same type. |
|
65
|
|
|
|
|
|
|
|
|
66
|
|
|
|
|
|
|
=cut |
|
67
|
|
|
|
|
|
|
|
|
68
|
1
|
|
|
1
|
|
81
|
use constant TCB_EPOCH_MJD => Math::BigRat->new("48469.0003725"); |
|
|
1
|
|
|
|
|
2
|
|
|
|
1
|
|
|
|
|
7
|
|
|
69
|
|
|
|
|
|
|
|
|
70
|
|
|
|
|
|
|
sub tcb_instant_to_mjd($) { |
|
71
|
4
|
|
|
4
|
1
|
5117
|
my($instant) = @_; |
|
72
|
4
|
|
|
|
|
15
|
return TCB_EPOCH_MJD + ($instant / 86400); |
|
73
|
|
|
|
|
|
|
} |
|
74
|
|
|
|
|
|
|
|
|
75
|
|
|
|
|
|
|
=item tcb_mjd_to_instant(MJD) |
|
76
|
|
|
|
|
|
|
|
|
77
|
|
|
|
|
|
|
Converts from a Modified Julian Date, interpreted in the manner |
|
78
|
|
|
|
|
|
|
conventional for TCB, to a count of seconds. The input must be a |
|
79
|
|
|
|
|
|
|
C object, and the result is the same type. |
|
80
|
|
|
|
|
|
|
|
|
81
|
|
|
|
|
|
|
=cut |
|
82
|
|
|
|
|
|
|
|
|
83
|
|
|
|
|
|
|
sub tcb_mjd_to_instant($) { |
|
84
|
4
|
|
|
4
|
1
|
4409
|
my($mjd) = @_; |
|
85
|
4
|
|
|
|
|
13
|
return ($mjd - TCB_EPOCH_MJD) * 86400; |
|
86
|
|
|
|
|
|
|
} |
|
87
|
|
|
|
|
|
|
|
|
88
|
|
|
|
|
|
|
=back |
|
89
|
|
|
|
|
|
|
|
|
90
|
|
|
|
|
|
|
=head1 SEE ALSO |
|
91
|
|
|
|
|
|
|
|
|
92
|
|
|
|
|
|
|
L, |
|
93
|
|
|
|
|
|
|
L, |
|
94
|
|
|
|
|
|
|
L |
|
95
|
|
|
|
|
|
|
|
|
96
|
|
|
|
|
|
|
=head1 AUTHOR |
|
97
|
|
|
|
|
|
|
|
|
98
|
|
|
|
|
|
|
Andrew Main (Zefram) |
|
99
|
|
|
|
|
|
|
|
|
100
|
|
|
|
|
|
|
=head1 COPYRIGHT |
|
101
|
|
|
|
|
|
|
|
|
102
|
|
|
|
|
|
|
Copyright (C) 2006, 2010, 2012 Andrew Main (Zefram) |
|
103
|
|
|
|
|
|
|
|
|
104
|
|
|
|
|
|
|
=head1 LICENSE |
|
105
|
|
|
|
|
|
|
|
|
106
|
|
|
|
|
|
|
This module is free software; you can redistribute it and/or modify it |
|
107
|
|
|
|
|
|
|
under the same terms as Perl itself. |
|
108
|
|
|
|
|
|
|
|
|
109
|
|
|
|
|
|
|
=cut |
|
110
|
|
|
|
|
|
|
|
|
111
|
|
|
|
|
|
|
1; |