%N w2wsco %F $MIR/src/subs/cosubs.for
%D Convert NEBK style coordinate for a single axis
%P nebk
%: coordinates
%B
      subroutine w2wsco (lun, iax, typei, win, typeo, wout)

      integer   lun, iax
      double precision win, wout
      character typei*(*), typeo*(*)
  ---------------------------------------------------------------------
  Convert one NEBK style coordinate with the COCVT routines.
  Coordinates for the other axes are assumed to be at the
  reference pixel

  Input
    lun     Handle of open file
    iax     Axis
    typei   Input coordinate type, Should be from list
               'hms',    'dms',    'arcsec', 'arcmin', 'absdeg',
               'reldeg', 'abspix', 'relpix', 'absghz', 'relghz',
               'abskms', 'relkms', 'absnat', 'relnat', 'none'
    win     Coordinate to be converted
               'hms', 'dms' in radians
               '*  deg'     in degrees
               'arcsec'     in arcsec
               'arcmin'     in arcmin
               '*  pix'     in pixels
               '*  ghz'     in GHz
               '*  kms'     in km/s
               '*  nat'     in natural axis coordinates
    typeo   Output coordinate type from above list
  Output
    wout    Converted output coordinate
