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