position method
Geocentric position of the observer in km, in the frame of the true
equator and equinox of date at jdTT.
Implementation
Vector3 position(double jdTT, {double? deltaTSeconds}) {
final theta = localSiderealTime(jdTT, deltaTSeconds: deltaTSeconds);
return Vector3(rhoCos * math.cos(theta), rhoCos * math.sin(theta), rhoSin) *
earthEquatorialRadius;
}