import numpy as np # Constants (in SI units) c = 299792458.0 # Speed of light (m/s) rE = 6.371E6 # Radius of Earth (m) rOrb = 1.496E11 # Radius of Earth's orbit (m) Td = 86164.092 # Sidereal day length (s) TYr = 366.2564 * Td # Seconds in a sidereal year tilt = 23.43928 # Degrees which Earth tilts relative orbit sampR= 60 # Position of the observatory detName = "L1" if detName == 'H1': lat = 46.4541048 # Degrees lng = -119.4175379 # Degrees, probably in the wrong system beta = 216.3 # Degrees from N/E detector, Angles from google maps elif detName == 'L1': # Note: to align to H1 data, use t0 = 3500 lat = 30.5021 # Degrees lng = 90.7479 # Degrees, not really useful; Probably in wrong system beta = 198.0 # Degrees from N/E detector, Angles from google maps else: raise Exception("Detector not specified.") # Data for the source (here: Crab Nebula) pulsName = 'J0534+2200' alp = 5.5755389 # Right ascension is in hours dlt = 22.0145 # Declination is in degrees iota = 61.3 # Degrees psi = 124.0 # Degrees phi0 = 0 # Phase angle of wave nu = 29.946923 # Rotational velocity of body (Hz) nuD = -3.77535E-10 # Rotational acceleration of body (Hz/s) nuDD = 1.1147E-20 # Rotational jerk of body (Hz/s^2) # Convert to radians alp *= np.pi / 12 dlt *= np.pi / 180 lat *= np.pi / 180 lng *= np.pi / 180 beta *= np.pi / 180 iota *= np.pi / 180 psi *= np.pi / 180 tilt *= np.pi / 180