import cosmoplotian
from matplotlib import pyplot as plt
ax = plt.axes(projection='astro hours mollweide')
ax.grid()