ISOCHRONES on CoCalc in Python 3

import numpy as np
import matplotlib.pyplot as plt
import sep, astroalign
import ccdproc

%matplotlib inline

import sys
sys.executable

sys.version

import isochrones
isochrones

isochrones.__version__

from isochrones import get_ichrone
iso = get_ichrone('mist', bands=['Bessell_B', 'Bessell_V'])

ages = np.linspace(1.4e9, 2e9, 5)

model = iso.isochrone(9.0)

model

model_bvs = []
model_vs = []
labels = []
for age in ages:
    model = iso.isochrone(np.log10(age))
    model_b = model.Bessell_B_mag
    model_v = model.Bessell_V_mag
    model_bv = model_b - model_v
    model_vs.append(model_v)
    model_bvs.append(model_bv)
    labels.append('{:.2f} Gyr'.format(age/1e9))

fig, axis = plt.subplots(figsize=(8,6))
axis.set_xlabel('B-V')
axis.set_ylabel('V')
for label, v, bv in zip(labels, model_vs, model_bvs):
    axis.plot(bv, v, label=label)
axis.legend()
axis.set_ylim(-2,4)
axis.invert_yaxis()

plt.show()

import astroquery
astroquery.__version__

import astroquery
astroquery.__version__

## uses too much memory?!

from isochrones.mist import MIST_Isochrone
mist = MIST_Isochrone()
mist.radius(1.0, 9.7, 0.0)

! ls -lh /ext/data/isochrones/mist/

! ls -lh /ext/data/isochrones/BC/mist/