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) plt.draw() plt.show()
# # 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) plt.draw() plt.show()
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() 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) plt.draw() plt.show()
# # 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'] print utc, tdb, rf_tdb, tdb1-0.5-int(tdb1-0.5)+tdb2-rf_tdb if False: tdb = np.linspace(0.,1.,5)