Compendio de Trabajos de Investigación CNDG – Biblioteca

Instituto Geofísico del Perú. V. 7 (2006) p. 165 - 174

CÁLCULO DE LAS COORDENADAS HELIOGRÁFICAS MEDIANTE UN

PROGRAMA EN FORTRAN

ISABEL LIZ CASTRO MERINO

Facultad de Ciencias Naturales y Matemáticas
Universidad Nacional del Callao
[email protected]

Prácticas dirigidas por: Dr. Mutsumi Ishitsuka Komaki

Observatorio de Ancón – Área de Astronomía

INTRODUCCION diámetro de abertura d = 150mm, y posee

Hoy en día, son muchos los astrónomos una montura ecuatorial y un motor de
amateurs y profesionales que observan el seguimiento (Lopez, 2005; Quispe, 2005).
Sol, mediante un telescopio. Por modesto El procedimiento observacional para
que este sea el observador que se inicia en determinar la posición de una mancha
este campo de investigación, utiliza como solar es: primero se determina la línea de
mínimo un refractor de 60 mm de declinación Este-Oeste terrestre. Para ello,
diámetro mejor si se tiene una montura se proyecta la imagen del Sol sobre una
ecuatorial y un motor de seguimiento. pantalla, luego señala la posición de la
Como el observador conoce el manejo del mancha en un extremo del campo de
telescopio y los fundamentos teóricos del visión y deteniendo el seguimiento del
estudio del Sol, los métodos de telescopio, esperar que por el movimiento
observación que utiliza son: Proyección de de rotación de la Tierra la sitúe al otro
Imagen, la utilización de filtros especiales extremo del campo. Una paralela a esta
situados delante del objetivo, o mediante recta que pase por el centro del Sol
telescopios reflectores especiales (sin determina la posición E-W y su
aluminar). perpendicular la posición del polo Norte
El área de Astronomía del Observatorio de terrestre.
Ancón, en convenio con la Universidad
San Luis Gonzaga de Ica, realiza las COORDENADAS HELIOGRÁFICAS
observaciones de manchas solares Para fijar la posición de un punto sobre la
utilizando el método de Proyección de superficie terrestre utilizamos dos
Imagen, y son realizados por estudiantes coordenadas: la latitud y la longitud. La
(observadores del sol), los mismos que latitud se mide desde el ecuador hacia los
vienen adquiriendo experiencia. El polos y desde 0º hasta 90º. Para la
instrumento que utilizan es un Telescopio longitud, se fija arbitrariamente un
Refractor de distancia focal f = 1050mm y meridiano origen (Greenwich) y, a partir
I. Castro

de él, se mide hacia el Este y el Oeste, Cefeo situado aproximadamente en 23h

desde 0º hasta 180º. Análogamente, en el +77º.1
Sol se definen una longitud y una latitud
heliográfica. La combinación de ambos ángulos
provoca, a lo largo del año, una desviación
La latitud se determina igual que en del eje solar respecto a la dirección Norte-
nuestro planeta: desde el ecuador solar Sur y una inclinación del ecuador respecto
hacia los polos y desde 0º hasta ±90º. Para a la visual (Figura 1).
definir la longitud se hace preciso fijar un
meridiano origen, pero como no hay
ninguno especialmente privilegiado, se
elige, por convenio, el meridiano
correspondiente al nodo ascendente del
ecuador solar sobre la eclíptica el 1 de
enero de 1854 a las 12 h. La longitud se
mide desde 360º hasta 0º hacia el Oeste, es
decir, en el sentido de la rotación solar.
Así, conociendo la latitud y la longitud
heliográficas podemos determinar la
Figura 1. Ángulos P, B0, L0.2
posición de cualquier punto sobre la
superficie solar. Este aspecto tan distinto que presenta el
Sol influye considerablemente a la hora de
Sin embargo, cabe resaltar que en el Sol calcular las coordenadas de las manchas
no hay estructuras estáticas y que por solares, pues lo que se mide en principio
tanto, lo normal es que las coordenadas de es la posición de esa mancha sobre el disco
una mancha solar vayan variando con el aparente. Para poder realizar
tiempo. Si los ejes de rotación del Sol y la posteriormente las debidas correcciones,
Tierra fuesen paralelos entre sí y se definen dos ángulos que caracterizan a
perpendiculares a la eclíptica, el eje de la posición del sistema de coordenadas
rotación solar siempre se vería en sobre el disco solar. El primero de estos
dirección Norte-Sur y el ecuador solar ángulos se le asigna la letra P y es el
sería una línea recta que pasaría por el
centro del disco aparente. El eje terrestre
1
está inclinado 23º26’ y apunta hacia la Parahelio (Página de Observación Solar)
http://www.astrocantabria.org/parhelio/doccoord.ht
estrella Polar, mientras que el eje de ml
2
Página de la Asociación Larense de Astronomía,
rotación solar está inclinado 7º15’ y se ALDA
http://www.tayabeixo.org/efemerides/efem_fis_sol.h
dirige hacia un punto de la constelación de tm

