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