orbitalPhaseAtStart again back to its meaning of initial mean anomaly
function to calculate the true anomaly where needed