'font.size': 22,
'font.family': 'txfonts',
'text.latex.unicode': True,
}
plt.rcParams.update(params)
sizel = 17
#####################################################################
obj = 'Pasiphae'
a = ['OPD', 'OHP', 'ESO']
b = [874, 511, 809]
c = read2('{}.dat'.format(obj),
coord_col=[0, 1, 2, 3, 4, 5],
time_col=[8],
time_fmt='jd')
coord = SkyCoord(c['coords'][0], c['coords'][1], unit=(u.hourangle, u.deg))
time = Time(c['times'], format='jd', scale='utc')
code = np.loadtxt('{}.dat'.format(obj), usecols=[12], ndmin=1, unpack=True)
for i in [0, 1, 2]:
d = read2('{}-jpl.eph'.format(b[i]),
coord_col=[4, 5, 6, 7, 8, 9],
skiprows=3)
coordjpl = SkyCoord(d['coords'][0],
d['coords'][1],
unit=(u.hourangle, u.deg))
anom = np.loadtxt('{}-jpl.eph'.format(b[i]),
skiprows=3,
usecols=[37],
'text.latex.unicode': True,
}
plt.rcParams.update(params)
sizel = 17
#####################################################################
a = 'Nereida_eme.eph'
b = 'Nereida_jpl.eph'
m = 'Triton_zhang.eph'
c = read2(a,
time_col=[0, 1, 2, 3, 4, 5],
time_fmt='iso',
coord_col=[6, 7],
skiprows=8)
#c = read2(a, time_col=[2], time_fmt='jd', coord_col=[4,5,6,7,8,9], skiprows=3)
d = read2(b,
time_col=[2],
time_fmt='jd',
coord_col=[4, 5, 6, 7, 8, 9],
skiprows=3)
n = read2(m,
time_col=[0, 1, 2, 3, 4, 5],
time_fmt='iso',
coord_col=[6, 7],
skiprows=8)
#e = read2('Nereid.dat', time_col=[8], time_fmt='jd', coord_col=[0,1,2,3,4,5])
params = {'text.usetex' : True,
'font.size' : 22,
'font.family' : 'txfonts',
'text.latex.unicode': True,
}
plt.rcParams.update(params)
sizel = 17
#####################################################################
obj = 'Elara'
a = ['OPD', 'OHP', 'ESO']
b = [874, 511, 809]
c = read2('{}.dat'.format(obj), coord_col=[0, 1, 2, 3, 4, 5], time_col=[8], time_fmt='jd')
coord = SkyCoord(c['coords'][0], c['coords'][1], unit=(u.hourangle, u.deg))
time = Time(c['times'], format='jd', scale='utc')
code = np.loadtxt('{}.dat'.format(obj), usecols=[12], ndmin=1, unpack=True)
for i in [0, 1, 2]:
d = read2('{}-jpl.eph'.format(b[i]), coord_col=[4, 5, 6, 7, 8, 9], skiprows=3)
coordjpl = SkyCoord(d['coords'][0], d['coords'][1], unit=(u.hourangle, u.deg))
anom = np.loadtxt('{}-jpl.eph'.format(b[i]), skiprows=3, usecols=[37], ndmin=1, unpack=True)
e = read2('{}-ste.eph'.format(b[i]), coord_col=[4, 5, 6, 7, 8, 9], skiprows=3)
coordlau = SkyCoord(e['coords'][0], e['coords'][1], unit=(u.hourangle, u.deg))
f = read2('{}-eme.eph'.format(b[i]), coord_col=[1, 2], skiprows=8)
coordeme = SkyCoord(f['coords'][0], f['coords'][1], unit=(u.hourangle, u.deg))
difjpldec = coord[code == b[i]].dec - coordjpl.dec
difjplra = (coord[code == b[i]].ra - coordjpl.ra)*np.cos(coordjpl.dec)
f = open('{}_jpl.txt'.format(obj), 'a')
'font.size' : 22,
'font.family' : 'txfonts',
'text.latex.unicode': True,
}
plt.rcParams.update(params)
sizel = 17
#####################################################################
a = 'Phoebe_eme.eph'
b = 'Phoebe_JPL.eph'
m = 'Phoebe_Josselin.eph'
c = read2(a, time_col=[0,1,2,3,4,5], time_fmt='iso', coord_col=[6,7], skiprows=8)
#c = read2(a, time_col=[2], time_fmt='jd', coord_col=[4,5,6,7,8,9], skiprows=3)
d = read2(b, time_col=[2], time_fmt='jd', coord_col=[4,5,6,7,8,9], skiprows=3)
n = read2(m, time_col=[2], time_fmt='jd', coord_col=[4,5,6,7,8,9], skiprows=3)
#e = read2('Nereid.dat', time_col=[8], time_fmt='jd', coord_col=[0,1,2,3,4,5])
#f = read2('Nereida_LNA.eph', time_col=[0], time_fmt='jd', coord_col=[1,2,3,4,5,6], skiprows=8)
eme_time = Time(c['times'], format='iso', scale='utc')
eme_coord = SkyCoord(c['coords'].transpose(), unit=(u.hourangle, u.deg))
jpl_time = Time(d['times'], format='jd', scale='utc')
jpl_coord = SkyCoord(d['coords'].transpose(), unit=(u.hourangle, u.deg))
ph15_time = Time(n['times'], format='jd', scale='utc')
ph15_coord = SkyCoord(n['coords'].transpose(), unit=(u.hourangle, u.deg))