Astroplan in CoCalc

import astroplan
astroplan.__version__

from astropy.coordinates import SkyCoord
from astroplan import FixedTarget

altair = FixedTarget.from_name('Altair')
from astroplan import Observer
observer = Observer.at_site('subaru')

from astroplan import Observer, FixedTarget
from astropy.time import Time
time_range = Time(["2015-08-01 06:00", "2015-08-01 12:00"])

%%writefile targets.txt
# name ra_degrees dec_degrees
Polaris 37.95456067 89.26410897
Vega 279.234734787 38.783688956
Albireo 292.68033548 27.959680072
Algol 47.042218553 40.955646675
Rigel 78.634467067 -8.201638365
Regulus 152.092962438 11.967208776

# Read in the table of targets
from astropy.table import Table
target_table = Table.read('targets.txt', format='ascii')

# Create astroplan.FixedTarget objects for each one in the table
from astropy.coordinates import SkyCoord
import astropy.units as u
targets = [FixedTarget(coord=SkyCoord(ra=ra*u.deg, dec=dec*u.deg), name=name)
           for name, ra, dec in target_table]
targets

from astroplan import (AltitudeConstraint, AirmassConstraint,
                       AtNightConstraint)
constraints = [AltitudeConstraint(10*u.deg, 80*u.deg),
               AirmassConstraint(5), AtNightConstraint.twilight_civil()]

from astroplan import is_observable, is_always_observable, months_observable
# Are targets *ever* observable in the time range?
ever_observable = is_observable(constraints, subaru, targets, time_range=time_range)

# Are targets *always* observable in the time range?
always_observable = is_always_observable(constraints, subaru, targets, time_range=time_range)

# During what months are the targets ever observable?
best_months = months_observable(constraints, subaru, targets)

from astroplan import download_IERS_A
download_IERS_A()

from astropy.table import Table
import numpy as np
observability_table = Table()
observability_table['targets'] = [target.name for target in targets]
observability_table['ever_observable'] = ever_observable
observability_table['always_observable'] = always_observable
print(observability_table)

from astroplan import observability_table
table = observability_table(constraints, subaru, targets, time_range=time_range)
print(table)

import matplotlib.pyplot as plt
from matplotlib import cm
from astroplan.plots import plot_sky

from astroplan import time_grid_from_range
time_range = Time(["2015-08-01 06:00", "2015-08-01 12:00"])
time_grid = time_grid_from_range(time_range)

plt.figure(figsize=(10,10))
cmap = cm.Set1             # Cycle through this colormap

for i, target in enumerate(targets):
    ax = plot_sky(target, subaru, time_grid,
                  style_kwargs=dict(color=cmap(float(i)/len(targets)),
                                    label=target.name))

legend = ax.legend(loc='lower center')
legend.get_frame().set_facecolor('w')
plt.show()