%N getfreq %F $MIR/src/subs/getfreq.for
%D Return frequency & increment for given image spectral pixel
%P nebk
%: utilities
%B
        subroutine getfreq (tin, pix, ifrq, freq, finc, ierr)

        integer   tin, ifrq, ierr
        real      pix
        double precision freq, finc
  ---------------------------------------------------------------------
  Return the frequency of a spectral pixel, regardless of spectral axis
  type, or which axis is the spectral axis.  Works for images only.

  Inputs:
    tin    i   Handle of dataset
    pix    r   Pixel of interest

  Output:
    ifrq   i   The spectral axis number.
    freq   d   Frequency in GHz.
    finc   d   Frequency increment in GHz.
    ierr   i   0 -> OK,
               1 -> no spectral axis,
 History:
    Refer to the RCS log, v1.1 includes prior revision information.

 $Id$
