CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
Avatar for Testing 18.04.

Real-time collaboration for Jupyter Notebooks, Linux Terminals, LaTeX, VS Code, R IDE, and more,
all in one place. Commercial Alternative to JupyterHub.

| Download
Project: Testing 18.04
Views: 836
Image: ubuntu1804-dev
Tweet
Share
Share
Kernel: Python 3 (system-wide)

ISOCHRONES on CoCalc in Python 3

https://isochrones.readthedocs.io/en/latest/

Data globally stored in $ISOCHRONES

command line:

  1. unset DISPLAY

  2. import matplotlib; matplotlib.use('agg')

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

%matplotlib inline

import sys
sys.executable
'/usr/bin/python3'
sys.version
'3.6.9 (default, Jul 17 2020, 12:50:27) \n[GCC 8.4.0]'
import isochrones
isochrones
PyMultiNest not imported. MultiNest fits will not work.
<module 'isochrones' from '/usr/local/lib/python3.6/dist-packages/isochrones/__init__.py'>
isochrones.__version__
'2.1'
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()
Image in a Jupyter notebook
import astroquery
astroquery.__version__
'0.4'

dependencies section

import astroquery
astroquery.__version__
'0.4'
## uses too much memory?!

from isochrones.mist import MIST_Isochrone
mist = MIST_Isochrone()
mist.radius(1.0, 9.7, 0.0)
array(nan)
! ls -lh /ext/data/isochrones/mist/
total 1.2G drwxr-xr-x 1 salvus salvus 1.4K Oct 14 09:05 MIST_v1.2_vvcrit0.4_full_isos -rw-r--r-- 1 salvus salvus 661M Sep 29 11:40 MIST_v1.2_vvcrit0.4_full_isos.txz -rw-r--r-- 1 salvus salvus 18M Oct 28 10:48 dm_deep_v1.2_vvcrit0.4_full_isos.h5 -rw-r--r-- 1 salvus salvus 336M Oct 28 10:48 full_grid_v1.2_vvcrit0.4_full_isos.npz -rw-r--r-- 1 salvus salvus 177M Oct 8 18:02 mist_v1.2_vvcrit0.4_full_isos.h5 drwxr-xr-x 1 salvus salvus 1.3K Oct 13 23:25 tracks 
! ls -lh /ext/data/isochrones/BC/mist/
total 337M -rw-r--r-- 1 salvus salvus 19M Oct 8 18:04 SDSSugriz.h5 -rw-r--r-- 1 salvus salvus 1.8M Oct 8 18:04 SDSSugriz.txz -rw-r--r-- 1 salvus salvus 71M Oct 8 18:04 UBVRIplus.h5 -rw-r--r-- 1 salvus salvus 8.5M Oct 8 18:04 UBVRIplus.txz -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm025.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm025.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm050.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm050.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm075.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm075.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm100.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm100.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm125.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm125.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm150.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm150.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm175.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm175.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm200.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm200.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm225.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm225.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm250.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm250.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm275.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm275.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm300.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm300.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm350.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm350.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehm400.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehm400.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehp000.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehp000.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 16 2016 fehp025.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehp025.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 29 2016 fehp050.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehp050.UBVRIplus -rw-rw-r-- 1 salvus salvus 3.0M Apr 30 2016 fehp075.SDSSugriz -rw-rw-r-- 1 salvus salvus 11M Feb 5 2019 fehp075.UBVRIplus