166
 Cálculo de las coordenadas heliográficas mediante un programa en Fortran

ángulo de posición del extremo Norte del lo largo del año desde –23º.32 hasta
eje de rotación medido desde el punto +23º.32.
Norte del disco, positivo hacia el Este y En la Figura 2 aparece representado P con
negativo hacia el Oeste. P puede variar a sus máximos valores:

-23°.32
- 23°.32

Figura 2. Representación de P.

El segundo ángulo se designa con la letra Así, por ejemplo, B0 = –5º significa que el
B0 y es la latitud heliográfica del centro del ecuador pasa 5º al Norte del centro del
disco. B0 puede variar desde –7º.25 hasta disco (Figura 3).3
+7º.25 y da idea de la inclinación del
ecuador respecto a la visual (Choque, Para determinar las coordenadas de una
2003). mancha, todavía se necesita un tercer
ángulo: la longitud del meridiano central,
Se utiliza para la proyección la plantilla que se denota por L0. Puede variar desde
correspondiente al valor de "B0" según las 0º hasta 360º y, tal como se define la
efemérides del día de la observación para longitud heliográfica, va disminuyendo día
ello se utilizaran el Programa de cálculo de a día a medida que el Sol gira.
coordenadas heliográficas diseñado para Para obtener ángulos P, B0, L0, mediante el
este fin. programa del cálculo de coordenadas
La posición de la plantilla será aquella heliográficas, es necesario calcular
cuya " B0= " (según efemérides) previamente el día juliano.
corresponda a la indicada en la parte
superior derecha.
3
Métodos de Observación del Sol – Página de
Heliofísica
http://idd008cq.eresmas.net/Pagina12.htm

167
I. Castro

Figura 3. Plantilla para B0 = 5 y -5

DIA JULIANO Los días julianos cubren un ciclo de 7980

Para fechar fenómenos astronómicos o años contados a partir del mediodía del 1
históricos lejanos es difícil considerar los de Enero de 4713 aC y es un calendario en
cambios que ha habido en el calendario. que solo se cuentan los días en forma
Como el día Juliano es una cuenta consecutiva, sin meses ni años; y que
sucesiva de días, los astrónomos lo utilizan terminará el 31 de Diciembre de 3267.
para determinar períodos de tiempo entre Este ciclo sale de multiplicar tres ciclos
dos fechas. Así, un cálculo inicialmente menores: uno de 28 años denominado
complejo, se resuelve con una sencilla solar, otro de 19 años que incorpora las
resta de números. fases lunares y uno de 15 años
El inicio del día Juliano se determinó denominado de indicción.
buscando la fecha en la cual coincidía el El día Juliano, a diferencia de la fecha
primer día de cada año. En 1582 Joseph civil no cambia a la medianoche, sino al
Justus Scaliger encontró que ese día era el mediodía, siendo sus horas, minutos y
1 de enero del año 4.713 aC, y lo llamo segundos de igual duración. De aquí que
día juliano en honor a su padre, Julius para la medianoche de un día contado
Cesar Scaliger (a diferencia del calendario como fecha civil, ha transcurrido medio
juliano inventado por el emperador día en fecha Juliana.
romano Julio Cesar).

168
 Cálculo de las coordenadas heliográficas mediante un programa en Fortran

Ejemplo: El día 27 de febrero del año 2003 determinar la confiabilidad de la

45
es el inicio del día Juliano 2452697.5. información obtenida.

RUTINA PRINCIPAL DEL

PROCESAMIENTO DEL MENU PRINCIPAL

PROGRAMA:
El procesamiento se realiza mediante el ELECCION (OPC)

ingreso de los datos y el formato utilizado

es numérico: Día (DD), Mes (MM), Año Valor 0 Valor 1
OPC SUBRUTINA
JULIANO
(YY), luego se ingresa la Hora (XE1),
Minuto (XE2). Valor 2

En primer lugar, se convierte los SUBRUTINA PBL

parámetros de fecha y hora y luego se

procede a resolverlos generando la
END
información del día juliano (jconver),
Figura 4. Estructura del Programa Principal de
después de este cálculo y se procede a Cálculo de Coordenadas Heliográficas (Menú)
calcular P, Bo, Lo y la oblicuidad.
El código contempla, además una función
alternativa. Calculo solo del día juliano.
La Estructura y la visualización del
programa principal y de las subrutinas de
coordenadas heliográficas y día juliano
son las siguientes: (Figuras 4, 5 y 6)

CONTROL DE CALIDAD DE LOS

Figura 5. Visualización del Programa de Cálculo
DATOS DEL PROGRAMA de Coordenadas Heliográficas (Menú)
Los datos obtenidos por el programa
pueden contener errores uno puede ser por
El análisis que se realizo para el control de
el método de procesamiento. Es por ello,
calidad de los datos fue comparar los datos
que cualquier dato obtenido requiere pasar
obtenidos con el programa MICA y el
por un control de calidad, a fin de
programa realizado (Figura 7). Se
considero que el programa MICA seria el
4
El calendario y el día Juliano, Página del más apropiado para la verificación.
Instituto Nacional de Astrofísica, Óptica y
Electrónica.
http://www.inaoep.mx/~rincon/juliano.html

