def test_works(self) :
UT = '2000-01-01T12:00:00.00'
el = 90.
az = 50.5
ra, dec = utils.elaz2radecGBT(el, az, UT)
self.assertAlmostEqual(dec, 38.43312, 1) #GBT Latitude
el = 38.43312
az = 0.
ra, dec = utils.elaz2radecGBT(el, az, UT)
self.assertAlmostEqual(dec, 90, 1)
el = 90 - 38.43312
az = 180.
ra, dec = utils.elaz2radecGBT(el, az, UT)
self.assertAlmostEqual(dec, 0, 1)
def calc_pointing(self) :
"""Calculates the telescope pointing.
Should only be used if in alt/az mode.
At every time the Ra and Dec of the telescope time is calculated.
These are stored as attributes (not fields) named ra and dec. This
requires the fields 'CRVAL3', 'CRVAL2' and 'DATE-OBS' to be set.
"""
self.ra = sp.zeros(self.dims[0])
self.dec = sp.zeros(self.dims[0])
for ii in range(self.dims[0]) :
self.ra[ii], self.dec[ii] = utils.elaz2radecGBT(
self.field['CRVAL3'][ii],
self.field['CRVAL2'][ii],
self.field['DATE-OBS'][ii])
