-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathinputs.py
42 lines (36 loc) · 1.75 KB
/
inputs.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
"""
Some inputs pointing to data locations and other variables that may change.
"""
import numpy as np
# current parameters for the model and their order
labels = ['$P$ (days)', '$t_{tran}$ (days)', '$e\cos\omega$', '$e\sin\omega$',
'$b$', '$M_{2,init}$', '$M_2$', '$M_1$', '[Fe/H]', 'Age (Gyr)',
'Distance (pc)', '$\sigma_{sys}$', 'h (pc)', '$A_\lambda$ scale']
# BJD - timeoffset for all light curves
timeoffset = 55000.
# location of the light curve fits files
keplerdata = './lightcurve/'
# file with a list of times to ignore
baddata = './ignorelist.txt'
# directory containing the PARSEC isochrones used
isodir = './PARSECv1.1/'
# what's in the isochrone and what its column index is
inds = {'feh': 0, 'age': 1, 'M': 2, 'Mact': 3, 'lum': 4, 'teff': 5, 'logg': 6,
'mbol': 7, 'Kp': 8, 'g': 9, 'r': 10, 'i': 11, 'z': 12, 'D51': 13,
'J': 14, 'H': 15, 'Ks': 16, 'int_IMF1': 17, '3.6': 18, '4.5': 19,
'5.8': 20, '8.0': 21, '24': 22, '70': 23, '160': 24, 'W1': 25,
'W2': 26, 'W3': 27, 'W4': 28, 'int_IMF2': 29, 'Rad': 30}
# magnitude names, values, errors, and extinction of the system in other
# filters that we want to simulate in the isochrones
magname = ['g', 'r', 'i', 'z', 'J', 'H', 'Ks', 'W1', 'W2']
magobs = np.array([15.245, 14.617, 14.435, 14.330, 13.447, 13.039, 12.977,
12.959, 13.020])
magerr = np.array([0.003, 0.004, 0.003, 0.004, 0.024, 0.022, 0.032, 0.026,
0.029])
# extinction measure in each filter
maglam = np.array([0.359, 0.248, 0.184, 0.137, 0.078, 0.050, 0.034, 0.021,
0.016])
# white dwarf cooling models
wdfiles = './wdmodels/Table_Mass*'
# what's in the model and what its index is
wdinds = {'g': 13, 'r': 14, 'i': 15, 'teff': 0, 'logg': 1, 'age': 26}