%N zedvhat %F $MIR/src/subs/zed.for
%D Calculate Zeeman estimate of a true V spectrum
%P nebk
%: zeeman
%B
      subroutine ZedVHat(mode,IHat,m,n,a,b,VHat)

      character mode*(*)
      integer m,n
      real IHat(m,n), VHat(M,n)
      real a,b
  ---------------------------------------------------------------------
  This returns an estimate of the true V spectrum, given a, b and IHat

  Inputs:
    mode       Flags to the algorithm. The only recognized flag is '2',
               for two sided derivative.
    IHat       Estimate of true noiseless spectrum
    m          Number of channels.
    n          Number of spectra.
    a,b        Alpha and beta parameters of interest.

  Output:
    VHat       An estimate of the true value of V.  Note that the first
               and possibly last (depending on derivative type) is NOT
               filled.
