Пример #1
0
ETNEARPRG = @2013-10-09T00:00:00

SPACECRAFT = 'JUNO'
TARGET = 'EARTH'

\begintext
"""

import SpiceyPy as spice

dpr = spice.dpr()

spice.furnsh('.'.join(__file__.split('.')[:-1]+['py']))

target = spice.gcpool('TARGET',0,1,99)[0]
spacecraft = spice.gcpool('SPACECRAFT',0,1,99)[0]
etnearca = spice.gdpool('ETNEARPRG',0,1)[0]

targetID = spice.bods2c(target)

et = etnearca
state,lt = spice.spkezr(spacecraft,et,'J2000','NONE',target)

dt = spice.vdot(state[:3],state[3:]) / spice.vdot(state[3:],state[3:])
oldet = et
et = et - dt
while abs(et-oldet)>0.:
  state,lt = spice.spkezr(spacecraft,et,'J2000','NONE',target)
  dt = spice.vdot(state[:3],state[3:]) / spice.vdot(state[3:],state[3:])
  oldet = et