Add method to multiply facet image data with scalar value for each spectral term

This functionality will come in handy to compute flat noise/flat gain (facet) images in wsclean.