%N coradec %F $MIR/src/subs/co.for
%D Create a simple RA/DEC coordinate system.
%P rjs
%: coordinates
%B
      subroutine coRaDec(lu,proj,ra0,dec0)

      integer lu
      character proj*(*)
      double precision ra0,dec0
  ---------------------------------------------------------------------
  Create a simple RA/DEC coordinate system.

  Input:
    proj       Projection code: NCP, SIN, etc.
    ra0,dec0   RA,DEC of the reference point.
  Output:
    lu         Handle of the output coordinate object.
