def gettraj(t, t0, u0, te, piex, piey, alpha, delta, t0par, qlat, qlong, hjds, ras, decs, diss, SPITZ): tau = (t - t0) / te beta = u0 qn, qe = getpar.geta(t, alpha, delta, t0par) if SPITZ == True: qnp, qep, qr = getpar.gets(t, hjds, ras, decs, diss, alpha, delta) else: qnp, qep = getpar.gett(t, qlat, qlong, alpha, delta) qn += qnp qe += qep dtau = piex * qn + piey * qe dbeta = piex * qe - piey * qn return tau + dtau, beta + dbeta
def getqnqe(dates, alpha, delta, t0par, qlat, qlong, SPITZ): qns, qes = [], [] qnqe = [] for t in dates: qn, qe = getpar.geta(t, alpha, delta, t0par) if SPITZ == True: qnp, qep, qr = getpar.gets(t, hjds, ras, decs, diss, alpha, delta) else: qnp, qep = getpar.gett(t, qlat, qlong, alpha, delta) qn += qnp qe += qep qnqe.append([qn, qe]) qnqe = np.array(qnqe).T return qnqe