;+
; NAME:
;     oc_ephicy_et.Pro
;
; PURPOSE:
;     Returns a structure with the following information:
;     {et, ra, dec, radot, decdot, ephinfo} where ephinfo
;     contains the observer range and observer range rate.
;
; DESCRIPTION:
;     Calculates the ra, dec, radot, decdot, observer range
;     and range rate using the NAIF ICY routines.
;
; CALLING SEQUENCE:
;     eph = oc_ephicy(target, observer, time, dt, neph)
;
; INPUTS:
;     target - the target for the ephemeris, i.e. "Pluto"
;     observer - the observer location, i.e. "Earth"
;     time - either a list of UTC strings or a list of et
; OUTPUTS:
;      eph - an ephemeris structure (see single_eph.pro)
;
; REVISON HISTORY:
;       11-May-2006 LAY SwRI
;       02-July-2006 LAY SwRI add kernel keyword to single_eph  to avoid extra calls
;                             to naif_laodedkernels
;       21-Sept-2006 LAY SwRI return correct neph
;-

function oc_ephicy_et, target, observer, et, neph, $
                    verbose = verbose, abcorr = abcorr, frame=frame, $
                    lonlatalt=lonlatalt

       naifinit

; Number of steps
       neph = n_elements(et)

       eph1 = single_eph( target, et[0], neph1, observer=observer, $
                          verbose = verbose, abcorr = abcorr, frame=frame, $
                          lonlatalt=lonlatalt)
       kernels = eph1.naif.kernels
       eph = replicate(eph1, neph)

       FOR i=1L, neph - 1 DO BEGIN
           eph1 = single_eph( target, et[i], neph1, observer=observer,  $
                          verbose = verbose, abcorr = abcorr, frame=frame, $
                          lonlatalt=lonlatalt, kernels=kernels)
           eph[i] = eph1
       ENDFOR

       return, eph
       
 END