In [7]:

    

import numpy as np
import matplotlib.pyplot as plt
from astropy.io import ascii


    

In [5]:

    

masses = [0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 1.000, 1.100, 1.200, 1.300, 1.400]
Ts = [3634., 3802., 3955., 4078., 4192., 4290., 4377., 4456., 4529., 4596., 4658.]
ll = [-0.32, -0.18, -0.07, 0.04, 0.13, 0.21, 0.28, 0.35, 0.41, 0.47, 0.52]


    

In [6]:

    

fig, ax = plt.subplots(nrows=2, sharex=True)
ax[0].plot(masses, Ts)
ax[0].plot(masses, Ts, "o")
ax[1].plot(masses, ll)
ax[1].plot(masses, ll, "o")
fig.savefig("linearOrLog.png")


    

In [8]:

    

data = ascii.read("/home/ian/Grad/Research/Stars/ScottiePippen/data/Dartmouth/PMS/fehp00afep0/m125fehp00afep0.jc2mass", names=["age", "LTeff", "logg", "LL", "U", "B", "V", "R", "I", "J", "H", "Ks"])


    

In [12]:

    

age = 1e-6 * data["age"] # [Myr]

ind = (age < 100)

age = age[ind]
ll = data["LL"][ind] # [log(L_sun)]
tt = 10**data["LTeff"][ind] #[K]


    

In [13]:

    

fig, ax = plt.subplots(nrows=2, sharex=True)
ax[0].plot(age, tt)
ax[0].plot(age, tt, "o")
ax[1].plot(age, ll)
ax[1].plot(age, ll, "o")
fig.savefig("linearOrLog_age.png")


    

In [ ]: