Now using projection of change vector onto position radius vector via dot product to set phi,theta components of infinitesimal rotation increment.

Have retained debug code, will remove on next commit.
Removed debug code.