Об ОЧЕПЯТКАХ в "Астрономия с ПК", Монтенбрук ИМЕННО 1993г. с прогами на ПАСКАЛЕ

[Александр2020]

В "Астрономия с ПК", Монтенбрук ИМЕННО 1993г. с программами на ПАСКАЛЕ обнаружены ОЧЕПЯТКИ.

стр 47.   В шапке программы "MINI_MOON":
ДОЛЖНО БЫТЬ: Т=(JD-2451545)/36525,
а НЕ (JD-2452545)


стр 48.   В выражении для DL:=,
2-я строчка, записан член: "-668*sin(LS)"
ДОЛЖНО БЫТЬ: "-668*sin(L)"

[Павел Васильев]

стр 48.   В выражении для DL:=,
2-я строчка, записан член: "-668*sin(LS)"
ДОЛЖНО БЫТЬ: "-668*sin(L)"

В кодах на паскале за 1993 год и на си за 2002 год одинаково
// Perturbations in longitude and latitude
  dL = +22640*sin(l) - 4586*sin(l-2*D) + 2370*sin(2*D) +  769*sin(2*l)
       -668*sin(ls) - 412*sin(2*F) - 212*sin(2*l-2*D) - 206*sin(l+ls-2*D)
       +192*sin(l+2*D) -...
Обоснуйте, что это должно быть исправлено.

[Александр2020]

1. Сравните с текстом в указанной книге на стр. 121
2. Сравните с текстом (с переменными, НО НЕ ЧИСЛОВЫМИ КОЭФФ.!!!) в книге Ж.Меёс "Астроном ф-лы для калькуляторов" на стр.120

[Александр2020]

Ранее Pavel Vassiliev любезно дал ссылку на библиотеку программ на Паскале к книге "Монтенбрук, Астрономия с ПК" (1993, ПАСКАЛЬ!!!)
(https://gitflic.ru/project/geoblock/apclib/file?file=source%2Fapc&branch=main)

В книге и в библиотеке обнаружено, помимо вышеописанных очепяток, следующее:

---------------------------------

В книге используется REAL переменные.
В библиотеке они же -- DOUBLE.

(к слову - лично я всё переделываю в EXTENDED).

---------------------------------

стр.22
4-я строка снизу есть оператор "B:=0;",
а в unit Apc.Time он отсутствует.

---------------------------------

стр.141
процедура GMS(),
она приведена на стр. 21,
а в библиотеке Apc.Mathem она фигурирует как DMS().

Аналогично: процедура GGG() приведена на стр. 20
в Apc.Mathem она фигурирует как DDD().

---------------------------------

стр.142
10-я строка сверху.
Коэффициент в операторе д.б
"MODJD:=T*36525+54544.5;"
а не:
"MODJD:=T*36525+51544.5;"

---------------------------------

пуск программы про Луну (стр. 140 - 142) осуществить не удалось.
Возможно из-за вышеуказанных очепяток. Однако сложилось впечатление, что затруднения - из-за ошибок реализации п/п, типа GMS() (см. выше).

Т.к. передо мной стояла задача практического плана (Луна), а не отладка программ - плюнул и написал свой вариант запуска.
Счёт по программе MOONEQU() (стр. 133) в режиме ручного задания даты показал ПОЛНОЕ совпадение с тестовой таблицей на стр. 143.

ОДНАКО !!!
==========

Оказалось, что, например, на 07-10-2023 00:00 UT имеем следующие значения из 3-источников:

                                                    прям.восхождение     склонение
                                                                  h  m  s            °  '  "
1. Результат по программе:           7:22:08.6     +27:23:15.0
2. По АК Козловского на 2023г.:   7:23:01.8     +26:44:06
3. по AstroPixels.com :                 07:22:11.30   +27:23:10.2  (Geocentric !!! Ephemeris for Moon. http://astropixels.com/ephemeris/moon/moon2023.html)


Т.о. расхождение с АК Козловского порядка полу градуса. Причём, я, ранее, уже замечал такого рода расхождение между Козловским и AstroPixels.com, но не мог сделать вывод т.к. не было 3-го источника.

Т.к. расхождение моего счёта по книге и AstroPixels.com минимальное (тем более, там Geocentric Ephemeris),
теперь делаю вывод, что Козловский считает "приблизительно" и это надо учитывать.

Т.о., точность MOON() - я принимаю (гарантировано) в 5-10 секунд (НО НЕ в 1 секунду, как указано в её шапке) и использую её результаты, как истинные.

...................

Т.к. меня интересовала только программа счёта Луны - прочие программы не смотрел и не тестировал.

[gasha]

3. по AstroPixels.com :                 07:22:11.30   +27:23:10.2  (Geocentric !!! Ephemeris for Moon. http://astropixels.com/ephemeris/moon/moon2023.html)

Сайт ИПА РАН даёт такой же результат:

Time scale            : UTC/GMT
Object                : Moon
Relative to           : Observer
Planetary theory      : DE405 (E.M. Standish, JPL Planetary and Lunar Ephemerides, DE405/LE405, JPL Interoffice Memorandum 312.F-98-048, 1998)
Coordinates           : Equatorial
Equator and equinox   : True of the date
Precession model      : IAU2000/2006
Nutation model        : IAU2000/2006
EOP                   : Not applied
Observer              : Geocentric
Light time correction : yes
Relativistic effects  : yes

          Date               Alpha            Delta       Distance(km)
YYYY MM DD HH MM SS.sss HH MM SS.ssssss  DD MM SS.ssssss
2023 10  7  0  0  0.000  7 22 11.296178  27 23 10.185180     397906.933

У Монтенбрука урезанная и устаревшая модель ELP-82. В АК тоже многое завязано на Монтенбруке.
По разнице указанных вами координат ошибки даже не в учёте/не учёте прецессии/нутации. А исходно неверные геоцентрические координаты.
Для восхода/захода Луны подойдёт. Но не более.