;
; Returns the position of Venus (in SOHO coordinates) as a function of time
; during the June 8, 2004 coronal transit.
;
; Positions are returned as an [N,2] array, where N is the number of times
; given in the input array (also works for scalar input, which will return a
; [1,2] array).
;
; Needs: venus_ephem.pro, and the ephemeris file SOHO_venus_snp.txt
;
; Use: print,venus_point('2003/05/07 10:34:00')
;
; Input times can be anything understood by the anytim2tai function.
;
FUNCTION venus_point,t
venus_ephem,x,y,d,time,timestr=timestr
i = (where(time EQ 0))[0]
t = anytim2tai(t)
t0 = anytim2tai(timestr[i])
delta_t_min = (t-t0)/60.0
xt = interpol(x,time,delta_t_min,/quad)
yt = interpol(y,time,delta_t_min,/quad)
return,[[xt],[yt]]
END