Exemplo n.º 1
0
        /// <summary>
        /// Положение и скорость невращающейся Земли.
        /// </summary>
        ///
        /// <param name="julianDate"> Юлианская дата.</param>
        ///
        /// <returns>
        /// Положение и скорость невращающейся Земли.
        /// </returns>
        private static (Vector Vector, Vector Velocity, Basis Basis) EarthParams(double julianDate)
        {
            var pvh = new double[2, 3];
            var pvb = new double[2, 3];

            IAUSOFA.iauEpv00(Date.J2000, julianDate - Date.J2000, pvh, pvb);
            var vector   = new Vector(pvh[0, 0], pvh[0, 1], pvh[0, 2]);
            var velocity = new Vector(pvh[1, 0], pvh[1, 1], pvh[1, 2]);

            return(vector * IAUSOFA.DAU, velocity *IAUSOFA.DAU / Date.JD_TO_SECOND, new Basis());
        }