vel_model seems to work! :)  Has error due to the innaccuracies in acos() as well as floating point precision.