;
; 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