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