%N colmn %F $MIR/src/subs/co.for
%D Convert celestial coordinates to direction cosines.
%P rjs
%: coordinates
%B
      subroutine coLMN(lu,in,x1,lmn)

      integer lu
      character in*(*)
      double precision x1(*), lmn(3)
  ---------------------------------------------------------------------
  Get the direction cosines corresponding to a particular coordinate.

  Input:
    lu         Handle of the coordinate object.
    in         As with coCvt
    x1         As with coCvt
  Output:
    lmn        The direction cosines (with respect to the reference
               position) of the celestial coordinate given by x1.
