%N coprjset %F $MIR/src/subs/co.for
%D Set celestial projection type.
%P rjs
%: coordinates
%B
      subroutine coPrjSet(lu, code, npv, pv)

      integer lu, npv
      double precision pv(30)
      character code*(*)
  ---------------------------------------------------------------------
  Set the coordinate system celestial projection type.

  Input:
    lu         Handle of the coordinate system.
    code       Projection code: NCP, SIN, etc. or blank, which on
               interpretation will be translated to CAR.  Can also be
               '-' to leave it alone in case all you want to do is set
               projection parameters.
    npv        Number of elements in pv.
    pv         Array of projection parameters.
