Use better interpolation in generating sky1D

The current code to generate sky1d uses linear interpolation. We need to do better.

Because of dithering of wavelengths relative to the pixel grid we should not have problems with under sampling. A least-squares spline might be a way to get something running quickly