Exemple #1
0
        /// Function to convert from UTC to sidereal time
        //  @param t         Epoch
        //
        //  @return sidereal time in hours.
        public static double UTC2SID(Time time)
        {
            // Hours of day (decimal)
            double h = time.SecondOfDay / 3600.0;

            // Fraction of day
            double frofday = time.SecondOfDay / 86400.0;

            // Compute Julian Day, including decimals
            double jd = time.Jd;

            // Temporal value, in centuries
            double tt = (jd - 2451545.0) / 36525.0;

            double sid = (24110.54841 + tt * ((8640184.812866) +
                                              tt * ((0.093104) - (6.2e-6 * tt))));

            sid = sid / 3600.0 + h;
            sid = MiscMath.Fmod(sid, 24.0);

            if (sid < 0.0)
            {
                sid += 24.0;
            }
            return(sid);
        }