Как расчитываются...-III (Затмения)

[an1440]

Та-а-а-к. Результаты стали гораздо ближе к описанным в "АнтиМорозинге". В частности,
-584/5/28 - Бозанне. Полное затмение с широкой полосой, наблюдалось на Панамском перешейке, во Франции, средней Италии, Греции и в Малой Азии . - у меня в Риме было 35% - стало 98%.
=================================
Сейчас я просчитываю все новолуния подряд. Но, например, после очередного новолуния ясно, что ожидать следующего следует не ранее, чем через 29 дней. А насчет затмений? Есть ли какой-то период, в течение которого после очередного затмения бесполезно ожидать следующего?
=================================
У меня в программе остались 2 дурацкие рекурсивные процедуры, занимающие 90% времени ее работы. Это вычисление JD по дате и вычисление арккосинуса - он используется при вычислении zенитного расстояния. Видимо, "фирменного" арккосинуса и в паскале нет - я смотрю, Вы его нигде не применяете. Тогда наверно для z должна быть другая формула, позволяющая без него обойтись?

[Анатолий Волчков]

Сейчас уезжаю, вернусь часам у 13, тогда и дам процедуру вычисления JD.

Про ArcCos. Я поступаю иначе: по Cos нахожу Sin и далее ArcTg. Дополнительные вычисления лишь сумма квадратов да корень.

Про периодичность затмений: очевидно, что затмение бывает, когда склонение Луны при соединении с Солнцем близко к склонению Солнца. Так что очередное затмение может произойти не ранее тех же 28-29 суток. Причем в интервале +- половина этого периода около Солнечных затмений могут происходить лунные, так как за пару недель Луна в полнолуние недалеко отклоняется от линии Солнце-Луна.

[Анатолий Волчков]

Типы переменных назначены исходя их конкретных условий работы компилятора.

Procedure JD8 (     IY : LONGINT ;
                      IM : INTEGER ;
              D : DOUBLE  ;
         Var JD : Double    ) ;

VAR
  JN, IYY, IMM, IA, IB, IS : LONGINT ;
  GRDAT : DOUBLE  ;

BEGIN

   IF ( IM > 2 )  THEN
     BEGIN
       IYY := IY;
       IMM := IM + 1;
     END
    ELSE
     BEGIN
       IYY := IY - 1;
       IMM := IM + 13;
     END;

   IF ( IYY = 0 ) THEN IS := 1
     ELSE
       BEGIN
      IS  := TRUNC ( IYY / ABS ( IYY ) ) ;
      IYY := ABS ( IYY );
       END;

   GRDAT := IYY * 10000.0 + IMM * 100.0 + D;

   IF ( IS > 0 )  THEN
     BEGIN
       IB := 0;

       IF ( GRDAT >= 15821015.0 )  THEN
         BEGIN
          IA := TRUNC ( IYY / 100 ) * IS;
        IB := 2 - IA + TRUNC ( ABS ( IA ) / 4 ) * IS;
         END;

       JD := 1720994.5 + D +
           TRUNC ( 365.25   * IYY ) +
           TRUNC (  30.6001 * IMM ) + IB;
     END
   ELSE
       JD := 1720994.5 + D -
           TRUNC ( 365.25   * IYY + 0.75 ) +
           TRUNC (  30.6001 * IMM ) ;

 END;

[an1440]

Дополнительные вычисления лишь сумма квадратов да корень.[/i]

Сумма квадратов чего?
Я понял так, что sin(x)=sqr(1-cos(x)^2) - по Пифагору так вроде получается ?
Правда, так обрабатывается только от 0 до 180°, но имея только косинус, 3-4 квадранты
все равно никак не отличишь от 1-2.

После всех уточнений и исправлений откопал еще одно затмение для Москвы:
22 июля 2009 в течение 3 минут начиная с UT=1:42  3% на высоте 2.5°
Действительно такое планируется?

[Анатолий Волчков]

Про сумму квадратов немного поспешил (имел в виду равенство суммы 1). Утро и сильно спешил.

По косинусу обычно отыскивают угол, величина которого не бывает больше 180 градусов, например, угол между векторами, зенитное расстояние.

В остальных случаях нужно искать дополнительную информацию о четверти.

Затмение в Москве не видимо. Видно в Китае, на Тихом Океане.

[an1440]

Ну наконец-то! Ожил форум.

