Create a function that takes the raw datacube and applies a known smile shift across spatial columns.

Use cubic spline interpolation to induce a subpixel shift:

from scipy.interpolate import CubicSpline