def update_azza(self): """ update_azza(): Update the AZ and ZA based on the current time with the vegas_status instance. """ (iptr, ang, stat) = s.sla_dafin(self['RA_STR'].replace(':', ' '), 1) self.update("RA", ang*15.0*RADTODEG) (iptr, ang, stat) = s.sla_dafin(self['DEC_STR'].replace(':', ' '), 1) self.update("DEC", ang*RADTODEG) MJD = astro.current_MJD() az, za = astro.radec_to_azza(self['RA'], self['DEC'], MJD, scope='GBT') self.update("AZ", az) self.update("ZA", za)
def update_azza(self): """ update_azza(): Update the AZ and ZA based on the current time with the nuppi_status instance. """ (iptr, ang, stat) = s.sla_dafin(self['RA_STR'].replace(':', ' '), 1) self.update("RA", ang * 15.0 * RADTODEG) (iptr, ang, stat) = s.sla_dafin(self['DEC_STR'].replace(':', ' '), 1) self.update("DEC", ang * RADTODEG) MJD = current_MJD() az, za = radec_to_azza(self['RA'], self['DEC'], MJD, scope='NRT') self.update("AZ", az) self.update("ZA", za)
def radvellist(self, sdate, rastring, decstring): # decode the date, convert to MJD ss = sdate.split('-') year = int(ss[0]) mont = int(ss[1]) day = int(ss[2]) mjdv = s.sla_cldj(year, mont, day) mjd0 = mjdv[0] rstat = mjdv[1] raval = s.sla_dafin(rastring, 1) # convert to radians dcval = s.sla_dafin(decstring, 1) rarad = raval[1] * 15. dcrad = dcval[1] # degrees radeg = rarad * 180.0 / math.pi decdeg = dcrad * 180.0 / math.pi # print "" # print ' %3s %7s %7s %7s %13s %13s %13s' % \ # ('UTHr', 'LST', 'Geo', 'Helio', 'Barycentric', 'LSRK', 'Galacto') # loop through 25 hours for i in range(0, 25): hour = i dfrac = hour / 24.0 mjd1 = mjd0 + dfrac rr = self.rvel(mjd1, radeg, decdeg) # ((mjdd, last*12.0/math.pi, Rgeo, totalhelio,totalbary, totallsrk,totalgal)) print(' %3d %7.3f %7.3f %7.3f %13.3f %13.3f %13.3f' % \ (hour,rr[1],rr[2],rr[3],rr[4],rr[5],rr[6]))