Added getAnglesByYear helper function to phase angle calculation