NAME: lsidtim PURPOSE: Compute local sidereal time at a given longitude and time. DESCRIPTION: This routine is based on the algorithms from p. 83 of "Astronomical Algorithms" by J. Meeus, 1st edition. CATEGORY: Astronomy CALLING SEQUENCE: lsidtim,jd,lon,sidtim [,UT=ut] INPUTS: jd - Julian Date (double precision), scalar or vector. lon - West longitude of observatory in radians (scalar). OPTIONAL INPUT PARAMETERS: KEYWORD INPUT PARAMETERS: UT - Time, in hours to add to JD to get the correct Universal Time. That the same as Universal Time minus the Local Time. OUTPUTS: lst - Local sidereal time for each of the input times (radians). KEYWORD OUTPUT PARAMETERS: COMMON BLOCKS: SIDE EFFECTS: RESTRICTIONS: PROCEDURE: MODIFICATION HISTORY: 94/05/05 - Written by Marc W. Buie, Lowell Observatory 2002/01/08, MWB, 2-d input on jd is now allowed. 2002/09/05, MWB, total rewrite to use modern formula