def testtpmstate(self):
"""tpm_state(x) => name of state."""
self.assertEqual(tpm.tpm_state(1).strip(), "Helio. mean FK4")
self.assertEqual(tpm.tpm_state(tpm.TPM_S20).strip(),
"Topo. obs. HA/Dec")
def testtpmstate(self):
"""tpm_state(x) => name of state."""
self.assertEqual(tpm.tpm_state(1).strip(), "Helio. mean FK4")
self.assertEqual(tpm.tpm_state(tpm.TPM_S20).strip(), "Topo. obs. HA/Dec")
tstate.utc = tpm.J2000
tstate.lon = tpm.d2r(-111.598333)
tstate.lat = tpm.d2r(31.956389)
tstate.alt = 2093.093
tstate.delta_ut = tpm.delta_UT(tstate.utc)
tpm.tpm_data(tstate, tpm.TPM_ALL)
v6 = tpm.V6S()
v6.r = 1e9
v6.alpha = ra
v6.delta = de
pvec[s1] = v6.s2c()
s2 = i
tpm.tpm(pvec, s1, s2, ep, eq, tstate)
v6 = pvec[s2].c2s()
ra1 = v6.alpha
de1 = v6.delta
ra1_d = tpm.r2d(ra1)
if ra1_d < 0.0 : ra1_d += 360.0
de1_d = tpm.r2d(de1)
if de1_d < 0.0 : de1_d += 360.0
s = "{0:02d}-{1:02d} {2:<17s} {3:s} {4:s} {5:8.4f} {6:8.4f}"
print s.format(s1, s2, tpm.tpm_state(s2),
tpm.fmt_alpha(ra1), tpm.fmt_delta(de1), ra1_d,
de1_d)
