Fixed a bug in readimfits which prevented it from dropping the frequency
and stokes axes from the wcs info.
find now multiplies out the beam from the flux so that the units are
just in flux.
Fixed a bug that caused the RA and Dec to be output only for the first
source found.
Added function match_source_lists to match the source lists produced by
the find function and create a multi-image source catalog.