%N zedihat %F $MIR/src/subs/zed.for
%D Calculate Zeeman estimate of a true I spectrum
%P nebk
%: zeeman
%B
      subroutine ZedIHat(mode,ispect,vspect,m,n,a,b,ihat)

      character mode*(*)
      integer m,n
      real ispect(m,n),vspect(m,n),ihat(m,n)
      real a,b
  ---------------------------------------------------------------------
  This returns an estimate of the true I spectrum, given a and b.

  Inputs:
    mode       Flags to the algorithm. The only recognized flag is '2',
               for two sided derivative.
    ispect,vspect Observed I and V spectra.
    m          Number of channels.
    n          Number of spectra.
    a,b        Alpha and beta parameters of interest.

  Output:
    ihat       An estimate of the true value of I.
