Ошибка или нет ? У меня или в StarCalc? Помогите разобраться...

[βλακασ]

Собственно, написал подобие программы для вычисления восходов/закатов Солнца и Луны. Точность, для оценки которой использовал StarCalc, меня устроила
вполне (несоответствие не более 1 мин.) Но вот только не понял я одного ньюанса :
Если взять в качестве места наблюдения город Норильск (координаты 69,33 СШ, 88,22 ВД) и даты близкие к окончанию полярного дня 24 - 28 июля 2016г. то в StarCalc-е 24 июля Солнце не заходит, а 25 имеет место восход в 01:46 и потом закат в 00:42. Не совсем понимаю как может полярный день оканчиваться восходом (ведь Солнце уже светит), к тому же наступающим после заката, но если эти моменты принять как происходящие в течение одних суток, то долгота дня получается очень короткой (чуть больше часа, - и это сразу по окончании полярного дня. Разве так бывает ?)
В моей же программе, как я догадался, закат происходит уже 24 июля в 00:43, а следующий восход будет в 01:45 на следущий день 25 июля. А затем, следующий закат будет уже в 00:27 следующих суток (то есть 26 июля). Вследствие чего, долгота дня всегда в этот период (сразу после полярного дня) - довольно велика - более 22 часов. Вопрос, собственно, в том, прав ли я ? И следует ли мне исправить представление результатов в программе - если да, то каким образом ?
Извините, если невнятно изложил, но если погонять генератор эфемерид StarCalc-а с координатами Норильска в окрестностях этих дат (около окончания полярного дня) и для сравнения воспользоваться  моей программкой - всё будет понятно и так. Программку (файл .exe) выкладываю . (написана в Visual Studio на языке VB).  Вводим дату и выбираем город (координаты при этом сами установятся) - нажимаем на "Пощщитать"

[gasha]

$ ./Sunset

      SUNSET: solar and lunar rising and setting times 
        (c) 1999 Oliver Montenbruck, Thomas Pfleger     

 First date (yyyy mm dd)               ... 2016 07 20

 Observing site:  east longitude [deg] ... 88.22
                  latitude [deg]       ... 69.33
                  local time - UT [h]  ... 7
 Twilight definition (c, n or a)       ... с

    Date             Moon              Sun             Twilight
                   rise/set          rise/set        beginning/end

 2016/07/20     22:59    03:22    always visible     always bright
 2016/07/21     23:06    04:57    always visible     always bright
 2016/07/22     23:11    06:38    always visible     always bright
 2016/07/23     23:13    08:21    always visible     always bright
 2016/07/24     23:15    10:05    always visible     always bright
 2016/07/25     23:17    11:49    01:46    00:42     always bright
 2016/07/26     23:20    13:33    02:01    00:27     always bright
 2016/07/27     23:25    15:19    02:13    00:15     always bright
 2016/07/28     23:33    17:04    02:22    00:05     always bright
 2016/07/29     23:48    18:46    02:31    23:49     always bright

 all times in local time (=UT+7h)

https://astronomy.ru/forum/index.php/topic,105562.msg3563594.html#msg3563594

У вас там, скорей всего, ошибка в формулах:

(кликните для показа/скрытия)

