;
; Returns the position of Mercury (in SOHO coordinates) as a function of time
; during the May 7, 2003 disk 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: merc_ephem.pro, and the ephemeris file merc_snp_L_5s_V041503.txt
;
; Use: print,merc_point('2003/05/07 10:34:00')
;
; Input times can be anything understood by the anytim2tai function.
;

FUNCTION merc_point,t
  
  merc_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)
  yt = interpol(y,time,delta_t_min)
  
  return,[[xt],[yt]]
  
END