Better UT1

JuliaAstro · Aug 15, 2018 · ef49c84 · ef49c84
1 parent 07aff95
commit ef49c84
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/src/offsets.jl b/src/offsets.jl
@@ -14,8 +14,8 @@ tai_offset(::GeocentricCoordinateTime, ep) = tai_offset(TT, ep) + LG_RATE * get(
 tai_offset(::BarycentricCoordinateTime, ep) = tai_offset(TDB, ep) + LB_RATE * get(ep - EPOCH_77)
 
 function tai_offset(::UniversalTime, ep)
-    # FIXME: This is not correct
-    tai_offset(UTC, ep) + getΔUT1(julian(ep))
+    jd = get(days(ep - UTCEpoch(AstroDates.JULIAN_EPOCH, AstroDates.H12)))
+    tai_offset(UTC, ep) + getΔUT1(jd)
 end
 
 function tai_offset(::CoordinatedUniversalTime, ep)

diff --git a/test/epochs.jl b/test/epochs.jl
@@ -24,7 +24,8 @@
         @test utc.epoch == 587513008
         @test utc.offset == 0.5512474363782758
         @test ut1.epoch == 587513008
-        @test_broken ut1.offset == 0.47200339587608653
+        # EarthOrientation.jl uses a different interpolation scheme than Orekit
+        @test ut1.offset ≈ 0.48504859616502927 atol=1e-6
         @test tdb.epoch == 587512939
         @test tdb.offset ≈ 0.3682890196414874 atol=1e-14
         @test tcb.epoch == 587512919