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