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