rca, dca, eo = erfa.atic13 ( ri, di, tt1, tt2)
#
reprd ( "CIRS -> astrometric:", rca, dca )
#ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13 ( rca, dca, 0.0, 0.0, 0.0, 0.0, tt1, tt2)
reprd ( "astrometric -> CIRS:", ri, di );
# Apparent place.
ra = erfa.anp ( ri - eo )
da = di
reprd ( "geocentric apparent:", ra, da )
# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(ri, di, utc1, utc2, dut1,
elong, phi, hm, xp, yp,
0.0, 0.0, 0.0, 0.0)
reprd( "CIRS -> topocentric:", rot, dot )
# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1,
elong, phi, hm, xp, yp,
phpa, tc, rh, wl)
reprd( "CIRS -> observed:", rob, dob )
# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13 ( rc, dc, pr, pd, px, rv,
utc1, utc2, dut1,
elong, phi, hm, xp, yp,
phpa, tc, rh, wl)
reprd ( "ICRS -> observed:", rob, dob )
rca, dca, eo = erfa.atic13(ri, di, tt1, tt2)
#
reprd("CIRS -> astrometric:", rca, dca)
# ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rca, dca, np.array([0.0]), np.array([0.0]), np.array([0.0]), np.array([0.0]), tt1, tt2)
reprd("astrometric -> CIRS:", ri, di)
# Apparent place.
ra = erfa.anp(ri - eo)
da = di
reprd("geocentric apparent:", ra, da)
# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(
ri, di, utc1, utc2, dut1, elong, phi, hm, xp, yp, np.array([0.0]), np.array([0.0]), np.array([0.0]), np.array([0.0])
)
reprd("CIRS -> topocentric:", rot, dot)
# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm, xp, yp, phpa, tc, rh, wl)
reprd("CIRS -> observed:", rob, dob)
# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13(
rc, dc, pr, pd, px, rv, utc1, utc2, dut1, elong, phi, hm, xp, yp, phpa, tc, rh, wl
)
reprd("ICRS -> observed:", rob, dob)
# ICRS to CIRS using some user-supplied parameters.
# SOFA heliocentric Earth ephemeris.
# CIRS to ICRS (astrometric).
rca, dca, eo = erfa.atic13(ri, di, tt1, tt2)
#
reprd("CIRS -> astrometric:", rca, dca)
#ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rca, dca, 0.0, 0.0, 0.0, 0.0, tt1, tt2)
reprd("astrometric -> CIRS:", ri, di)
# Apparent place.
ra = erfa.anp(ri - eo)
da = di
reprd("geocentric apparent:", ra, da)
# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm,
xp, yp, 0.0, 0.0, 0.0, 0.0)
reprd("CIRS -> topocentric:", rot, dot)
# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm,
xp, yp, phpa, tc, rh, wl)
reprd("CIRS -> observed:", rob, dob)
# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13(rc, dc, pr, pd, px, rv, utc1, utc2,
dut1, elong, phi, hm, xp, yp, phpa,
tc, rh, wl)
reprd("ICRS -> observed:", rob, dob)
# ICRS to CIRS using some user-supplied parameters.
# SOFA heliocentric Earth ephemeris.