5
Día Juliano, Observatorio Astronómico “ Ing.
Ángel Di Palma”,
http://oamm.50megs.com/maweb6.html

169
I. Castro

SUBRUTINA SUBRUTINA JULIANO

PBL

FECHA Y HORA
JULIANO

FASE DE EJECUCION
CONVERSION
PARÁMETROS
JCONVER
Día juliano Día juliano Día juliano
FASE DE EJECUCION
Po, Bo, OBLICUIDAD

FASE DE EJECUCION RETORNAR

Lo

Po, Bo, Lo, Po, Bo, Lo, Po, Bo, Lo,

Oblicuidad, Oblicuidad, día juliano Oblicuidad,
día juliano día juliano

RETORNAR

Figura 6. Estructura de la Subrutina PBL y Juliano

La base de datos analizados corresponde RESULTADOS DEL PROGRAMA

desde el 01 Enero de 1900 al 01 de Enero A continuación se exponen los resultados
2050 a horas 05:00 UT (intervalo tomado obtenidos al ejecutar el programa con los
cada 10 años). Estos cálculos se realizaron parámetros ingresados.
para el día Juliano. Mientras que para los Todos los datos son almacenados
parámetros P, B0 y L0, los datos obtenidos temporalmente en memoria de la PC en un
son del 01 de Enero de 1990 hasta el 2005 archivo de formato texto llamado
(intervalo de 1 año) a horas 05:00UT. “ResulCH.txt”, también estos resultados
Comparando los parámetros P, B0 y L0, son impresos por intermedio de una
observamos que el error es: -0.01° a 0.01° impresora instalada en el puerto LPT1 y
para B0 y L0, mientras que para P el error son visualizados también en pantalla. Los
es 0.0° datos de salida corresponden a grados,
minutos tanto en la latitud como de la
longitud (Figuras 8 y 9).

170
 Cálculo de las coordenadas heliográficas mediante un programa en Fortran

Comparacion entre el Programa MICA

y el Programa de Ancón

2470000

2460000

2450000
Dia Juliano

2440000

2430000

2420000

2410000
1890 1905 1920 1935 1950 1965 1980 1995 2010 2025 2040 2055

Años MICA
JulAncón

Comparación entre el Programa MICA

y el Programa de Ancón para P

3
P

-3

-6

-9

1988 1990 1992 1994 1996 1998 2000 2002 2004 2006
Tiempo P1 MICA
P2 Ancón

Figura 7. Comparación de los datos obtenidos con el programa MICA y el Programa desarrollado en Ancón.

171
I. Castro

Figura 8. Visualización de los Resultados de la SubrutinaPBL en el Archivo ResulCH.txt.

Figura 9. Visualización del Resultado de la Subrutina juliano, en pantalla.

CONCLUSIÓN programa tanto en la presentación como

En el desarrollo del siguiente estudio, se completar nuevas funciones.
ha llegado a las siguientes conclusiones:
Este programa consiste en calcular algunas AGRADECIMIENTOS
efemérides solares, como el caso de las Mi agradecimiento por el apoyo del
coordenadas heliográficas. Instituto Geofísico del Perú, mediante el
Con el desarrollo de este software uso de sus instalaciones. Al apoyo dado
obtenemos una herramienta muy por el Director de la Estación Telemétrica
importante que con el tiempo se podrá y Satelital de Ancón, Dr. Mutsumi
implementar valiéndonos de su código Ishitsuka, y al MSc. Hugo Trigoso, como
fuente, y se realizaría cambios en el también al Dr. José Ishitsuka del
Observatorio Nacional de Astronomía del

172
 Cálculo de las coordenadas heliográficas mediante un programa en Fortran

Japón por el asesoramiento en la Quispe, A. (2005): Descripción del

elaboración del presente informe. Método de Proyección para el Conteo de
A Dios, a mi familia, como a todo el Manchas Solares en la Estación Solar
personal de Ancón por su amistad y la UNICA-IGP. Compendio de Trabajos de
colaboración que me brindan. Investigación Realizados por Estudiantes
durante el año 2004. Vol.6.
BIBLIOGRAFIA
U.S. Naval Observatory, (1966): The
Choque, E. (2003): Sistema de American Ephemeris and nautical
coordenadas para datos de Imagen solar. Almanac for the year 1968, U.S.
Compendio de trabajos de investigación Government Printing Office, Whasington,
realizados por estudiantes durante el año D.C.
2002. Vol.4, p.109-118.

Garcia de la Rosa, I., (1989): Apuntes de

Física Solar; Instituto Astrofísico de
Canarias.

López, C. (2005): Procesamiento y

Comparación de los Datos de Manchas
Solares Registrados en el Observatorio de
La Universidad Nacional de Ica.
Compendio de trabajos de investigación
realizados por estudiantes durante el año
2004. Vol.6.

173
I. Castro

174