Example #1
0
 private void PrintScr(PosData posData)
 {
     lblAlti.Text    = Convert.ToString(posData.Elevation);
     lblAzi.Text     = Convert.ToString(posData.Azimuth);
     lblDays.Text    = Convert.ToString(posData.Days);
     lblDec.Text     = Convert.ToString(posData.Declination);
     lblEot.Text     = Convert.ToString(posData.EoT);
     lblHra.Text     = Convert.ToString(posData.HRA);
     lblRise.Text    = Convert.ToString(posData.Sunrise);
     lblSet.Text     = Convert.ToString(posData.Sunset);
     lblTimeCor.Text = Convert.ToString(posData.TC);
 }
Example #2
0
        private void btnCalculate_Click(object sender, EventArgs e)
        {
            try
            {
                PosData posData = new PosData();

                posData.Latitude   = Convert.ToDouble(txtLat.Text);
                posData.Longtitude = Convert.ToDouble(txtLong.Text);

                posData.DtGMT = Convert.ToInt32(txtGmt.Text);
                posData.Date  = datePicker.Value;

                Calculate(posData);
                PrintScr(posData);
            }
            catch (Exception ex)
            {
                MessageBox.Show(ex.Message);
            }
        }
Example #3
0
        private void Calculate(PosData posData)
        {
            int year  = posData.Date.Year;
            int month = posData.Date.Month;
            int day   = posData.Date.Day;

            double hour   = Convert.ToDouble(posData.Date.Hour);
            double minute = Convert.ToDouble(posData.Date.Minute);

            if (CheckLeap(year))
            {
                months[1] = 29;
            }

            for (int i = 0; i < month - 1; i++)
            {
                posData.Days = posData.Days + months[i];
            }

            posData.Days = posData.Days + day;

            posData.LSTM = 15 * posData.DtGMT;
            posData.B    = Math.Round((posData.Days - 81) * (360d / 365d), 4, MidpointRounding.AwayFromZero);

            double eot_1 = 0;
            double eot_2 = 0;
            double eot_3 = 0;

            double degRad = posData.B * rad;

            eot_1 = Math.Round((9.87 * 2 * Math.Sin(degRad) * Math.Cos(degRad)), 4, MidpointRounding.AwayFromZero);
            eot_2 = Math.Round((-7.53 * Math.Cos(degRad)), 4, MidpointRounding.AwayFromZero);
            eot_3 = Math.Round((-1.5 * Math.Sin(degRad)), 4, MidpointRounding.AwayFromZero);

            posData.EoT = Math.Round((eot_1 + eot_2 + eot_3), 2, MidpointRounding.AwayFromZero);
            posData.TC  = Math.Round((4 * (posData.Longtitude - posData.LSTM) + posData.EoT), 2, MidpointRounding.AwayFromZero);

            posData.LT  = Math.Round((hour + (minute / 60d)), 4, MidpointRounding.AwayFromZero);
            posData.LST = Math.Round((posData.LT + (posData.TC / 60d)), 4, MidpointRounding.AwayFromZero);
            posData.HRA = Math.Round((15 * (posData.LST - 12d)), 2, MidpointRounding.AwayFromZero);

            posData.Declination = Math.Round((23.45d * Math.Sin(degRad)), 2, MidpointRounding.AwayFromZero);

            double degRadDec = posData.Declination * rad;
            double degRadLat = posData.Latitude * rad;
            double degRadHra = posData.HRA * rad;

            double elev_1 = 0;
            double elev_2 = 0;
            double elev   = 0;

            elev_1 = Math.Round((Math.Sin(degRadDec) * Math.Sin(degRadLat)), 4, MidpointRounding.AwayFromZero);
            elev_2 = Math.Round((Math.Cos(degRadDec) * Math.Cos(degRadLat) * Math.Cos(degRadHra)), 4, MidpointRounding.AwayFromZero);
            elev   = elev_1 + elev_2;

            posData.Elevation = Math.Round(((Math.Asin(elev) * 180d) / Math.PI), 2, MidpointRounding.AwayFromZero);


            double azi_1      = 0;
            double azi_2      = 0;
            double azi_3      = 0;
            double azi        = 0;
            double degRadElev = posData.Elevation * rad;

            azi_1 = Math.Round((Math.Sin(degRadDec) * Math.Cos(degRadLat)), 4, MidpointRounding.AwayFromZero);
            azi_2 = Math.Round((-1d * Math.Cos(degRadDec) * Math.Sin(degRadLat) * Math.Cos(degRadHra)), 4, MidpointRounding.AwayFromZero);
            azi_3 = azi_1 + azi_2;
            azi   = Math.Round((azi_3 / Math.Cos(degRadElev)), 4, MidpointRounding.AwayFromZero);

            posData.Azimuth = Math.Round((360 - (Math.Acos(azi) * 180d) / Math.PI), 2, MidpointRounding.AwayFromZero);

            //MessageBox.Show(posData.Days.ToString());
        }