%N azel %F $MIR/src/subs/ephem.for
%D Calculate azimuth and elevation from ra/dec.
%P rjs
%: utilities
%B
        subroutine azel(obsra,obsdec,lst,latitude,az,el)

        implicit none
        double precision obsra,obsdec,lst,latitude
        double precision az,el

  This computes the azimuth and elevation.

  Input:
    obsra )    Apparent RA and DEC of the source of interest (radians).
    obsdec)
    lst        Local sidereal time (radians).
    latitude   Observatory geodetic latitude (radians).

  Output:
    az,el      Azimuth and elevation angle (radians).
