%N coaxtype %F $MIR/src/subs/co.for
%D Get information about the world coordinates on an axis.
%P mrc
%: coordinates
%B
      subroutine coAxType(lu,iax,axtype,wtype,units)

      integer   lu, iax
      character axtype*(*), wtype*(*), units*(*)
  ---------------------------------------------------------------------
  Get information about the world coordinates on an axis.  Provides a
  convenient interface to coCtype, returning only the information that
  should be required by the outside world.

  Input:
    lu         Handle of the coordinate object.
  In/out:
    iax        The axis number.  If zero, then it will be set from
               axtype.
    axtype     Generic world coordinate type:
                 'longitude'  for RA,  GLON and ELON axes.
                 'latitude'   for DEC, GLAT and ELAT axes.
                 'spectral'   for spectral axes.
               Not reset if iax was zero on input.
  Output:
    wtype      World coordinate type, i.e. ctype with the algorithm
               code stripped off, e.g. 'RA---NCP' yields 'RA', and
               'VOPT-F2W' yields 'VOPT'.
    units      Miriad internal units for this world coordinate in a
               form suitable for reporting, e.g. 'VOPT-F2W' would
               yield 'km/s'.  May be used to identify angular,
               frequency, velocity, etc. axes.