$ wine ./Календарь.exe
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style                                                                                                                                   
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
fixme:gdiplus:GdipGetLineSpacing ignoring style
System.ArithmeticException: NAN
  at System.Math.Sign (Double value) [0x00000] in <filename unknown>:0
  at Microsoft.VisualBasic.Conversion.Fix (Double Number) [0x00000] in <filename unknown>:0
  at Колендаръ.Form1.Button1_Click (System.Object sender, System.EventArgs e) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Control.OnClick (System.EventArgs e) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Button.OnClick (System.EventArgs e) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.ButtonBase.OnMouseUp (System.Windows.Forms.MouseEventArgs mevent) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Button.OnMouseUp (System.Windows.Forms.MouseEventArgs mevent) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Control.WmLButtonUp (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Control.WndProc (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.ButtonBase.WndProc (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Button.WndProc (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Control+ControlWindowTarget.OnMessage (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.Control+ControlNativeWindow.WndProc (System.Windows.Forms.Message& m) [0x00000] in <filename unknown>:0
  at System.Windows.Forms.NativeWindow.WndProc (IntPtr hWnd, Msg msg, IntPtr wParam, IntPtr lParam) [0x00000] in <filename unknown>:0

P.S. Я понимаю, что тренируетесь/учитесь, но ЗАЧЕМ такую страшную форму делать?))

[βλακασ]

скорей всего, ошибка в формулах

может быть, конечно..... но для других (не полярных) широт - всё прекрасно работает. Да и сами формулы не особо сложные, брал их из книги Ж Меёса "Астрономические формулы для калькуляторов".  Подозреваю, что нужно как-то оговорить условиями момент принятия результирующей переменной значения NaN.
из приведённой Вами ссылки, практически ничего не понял. Я кроме Бейсика, который учил в конце восьмидесятых годов прошлого века, с программированием не знаком. потому и пишу "написал подобие программы". Слишком нескромно, назвать меня программистом, хоть и начинающим; так, скорее хобби. А с другой стороны - точность получилась высокая. Вот под спойлером  часть кода, которая отвечает за рассчёт солнца с комментариями - что к чему:

(кликните для показа/скрытия)

REM рассчёт данных Солнца Начальная эпоха J1900
        ' здесь P - юлианская дата в полдень по гринвичу, JD - юлианская дата в полдень по местному поясному времени (Rad - перевод в радианы, Deg - перевод в градусы )
        J1 = ((P - Lon / 15 / 24) - 2415020.0) / 36525 ' юлианские столетия от начальной эпохи до местного истинного полудня (кульминации) заданной даты (правда пока без уравнения времени)
        J = J1 ' задание момента юлианских столетий для вычисления координат и времени
        SunCalc(J) 'расчет моментов восхода и заката от юлианского полдня (результаты Tr и Ts)
        Tr1 = Tr ' сохранение значения времени восхода для последующего уточнения
        Ts1 = Ts ' сохранение значения времени заката для последующего уточнения
        δk = δ ' сохранение значения склонения на момент кульминации
        ' уточнение восхода (то есть пересчёт момента восхода, исходя из предварительно полученного значения момента)
        J = ((JD - 2415020.0) / 36525) + ((Tr1 - 12) / 876600) ' берем значение юлианских столетий от начальной эпохи до полученного момента местного восхода
        SunCalc(J) 'уточненный расчёт для юлианского восхода
        δr = δ ' сохранение склонения момента восхода
        Trr = Tr ' сохранение значения уточнённого времени восхода под вывод результата (время заката тоже рассчитывается в подпрограмме, но результат здесь игнорируется)
        BR = Fix(Trr) ' и преобразование его в форму часы-минуты с выводом результата в текст-боксе
        NR = Math.Round((Trr - Fix(Trr)) * 60)
        If NR = 60 Then
            BR = BR + 1
            NR = "00"
        End If
        TextBox8.Text = Format(BR, "0") + " ч.  :  " + Format(NR, "00") + " мин."
        ' уточнение заката
        J = ((JD - 2415020.0) / 36525) + ((Ts1 - 12) / 876600) ' берем значение юлианских столетий от начальной эпохи до момента местного заката
        SunCalc(J) 'уточненный расчёт для юлианского заката
        δs = δ ' сохранение склонения момента заката
        Tss = Ts ' сохранение значения уточнённого времени заката под вывод результата (время восхода тоже рассчитывается в подпрограмме, но результат здесь игнорируется)
        BS = Fix(Tss) ' и преобразование его в форму часы-минуты с выводом результата в текст-боксе
        NS = Math.Round((Tss - Fix(Tss)) * 60)
        If NS = 60 Then
            BS = BS + 1
            NS = "00"
        End If
        TextBox10.Text = Format(BS, "0") + " ч.  :  " + Format(NS, "00") + " мин."
        If BS >= 24 Then ' в приполярных широтах так случается иногда почему-то ...... вот в чём вопрос: кто виноват и чё делать ?
            TextBox10.Text = Format((BS - 24), "00") + " ч. : " + Format(NS, "00") + "мин. След. сут."
        End If
        ' долгота дня
        U = Tss - Trr
        HU = Fix(U) ' и преобразование долготы дня в форму часы-минуты с выводом результата в текст-боксе
        MU = Math.Round((U - Fix(U)) * 60)
        If MU = 60 Then
            HU = HU + 1
            MU = "00"
        End If
        TextBox12.Text = Format(HU, "0") + " ч.  :  " + Format(MU, "00") + " мин."
        ' вычисление высоты Солнца и различные условия для полярных дней и ночей
        W = (Math.PI / 2 - Lat + δk) * Deg  'высота солнца при склонении для момента кульминации (истинного полдня) без учёта рефракции, а ведь она значительна у горизонта
        If W > 90 Then W = 180 - W ' высота над горизонтом в южном полушарии измеряется от севера
        TextBox7.Text = Math.Round(W, 1) ' выведение в конечном результате высоты кульминации в текст-боксе


' вышеупомянутое вычисление обращается к подпрограмме SunCalc(J), которая возвращает момент восхода и захода солнца в десятичных долях часа по поясному времени . Аргументом
' служит переменная J - момент времени, на который производится рассчёт, выраженный в юлианских столетиях. Вот её текст:

Private Sub SunCalc(ByVal J As Double) ' подпрограмма вычисляет из конкретного момента времени (в юлианских столетиях) моменты восхода и заката по местному поясному времени в долях часа
        ' На входе (J) - время в юлианских столетиях от начальной даты J1900 до указанного момнета. На выходе (Tr, Ts) - моменты восхода и заката по местному поясному времени
        ' взято из книги Ж. Меёса
        Lat = TextBox3.Text * Rad ' широта места в радианах
        Lon = TextBox4.Text ' долгота места в градусах
        Gr = ComboBox3.Text ' разница в часах между местом и Гринвичем
        L = 279.69668 + (36000.76892 * J) + (0.0003025 * (J ^ 2)) ' средняя долгота солнца в градусах
        M = 358.475833 + (35999.04975 * J) - (0.00015 * (J ^ 2)) - (0.0000033 * (J ^ 3)) ' средняя аномалия солнца в градусах
        M = M Mod 360 ' введение в диапазон от 0 до 360 градусов
        While M < 0
            M = M + 360
        End While
        M = M * Rad ' выражение средней аномалии в радианах
        C = (1.91946 - (0.004789 * J) - (0.000014 * (J ^ 2))) * Math.Sin(M) + (0.020094 - (0.0001 * J)) * Math.Sin(2 * M) + 0.000293 * Math.Sin(3 * M) 'уравнение центра в градусах
        Ω = 259.183275 - (1934.142 * J) + (0.002078 * (J ^ 2)) + (0.0000022 * (J ^ 3)) ' долгота восходящего узла орбиты луны в градусах
        Ω = Ω Mod 360 ' введение в диапазон от 0 до 360 градусов
        While Ω < 0
            Ω = Ω + 360
        End While
        Ω = Ω * Rad ' выражение долготы восходящего узла в радианах
        λ = L + C 'истинная эклиптическая долгота солнца в градусах
        λ = λ - 0.00569 - (0.00479 * Math.Sin(Ω)) ' видимая долгота солнца (введение поправок на аберрацию и нутацию)
        Fa = (153.23 + 22518.7541 * J) * Rad ' поправки на планетные возмущения по Меёсу.
        Fb = (216.57 + 45037.5082 * J) * Rad
        Fc = (312.69 + 32964.3577 * J) * Rad
        Fd = (350.74 + 445267.1142 * J - 0.00144 * (J) ^ 2) * Rad
        Fe = (231.19 + 20.2 * J) * Rad
        F = 0.00134 * Math.Cos(Fa) + 0.00154 * Math.Cos(Fb) + 0.002 * Math.Cos(Fc) + 0.00179 * Math.Sin(Fd) + 0.00178 * Math.Sin(Fe) ' сумма планетных возмущений
        λ = λ + F ' видимая долгота исправленная за планетные возмущения
        λ = λ Mod 360 ' введение в диапазон от 0 до 360 градусов
        While λ < 0
            λ = λ + 360
        End While
        λ = λ * Rad ' эклиптическая долгота в радианах
        ε = 23.452294 - (0.0130125 * J) - (0.00000164 * (J ^ 2)) + (0.000000503 * (J ^ 3)) + (0.00256 * Math.Cos(Ω)) 'наклон эклиптики в градусах
        ε = ε * Rad  'наклон эклиптики в радианах
        Ex = 0.01675104 - 0.0000418 * J - 0.000000126 * (J) ^ 2 ' эксцентриситет орбиты земли
        δ = Math.Asin(Math.Sin(ε) * Math.Sin(λ)) 'склонение в радианах, исходя из широты солнца равной 0. Широта Солнца всегда равна 0, потому что оно лежит в плоскости орбиты Земли.
        H = Math.Acos((-(0.0147) - Math.Sin(Lat) * Math.Sin(δ)) / (Math.Cos(Lat) * Math.Cos(δ))) ' часовой угол для восхода и захода с поправкой на рефракцию в радианах и размер диска светила в радианах
        H = (H * Deg) / 15 ' часовой угол в долях часа
        ' вычисление моментов времени
        Ygr = (Math.Tan(ε / 2)) ^ 2 ' коэффициент для формулы Смарта (уравнения времени)
        η = (Ygr * Math.Sin(2 * L * Rad) - 2 * Ex * Math.Sin(M) + 4 * Ex * Ygr * Math.Sin(M) * Math.Cos(2 * L * Rad) - 0.5 * Math.Sin(4 * L * Rad) * (Ygr) ^ 2 - 1.25 * Math.Sin(2 * M) * (Ex) ^ 2) ' уравнение времени в радианах (формула Смарта) по Меёсу.
        η = (η * Deg) / 15 ' уравнение времени в долях часа. По отношению к принятому в России понятию оно тут с противоположным знаком.
        Tr = 12 - H - η - (Lon / 15) + Gr ' местное поясное время восхода (результат)
        Ts = 12 + H - η - (Lon / 15) + Gr ' местное поясное время заката (результат)
    End Sub

и ещё - а что это у Вас под спойлером ? Некорректности в моем exe-шнике ?  а как Вы их получили ?
PS в любом случае - Спасибо за ответ !

[xd]

Есть граница, разделяющая заходящие и незаходящие объекты. Для незаходящих объектов зенитное расстояние в нижней кульминации менее 90°. Высота/зенитное расстояние в кульминации зависит от широты на Земле и склонения светила, а также от атмосферной рефракции и горизонтального параллакса. Для Луны и Солнца склонение меняется в течение дня. Для приполярных широт, когда продолжительность суток короткая, нельзя брать склонение от балды. Если для расёта захода Вы возьмёте склонение на предыдущий полдень, а для восхода - на следующий, то ошибка будет большой, а для совсем приграничных случаев можно и вовсе получить абсурдные варианты. А вот брать склонение в момент нижней кульминации можно, и если  сделать ещё одну итерацию, уточняющую склонение в момент захода, результат получится очень хороший.

[βλακασ]

Полностью согласен! Только я ничего от балды и не брал. Тем более не задавался моментами из предыдущей или последующей даты. Такие моменты у меня вылезают в результате вычислений, вблизи начала и окончания полярного дня. Я  рассуждал так:
Создал подпрограмму (функцию, назвал её SunCalc) которая вычисляет моменты восхода Tr и заката Ts солнца в зависимости от произвольного момента времени выраженного в юлианских столетиях J.
То есть - задаёшь ей в качестве входных данных определённый момент времени J и на выходе получаешь готовые моменты восхода и заката, выраженные в десятичных долях часа Tr и Ts.
Затем, вычисляю юлианскую дату, на момент Гринвичского полудня (переменная P, выраженная в сутках) той самой даты, которая используется для рассчёта (не предыдущая и не последующая, - а именно та, которая введена в программу). Чтобы определить момент полдня по среднему солнечному времени для заданной местности из момента полдня по Гринвичу вычитаю долготу места, выраженную в долях суток.  Ведь местный средний полдень на восточной долготе наступит на Lon/15/24 суток раньше Гринвича. (Lon - это долгота в градусах). Для западной долготы - долгота со знаком минус.
Отсюда выражение для юлианской даты среднего местного полдня  P - Lon / 15 / 24. Дальше находим переменную J , то есть находим количество юлианских столетих, между местным средним полднем указанной даты и началом эпохи 1900.
J = ((P - Lon / 15 / 24) - 2415020.0) / 36525. Именно этот момент и использую для первого приближения. Момент среднего местного полдня на указанную дату.
Теперь, получив в первом приближении два значения (момент восхода Tr и момент заката Ts), уточняем их. Для этого находим количество юлианских столетий от полученного момента восхода до начала эпохи и повторяем вычисления. Также находим и количество юлианских столетий от полученного момента заката до начала эпохи и также повторяем рассчёты.
Поскольку Tr и Ts выражены в долях часа поясного времени, в первом случае, для восхода (здесь Gr - разница между местным поясным и Гринвичским временем, положит. для вост. долготы)
J = ((P - Gr/24 - 2415020.0) / 36525) + ((Tr - 12) / 876600)   
тут ((P - Gr/24 - 2415020.0) / 36525)  - количество юлианских столетий от полдня местного поясного времени указанной даты до начала эпохи
((Tr - 12) / 876600)  - время от момента полдня  до рассвета, как правило число отрицательное, выраженное в юлианских столетиях (876600 это число часов в столетии)
Аналогично, для заката
J = ((P - Gr/24 - 2415020.0) / 36525) + ((Ts - 12) / 876600)
Таким образом, первое приближение берётся на средний местный полдень, а второе - на непосредственно полученные моменты восхода и заката указанной даты.  по-моему всё правильно....
но если нет - как тогда вычислить количество юлианских столетий от моментов восхода и заката, выраженных в часах по местному поясному времени на указанную дату до начала эпохи J1900. ? 

[βλακασ]

полярный день заканчивается закатом, поэтому сначала поместил колонку заката, а потом восхода. Моя программа выдаёт довольно точные результаты, если не считать того глюка, что закат, приходящийся на начало суток она выдаёт при расчёте предыдущей даты (при вычислениях , например, на дату 26.07.2016 результат для заката выдаёт " в 00:16 следующих суток, то есть 27.07.2016").
И ещё: в старкальке нету двух закатов 28 июля .....получается закатится в 00:05 потом через два часа взойдёт в 02:23, потом не закатываясь ни разу, на следующий день в 02:32 опять восход и следующий закат будет в 23:48  ? чес.слово - что то тут
не так.....

программа стар кальк  (выдаёт сразу в таблице)                   моя программа (таблицу заполняю вручную, последовательно вводя даты)
Норильск, период от 24 по 29 июля 2016 г.                          Норильск, период от 24 по 29 июля 2016 г.

   дата            закат         восход                                                дата                закат            восход
24.07.2016       нету           нету                                              24.07.2016            нету               нету
25.07.2016      00:42          01:46                                             25.07.2016            00:43            01:45     
26.07.2016      00:27          02:01                                             26.07.2016            00:27            02:01           
27.07.2016      00:15          02:13                                             27.07.2016            00:16            02:12     
28.07.2016      00:05          02:23                                             28.07.2016            00:06            02:22      снова закат  23:57
29.07.2016      23:48          02:32                                             29.07.2016            23:49            02:31     
30.07.2016      23:41          02:40                                             30.07.2016            23:41            02:40

[βλακασ]

И всё-таки в приполярных областях бывают дни, когда в течение суток наблюдается восход и два заката !
Полностью переделал алгоритм расчёта восходов и закатов Солнца. Теперь программа с шагом 12 мин. вычисляет высоту Солнца над горизонтом в течение суток и проверяет моменты, когда высота меняет знак (Солнце пересекает горизонт). Момент восхода/заката вычисляется интерполяцией. Для Луны - аналогично. Во вложении - исправленная версия и код программы (в виде текстового документа с комментариями). Точность - на мой взгляд, очень даже ничего....
Буду признателен за отзывы (в т. ч. критические) по существу !  Поскольку нисколько не программист.... 

[1valdis]

восход и два заката !

Если Солнце второй раз заходит, оно должно и второй раз взойти. Что-то здесь нечисто.

[βλακασ]

Если Солнце второй раз заходит, оно должно и второй раз взойти

Оно взойдёт непременно, НО уже в течение СЛЕДУЮЩИХ суток. Как например, упомянуто выше, в Норильске 28 июля 2016г. на момент начала суток (в полночь) солнце будет светить. Первый закат в 0.05 , потом восход в 02.23 и снова закат в 23.57 (за три минуты до окончания суток).
Таким образом каждый раз восход и закат чередуются как надо, но в течение конкретных суток (интервала от 00.00 до 24.00 часов) иногда случается два заката и один восход. (Второй восход будет, но уже на следующие сутки.)
StarCalc при этом показывает первый закат 0ч.05м., не отображая нигде закат в 23ч.57м., а моя программа наоборот указывает момент заката 23ч.57м., и ни слова о том, что будет закат в 0ч.05м. Бесспорно, редчайший случай, но факт.