I would be dangerous if I could count...

Shoot couple of rays after interpolation, results in more accurate output especially for long phases, like PKPPKP, especially for paths where small errors in ray parameter due to interpolation result in rays not arriving at correct distance. This does make calculations slightly slower. Also allow receiver at depth and new phases Ped and Sed for "exclusively downgoing" for cases where receiver is below source. Also usual bug fixes...