In [1]:

    

import pandas
import numpy as np
import isochrones.isochrones as iso
from isochrones import dartmouth


    

In [2]:

    

dar = dartmouth.Dartmouth_Isochrone()


    

In [3]:

    

ev = dar.evtrack(1.5,dage=0.02)
plot((ev['Teff']),ev['logL'],'.')


    

    
Out[3]:





[<matplotlib.lines.Line2D at 0x10e2a4e50>]






    

In [27]:

    

import isochrones.starmodel as sm
reload(iso)
reload(dartmouth)
dar = dartmouth.Dartmouth_Isochrone()
mod = sm.StarModel(dar,Teff=(5783,64),logg=(4.36,0.07),feh=(-0.06,0.1))#,
                    #H=(9.92,0.04),J=(10.23,0.04),K=(9.85,0.04))


    

In [28]:

    

mod.fit_mcmc()


    

In [29]:

    

mod.plot_samples('radius')
mod.plot_samples('mass')


    

In [ ]: