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