def test_observer_test():
# Setup observatory location - in this case, Parkes Australia
(latitude, longitude, elevation) = ('-32.998370', '148.263659', 100)
ov = GSMObserver2016()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(1400)
d = ov.view(logged=True)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(1400)
d = ov.view(logged=True)
plt.show()
def test_gsm_observer():
""" Test GSMObserver() is working
"""
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.view(logged=True)
ov.view_observed_gsm(logged=True)
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver2016()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.view(logged=True)
ov.view_observed_gsm(logged=True)
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = LFSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.view(logged=True)
ov.view_observed_gsm(logged=True)
def test_observed_mollview():
""" Generate animated maps of observing coverage over 24 hours """
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
obs = []
if not os.path.exists('generated_sky'):
os.mkdir('generated_sky')
for ii in range(0, 24, 4):
ov.date = datetime(2000, 1, 1, ii, 0)
ov.generate(50)
sky = ov.view_observed_gsm(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/galactic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'E'], min=9, max=20)
plt.savefig('generated_sky/ecliptic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'C'], min=9, max=20)
plt.savefig('generated_sky/equatorial-%02d.png' % ii)
plt.close()
ov.view(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/ortho-%02d.png' % ii)
plt.close()
print(ii)
os.system('convert -delay 20 generated_sky/ortho-*.png ortho.gif')
os.system('convert -delay 20 generated_sky/galactic-*.png galactic.gif')
os.system('convert -delay 20 generated_sky/ecliptic-*.png ecliptic.gif')
os.system(
'convert -delay 20 generated_sky/equatorial-*.png equatorial.gif')
def test_generate_with_and_without_args():
""" Test generating without frequency argument """
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.date = datetime(2000, 1, 1, 22, 0)
ov.generate()
ov.generate(51)
ov.generate(52)
ov.generate()
now = Time(datetime.now())
ov.generate(obstime=now)
ov.generate(obstime=now)
ov.generate(obstime=now, freq=53)
ov.generate(obstime=now, freq=52)