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])