private static void FillInPhenomenaDetails(ref AASGalileanMoonDetail detail) { double Y1 = 1.071374 * detail.ApparentRectangularCoordinates.Y; double r = Y1 * Y1 + detail.ApparentRectangularCoordinates.X * detail.ApparentRectangularCoordinates.X; if (r < 1) { if (detail.ApparentRectangularCoordinates.Z < 0) { //Satellite nearer to Earth than Jupiter, so it must be a transit not an occultation detail.bInTransit = true; detail.bInOccultation = false; } else { detail.bInTransit = false; detail.bInOccultation = true; } } else { detail.bInTransit = false; detail.bInOccultation = false; } }
public AASGalileanMoonsDetails() { Satellite1 = new AASGalileanMoonDetail(); Satellite2 = new AASGalileanMoonDetail(); Satellite3 = new AASGalileanMoonDetail(); Satellite4 = new AASGalileanMoonDetail(); }
public static AASGalileanMoonsDetails Calculate(double JD, bool bHighPrecision) { //Calculate the position of the Sun double sunlong = AASSun.GeometricEclipticLongitude(JD, bHighPrecision); double sunlongrad = AASCoordinateTransformation.DegreesToRadians(sunlong); double beta = AASSun.GeometricEclipticLatitude(JD, bHighPrecision); double betarad = AASCoordinateTransformation.DegreesToRadians(beta); double R = AASEarth.RadiusVector(JD, bHighPrecision); //Calculate the the light travel time from Jupiter to the Earth double DELTA = 5; double PreviousEarthLightTravelTime = 0; double EarthLightTravelTime = AASElliptical.DistanceToLightTime(DELTA); double JD1 = JD - EarthLightTravelTime; bool bIterate = true; double x; double y; double z; double l; double lrad; double b; double brad; double r; while (bIterate) { //Calculate the position of Jupiter l = AASJupiter.EclipticLongitude(JD1, bHighPrecision); lrad = AASCoordinateTransformation.DegreesToRadians(l); b = AASJupiter.EclipticLatitude(JD1, bHighPrecision); brad = AASCoordinateTransformation.DegreesToRadians(b); r = AASJupiter.RadiusVector(JD1, bHighPrecision); x = r * Math.Cos(brad) * Math.Cos(lrad) + R * Math.Cos(sunlongrad); y = r * Math.Cos(brad) * Math.Sin(lrad) + R * Math.Sin(sunlongrad); z = r * Math.Sin(brad) + R * Math.Sin(betarad); DELTA = Math.Sqrt(x * x + y * y + z * z); EarthLightTravelTime = AASElliptical.DistanceToLightTime(DELTA); //Prepare for the next loop around bIterate = (Math.Abs(EarthLightTravelTime - PreviousEarthLightTravelTime) > 2e-6); //2e-6 corresponds to 0.17 of a second if (bIterate) { JD1 = JD - EarthLightTravelTime; PreviousEarthLightTravelTime = EarthLightTravelTime; } } //Calculate the details as seen from the earth AASGalileanMoonsDetails details1 = CalculateHelper(JD, sunlongrad, betarad, R, bHighPrecision); AASGalileanMoonDetail details1Satellite1 = details1.Satellite1; AASGalileanMoonDetail details1Satellite2 = details1.Satellite2; AASGalileanMoonDetail details1Satellite3 = details1.Satellite3; AASGalileanMoonDetail details1Satellite4 = details1.Satellite4; FillInPhenomenaDetails(ref details1Satellite1); FillInPhenomenaDetails(ref details1Satellite2); FillInPhenomenaDetails(ref details1Satellite3); FillInPhenomenaDetails(ref details1Satellite4); //Calculate the the light travel time from Jupiter to the Sun JD1 = JD - EarthLightTravelTime; l = AASJupiter.EclipticLongitude(JD1, bHighPrecision); lrad = AASCoordinateTransformation.DegreesToRadians(l); b = AASJupiter.EclipticLatitude(JD1, bHighPrecision); brad = AASCoordinateTransformation.DegreesToRadians(b); r = AASJupiter.RadiusVector(JD1, bHighPrecision); x = r * Math.Cos(brad) * Math.Cos(lrad); y = r * Math.Cos(brad) * Math.Sin(lrad); z = r * Math.Sin(brad); DELTA = Math.Sqrt(x * x + y * y + z * z); double SunLightTravelTime = AASElliptical.DistanceToLightTime(DELTA); //Calculate the details as seen from the Sun AASGalileanMoonsDetails details2 = CalculateHelper(JD + SunLightTravelTime - EarthLightTravelTime, sunlongrad, betarad, 0, bHighPrecision); AASGalileanMoonDetail details2Satellite1 = details2.Satellite1; AASGalileanMoonDetail details2Satellite2 = details2.Satellite2; AASGalileanMoonDetail details2Satellite3 = details2.Satellite3; AASGalileanMoonDetail details2Satellite4 = details2.Satellite4; FillInPhenomenaDetails(ref details2Satellite1); FillInPhenomenaDetails(ref details2Satellite2); FillInPhenomenaDetails(ref details2Satellite3); FillInPhenomenaDetails(ref details2Satellite4); //Finally transfer the required values from details2 to details1 details1.Satellite1.bInEclipse = details2.Satellite1.bInOccultation; details1.Satellite2.bInEclipse = details2.Satellite2.bInOccultation; details1.Satellite3.bInEclipse = details2.Satellite3.bInOccultation; details1.Satellite4.bInEclipse = details2.Satellite4.bInOccultation; details1.Satellite1.bInShadowTransit = details2.Satellite1.bInTransit; details1.Satellite2.bInShadowTransit = details2.Satellite2.bInTransit; details1.Satellite3.bInShadowTransit = details2.Satellite3.bInTransit; details1.Satellite4.bInShadowTransit = details2.Satellite4.bInTransit; return(details1); }