def radvel(self, sdate, time, radeg, decdeg) : # 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] uttime,tf1 = self.sex2deg(time) # time in hours #radeg,tf1 = sex2deg(rastring) #if tf1 == 0 : radeg = radeg*15.0 #decdeg,tf1 = sex2deg(decstring) # print '<br>RA=%8.4f deg, DEC=%8.4f deg' % (radeg,decdeg) mjdt = mjd0 + uttime/24.0 # print 'Using ', sdate, mjdt, uttime, radeg,decdeg rr = self.rvel( mjdt, radeg, decdeg) # ((mjdd, last*12.0/math.pi, Rgeo, totalhelio,totalbary, totallsrk,totalgal)) return rr
def current_MJD(): """ current_MJD(): Return the current MJD accurate to ~1 sec. """ YY, MM, DD, hh, mm, ss, wday, yday, isdst = time.gmtime() mjd_f, J = s.sla_dtf2d(hh, mm, ss) mjd_i, J = s.sla_cldj(YY, MM, DD) return mjd_i + mjd_f
def datetime2mjd_utc(d): """Converts a passed datetime object in UTC to the equivalent Modified Julian Date (MJD), which is returned""" # Compute MJD for UTC (mjd, status) = S.sla_cldj(d.year, d.month, d.day) if status != 0: return None (fday, status) = S.sla_dtf2d(d.hour, d.minute, d.second + (d.microsecond / 1e6)) if status != 0: return None mjd_utc = mjd + fday return mjd_utc
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])