spline fit and evaluate errors
apply just a smoothing kernel [0.25,0.5,0.25] to the flat that suffer from artificial ringing
error propagation of sky model stat. uncertainty (only diagonal but very small anyway)
clipped mean and propagate calibration errors