Емельянов Эдуард Владимирович (eddy_em) wrote,
Емельянов Эдуард Владимирович
eddy_em

Categories:

Сравнение астрометрических библиотек - 2

Как я уже писал, получается, что разные библиотеки для вычисления apparent place дают разные данные. С расхождением в минуты!
Один из авторов любезно поделился со мной статьей, где таких расхождений нет: там говорится о 5-6mas, а не минутах!!!

Берем все тот же comparelibNOVA-NOVAS-SOFA-ERFA. Вводим туда координаты БТА, сегодняшнюю поправку DUT1 и смотрим, что поменяется, если в первом случае атмосферное давление сделать равным 100кПа, а во втором — 78кПа (обычное давление на горе). Давление как на уровне моря:
                                                          d J2000             d NOVAS          d SOFA Atco13
                  ra deg           ,dec deg          ra''      ,dec''    ra''      ,dec''    ra''      ,dec''
J2000             345.986294056000,28.150389198100
libNOVA equ_prec2 346.228458943796,28.258597189533   871.794,  389.549    -0.009,    0.039   324.298, -266.516
NOVAS precession  346.228461429350,28.258586387865   871.803,  389.510                       324.306, -266.555
SOFA Atci13       346.221443558472,28.258637867687   846.538,  389.695   -25.264,    0.185   299.042, -266.370
SOFA Atco13       346.138376291446,28.332629488577   547.496,  656.065  -324.306,  266.555
SOFA Atco13 NO r  346.221423118462,28.258608573980   846.465,  389.590   -25.338,    0.080   298.969, -266.475
SOFA Atco13 NO eo 345.886031202046,28.332629488577  -360.946,  656.065 -1232.749,  266.555  -908.442,    0.000
ERFA Atco13       346.138376291446,28.332629488577   547.496,  656.065  -324.306,  266.555     0.000,    0.000

Давление как на БТА:
                                                          d J2000             d NOVAS          d SOFA Atco13
                  ra deg           ,dec deg          ra''      ,dec''    ra''      ,dec''    ra''      ,dec''
J2000             345.986294056000,28.150389198100
libNOVA equ_prec2 346.228458943945,28.258597189600   871.794,  389.549    -0.009,    0.039   258.673, -208.067
NOVAS precession  346.228461429500,28.258586387931   871.803,  389.510                       258.682, -208.106
SOFA Atci13       346.221443558383,28.258637867562   846.538,  389.695   -25.264,    0.185   233.418, -207.921
SOFA Atco13       346.156605302540,28.316393568600   613.120,  597.616  -258.682,  208.106
SOFA Atco13 NO r  346.221423120371,28.258608573613   846.465,  389.590   -25.338,    0.080   233.344, -208.026
SOFA Atco13 NO eo 345.904260212808,28.316393568600  -295.322,  597.616 -1167.124,  208.106  -908.442,    0.000
ERFA Atco13       346.156605302540,28.316393568600   613.120,  597.616  -258.682,  208.106     0.000,    0.000

Запустил однострочником, чтобы как можно меньше было расхождение по JD.
В принципе, судя по коду, сравнение делается некорректно. В libnova вызывается функция ln_get_equ_prec2, никакого отношения к apparent place не имеющая, да и julian date не вычисляется нормально. В NOVAS julian date вычисляется, исходя из поправок для TT/TAI, не учитывая DUT1, но все равно идет просто пересчет координат на новую эпоху, без apparent place.

В общем, сравнивать их надо корректно. Но, судя по тому, что NOVA не использует DUT1, ее результаты таки должны быть менее правильными. Хотя, понятное дело: даже для БТА ошибка вычисления координат в ±1'' — мелочь по сравнению с ошибками наведения и сопровождения.
Tags: c, астрофизика
Subscribe

  • А как народ книжки покупает?

    Вот хочу почитать относительно свеженького Гонсалеса-Вудса. Читалки формата А4 у меня, к сожалению, нет. Смотрю на бумажную версию: 113 баксов(!!!)…

  • M$ teams…

    Начал с сегодняшнего дня студентам ЮФУ удаленно лекции читать. У них все завязано на различные корпорации зла. И базовая работа - через teams. ОК,…

  • Почему systemd — дерьмо

    Уже давно на эту статейку натыкался, но все забывал в "закладки" добавить. Вот, добавляю: "systemd — отстой". Советую эту статейку почитать…

promo eddy_em september 3, 12:13 8
Buy for 10 tokens
Уже больше полугода занимаюсь разработкой, вот, наконец-то в мастерских взялись за меня и начали выдавать первые детали. Сегодня сделал тестовую сборку (как обычно, местами пришлось "доработать напильником"): Пока прибор без названия (да и как-то не лезет в голову ничего, у меня нет…
  • Post a new comment

    Error

    Anonymous comments are disabled in this journal

    default userpic
  • 0 comments