Thanks, but unfortunately this doesn’t really do anything if the desiredTimeStep < deltaTime apart from forcing the velocity to only update in multiples of desiredTimeStep. It still only updates at variable deltaTime intervals, with engine physics still using the velocity in between those updates, thus giving different results for different deltaTime values.
I posted an image as a comment to my OP which illustrates what I think the problem is, I hope.