# if True:
# # try SOFA routines (but without UTC -> UT1)
# import sidereal
# # SHOULD TRANSFER TO UT1!!
# gmst = sidereal.gmst82(mjd.utc.jd1, mjd,utc.jd2)
if False:
# check with Fisher's ephemeris
import rf_ephem
rf_ephem.set_ephemeris_dir('/data/mhvk/packages/jplephem', 'DEc421')
rf_ephem.set_observer_coordinates(*xyz_gmrt)
rf_delay = rf_ephem.pulse_delay(
eph1957.evaluate('RAJ',mjd.tdb.mjd[0])/15.,
eph1957.evaluate('DECJ',mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0]-int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1]-mjd.utc.mjd[0])*24.*3600.)['delay']
rf_rv = rf_ephem.doppler_fraction(
eph1957.evaluate('RAJ',mjd.tdb.mjd[0])/15.,
eph1957.evaluate('DECJ',mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0]-int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1]-mjd.utc.mjd[0])*24.*3600.)['frac']
import matplotlib.pylab as plt
plt.ion()
plt.plot(mjd.utc.mjd, delay-rf_delay-d_orb)
plt.plot(mjd.utc.mjd, (rv-rf_rv-v_orb)*c.to(u.km/u.s).value)
# if True:
# # try SOFA routines (but without UTC -> UT1)
# import sidereal
# # SHOULD TRANSFER TO UT1!!
# gmst = sidereal.gmst82(mjd.utc.jd1, mjd,utc.jd2)
if False:
# check with Fisher's ephemeris
import rf_ephem
rf_ephem.set_ephemeris_dir('/data/mhvk/packages/jplephem', 'DEc421')
rf_ephem.set_observer_coordinates(*xyz_gmrt)
rf_delay = rf_ephem.pulse_delay(
eph1957.evaluate('RAJ',mjd.tdb.mjd[0])/15.,
eph1957.evaluate('DECJ',mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0]-int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1]-mjd.utc.mjd[0])*24.*3600.)['delay']
rf_rv = rf_ephem.doppler_fraction(
eph1957.evaluate('RAJ',mjd.tdb.mjd[0])/15.,
eph1957.evaluate('DECJ',mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0]-int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1]-mjd.utc.mjd[0])*24.*3600.)['frac']
import matplotlib.pylab as plt
plt.ion()
plt.plot(mjd.utc.mjd, delay-rf_delay-d_orb)
plt.plot(mjd.utc.mjd, (rv-rf_rv-v_orb)*c.to(u.km/u.s).value)
# if True:
# # try SOFA routines (but without UTC -> UT1)
# import sidereal
# # SHOULD TRANSFER TO UT1!!
# gmst = sidereal.gmst82(mjd.utc.jd1, mjd,utc.jd2)
if False:
# check with Fisher's ephemeris
import rf_ephem
rf_ephem.set_ephemeris_dir("/data/mhvk/packages/jplephem", "DEc421")
rf_ephem.set_observer_coordinates(*xyz_gmrt)
rf_delay = rf_ephem.pulse_delay(
eph1957.evaluate("RAJ", mjd.tdb.mjd[0]) / 15.0,
eph1957.evaluate("DECJ", mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0] - int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1] - mjd.utc.mjd[0]) * 24.0 * 3600.0,
)["delay"]
rf_rv = rf_ephem.doppler_fraction(
eph1957.evaluate("RAJ", mjd.tdb.mjd[0]) / 15.0,
eph1957.evaluate("DECJ", mjd.tdb.mjd[0]),
int(mjd.utc.mjd[0]),
mjd.utc.mjd[0] - int(mjd.utc.mjd[0]),
len(mjd),
(mjd.utc.mjd[1] - mjd.utc.mjd[0]) * 24.0 * 3600.0,
)["frac"]
import matplotlib.pylab as plt
plt.ion()
# if True:
# # try SOFA routines (but without UTC -> UT1)
# import sidereal
# # SHOULD TRANSFER TO UT1!!
# gmst = sidereal.gmst82(mjd.utc.jd1, mjd,utc.jd2)
if False:
# check with Fisher's ephemeris
import rf_ephem
rf_ephem.set_ephemeris_dir('/data/mhvk/packages/jplephem', 'DEc421')
rf_ephem.set_observer_coordinates(*xyz_gmrt)
rf_delay = []
rf_rv = []
for m in t.utc.mjd:
rf_delay += rf_ephem.pulse_delay(
eph1957.evaluate('RAJ',m)/15., eph1957.evaluate('DECJ',m),
int(m), m-int(m), 1, 0.)['delay']
rf_rv += rf_ephem.doppler_fraction(
eph1957.evaluate('RAJ',m)/15., eph1957.evaluate('DECJ',m),
int(m), m-int(m), 1, 0.)['frac']
import matplotlib.pylab as plt
plt.ion()
#plt.plot(t.utc.mjd, delay-rf_delay-d_orb)
plt.plot(t.utc.mjd, d_earth-rf_delay)
#plt.plot(t.utc.mjd, (rv-rf_rv-v_orb)*c.to(u.km/u.s).value)
plt.draw()
if False:
for utc, tdb1, tdb2 in zip(t.utc.mjd, t.tdb.jd1, t.tdb.jd2):
rf_tdb = rf_ephem.utc_to_tdb(int(utc), utc-int(utc))['tdb']
