Efemérides Eclipse de sol 4 de enero de 2011

Hace unos días ya anuncié que mañana 4 de enero, sería visible un eclipse parcial de sol desde España. Existen diferentes tablas con las efemérides para distintas ciudades, pero muchas veces nos interesa conocerlas para otras posiciones. A continuación he escrito un programa para calcular los instantes de inicio, máximo y fin del eclipse, así como la magnitud del mismo. Aprovechando las posibilidades que ofrecen los Sistemas de Información Geográfica, he integrado los cálculos tradicionales con la obtención de la posición a partir de un mapa.

No hace falta introducir a mano las posiciones. Para utilizarlo hay que elegir el lugar sobre el mapa, haciendo zoom sobre nuestra posición para acercarnos y pulsar sobre “Calcular efemérides para este punto”. Se mostrará la latitud y longitud del lugar, y a partir de ellos, la altitud y las efemérides del eclipse. Todas las horas están expresadas en Tiempo Universal. Para la hora civil hay que sumar 1 hora en invierno.

Para los cálculos he utilizado los siguientes Elementos Besselianos, tomados del Five Millenium Catalog of Solar Eclipses

Polynomial Besselian Elements for:   2011 Jan 04    9.000 TDT  (=t0)

  n        x          y         d          l1         l2          μ

  0  -0.1406310  1.0558220 -22.7412205  0.5635920  0.0173670 313.811188
  1   0.5162760  0.1051387  0.0040620  0.0001098  0.0001092  14.996630
  2  -0.0000418  0.0001064  0.0000060 -0.0000108 -0.0000107   0.000000
  3  -0.0000065 -0.0000015  0.0000000  0.0000000  0.0000000   0.000000

                tan f1 = 0.0047557        tan f2 = 0.0047320

El algoritmo para la determinación de los contactos lo he adaptado de Explanatory Supplement to the Astronomical Almanac.