%N jullst %F $MIR/src/subs/ephem.for
%D Convert Julian date to local mean sidereal time.
%P rjs
%: utilities
%B
        subroutine Jullst(jday,long,lst)

        implicit none
        double precision jday,lst,long

  Convert from Julian day to local mean sidereal time.

  Reference: Explanatory Supplement to the Astronomical Almanac, p50-52.
  Accuracy appears to be 0.01 sec of time.

  Input:
    jday       Julian day of interest.
    long       Observatory longitude (radians). East of Greenwich is
               positive.
  Output:
    lst        Local mean sidereal time (radians), in the range [0,2*pi].
