dc = erfa.af2a(-12,31,10.3965)
##print('rc, dc', rc,dc)
#
reprd("ICRS, epoch J2000.0:", rc, dc )
# Proper motion: RA/Dec derivatives, epoch J2000.0.
pr = math.atan2(-354.45e-3 * erfa.DAS2R, math.cos(dc))
pd = 595.35e-3 * erfa.DAS2R
#print('pr, pd: ', pr, pd)
# Parallax (arcsec) and recession speed (km/s).
px = 164.99e-3
rv = 0.0
# ICRS to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rc, dc, pr, pd, px, rv, tt1, tt2)
#print('ri, di', ri, di)
#
reprd ( "catalog -> CIRS:", ri, di )
# 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 )
dc = erfa.af2a(np.array([[-12, 31, 10.3965]]))
##print('rc, dc', rc,dc)
#
reprd("ICRS, epoch J2000.0:", rc, dc)
# Proper motion: RA/Dec derivatives, epoch J2000.0.
pr = np.array([math.atan2(-354.45e-3 * erfa.DAS2R, math.cos(dc))])
pd = np.array([595.35e-3 * erfa.DAS2R])
# print('pr, pd: ', pr, pd)
# Parallax (arcsec) and recession speed (km/s).
px = np.array([164.99e-3])
rv = np.array([0.0])
# ICRS to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rc, dc, pr, pd, px, rv, tt1, tt2)
# print('ri, di', ri, di)
#
reprd("catalog -> CIRS:", ri, di)
# 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, 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)
dc = erfa.af2a(-12, 31, 10.3965)
##print('rc, dc', rc,dc)
#
reprd("ICRS, epoch J2000.0:", rc, dc)
# Proper motion: RA/Dec derivatives, epoch J2000.0.
pr = math.atan2(-354.45e-3 * erfa.DAS2R, math.cos(dc))
pd = 595.35e-3 * erfa.DAS2R
#print('pr, pd: ', pr, pd)
# Parallax (arcsec) and recession speed (km/s).
px = 164.99e-3
rv = 0.0
# ICRS to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rc, dc, pr, pd, px, rv, tt1, tt2)
#print('ri, di', ri, di)
#
reprd("catalog -> CIRS:", ri, di)
# 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)
