%N calstoke %F $MIR/src/subs/calstoke.for
%D Flux characteristics of selected calibrators.
%P nebk,
%: utilities
%B
      subroutine calstoke(source,stokes,freq,flux,nchan,ierr)

      implicit none
      character source*(*),stokes*1
      integer ierr,nchan
      double precision freq(nchan)
      real flux(nchan)

     Return IQUV for known calibrator sources.  Spectra generally
     described by third or fourth order polynomials in either the
     log(f)-log(S), f-S or log(f)-S domains.

   Input:
     source    Source name.
     stokes    One of IQUV.
     freq      Frequency (GHz)
     nchan     Number of channels.
   Output:
     flux      The source flux for those frequencies for that Stokes
               parameter.
     ierr      Status return: 0 = All OK.
                              1 = OK, but extrapolation performed.
                              2 = No match found.
