from ligo.skymap.tool.ligo_skymap_plot_coherence import main
from astropy.utils.data import download_file
filename = download_file('https://gracedb.ligo.org/api/superevents/S230707ai/files/bayestar.multiorder.fits,2', cache=True)
main([filename])