Draft: Update interpolation function/method

interp1d is depcrecated and there is concern that passing y-values that are complex valued does not work properly. Instead use numpy.interp and interpolate real and imaginary terms separately before combining

This should hopefully address a potential bug or weird behaviour in the uncertainty budgets for H1 because of the way a narrow resonant feature in the suspension model together with interpolation is calculated.

Closes #215