קאווע שטיבל

קוד: וועל אויס אלע

        private int Calculate(Direction direction)
        {
            /* doy (N) */
            int N = date.DayOfYear;

            /* appr. time (t) */
            double lngHour = longitude / 15.0;

            double t;

            if (direction == Direction.Sunrise)
                t = N + ((6.0 - lngHour) / 24.0);
            else
                t = N + ((18.0 - lngHour) / 24.0);

            /* mean anomaly (M) */
            double M = (0.9856 * t) - 3.289;

            /* true longitude (L) */
            double L = M + (1.916 * Math.Sin(Deg2Rad(M))) + (0.020 * Math.Sin(Deg2Rad(2 * M))) + 282.634;
            L = FixValue(L, 0, 360);

            /* right asc (RA) */
            double RA = Rad2Deg(Math.Atan(0.91764 * Math.Tan(Deg2Rad(L))));
            RA = FixValue(RA, 0, 360);

            /* adjust quadrant of RA */
            double Lquadrant = (Math.Floor(L / 90.0)) * 90.0;
            double RAquadrant = (Math.Floor(RA / 90.0)) * 90.0;
            RA = RA + (Lquadrant - RAquadrant);

            RA = RA / 15.0;

            /* sin cos DEC (sinDec / cosDec) */
            double sinDec = 0.39782 * Math.Sin(Deg2Rad(L));
            double cosDec = Math.Cos(Math.Asin(sinDec));

            /* local hour angle (cosH) */
            double cosH = (Math.Cos(Deg2Rad((double)zenith / 1000.0f)) - (sinDec * Math.Sin(Deg2Rad(latitude)))) / (cosDec * Math.Cos(Deg2Rad(latitude)));

            /* local hour (H) */
            double H;

            if (direction == Direction.Sunrise)
                H = 360.0 - Rad2Deg(Math.Acos(cosH));
            else
                H = Rad2Deg(Math.Acos(cosH));

            H = H / 15.0;

            /* time (T) */
            double T = H + RA - (0.06571 * t) - 6.622;

            /* universal time (T) */
            double UT = T - lngHour;

            UT += utcOffset;  // local UTC offset

            if (daylightChanges != null)
                if ((date > daylightChanges.Start) && (date < daylightChanges.End))
                    UT += daylightChanges.Delta.TotalHours;

            UT = FixValue(UT, 0, 24);

            return (int)Math.Round(UT * 3600);  // Convert to seconds
        }