ApJdataFrames Liu et al. 2010

Title: Discovery of a Highly Unequal-mass Binary T Dwarf with Keck Laser Guide Star Adaptive Optics: A Coevality Test of Substellar Theoretical Models and Effective Temperatures
Authors: Liu, Dupuy, Leggett

Data is from this paper:
http://iopscience.iop.org/article/10.1088/0004-637X/722/1/311/


In [1]:
import warnings
warnings.filterwarnings("ignore")

In [2]:
from astropy.io import ascii

In [3]:
import pandas as pd

Table 1 - Target Information for Ophiuchus Sources


In [4]:
#! mkdir ../data/Liu2010
#! curl http://iopscience.iop.org/0004-637X/722/1/311/suppdata/apj336343t6_ascii.txt > ../data/Liu2010/apj336343t6_ascii.txt

In [5]:
tbl6 = pd.read_csv("../data/Liu2010/apj336343t6_ascii.txt",
                   sep="\t", na_values=" ... ", skiprows=[0,1,2], skipfooter=1, usecols=range(9))
tbl6.head()


Out[5]:
Bandpass c_0 c_1 c_2 c_3 c_4 c_5 c_6 rms (mag)
0 J_MKO 18.904480 -8.053993 1.491738 -0.136761 0.006541 -0.000156 1.462266e-06 0.14
1 H_MKO 13.667090 -5.426683 1.027185 -0.096185 0.004734 -0.000117 1.148133e-06 0.07
2 K_MKO 5.795845 -1.471358 0.286801 -0.026475 0.001283 -0.000032 3.159780e-07 0.08

We want the J-Band coefficients, in order of highest order coefficient to lowest order coefficient.


In [6]:
J_s = tbl6.loc[0]

In [7]:
coeffs = J_s[["c_"+str(i) for i in range(6, -1, -1)]].values

In [8]:
func = np.poly1d(coeffs)

In [9]:
print(func)


           6             5            4          3         2
1.462e-06 x - 0.0001559 x + 0.006541 x - 0.1368 x + 1.492 x - 8.054 x + 18.9

The end