Портит всю картину рефракция.
Ее приходится учитывать по формуле r=0.016111°*tg(z), а так как эта формула неприменима вблизи горизонта, то при высоте менее 1.58° приходится переходить на постоянную величину 0.583°.
В результате функция имеет излом в этой точке, излом же получается и на траектории видимого движения Солнца и Луны. При расчете затмения в непосредственной близости от горизонта нарушантся монотонность возрастания фазы (провал), а так же возникают ложные зоны видимости, окруженные зоной невидимости...
При б`ольших же высотах расчет получается гораздо лучше, в частности, упомянутое мной 22/07/2009 действительно дает центральную линию через Китай в Тихий океан. Москва же вероятно находится где-то вблизи от границы видимости частной фазы в восходе.
==================================
Процедура JulianDay почему-то пропускает 5-14 сентября 1582 года - то есть JD(15,09,1582)=JD(05,09,1582)=2299130.5
Обратная же процедура пропкскает 5-14 октября, как и приказывали их благородие Григорий XIII.

[Анатолий Волчков]

Для учета рефракции у горизонта (h<15 град) применяют иные формулы. Только вот смысла особого это не имеет, так как в этом случае сильно влияют уклонения реального строения атмосферы от модели, принятой при построении формул. Просто наблюдения организуют так, чтобы не требовалось измерять точные значения высоты светила вблизи горизонта.

Относительно провала JD мне не понятно. Я, правда, не возился с конкретными значениями JD в тот период, но программа, вычисляющая моменты восходов-заходов, не дает никаких скачков в этот период. Если бы JD повторялись через 10 суток, то и моменты тоже должны были бы повторяться, ибо аргумент для вычисления положений светил – JD. Может быть, у Вас переменные назначены не надлежащего размера и происходит переполнение разрядной сетки?

[an1440]

Да нет - не провал JD - JD как раз идут по-порядку. Но по Григорию-XIII не должно существовать в природе дней с 5 по 14 октября 1582 (т.е. после 4 октября было 15).
Обратная процедура так и делает. А прямая - почему-то этот пропуск КАЛЕНДАРНЫХ ДАТ делает не в октябре, а в сентябре.

[an1440]

Проц. JD8 работает нормально, если вместо
GRDAT := IYY * 10000.0 + IMM * 100.0 + D
переписать
GRDAT := IY * 10000.0 + IM * 100.0 + D
Или есть какие возражения?

Опять же никак не пойму чересчур сложную конструкцию:
IF ( IYY = 0 ) THEN IS := 1
    ELSE
      BEGIN
     IS  := TRUNC ( IYY / ABS ( IYY ) ) ;
     IYY := ABS ( IYY );
      END;

По-моему, она равноценна
IS=1
IF IYY<0 THEN IS=-1
IYY = ABS ( IYY )

Или тоже есть какие-то скрытые (от меня пока) соображения?

[Анатолий Волчков]

Процедура сделана по формулам, приведенным в Меесе. Разбираться еще раз в тонкостях операторов вряд ли стоит. Все усложнения имели определенные причины, которые выявлялись при отладке. Повторять процесс отладки не совсем хочется.

В конце концов, Вы вправе изменять программу по своему усмотрению, проверяя правильность результата. Я в работе не замечал ошибок по причине неверности этой процедуры.

[Nazar]

Red Shift 4 и нет проблем.

[an1440]

Red Shift 4 и нет проблем.

Шифт-то он, конечно, ред, но вопрос не в том, чтобы получить готовый результат, а в том, чтобы иметь возможность получить его самостоятельно, и при необходимости вставить куда угодно будет (в программу).
Как, к примеру, перевести часы компъютера на летнее время?
Скажете, нет проблем - ВИНДУС сделает!
А компъютер - 286 и без жесткого диска...
Вот и возникает потребность в подобных процедурах.

[Barnard 33]

Я вот по данным RedShift составил каталог затмений на
www.jenix.chat.ru/astro.htm
Вычисляю затмения. Смотрю с Солнца потом с земли, где живу. Получились такие затмения:
            
Дата   Время   Фаза   Место   Высота
31.05.2003   8-20 до 10-30   0,5   луна выше   
03.10.2005   15-50 до 16-30   0,05   луна ниже   
29.03.2006   16-30 до 18-30   0,85   луна ниже,слева   17-35
19.03.2007   7-10 до 8-20   0,85   луна выше   с восхода
01.08.2008   15-20 до 17-40   0,7   луна выше   
22.07.2009   7-00 до 7-30   0,05   луна ниже   с восхода
15.01.2010   12-10 до 13-20   0,05   луна ниже   

[an1440]

Я вот по данным RedShift составил каталог затмений ...
... где живу.

Что тут сказать? - хорошо живете! Каждый год по затмению! Да еще какие!
Только по моей "самопальной" программе не просчиталось 2-е из перечисленных Вами, а предпоследнее - в настолько малой фазе, что то ли будет - то ли нет.
Кстати, у Вас 2 последних - по 0.05, а у меня 0(чуть больше) и 0.03 - то есть разные...

[Анатолий Волчков]

2 октября 2005 г будет кольцеобразное затмение Солнца. Полоса кольцевого затмения пройдет от Испании через Северную Африку и уйдет в Индийский океан севернее Мадагаскара.

22 июля 2009 г полоса полного затмения пройдет по Китаю и уйдет в Тихий океан. Половину диаметра Солнца закроет во Владивостоке и почти столько же на южном Сахалине.

Что касается расхождений в топоцентрических условиях, так вообще непонятно, как сравнивать: Barnard не указал своего места жительства, нет этих сведений и в данных регистрации. Вот почему я постоянно настаивал на обязательном указании места жительства при регистрации. Это же астрономия, условия наблюдений сильно зависят от места жительства.

[Barnard 33]

Ещё раз проверил по RedShift глядя из Челябинска
Подтверждаю 3 октября 2005 года солнечного затмения не будет - поторопился
22 июля 2009 года Луна тока немного коснётся диска Солнца - затмением назвать трудно.
Ошибка была следующая - по умолчанию установлены размеры Солнца и Луны увеличенными в 3 раза, чтобы легче смотреть на звёздном небе. При вычислениях вернуть забыл - буду пересматривать.  

[Анатолий Волчков]

Только если уж вычислять затмения, так для всей Земли, то есть все затмения, которые видны хоть в одной точке земной поверхности. Смысла публиковать таблицу затмений, видимых в Челябинске, вроде совсем мало. Ведь люди в Мурманске или Владивостоке не смогут воспользоваться этими данными, да и данные для их местности будут неполными. Нужен весь список по Земле.

[Barnard 33]

Хе-хе...
Точно - 22 июля в Москве не видно затмения, а в Челябе (2000 км на восток) Луна минуты на 3-5 касается Солнца в UT=1-20.

[Анатолий Волчков]

Вот примерный список затмений на ближайшие 10 лет. Делал быстро, мог и пропустить кое-что.

Sun Eclipse  2002  4 dec
Sun Eclipse  2003 31 mau
Sun Eclipse  2003 23 nov
Sun Eclipse  2004 19 apr
Sun Eclipse  2004 14 oct
Sun Eclipse  2005  8 apr
Sun Eclipse  2005  3 oct
Sun Eclipse  2006 29 mar
Sun Eclipse  2006 22 sep
Sun Eclipse  2007 19 mar
Sun Eclipse  2007 11 sep
Sun Eclipse  2008  7 feb
Sun Eclipse  2008  1 aug
Sun Eclipse  2009 26 jan
Sun Eclipse  2009 22 jul
Sun Eclipse  2010 15 jan
Sun Eclipse  2010 11 jul
Sun Eclipse  2011  4 jan
Sun Eclipse  2011  1 jun
Sun Eclipse  2011 25 nov

Moon Eclipse  2002 20 nov
Moon Eclipse  2003 16 may
Moon Eclipse  2003  9 nov
Moon Eclipse  2004  4 may
Moon Eclipse  2005 24 apr
Moon Eclipse  2005 17 oct
Moon Eclipse  2006 14 mar
Moon Eclipse  2006  7 sep
Moon Eclipse  2007  3 mar
Moon Eclipse  2007 28 aug
Moon Eclipse  2008 21 feb
Moon Eclipse  2008 16 aug
Moon Eclipse  2009  9 feb
Moon Eclipse  2009  7 jul
Moon Eclipse  2009  5 aug
Moon Eclipse  2009 31 dec
Moon Eclipse  2010 26 jun
Moon Eclipse  2010 21 dec
Moon Eclipse  2011 15 jun
Moon Eclipse  2011 10 dec

[Barnard 33]

Ну вот прогнал RedShift. Конечно не указываю фазу и тип. Ну когда близко будет тогда и посчитаем.
Затмения Солнца на 2002-2010 год для Земли.


[вложение удалено с форума]