2 BEGIN { plan tests => 12 }
3 # Test the calculation of (modified) Julian date
6 # First a lookup table of epoch and MJD
7 # Use 3 sig fig in MJD (hence the use of strings)
8 # This will not work on systems that use a different reference
9 # epoch to unix time. To be more general we should use strptime
10 # to parse the reference date.
12 951827696 => '51603.524', # 2000-02-29T12:34:56UT
13 1000011 => '40598.574', # 1970-01-12T13:46:51UT
14 1021605703 => '52411.140', # 2002-05-17T03:21:43UT
15 1121605703 => '53568.547', # 2005-07-17T13:08:23UT
16 1011590000 => '52295.218', # 2002-01-21T05:13:20UT
17 1011605703 => '52295.399', # 2002-01-21T09:35:03
20 # Now loop over each MJD
21 for my $time (keys %mjd) {
23 # First check using GMT
24 my $tp = gmtime( $time );
25 ok(sprintf("%.3f",$tp->mjd),$mjd{$time});
27 # Now localtime should give the same answer for MJD
28 # since MJD is always referred to as UT
29 $tp = localtime( $time );
30 ok(sprintf("%.3f",$tp->mjd),$mjd{$time});