coords = u[:, (None), (None)] * np.linspace(0, 9, 10)[(None), :, (None)] + v[:,
    (None), (None)] * np.linspace(0, 2.5, 10)[(None), (None), :]
Fast interpolation of regularly sampled 3D data with different intervals in x,y, and z