# ----------------------------- MAIN ----------------------------- # ====================================================================== # Read the catalog # ====================================================================== # Read Gaia DR2 IERS quasar data gdr2 = read_dr2_iers() # Read ICRF1 catalog icrf1 = read_icrf1() # Read ICRF2 catalog icrf2 = read_icrf2() # Read ICRF3 S/X catalog icrf3sx = read_icrf3(wv="sx") # Read ICRF3 K catalog icrf3k = read_icrf3(wv="k") # Read ICRF3 S/X catalog icrf3xka = read_icrf3(wv="xka") # ====================================================================== # Calculate the position difference # ====================================================================== # ICRF1 - Gaia DR2 icrf1_gdr2 = radio_cat_diff_calc(icrf1, gdr2, "iers_name", ["icrf1", "gdr2"]) # ICRF2 - Gaia DR2 icrf2_gdr2 = radio_cat_diff_calc(icrf2, gdr2, "iers_name", ["icrf2", "gdr2"])
dra_ddec_cov / dra_err / ddec_err)
sep_table = Table([com_sou["iers_name"], sep[0], sep[-1]],
names=["iers_name", "ang_sep", "nor_sep"])
return sep_table
# ------------------------------ MAIN -- -----------------------------
# Read the catalogs.
# ICRF2 catalog
icrf2 = read_icrf2()
# ICRF3 S/X catalog
icrf3sx = read_icrf3(wv="sx")
# Gaia DR2 aux_iers catalog
gdr2 = read_dr2_iers()
# Calculate the positional differences.
# ICRF2 - Gaia DR2
icrf2_gdr2 = radio_cat_diff_calc(icrf2, gdr2, "iers_name")
# ICRF3 SX - Gaia DR2
icrf3sx_gdr2 = radio_cat_diff_calc(icrf3sx, gdr2, "iers_name")
# VSH parameters
vshpar_icrf2 = np.genfromtxt("../logs/icrf2_gaiadr2_vsh02.log",
usecols=(1, ),
skip_header=1)