plot_sky scatter labels for single target in time

[HealthyPear]

When doing something like

from astroplan.plots import plot_sky
from astroplan import time_grid_from_range
import astropy.units as u

start_time = Time("2024-01-01 00:00:00")
time_grid = time_grid_from_range([start_time, start_time+30*u.day*12], time_resolution=30*u.day)

plot_sky(target_source, observer, time_grid, ax=None, style_kwargs=None, north_to_east_ccw=False, grid=True, az_label_offset=0*u.deg, warn_below_horizon=False, style_sheet=None)

the timing information is missing from the plot (in this case it would be a small number from 1 to 12 months, but in the case of an arbitrary amount of time it would be a finer array)

I did not find a simple way from astroplan API to do it, so I would propose to do it in the form of a (optional) colorbar which with a time scalar mappable maps time to color.

[HealthyPear]

Something like this for e.,g. the culmination of a source very month for a year

I can open a PR if you like it.

[HealthyPear]

in case of multiple targets pone can always use different markers