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.





