# Apply IERS corrections.
x += dx
y += dy
# SOFA CIO locator. */
s = erfa.s06(tt1, tt2, x, y)
# Populate the context.
astrom = erfa.apci(tt1, tt2, pvb, pvh[0], x, y, s)
# Carry out the transformation and report the results.
ri, di = erfa.atciq(rc, dc, pr, pd, px, rv, *astrom)
reprd ( "ICRS -> CIRS (JPL, IERS):", ri, di )
# The same but with Saturn then Jupiter then Sun light deflection.
b0 = erfa.LDBODY((0.00028574, 3e-10,
((-7.8101442680818964, -5.6095668114887358, -1.9807981923749924),
(0.0030723248971152, -0.0040699547707598, -0.0018133584165345))))
b1 = erfa.LDBODY((0.00095435, 3e-9,
((0.7380987962351833, .6365869247538951, 1.9693136030111202),
(-0.0075581692172088, 0.0012691372216750, 0.0007279990012801))))
b2 = erfa.LDBODY((1.0, 6e-6,
((-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
(0.0000062923521264, -0.0000003308883872, -0.0000002964866231))))
b= [b0, b1, b2]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, b, *astrom)
reprd ( "ICRS -> CIRS (+ planets):", ri, di )
# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, b, *astrom)
reprd ( "CIRS -> astrometric:", rca, dca )
(
0.00095435,
3e-9,
np.array(
(
(0.7380987962351833, 0.6365869247538951, 1.9693136030111202),
(-0.0075581692172088, 0.0012691372216750, 0.0007279990012801),
)
),
)
)
b2 = erfa.LDBODY(
(
1.0,
6e-6,
np.array(
(
(-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
(0.0000062923521264, -0.0000003308883872, -0.0000002964866231),
)
),
)
)
b = [[b0, b1, b2]]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, astrom, b)
reprd("ICRS -> CIRS (+ planets):", ri, di)
# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, astrom, b)
reprd("CIRS -> astrometric:", rca, dca)
# SOFA CIO locator. */
s = erfa.s06(tt1, tt2, x, y)
# Populate the context.
astrom = erfa.apci(tt1, tt2, pvb, pvh[0], x, y, s)
# Carry out the transformation and report the results.
ri, di = erfa.atciq(rc, dc, pr, pd, px, rv, *astrom)
reprd("ICRS -> CIRS (JPL, IERS):", ri, di)
# The same but with Saturn then Jupiter then Sun light deflection.
b0 = erfa.LDBODY(
(0.00028574, 3e-10,
((-7.8101442680818964, -5.6095668114887358, -1.9807981923749924),
(0.0030723248971152, -0.0040699547707598, -0.0018133584165345))))
b1 = erfa.LDBODY(
(0.00095435, 3e-9,
((0.7380987962351833, .6365869247538951, 1.9693136030111202),
(-0.0075581692172088, 0.0012691372216750, 0.0007279990012801))))
b2 = erfa.LDBODY(
(1.0, 6e-6,
((-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
(0.0000062923521264, -0.0000003308883872, -0.0000002964866231))))
b = [b0, b1, b2]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, b, *astrom)
reprd("ICRS -> CIRS (+ planets):", ri, di)
# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, b, *astrom)
reprd("CIRS -> astrometric:", rca, dca)