Astronomia basica

Consumo responsable y trabajo etico

2013-11-05T13:18:00.001+01:00

Trabajadora herida en el desplome del edificio Rana Plaza en Bangladesh,
donde se ubicaban los talleres de las firmas internacionales

Por un consumo responsable, apoyemos el trabajo ético

Las víctimas del derrumbe de los talleres textiles en Bangladesh de las empresas multinacionales, producido el 24 de abril del 2013 a las 09:00 am, aún no han cobrado las miserables indemnizaciones, 6 meses después del desastre que acabo con la vida de casi 1200 trabajadores y dejo mas de 2000 heridos muy graves.

Fuente: ONG Solidarity Center, basado en entrevistas a 900 supervivientes y familiares de víctimas / HEBER LONGAS / EL PAÍS (FOTOGRAFÍA: ANDREW BIRAJ (REUTERS).)

Tan solo Primark ha prometido dar (al parecer aún no ha pagado nada, pero se cree que está en ello) una ridícula ayuda, de 140 dolares a los afectados, las demás empresas racanean miseráblemente para evitar pagar las indemnizaciones, a pesar de lo exiguo de su monto.

Los heridos, cerca de 2000 quedaron con secuelas tan graves que en su mayoría no podrán trabajar mas, quedando inválidos de por vida y sin pensión de invalidez ni nada que se le parezca, pues las leyes de ese pais protegen a las empresas internacionales en lugar de proteger a sus ciudadanos.

Asi, con impuestos pequeños y corrupcion alta, salen beneficiadas las empresas (occidentales y locales) y los politicos y funcionarios locales, dejando a los trabajadores sin casi ninguna cobertura social., las protestas se reprimen a palos o a tiros y los muertos se entierran sin demasiado ceremonial.

Edificio Rana Plaza antes de caer

Edificio Rana Plaza después del derrumbe

Los consumidores de las grandes marcas mencionadas no hemos reaccionado y seguimos consumiendo esas prendas de ropa manchadas con la sangre de tantos inocentes, impasibles ante tanto dolor e injusticia. Las autoridades locales han dado muestras de complicidad interesada con esas empresas y no les exigen casi ningún requisito para ofrecer un mínimo de seguridad laboral a sus trabajadores. Las únicas medidas que se están tomando son ser meras acciones cosméticas para evitar el escandalo y la reacción de los consumidores. Para ello en lugar de exigir medidas de seguridad y salarios dignos, se gastan grandes cantidades de dinero en campañas de publicidad en los medios internacionales de modo que no llegue la verdad a los consumidores, que en definitiva son los únicos que pueden obligar a estas empresas a trabajar de forma ética.

La explotación de los obreros por estas empresas es de tal magnitud, salarios de entre 28 y 75 euros al mes y ningún sistema de seguridad social, que parece incomprensible que nadie que tenga un mínimo de decencia pueda contribuir al enriquecimiento de esos piratas modernos disfrazados de honorables empresarios.

Víctimas del accidente han sido olvidadas por las empresas

Estas grandes firmas multinacionales pretenden quedar ajenas a estos crímenes, aparentando no tener responsabilidad en tantas muertes que se suceden tan continuamente. Estas grandes empresas multinacionales son los auténticos responsables de la miseria de muchos países, y sus clientes, es decir nosotros, somos colaboradores necesarios en la creación y en el mantenimiento de esa miseria global, de ese moderno genocidio, al que asistimos indiferentes como si no tuviéramos nada que ver en el asunto. Si al ver las prendas de ropa de estas empresas en sus escaparates y en sus expositores, viéramos impregnadas en ellas las manchas de sangre y las gotas sudor de tantos seres humanos explotados y envilecidos con salarios de hambre, quizás por pudor, o quizás por solidaridad, miraríamos a otro sitio y comprariamos una prenda quizás algo mas cara , quizás no tan de moda, pero confeccionada de forma justa y ética.

La interpretacion de la carta astral 3ª parte. Los significadores universales.

2013-07-23T18:34:00.003+02:00

Interpretacion de la Carta Astral, parte 3: Los significadores universales

Como continuación de lo dicho en las partes 2 y 1 seguiremos básicamente con alguna indicaciones de Adolfo Weiss, relativas de dos típicos errores
que se cometen habituálmente en la interpetración de la cartas astrales. Weiss sostiene, siguiendo en esto a Morín, que un planeta solo puede inducir una determinada influencia si está locálmente determinado para ello. Un planeta por si mismo, o accidentálmente, por su posicion en un signo, ejercería una influencia universal, que no se materializará individuálmente más que a través de la Casa en que se encuentre o rija dicho planeta. Para Weiss esta sería la forma más racional y segura de establecer un pronóstico, evitando que sigamos caminos equivocados nos advierte de ese error tan comúnmente aceptado.

Segun Weiss son erróneas las ideas difundidas por Claudio Ptolomeo en relación a los significadores universales: el ejemplo más característico de este error sería considerar al Sol, de forma general, como significador del padre en los natalicios diurnos y a Saturno en los nocturnos. Este axioma tan difundido, incluso entre astrólogos de renombre, nos llevaría a posiciones absurdas a causa, principálmente, de su excesiva generalidad. Weiss alerta del uso indiscriminado de los significadores y sus analogías citando a Morín: “Como se sabe, hasta ahora, ha sido admitido por todos los astrólogos que el Sol significaría en sí al padre, al conyugue, al soberano, a los nobles, la gloria, el honor, la vida, etc. ; la Luna, a la madre, las soberanas, el pueblo, etc. ; Mercurio, las facultades intelectuales, et., etc.; … Aunque los planetas poseen analogías … con las mencionadas personas o cosas ….no es lícito preferir de antemano una de estas analogías a las demás ...ya que el resultado que le corresponde en cada natividad particular nunca puede resultar de sus analogías, sino tan solo de su determinación local. … tenerlo simultáneamente (hablando del Sol) por el significador del padre,el conyugue, la vida, los honores etc., se llega a pronósticos idénticos a todas estas categorías, y por tanto a los resultados mas necios que pueda imaginarse. El mismo Cardano se ve obligado a confesar “ Ptolomeo ha introducido en la astrología otra confusión atribuyendo varios significados a un mismo significador, así ha declarado que la Luna significa el cuerpo, las propiedades morales, la vida, la madre, las criadas, las hermanas. Pués bién, ¿como debe estar condicionada la Luna en la figura natal de un hombre longevo, pero cuya esposa muere durante un parto, a quién la muerte priva de alguna de sus hijas, mientras que las demás le quedan a salvo, cuyas criadas le traicionan, cuya madre muere prematúramente, mientras el mismo goza de buena salud, pero posee una naturaleza cambiadiza y malos instintos.?”

“Si Ptolomeo tuviera razón afirmando que en un nacimiento diurno los juicios sobre el padre han de basarse en el estado cósmico del Sol, y en un nacimiento nocturno, en el de Saturno, le seguiría que el día que estuviera el Sol en el signo de Leo, y además en buén aspecto con Saturno o con Venus, no nacería en toda la Tierra nadie cuyo padre no fuera feliz y longevo, o en caso contrario, con un Sol mal colocado y afligido, que no fuese infeliz y de corta vida. Contémplese otra consecuencia absurda, ¡Todos cuantos hayan nacido durante ese periodo serían influenciados de manera igual en su destino por sus padres!. Tal teoría no es solo contraria a la experiencia sino que anula el significado de las Casas y haría inutil establecerlas. “

Se pueden encontrar otros ejemplos de los absurdos a que se llega con esta teoría de los significadores universales, como es el caso de la natividad del sabio Isaac Newton, en la que se encuentra Mercurio en su mayor debilidad , en exilio, en el signo de Sagitario, de donde se demuestra la inutilidad de tratar de sacar conclusiones respecto de la inteligencia basandose en la idea de Mercurio como significador universal de esa cualidad.

En resumen, a los planetas podemos atribuirles unas analogías generales, consecuencia de su naturaleza esencial, esas analogías se concretarán y valorarán según su estado cósmico, es decir a su posición por signo, aspectos, y demás características generales que nos daría la influencia general que podrían llegar a ejercer y, por último, las conclusiones concretas o particulares de cada natividad se podrían extraer concretando o individualizando esas posibles influencias por sus posiciones y regencias en cada Casa.

Tabla de Casas de Regiomontano. Programa en basic

2013-01-07T11:13:00.001+01:00

'regio.bas
'Este programa produce una Tabla de Casas por el sistema de Regiomontano que se graba en un 'fichero con el nombre Regiomon.doc. Este fichero puede abrirse con cualquier procesador de 'textos.

'Las instrucciones utilizadas en el codigo son de Basic y pueden ser facilmente traducibles a 'cualquier lenguaje de programacion.

'Se producen las longitudes de las Casas 10, 11, 12, 1, 2 y 3, las otras no son necesarias pues basta con añadirles 6 signos.

cls
pi=3.14159296535
cdr= pi/180
gecl=23:mecl=27:secl=15
obl=(gecl+mecl/60+secl/3600)*cdr

'print tan(105/cdr)
'input rrr

dim z(12)
dim p(12)
dim v(12)
dim x(12)
dim y(12)
dim l(12)
dim sic(12)
dim sic$(12)
dim gr(12)
dim mi(12)
dim po$(12)
open "regiomon.doc" for output as #1
for si=1 to 12
for grado=0 to 29
g=grado
gr=int(g)
l=(si-1)*30+g                           'l=longitud del medio cielo
ss=si:gosub signos:si$=ss$     'signo del MC en letra
ra=atn(tan(l*cdr)*cos(obl))/cdr 'ascension recta en grados
if ra<0 then ra=ra+360
if ra>360 then ra=ra-360

if l=0   then ra=0.0001
if l=90 then ra=90
if l>90 and l<270 then ra=ra+180
if l=180 then ra=180
if l=270 then ra=270
if l>270 then ra=ra+360
if l>360 then ra=ra-360
if l=360 then ra=360

rag=int(ra):ram=(ra-rag)*60 ' ascension recta en grados y minutos
st=ra*24/360                         ' st=hora sideral

' h.s. en horas,minutos y segundos
sth=int(st)
stmm=(st-sth)*60 :stm=int(stmm)
stss=(stmm-stm)*60:sts=cint(stss)

'for la=22 to 37
la=22
' como ejemplo se crean las Casas para la latitud 22N, basta con sustituir este valor para
' crear la tabla de otra latitud. Tambien puede introducirse un bucle para varias latitudes
' ver la linea anterior donde se crea un bucle de 22 a 37 grados de latitud.

lat=la*cdr

for n=1 to 12 ' bucle de las 12 casas

   z(n)=60+n*30 : if z(n)>360 or z(n)=360 then z(n)=z(n)-360
   p(n)=atn(tan(z(n)*cdr)*cos(lat))/cdr
   sz=sin(z(n))
   cz=cos(z(n))
   if sz>0 and cz<0 then l(n)=l(n)+180        ' cuadrante 2

   if sz<0 and cz<0 then l(n)=l(n)+180        ' cuadrante 3
   if sz<0 and cz>0 then l(n)=l(n)+360        ' cuadrante 4

   if n<7 then then p(n)=p(n)-180
   if p(n)<0 then p(n)=p(n)+360

v(n)=sin(p(n)*cdr) * sin(lat) * sin(obl)

x(n)=cos(ra*cdr+z(n)*cdr) * cos(lat) * cos(obl)- v(n)
y(n)=sin(ra*cdr+z(n)*cdr) * cos(lat)

t=y(n)/(x(n))
l(n)=atn(t)/cdr

if y(n)>0 and x(n)<0 then l(n)=l(n)+180
if y(n)<0 and x(n)<0 then l(n)=l(n)+180
if y(n)<0 and x(n)>0 then l(n)=l(n)+360

if l(n)<0 then l(n)=l(n)+360
if l(n)>360 then l(n)=l(n)-360

sic(n)=int(l(n)/30)+1 'signo de la Casa
ss=sic(n):gosub signos:sic$(n)=ss$
gc=l(n)-30*(sic(n)-1):gr(n)=int(gc)
m=(gc-gr(n))*60:mi(n)=cint(m)

next n

if g=0 then
   print #1, "Tabla de Casas de Regiomontano. Alba de Hermes"
   print #1, " C10 AR   HS"
   print #1, " s g gr mi h m s Lo Casa11 Casa12 Casa 1 Casa 2 Casa 3"
end if
print #1, using "\ \ ## ### ##,.# ## ## ## ## ##\\## ##\\## ##\\## ##\\## ##\\##";_
si$;g;rag;ram;sth;stm;sts;la;_
gr(11);sic$(11);mi(11);gr(12);sic$(12);mi(12);gr(1);_
sic$(1);mi(1);gr(2);sic$(2);mi(2);gr(3);sic$(3);mi(3)
' next la
next grado
next si
close #1
end

signos:
select case ss
case 1: ss$="ar"
case 2: ss$="ta"
case 3: ss$="ge"
case 4: ss$="cn"
case 5: ss$="le"
case 6: ss$="vi"
case 7: ss$="li"
case 8: ss$="es"
case 9: ss$="sa"
case 10:ss$="cp"
case 11:ss$="ac"
case 12:ss$="pi"
end select
return
sub s

while inkey$=""
wend
end sub

Dia de la semana en que naci

2013-01-06T16:04:00.000+01:00

Para saber el dia de la semana que fue el dia que naci, lo primero es ver que dia fue el 1 de enero de cada año en la tabla siguiente:

Bisiestos
En esta columna

1909 Viernes	1910 Sabado	1911 Domingo	1912 Lunes
1913 Miercoles	1914 Jueves	1915 Viernes	1916 Sabado
1917 Lunes	1918 Martes	1919 Miercoles	1920 Jueves
1921 Sabado	1922 Domingo	1923 Lunes	1924 Martes
1925 Jueves	1926 Viernes	1927 Sabado	1928 Domingo
1929 Martes	1930 Miercoles	1931 Jueves	1932 Viernes
1933 Domingo	1934 Lunes	1935 Martes	1936 Miercoles
1937 Viernes	1938 Sabado	1939 Domingo	1940 Lunes
1941 Miercoles	1942 Jueves	1943 Viernes	1944 Sabado
1945 Lunes	1946 Martes	1947 Miercoles	1948 Jueves
1949 Sabado	1950 Domingo	1951 Lunes	1952 Martes
1953 Jueves	1954 Viernes	1955 Sabado	1956 Domingo
1957 Martes	1958 Miercoles	1959 Jueves	1960 Viernes
1961 Domingo	1962 Lunes	1963 Martes	1964 Miercoles
1965 Viernes	1966 Sabado	1967 Domingo	1968 Lunes
1969 Miercoles	1970 Jueves	1971 Viernes	1972 Sabado
1973 Lunes	1974 Martes	1975 Miercoles	1976 Jueves
1977 Sabado	1978 Domingo	1979 Lunes	1980 Martes
1981 Jueves	1982 Viernes	1983 Sabado	1984 Domingo
1985 Martes	1986 Miercoles	1987 Jueves	1988 Viernes
1989 Domingo	1990 Lunes	1991 Martes	1992 Miercoles
1993 Viernes	1994 Sabado	1995 Domingo	1996 Lunes
1997 Miercoles	1998 Jueves	1999 Viernes	2000 Sabado
2001 Lunes	2002 Martes	2003 Miercoles	2004 Jueves
2005 Sabado	2006 Domingo	2007 Lunes	2008 Martes
2009 Jueves	2010 Viernes	2011 Sabado	2012 Domingo
2013 Martes	2014 Miercoles	2015 Jueves	2016 Viernes
2017 Domingo	2018 Lunes	2019 Martes	2020 Miercoles

Por ejemplo si naciste el año 1985, el dia 1 de enero de se año fue Martes

y despues contamos desde el martes tantos dias como se indique en la tabla siguiente segun el dia y el mes de nacimiento (uno mas si es bisiesto en las fechas posteriores al 28 de febrero):

1 8 15 22 29	2 9 16 23 30	3 10 17 24 31	4 11 18 25	5 12 19 26	6 13 20 27	7 14 21 28
0	1	2	3	4	5	6	Enero y Octubre
1	2	3	4	5	6	0	Mayo
2	3	4	5	6	0	1	Agosto
3	4	5	6	0	1	2	Febrero, Marzo y Noviembre
4	5	6	0	1	2	3	Junio
5	6	0	1	2	3	4	Septiembre y Diciembre
6	0	1	2	3	4	5	Abril y Julio

Por ejemplo si nací el 15 de Agosto, ese día será 2dias después del Martes, es decir Jueves.

Nociones elementales sobre los sistemas de Casas en las Cartas Astrales

2013-01-06T12:29:00.003+01:00

Los primeros horóscopos que se conocen no contenían mas que una relación de planetas con su posición en el zodiaco. A partir de un determinado momento empieza a ganar fuerza la idea de que las estrellas que primero aparecían en el horizonte oriental tendrían una importancia decisiva en los acontecimientos o en la vida de la persona para la cual se levantaba ese horóscopo. Así aparece la noción de signo naciente o Ascendente.

En la época romana, durante los primeros siglos de nuestra era, era el signo ascendente lo que se consideraba que representaba al nativo en lugar del signo solar (signo donde se posiciona el Sol) como actualmente, de ahí que pronto se tendió a considerar el resto de la bóveda celeste tal como la vemos desde nuestra posición, es decir la bóveda celeste que se apoya en nuestro horizonte.

En esa bóveda celeste, una vez aceptada su parte oriental o Ascendente como significativa parece lógico que se conceda una cierta importancia a su parte opuesta, es decir la parte de descenso de las estrellas, que es el opuesto al Ascendente y las partes que formar 90º con ellos, el Medio Cielo , y su opuesto el Inmun Coeli. Esta idea , que surgiría de forma intuitiva, fue ganando firmeza a medida que se observaran los resultados en su aplicación hasta el punto que en la antigüedad llega a ser aceptado unánimemente por todos los astrólogos. Una vez que se acepto el concepto de dividir la bóveda celeste y el horizonte en cuatro partes ya estaba abierto el camino a una mayor sofisticación o refinamiento de la idea y por analogía con la división del circulo zodiacal en 12 signos se dividió la bóveda celeste en 12 partes que se denominaron Casas, Domicilios, Mansiones, etc. Otras divisiones en numero de 8, 16, etc, se sabe que se había experimentado a lo largo de la historia, por ejemplo en un manuscrito esenio de las cuevas de Qumram hay horóscopos con división en 16 casas, Casas de "luz" y de "oscuridad".

El Horizonte: La idea de considerar el horizonte de cada lugar y las variaciones que en la bóveda celeste se producen a partir de ese horizonte lleva a un punto de gran interés para la astrología horoscopica, pues consigue un sistema de representar el cielo de una forma muy individualizada no solo a partir del momento sino también del lugar de nacimiento. De esta forma se consigue que dos nacimientos simultáneos pero en lugares diferentes tengan horóscopos diferentes y que reflejen con mas exactitud el cielo tan como es observado en ese momento y en ese lugar. Los antiguos que idearon este sistema de representación no conocían, debido al desarrollo de la astronomía, que la Tierra tiene dos movimientos simultáneos: el de traslación (1 año) y el de rotación (1 día), pero intuitivamente llegaron a una representación que ofrece una fusión de ambos movimientos.

El zodiaco con sus 12 signos es un marco de representación respecto del movimiento de traslación, mientras que el circulo de Casas al estar en relación con la variación del horizonte del lugar representa el movimiento de rotación de la tierra respeto de su eje.

Importancia del circulo de Casas: Si bien la representación de primer circulo, el circulo zodiacal, no ha sufrido variaciones desde su formulación, el segundo circulo ha sido objeto de numerosas conjeturas y perfeccionamientos. Esto ha hecho que surjan diversas formas o sistemas técnicos de realizar esa división, fruto de la reflexión de investigadores que han querido ver en sus intuiciones e investigaciones la solución definitiva a un problema sobre el cual aun no se ha dicho la ultima palabra.

Desde el principio al circulo de Casas se le ha asignado un papel mucho mas importante en el entramado de la interpretación del horóscopo. que al primer circulo, el circulo zodiacal. Al circulo de los signos se atribuye la importancia de modificar, por la posición del planeta en tal o cual signo, la naturaleza esencial del planeta y deducir de ello lo que llamamos su estado cósmico. Sin embargo al segundo circulo se la ha asignado la difícil tarea de describir la relación del nativo con todo el conjunto de circunstancias que acompañan su vida y que modifican su carácter y destino. Por la propia naturaleza amplia y subjetiva de esos conceptos no podía menos que esperarse que surgiera la controversia y el desacuerdo.

Principales sistemas de división de horóscopo: Podemos clasificar los principales sistemas de división de Casas, viendo cada uno de ellos en relación a cuando surge y en la medida de lo posible porque surge, esto seria una clasificacion histórica, pero resulta mas interesante la clasificacion que hace el astrólogo W. Holden de acuerdo con sus características astronómicas. A los tres grupos que hace este astrólogo hay que añadir el método cromático de Demetrio Santos. Este ultimo es un cuarto grupo pues no es posible clasificarlo con los otros pues se aparta radicalmente de los conceptos que han inspirado todos los otros sistemas.

De todo el conjunto de sistemas, muy extenso y en crecimiento, tan solo cuatro sistemas han ganado suficiente firmeza para que fueran considerados, cada cual en su época, el estándar aceptado de forma mas o menos universal, estos son el sistema de Casas-signo, en la antigüedad, el sistema de Alcabitus, durante toda la edad media, el sistema de Regiomontano en la época moderna y el sistema de Placidus hasta hoy. El resto de sistemas son aportaciones mas modernas a la teoría astrológica que son objeto de investigaciones actuales.

Este tema no esta acabado, es posible que tengan que verse todavía importantes novedades, como la aportada por el astrólogo español Santos, que abandona los sistemas geométricos para conectar con los nuevos conceptos de la astrofísica y la astrobiologia. Resumiendo lo conocido hasta hoy, los sistemas de división de Casas (Domificacion) los podemos clasificar en esta forma:

Métodos eclípticos
Métodos espaciales
Métodos horarios y
otros Métodos

Conceptos astronómicos: Para no hacer muy difícil la lectura de lo anterior haré un recordatorio de algunos conceptos necesarios para seguir con la exposición del tema:

COORDENADAS: la localización de un punto en una esfera se realiza con dos cantidades o coordenadas. Las coordenadas son siempre dos ángulos para determinar una posición. Son siempre dos, una sobre un plano de referencia y la otra en un plano perpendicular. Si se trata de la esfera terrestre, la tierra, se utilizan las coordenadas geográficas; si se trata de la esfera celeste se utilizan las coordenadas astronómicas.

PLANO DE REFERENCIA: es el plano donde se miden las coordenadas, se puede denominar PR.

Ángulo en el plano de referencia: es la medida del ángulo generalmente en grados (dividiendo un circulo completo entre 360), pero también puede medirse en horas (dividiendo un circulo completo entre 24).

Perpendicular al plano de referencia: es un plano perpendicular al plano de referencia, por ejemplo el plano vertical al horizonte.

POLOS: son los puntos donde corta el eje perpendicular al plano de referencia con la bóveda celeste

Círculos máximos: son círculos imaginarios, que se suponen en el cielo a gran distancia para medir ángulos entre estrellas.

Resumiendo lo dicho antes los sistemas de coordenadas utilizados por los astrónomos son los siguientes:

A) Coordenadas Geográficas: son los ángulos con los que situamos un punto sobre la superficie terrestre.

Ecuatoriales : son las coordenadas que toman como plano de referencia el ecuador.

Longitud, la coordenada sobre el plano de referencia, el ecuador terrestre. No confundir con la longitud astronómica que se mide sobre la eclíptica.

Latitud, : es la coordenada perpendicular al PR

Co-Latitud: 90 grados menos la Latitud

Polo Norte y Polo Sur ecuatorial , son sus polos.

Meridianos: circulos maximos que pasan por los polos.

Paralelos: Círculos sobre la superficie terrestre paralelos al Ecuador.

B) Coordenadas Astronómicas: son los ángulos con los que situamos un punto del cielo.

Horizontales: son las coordenadas que toman como plano de referencia el horizonte.

Azimut , la coordenada sobre el PR, el horizonte del observador.

Altitud: es la coordenada perpendicular al PR

Distancia zenital: es 90 grados menos la Altitud.

Zenit y Nadir , son sus polos

Primer Vertical

Meridianos: círculos formados por planos paralelos al horizonte

Ecuatoriales Horarias: son las coordenadas horarias que tomas como ángulo sobre el ecuador el ángulo horario, es decir la hora

Ascensión Recta, Ángulo Horario , es la coordenada sobre el PR

Declinación, es la coordenada perpendicular al PR.

Distancia polar: 90 grados menos la Declinacion.

Norte y Sur ecuatorial , son sus polos

Eclípticas : son las coordenadas que toman como plano de referencia el plano de la Eclíptica.

Longitud: es la coordenada sobre el PR, la Eclíptica.

Latitud, es la coordenada perpendicular al PR.

Elevación polar: 90 grados menos la latitud

Norte y Sur eclíptico: son sus polos.

C) y algunos otros conceptos:

Primer vertical : Circulo máximo perpendicular al plano horizontal, con centro en el punto de observación y orientado en dirección Este-Oeste

Ascensión oblicua : Ángulo o movimiento aparente de las estrellas en Círculos paralelos al ecuador.

Domificacion : División en Casas de la bóveda celeste.

DESCARGAR GRATIS PROGRAMAS DE CARTAS ASTRALES

2013-01-02T16:25:00.002+01:00

DESCARGAR GRATIS PROGRAMAS PARA LEVANTAR CARTAS ASTRALES

Existen en la red numerosos programas informáticos que puedes bajar a tu ordenador, gratuítamente, y crear cartas astrales o realizar cálculos de forma automática.
De entre todos ellos solo unos pocos son aconsejables, pués muchos son versiones de demostración, "demos", solo válidas por un námero de días determinado, otros son versiones reducidas del programa completo, que no son muy útiles. Al final solo quedan unos pocos que son auténticamente prácticos. De ellos recomendamos los siguientes:

Astrolog: Hay muchas versiones, la mas completa, con un atlas, amplia base de datos de natividades de famosos (desgraciádamente solo del ámbito anglosajón) y ficheros de texto de interpretaciones que pueden traducirse o modificarse.

El programa principal es este:

Astrolog32

Puede completarse con dos ficheros de datos que forman el atlas internacional y el atlas de Estados Unidos. También es aconsejable incluir el fichero con datos de natividades de personas famosas y las efemérides para calcular horóscopos antiguos.

En principio este programa se desarrolló para 16-bit por Walter Pullen, más tarde fué mejorado por Alois Treindl y Valentín Abramov. Antonio Ramos lo reprogramó para 32-bit y añadió las mejoras de la versión 32, el desarrollo original de este programa puede descargarse desde aquí Astrolog.
La versión original de Walter Pullen también puede descargarse desde la página personal del Sr. Pullen http://www.astrolog.org/
Desde aquí se puede obtener la version 5.4 que fué la última que se desarrolló por su creador. En esta página también se dan interesantes enlaces a otras páginas de astrología y artículos sobre sus concepciones espirituales.

No obstante, es mas completa la versión astrolog32 de Antonio Ramos, indicada antes.
Esta version se puede descargar diréctamente desde el sitio de su autor el Sr. Ramos http://www.southville.net/astrolog32/en/ . En este sitio se ofrece la posibilidad de enlace con otra versión más moderna, la versión 2 de este programa realizada por el Sr. Tomas Kubec, que esta disponible en varios idiomas (inglés, francés, italiano y checo) , pero no en español lamentáblemente.

El enlace que lleva diréctamente a la versión 2 es este . Aqui puede descargarse el programa y sus complementos: atlas, etc. Su autor, el Sr. Tomas Kubec, ha sustituido este programa por este otro StarFicher, que puede descargarse desde esta pagina. De los dos programas me sigue gustando más el anterior, es decir la versión 2 de astrolog, pues aunque StarFicher cambia totálmente la filosofía y entorno gráfico del programa, en mi opinión no lo mejora. No obstante es posible que otro usuario sepa encontrar algunas virtudes diferentes y apropiadas a ciertos casos. No debe menospreciarse este último programa con solo esta primera opinión ya que este programa está evolucionando e incorportando nuevas opciones y posibilidades que lo van mejorando en sus nuevas versiones.

Zet: Este programa puede ser el principal competidor de Astrolog, pués le supera en algunos aspectos: en primer lugar presenta también funciones de astronomía muy interesantes y tiene la posibilidad de ponerse en español (en versiones antiguas, no en las últimas). Esta posibilidad se refiere a las pantallas del programa, no a las interpretaciones que estan todas en ingles. No obstante como las interpretaciones nos llegan en ficheros de texto, pueden modificarse, personalizarse o traducirse. Puede encontrarse en:

ZET

Halloran: la versión de demostración de este programa es ideal para principiantes, pués tiene pocas, pero las imprescindibles, opciones y puede ponerse en español desde la instalación. Por eso yo siempre la recomiendo a quién me pregunta con que programa puede iniciarse. Puede bajarse desde:

HALLORAN

Astrowin: Esta web creada por Allen Edwall, de orientación rosacruz, ofrece varios programas, el mas interesante es el que tiene el nombre AstroWin, pero hay varios mas que pueden ser útiles para algunas aplicaciones especiales:

ASTROWIN

Whatwatch: Este interesante programa, realizado por el astrólogo alemán Robert Amlung, permite realizar sinastrías y cartas compuestas. También, en su ultima versión aporta una muy buena opción de astro-cartografía. Lo mas característico de este programa es la representación lineal del aspectario de cada carta astral. Puede bajarse diréctamente desde:

http://papertv.com/

Jagannatha Hora: Este es un programa excelente de astrología hindú, prácticamente tiene disponibles todas las funciones de la versión completa por lo que puede utilizarse para profundizar en el estudio de esta rama de la astrología. Se puede bajar diréctamente desde:

http://www.vedicastrologer.org/jh/

Canopus: Este programa ofrecido por la Clairvision School de Australia del Dr. Samuel Sagan tiene una serie de opciones que le hacen especialmente atractivo. Esta escuela ofrece cursos de astrología, publica los libros de Samuel Sagan y puedes descargarte música relajante. Con este programa puedes componer una presentación bastante personal e incorporar tus propios textos. Actualmente no se puede bajar directamente desde la pagina de sus creadores, pero ofrecen un servicio on-line para la confección de cartas astrales.

King of the Stars: Este programa permite realizar cartas astrales según los criterios de la astrología china. Se puede bajar desde:

http://www.delemme.com/etel1.htm

Esta web también ofrece otros programas interesantes, de numerología, Feng Shui, etc

CALCULO DEL DIA DE LA SEMANA DE UNA FECHA

2013-01-02T16:17:00.000+01:00

Para calcular el dia de la semana que corresponde a una fecha
cualquiera, se procede de la forma siguiente:

1º: se calcula el Dia Juliano (DJ) que corresponde a esa fecha:
Veamos un ejemplo: calcular el DJ corresponde al 25 de Enero de 1945 a las 0 horas

llamamos D al dia, asi D=25, M al Mes, asi M=1, Y al año, asi Y=1945

Primer paso: vemos si M es mayor o menor de 2:
si M es menor de 3, es decir si M es 1 o 2 (como es en nuestro ejemplo) quitamos 1 a Y y sumamos 12 a M: en nuestro ejemplo tendremos Y=1945 - 1 = 1944, M=1+12=13

Segundo paso: calculamos A y B
A es igual a Y dividido por 100 y quitamos los decimales. En nuestro ejemplo 1944 / 100 = 19,44 asi que A=19

B es igual a 2 - A + int(A/4)
En nuestro ejemplo sera: B = 2 - 19 +int(19 / 4) = 2 - 19 + 4 = -13

int quiere decir que de la operacion que se indica entre parentesis se desprecian los decimales.

Si la fecha fuera antigua (anterior a la reforma gregoriana) y esta dada segun el Calendario Juliano, B = 0

Tercer paso: Con estos datos calculamos el DJ

DJ = int(365,25 x (Y+4716)) + int(30,6001 x (M+1) + D + B - 1524,5

En nuestro caso tendremos:

DJ = int(365,25 x (1944 + 4716)) + int(30,6001 x (13 + 1) + 25 -13 - 1524,5 =

int(365,25 x 6660) + int(30,6001 x 14) +12 =

2.432.565 + 428 + 12 - 1524,5 = 2.433.005 -1524,5 = 2431480.5

2º Ahora, con el DJ calculamos el dia de la semana
sumamos al DJ la cantidad de 1,5, y el resultado se divide por 7. El resto de la division es el dia de la semana:

1 Lunes

2 Martes

3 Miercoles

4 Jueves

5 Viernes

6 Sabado

0 Domingo

Sumamos a DJ + 1,5 = 2431480.5 + 1.5 = 2431482

Dividimos el resultado por 7 : 2431482 / 7 = 347354.571428571

Multiplicamos por 7 la parte entera del resultado anterior: 7 x 347354 = 2431478

Restamos los dos resultados anteriores: 2431482 - 2431478 = 4

Como sale 4 el dia es Jueves

Determinacion de cuspides. Sistema de Morinus. Programa en Basic

2013-01-02T13:57:00.000+01:00

'morinus
cls

pi=3.14159296535
cdr= pi/180

obl=23.45*cdr

dim z(12)
dim p(12)
dim v(12)
dim x(12)
dim y(12)
dim l(12)
dim signo(12)
dim grado1(12)
dim minuto(12)

lat=43*cdr
st=226.02
for n=1 to 12
z(n)=st+60 + n*30
if z(n)>360 then z(n)=z(n)-360
l(n)=(atn(sin(z(n)*cdr) * cos(obl) / cos (z(n)*cdr)))/cdr
if n<4 then l(n)=l(n)+180
if n>3 and n<10 then l(n)=l(n)+270
if l(n)>360 then l(n)=l(n)-360
if l(n)<0 then l(n)=l(n)+180
signo(n)=int(l(n)/30)+1
grado=l(n)-30*(signo(n)-1)
grado1(n)=int(grado)
minuto(n)=int((grado - grado1(n))*60)
next

for n=1 to 12
print using "## \\ ###,.#### \\ ###,.#### ##(##)##"_
; n;" p";p(n);" l";_
l(n);grado1(n);signo(n);minuto(n)
next

Coordenadas geocentricas del sol. Programa en Basic

2013-01-02T13:56:00.000+01:00

'coordenadas del sol (basado en las formulas de Jean Meeus "Astronomical algorithms")
screen 9
cls
pi=4*atn(1)

ra=2*pi/360

'dj=2448543.5 ' dia juliano de ejemplo
dj=2448843.5

t=(dj-2451545)/36525

'Longitud geometrica media del sol
Lms=280.46646+36000.76983*t+0.0003032*t^2
print "Lms = ";Lms
call cuadrante (Lms)
print "Lms = ";Lms

'Anomalia media del sol
Ams=357.52911+35999.05029*t-0.0001537*t^2
call cuadrante (Ams)
call cuadrante (2*Ams)
call cuadrante (3*Ams)
'Ecuacion del centro del sol
Ecs=(1.914602-0.004817*t-0.000014*t^2)*sin(Ams*ra)+_
(0.019993-0.000101*t)*sin(2*Ams*ra)+_
(0.000289)*sin(3*Ams*ra)
call cuadrante(Ecs)
'Longitud verdadera del sol
Lvs=Lms+Ecs
call cuadrante(Lvs)
'correccion
u=125.04-1934.136*t
'Longitd aparente del sol
Las=Lvs-0.00569-0.00478*sin(u*ra)
call cuadrante(Las)

lo=int(Las)
si=int(lo/30)+1
gr=lo-(si-1)*30
if gr=30 then si=si+1:gr=0
if si>12 then si=si-12
mi=(las-lo)*60
m=int(mi)
s=int((mi-m)*60)

lo$=str$(lo)+"g :"+str$(m)+"m :"+str$(s)+"s"
lo1$=str$(gr)+"("+str$(si)+")"+str$(cint(mi))
print "lo.sol: ";Las#;lo$;" = ";
print lo1$

end

sub s
while inkey$="":wend
end sub

sub cuadrante(k#)
do
if k#<0 then k#=k#+360 else exit loop
loop
do
if k#>360 then k#=k#-360 else exit loop
loop
end sub

Longitud, ascension recta y declinacion del sol en un periodo determinado. Programa en Basic

2013-01-02T13:54:00.000+01:00

cls
'calculamos el valor de pi y el coeficiente ra que pasa grados
'a radianes
pi=4*atn(1):ra=2*pi/360

'damos el valor de la primera fecha de la serie que
'queremos calcular. a es el año, m es el mes y d es el dia.
'Si d es un numero fraccionario indica la fraccion decimal de ese dia
'a partir de las 00.00 h, es decir la media noche. Un dia con fraccion
'de 0.5 dara dias julianos enteros pues estos se miden desde las
'12.00 del medio dia.

print "dj = Dia juliano "
print "t = Tiempo T en centurias julianas desde J2000"
print "Lms = Longitud geometrica media del sol "
print "Ams = Anomalia media del sol"
print "Ecs = Ecuacion del centro del sol "
print "Lvs = Longitud verdadera del sol"
print "Las = Longitud aparente del sol"
print "eot = excentricidad de la orbita de la tierra"
print ""
print " dj t Lms Ams Ecs Lvs Las eot Avs rvs oe arc ds Las Las"
print "____________________________________________________________________"
'calculamos el dia juliano como se explico en otro post.
'a=2000:m=1:d=1.5
a=1992:m=10:d=13
an=a:me=m
if m<3 then
m=m+12
a=a-1
end if
A=int(a/100)
B=2-A+int(A/4)
dj = int(365.25 * (a + 4716)) + int(30.6001 * (m+1)) + d + B - 1524.5
'si queremos hacer la serie a partir de un dia juliano determinado
'hay que poner ese dato aqui.
'dj=2448976.5
'dj=2448908.5
'................................
'Hacemos un bucle a partir del dia inicial, en este caso para 20 dias mas.
for n=1 to 20
'damos los valores sucesivos
dj=dj+n-1

print dj;" ";
'..............
' t es el numero de periodos de 36525 dias desde la epoca J2000,
' es decir desde el 1.5 de enero de ese año.
t=(dj-2451545)/36525
print using ("###.#######",t);" ";
'..............
'despues calculamos la longitud media del sol con referencia
'al equinocio de cada fecha
Lms=280.46646+36000.76983*t+0.0003032*t^2
'print using ("###.#####",Lms);" ";
call cuadrante Lms
print using ("###.#####",Lms);" ";
'..............
'a continuaion se calcula la anomalia media del sol, que
'es igual a la anomalia media de la tierra.
Ams=357.52911+35999.05029*t-0.0001537*t^ 2
call cuadrante Ams
print using ("###.#####",Ams);" ";
Ams2=2*Ams
call cuadrante (Ams2)
Ams3=3*Ams
call cuadrante (Ams3)
'...................
'el siguiente paso es calcular la ecuacion del centro del sol
Ecs=(1.914602-0.004817*t-0.000014*t^2)*sin(Ams*ra)+_
(0.019993-0.000101*t)*sin(Ams2*ra)+_
(0.000289)*sin(Ams3*ra)
call cuadrante Ecs
print using ("###.#####",Ecs);" ";
'..................
'Sumando a la ecuacion anterior la longitud media calculada
'antes tenemos la longitud verdadera del sol
Lvs=Lms+Ecs
call cuadrante Lvs
print using ("###.#####",Lvs);" ";
'...................
'hacemos una correccion por la nutacion y la aberracion:
u=125.04-1934.136*t
call cuadrante u
'print "Correccion u = ";u
'y ya tenemos la longitud aparente del sol, referida al
'verdadero equinocio de la fecha considerada:
Las=Lvs-0.00569-0.00478*sin(u*ra)
call cuadrante Las
print using ("###.#####",Las);" ";
lo=int(Las)
si=int(lo/30)+1
gr=lo-(si-1)*30
if gr=30 then si=si+1:gr=0
if si>12 then si=si-12
mi=(Las-lo)*60
m=int(mi)
s=int((mi-m)*60)

lo$=str$(lo)+"g :"+str$(m)+"m :"+str$(s)+"s"
lo1$=str$(gr)+"("+str$(si)+")"+str$(int(mi))
'---------------
'A continuacion calculamos el radio vector, es decir la distancia
'entre el centro del sol y el centro de la tierra.
'(en unidades astronomicas)
'Para este calculo tenemos primero que conocer la excentricidad de la
' tierra y la anomalia verdadera del sol:
'La excentricidad de la orbita de la tierra es:
eot=0.016708634 - 0.000042037*t - 0.0000001267*t^2
print using ("##.########",eot);" ";
'................
'para saber la anomalia verdadera del sol basta con sumar
'a la anomalia media su ecuacion del centro:
Avs=Ecs+Ams
call cuadrante Avs
print using ("###.#####",Ams);" ";
'-------------------
' y el radio vector del sol sera:
rvs=1.000001018*(1-eot^2)/(1+eot*cos(Avs*ra))
'print 1.000001018*(1-eot^2)
print using ("##.########",rvs);" ";
'............
'Por ultimo, para calcular la ascension reta y la declinacion, debemos
'primero calcular la oblicuidad de la ecliptica:
oe=(23+26/60+21.448/3600)-46.8150/3600*t-0.00059/3600*t^2+0.001813/3600*t^3
print using ("##.########",oe);" ";
'..............
'y por tanto la ascension recta del sol sera:
oeu=oe+0.00256*cos(u*ra)
call cuadrante oeu
arc=atn(cos(ra*oeu)*sin(Las*ra)/cos(Las*ra))/ra+180
call cuadrante arc
print using ("###.#####",arc);" ";
'................
'y la declinacion
ds=asn(sin(ra*oeu)*sin(Las*ra))/ra
print using ("###.#####",ds);" ";
'..............
print using ("###",lo);"º ";using ("##",m);"m ";using ("##",s);"s ";" = ";lo1$

next n
end

'-----------------
sub s
while inkey$="":wend
end sub

sub cuadrante byref k
if k<0 then
do: k=k+360: loop until k>0
end if

if k>360 then
do: k=k-360: loop until k<360
end if
end sub
'Antonio Martinez

Calculo del Dia Juliano. Programacion

2013-01-02T13:52:00.000+01:00

'julia-da.bas
'--------------------------
cls
'd=19 :m=3 :a=1920 'fechas de ejemplo

'd=4.81:m=10:a=1957
'd=27.5:m=1 :a=333 'calendario juliano
'd=17.9:m=8 :a=-1001
'd=29 :m=2 :a=-1000
'd=1.5 :m=1 :a=-4712
d=17 :m=11:a=1858

'Este programa para calcular el día juliano,
'más completo que el dado antes en este blog
'permite la introducción de fechas negativas,
'hasta el día juliano 0.
'Las fechas deben introducirse en forma astronómica
'y no histórica, por ejemplo el año anterior al
'año 1, segun los historiadores es el año -1,
'sin embargo astronómicamente es el año 0.
'Si tenemos un año como el 453BC 453 antes de Cristo) 'astronómicamentes será el año -452

'lo primero que hacemos es esta conversion de los valores de m
' de a
if m<3 then a=a-1:m=m+12
x=int(a/100)

y=2-x+int(x/4) 'fechas dadas en calendario gregoriano
y=0 ' en fechas dadas en calendario juliano
' segun el caso hay que quitar una de las dos instrucciones
'anteriores

DJ#=int(365.25*(a+4716))+int(30.6001*(m+1))+d+y-1524.5
' formula de Jean Meeus en "Astronomical algorithms" 1998
if DJ#<=2299149.5 then DJ#=DJ#-y
' dia juliano del 4-Octubre-1582 fecha de
' inicio del calendario gregoriano.
' esta instruccion equivale a poner y=0
' si la fecha es juliana, pero evita
' modificar el programa en ese caso
print using "########,.##";DJ#

end

'NOTAS:
'Este programa esta escrito originalmente en Turbobasic,
'si se va a traducir a otro basic mas moderno, por
'ejemplo a Justbasic, basta con quitar las # de las
'variables.

'** En algunos lenguajes de programación,
'por ejemplo Fortram, la instrucción int aplicada
'a números negativos solo elimina la parte decimal
'dejando la misma parte entera, por ejemplo,
'en Fortram int(-5.12)=-5, en Basic int(-5.12)=-6
'es decir el entero inmediatamente inferior.
'En estas formulas se utiliza este sentido por lo
'que debe tenerse cuidado si se traduce el
'programa a otro lengaje distinto del BASIC

'** Si vemos, en algun texto astronómico,
'el Día Juliano Modificado "MJD", se
'refiere al número de días respecto del 17
'de noviembre de 1858 a las 00h UT,
'MJD = DJ-DJ(17/11/1858/0:00) = DJ-2,400,00.5

'** Julio Cesar unificó el calendario romano en el
'año -45, pero su forma definitiva no llego hasta
'el año +8. A pesar de esto las fechas anteriores
'a la reforma gregoriana, aunque sean anteriores
'al año -45, se suelen dar en calendario juliano
'en la practica astronómica.

'** Los días julianos comienzan el 1 de enero del
'año 4.713BC a las 12:00h del medio día, así que
'el día juliano 0 será el día 1.5 del mes 1 del
'año -4712

Longitud, latitud, distancia y paralaje de la Luna

2013-01-02T13:51:00.000+01:00

Fecha de ejemplo: año a = 1972 mes m = 3 dia d = 6
s=seno
asn(a)=arco cuyo seno vale a

Si d es un numero fraccionario indica la fraccion decimal de ese dia
a partir de las 00.00 h, es decir la media noche. Un dia con fraccion
de 0.5 dara dias julianos enteros pues estos se miden desde las 12.00 del medio dia.
si m fuese 1 o 2 hacemos m=m+12 y a=a-1

Valores auxiliares para el calculo de dj A=int(a/100)= 19 B=2-A+int(A/4) =-13
Dia Juliano dj=int(365.25*(a+4716))+int(30.6001*(m+1))+d+B-1524.5= 2441382.5 dias
t es el numero de periodos de 36525 dias desde la epoca J2000,
es decir desde el 1.5 de enero de ese año. t=(dj-2451545)/36525 = -0.2782341 centurias

longitud media de la luna referida al equinocio medio de la centuria t
lml=218.3164477+481267.88123421*t-0.0015786*t^2+(t^3)/538841-(t^4)/65194=233.187º
Elongacion media de la luna
Eml=297.8501921+445267.1114034*t-0.0018819*t^2+(t^3)/545868-(t^4)/113065= 249.3621º
Anomalia media del sol
Ams=357.5291092+35999.0502909*t-0.0001536*t^2+(t^3)/2449 = 61.3662 º
Anomalia media de la luna
Aml=134.9633964+477198.8675055*t+0.0087414*t^2+(t^3)/69699-(t^4)/14712= 201.9732º
Argumento de latitud de la luna
All=93.272095+483202.0175233*t-0.0036539*t^2-(t^3)/3526+(t^4)/86331 = 290 º
Argumento a1=119.75+ 131.849*t = 83.07º
" a2= 53.09+479264.290*t = 265.43º
" a3=313.45+481266.484*t = 328.71º
Adicion a long.(Venus y Jupiter) ad.lo=3958*sen(a1)+1962*sen(l)+318*sen(a2)=1970
Adicion a la latitud (a=a1-All l=lml-Aml)
ad.la=-2235*s(lml)+382*s(a3)+175*s(a)+175*s(a)+127*s(l)-115*s(lml+Aml) = 1664
Coeficiente de excentricidad de la Tierra e=1-0.002516*t-0.0000074*t^2 = 1.000699

Terminos periodicos principales
de longitud y distancia
co
1 2 3 4 5 6 7
_______________________________
1, 0, 0, 1, 0,6288774,-20905355
2, 2, 0,-1, 0,1274027, -3699111
3, 2, 0, 0, 0, 658314, -2955968
4, 0, 0, 2, 0, 213618, -569925
5, 0, 1, 0, 0,-185116, 48888
6, 0, 0, 0, 2,-114332, -3149
7, 2, 0,-2, 0, 58793, 246158
________________________________
Argumento del sumando a1=co(2)*Eml+co(3)*Ams+co(4)*Aml+co(5)*All = 201.9732
Sumando de longitud del termino periodico 1 l=co(6)*sin(a1) =-2353090
" distancia " " 1 r=co(7)*cos(" ) =19386766.3
a2 = 296.7511
" longitud " " 2 l=co(6)*sin(" ) =-1137668
" distancia " " 2 r=co(7)*cos(" ) =-1665029.01
a3 = 498.7243
" longitud " " 3 l=co(6)*sin(" ) = 434278
" distancia " " 3 r=co(7)*cos(" ) =2221541.41
a4 = 403.9464
" longitud " " 4 l=co(6)*sin(" ) = 148248
" distancia " " 4 r=co(7)*cos( ") =-410339.694
a5 = 61.3662
" longitud " " 5 l=co(6)*sin(" ) =-162590
" distancia " " 5 r=co(7)*cos(" ) = 23427.5791
a6 = 580
" longitud " " 6 l=co(6)*sin(" ) = 73491
" distancia " " 6 r=co(7)*cos(" ) = 2412.2742
a7 = 94.7778
" longitud " " 7 l=co(6)*sin(" ) = 58589
" distancia " " 7 r=co(7)*cos(" ) = -20503.3952
Suma de terminos periodicos de longitud (-1127527) = -2923699
" " " distancia (-16590875) = 19542553
" " " latitud (-3229126) = -4961476

longitud de la luna l.l =lml+su.l/1000000 = 230.2634
latitud de la luna la.l=su.la/1000000 = 355.0385
distancia a la luna d.l =385000.56+su.r/1000 = 404543.1 km
paralaje ecuatorial horizontal peh =asn(6378.14/d.l) = 0.9034

DIA JULIANO

2011-09-26T22:07:00.001+02:00

DÍA JULIANO:

La forma de referirnos a una fecha en astronomía no es mediante el año, més y día, sino mediante el Día Juliano.

El Día Juliano es el número de días desde el día que consideramos día 0 a las 12 horas (Medio Día), así pués el Día Juliano correspondiente a un día cualquiera a las 0 horas (media noche) será un número entero más la fracción 0.5 correspondiente a las 12 horas de diferencia. Con el Día Juliano podemos realizar cualquier cálculo astronómico, por ejemplo podemos calcular la Hora Sideral y con esta podemos calcular el Ascendente, también es posible el cálculo de las efemérides de cualquier planeta o satélite. En general todos los cálculos astronómicos precisan la definición de la fecha mediante el Día Juliano.

Para realizar el cómputo del Día Juliano se ha considerado mediante convenio internacional que las fechas comienzan a las 12 horas del 1 de Enero del año -4712, de esta forma nos aseguramos que cualquier fecha histórica tendrá un Día Juliano positivo. (No se debe confundir el Día Juliano con el Calendario Juliano o con una fecha juliana pues son cosas totálmente distintas).

La forma de cálculo que expongo en este post, basado en la fórmula indicada por el Anuario Astronómico de España, tiene limitaciones importantes, pues no permite calcular el Día Juliano para fechas antiguas, si se desea un cálculo más general, que incluya fechas tanto las fechas gregorianas como julianas (las del Calendario Juliano), debe utilizarse el programa que he publicado en este otro post, pulsando aqui, y que esta basado en las formulas mas precisas, de Jean Meeus.

Este cálculo utiliza dos tablas:

Tabla 1

Centuria DJC
____________________

1600... 2305447.5

1700.... 2341971.5

1800... 2378495.5

1900... 2415019.5

2000... 2451544.5

2100... 2488068.5
____________________

tabla 2
_______________
Mes..... DJM bis.
_______________

ENE....... 0 ...-1

FEB..... 31 ...-1

MAR.... 59 ....0

ABR..... 90 ....0

MAY.. 120 ....0

JUN.... 151 ....0

JUL.... 181 ....0
AGO... 212 ....0

SEP.... 243 ....0

OCT.... 273 ....0

NOV... 304 ....0

DIC..... 334 ....0

____________

Con las tablas anteriores podremos calcular cualquier Dia Juliano mediante siguiente formula:

DJ = DJc + DJf + DJm + bis + día + h/24

En la formula anterior tendremos que:
DJ = Dia Juliano.
DJc = Dia Juliano de la centuria.
Df = Numero de dias de la fraccion de centuria.
Dm = Numero de dias desde el 1 de enero hasta el dia 1 del mes de la fecha.
bis = indicador de año bisiesto.
dia = Numero de dia del mes de la fecha.
h = hora a partir de las 0 horas (sistema de 24 horas)

Por ejemplo: supongamos que queremos calcular el Dia Juliano correspondiente al 14 de Junio de 1997 a las 00:00 horas.

DJc=2415019.5 (valor de la tabla 1 correspondiente a 1900)
Df =365*97+INT(97/4)=35429
Dm = 151 (Valor de la tabla 2 para el mes de junio)
bis= 0 (Valor de la tabla 2 para el mes de junio)
dia = 14

h = 0 (0/24)

Notas:

a) En la ecuación INT(97/4) significa la parte entera de dividir las dos ultimas cifras del año por cuatro.

b) Si en lugar de ser las 0 horas, hubiera sido una hora cualquiera del dia, h seria esa hora dividida por 24.

DJ = 2415.019'5 + 35.429 + 151 + 0 + 14 + 0 = 2.450.613'5

Calculo aproximado de efemerides planetarias. Programa en Basic

2011-09-23T20:08:00.000+02:00

' EFEMERIDES PLANETARIAS
cls

dim b$(32)
c%=360
pi#=4*atn(1)
cd#=360/365.2422 ' grados recorridos en un dia medio
ra#=180/pi# ' factor para transformar radianes en grados multiplicando
' grados en radianes dividiendo
rr#=360/pi# ' el doble de ra#

for f=1 to 8
read a$(f),t(f),w(f),p(f),e(f),s(f),i(f),no(f)
next f

principio:
'input "Fecha aaaa,mm,dd", y,m,d
y=1984:m=7:d=12

' calculo del periodo de tiempo en dias desde el 1 de enero de 1980
if m<3 then k=m else k=0
j=int(365.25*(y-k))+int(30.4375*(m+1+12*k))+d-int(y/100)+_
int(int((y/100)/4))+1720996.5

ed%=j-2444238.5+1 ' periodo de dias desde el dia con datos iniciales
print ed%
' POSICION DE LA TIERRA
'ed%=1655
we#=98.83354 ' longitud del periodo
' 102.596403= longitud del perihelio
ee#=0.016718 ' excentricidad de la tierra

' cd#*ed% = grados que recorre en los dias del periodo considerado
' en su orbita, referido al
' 1.00004 = periodo orbital de la tierra en a¤os sidereos
z#=cd#*ed%/1.00004 :gosub rango:ne#=z#

'le sumamos la longitud del punto de partida y le restamos la del perihelio.
' y tenemos la distancia desde el punto a calcular hasta el perihelio
z#=ne#+we#-102.596403 :gosub rango:me#=z#

z#=ne#+rr#*ee#*sin(me#/ra#)+we#:gosub rango:le#=z#
ve#=le#-102.596403 ' anomalia verdadera
re#=(1-ee#^2)/(1+ee#*cos(ve#/ra#)) ' radio vector

print " EFEMERIDES PLANETARIAS"
print
print " Planeta Ascension Recta Declinacion"
d1$= "\ \ ### \\ ## \\ ## \\ ### \\ ## \\ ## \\ ## ## ##"

for n=1 to 8
z#=cd#*ed%/t(n) :gosub rango:np#=z#
z#=np#+w(n)-p(n) :gosub rango:mp#=z#
z#=np#+rr#*e(n)*sin(mp#/ra#)+w(n):gosub posit:lp#=z#
vp#=lp#-p(n)
rp#=s(n)*(1-e(n)^2)/(1+e(n)*cos(vp#/ra#))

pps=sin((lp#-no(n))/ra#)*sin(i(n)/ra#)
ppc=sqr(1-pps^2)
psi=ra#*atn(pps/ppc)
yx=lp#-no(n)
y=sin(yx/ra#)*cos(i(n)/ra#)
x#=cos(yx/ra#)
ct=ra#*atn(y/x#)
q=ct: gosub cuadrante: ct=q
l1=ct+no(n)
r1=rp#*cos(psi/ra#)

if n>2 then goto 430 ' Planetas interiores Mercurio y venus
ll=le#-l1
a=ra#*atn((r1*sin(ll/ra#))/(re#-r1*cos(ll/ra#)))
z#=180+le#+a:gosub posit:lam=z#
goto 470

430 ' planetas exteriores jupiter - pluton
ll=l1-le#
' longitud ecliptica
z#=l1+ra#*atn((re#*sin(ll/ra#)/(r1-re#*cos(ll/ra#)))):gosub posit:lam=z#
470 ' latitud ecliptica
bet=ra#*atn(r1*tan(psi/ra#)*sin((lam-l1)/ra#)/(re#*sin((l1-le#)/ra#)))

' coordenadas ecliptico-ecuatoriales
la=lam/ra#
ec=23.441884/ra#
be=bet/ra#
y=sin(la)*cos(ec)-tan(be)*sin(ec)
x#=cos(la)
rh#=ra#*atn((sin(la)*cos(ec)-tan(be)*sin(ec))/(cos(la)))

' va al cuadrante
q=rh# : gosub cuadrante: rh#=q
' formato en signos, grados y minutos
si=int(rh#/30)+1
grr=rh#-(si-1)*30
gr=int(grr)
mi=(grr-gr)*60
' ascension recta en formato en horas, m. y s.
rx=rh#/15
ry=int(rx)
mx=60*(rx-ry)
my=int(mx)
sx=cint(60*(mx-my))
' declinacion en formato en grados, m. y s.
pps#=sin(be)*cos(ec)+cos(be)*sin(ec)*sin(la)
ppc#=sqr(1-pps#^2)
dc=ra#*atn(pps#/ppc#)
sig=1:if dc<0 then sig=-1
dc1=abs(dc)
dh=int(dc1)
dm=60*(dc1-dh)
dmi=int(dm)
ds=cint(60*(dm-dmi))
dh=sig*dh
print using d1$;a$(n);ry;"h";my;"m";sx;"s";dh;"§";dmi;"m";ds;"s";si;gr;mi
next n

' calculo de longitud y latitud
for n=1 to 8

next

' dibuja los planetas en un atlas estelar
print "t para terminar , s para otra fecha"
a$=input$(1)
cls
if a$="t" then end else goto principio

end

' rango de 0-360 grados
rango:
z#=360*(z#/360-int(z#/360))
return

' numeros positivos
posit:
if z#>360 then z#=z#-360
if z#<0 then z#=z#+360
return

cuadrante:
if x#<0 then q=q+180:return
if y>0 then return
q=q+360
return

sub s
print "para seguir pulsa una tecla"
while inkey$="" :wend
end sub

' (1) (2) (3) (4) (5) (6) (7)
' t we# p ee# s i n
' periodo long. peri. l.perihelio excentr. semie ma. inclinac. nodo asc.
data "Mercu", 0.24085,231.2973 , 77.1442128,.2056306,0.3870986,7.0043579 ,48.0941733
data "Venus", 0.61521,355.73352, 131.2895792,.0067826,0.7233316,3.394435 ,76.4997524
data "Marte", 1.88089,126.30783, 335.6908166,.0933865,1.5236883,1.8498011 ,49.4032001
data "Jupit", 11.86224,146.966365, 14.0095493,.0484658,5.202561 ,1.3041819 ,100.2520175
data "Satur", 29.45771,165.322242, 92.6653974,.0556155,9.554747 ,2.4893741 ,113.4888341
data "Urano", 84.01247,228.070855, 172.7363288,.0463232,19.21814 ,0.7729895 ,73.8768642
data "Neptu",164.79558,260.3578998, 47.8672148,.0090021,30.10957 , 1.7716017,131.5606494
data "Pluto",250.9 ,209.439 ,222.972 ,.25387 ,39.78459 ,17.137 ,109.941
' tierra 1.00004 98.83354 102.596403 .016718 1 0 0
' (1) periodo orbital en a¤os sidereos
' (2) longitud el 1 de enero de 1980
' (3) longitud del perihelio
' (4) excentricidad de la orbita
' (5) valor medio del semieje mayor en relacion con el valor de la tierra
' (6) declinacion de la orbita a la ecliptica
' (7) nodo ascendente

Programa en Basic. Cuspides, Ascendente y Medio Cielo

2011-09-06T21:43:00.000+02:00

'casa10.bas

' Hora Sid Casa X CASAXI CASA XII Casa I CASA II CASA III
data "00:00:00"," 0( 1)"," 8.6( 2)","22.0( 3)","26( 4)10","12.4( 5)"," 2.4( 6)"
data "00:03:40"," 1( 1)"," 9.6( 2)","22.9( 3)","26( 4)51","13.1( 5)"," 3.4( 6)"
data "02:03:08"," 3( 2)","12.4( 3)","19.2( 4)","18( 5)12"," 6.6( 6)"," 0.7( 7)"
data "04:16:27"," 6( 3)","13.4( 4)","15.6( 5)","11( 6)39"," 3.3( 7)"," 1.3( 8)"
data "05:55:38","29( 3)"," 5.0( 5)"," 5.1( 6)","29( 6)13","23.9( 7)","23.2( 8)"
data "06:00:00"," 0( 4)"," 5.9( 5)"," 6.0( 6)"," 0( 7)00","24.0( 7)","24.1( 8)"
data "06:04:22"," 1( 4)"," 6.8( 5)"," 6.9( 6)"," 0( 7)47","24.9( 7)","25.0( 8)"
data "08:12:54"," 1( 5)"," 5.4( 6)"," 2.6( 7)","23( 7)31","20.1( 8)","23.2( 9)"
data "10:23:35"," 4( 6)"," 5.4( 7)","28.6( 7)","16( 8)30","16.3( 9)","24.4(10)"
data "11:56:20","29( 6)","26.8( 7)","16.9( 8)"," 3( 9)09"," 7.1(10)","20.4(11)"
data "12:00:00"," 0( 7)","27.6( 7)","17.6( 8)"," 3( 9)50"," 8.0(10)","21.4(11)"
data "12:03:40"," 1( 7)","28.3( 7)","18.3( 8)"," 4( 9)30"," 8.9(10)","22.6(11)"
data "14:30:21","10( 8)"," 0.7( 9)","17.7( 9)"," 4(10)55","23.6(11)","10.0( 1)"
data "16:46:16","13( 9)"," 0.7(10)","12.9(10)","20(11)09","19.4( 1)","22.0( 2)"
data "17:55:38","29( 9)","17.2(10)","11.2(11)","27(12)27","15.8( 2)","10.5( 3)"
data "18:00:00"," 0(10)","18.3(10)","12.7(11)"," 0( 1)00","17.3( 2)","11.7( 3)"
data "18:04:22"," 1(10)","19.5(10)","14.2(11)"," 2( 1)33","18.8( 2)","12.8( 3)"
data "20:41:41"," 8(11)"," 4.4(12)","17.5( 1)","12( 3)18"," 1.8( 4)","19.0( 4)"
data "22:49:53","11(12)","16.4( 1)"," 3.1( 3)","12( 4)43","28.7( 4)","16.8( 5)"
data "23:56:20","29(12)"," 7.4( 2)","21.1( 3)","25( 4)30","11.7( 5)"," 1.7( 6)"
data "24:00:00"," 0( 1)"," 8.6( 2)","22.0( 3)","26( 4)10","12.4( 5)"," 2.4( 6)"
'..................
for nn=1 to 21:read hs$,h10$,h11$,h12$,h1$,h2$,h3$: next
la$="51N00"
'..................
cls

e$="23:27:15"
pi#=4*atn(1)
ra#=2*pi#/360
'......................
hs#=val(left$(hs$,2))+val(mid$(hs$,4,2))/60+val(right$(hs$,2))/3600
if hs#=24 then hs#=0
la=val(left$(la$,2))+val(right$(la$,2))/60
la.c=la-11/60-29/(60*60) 'correccion de la latitud geografica

e#=val(left$(e$,2))+val(mid$(e$,4,2))/60+val(right$(e$,2))/3600

' ................

'Ascension Recta del Medio Cielo
armc=hs#*360/24

'Longitud del Medio Cielo
MC=atn(tan((armc)*ra#)/cos(e#*ra#))/ra#

if MC<0 then MC=MC+360
MC=MC+180
if MC>=360 then MC=MC-360

if hs#>=0 and hs#=<6 then MC=MC-180
if hs#>18 and hs#=<24 then MC=MC+180
'................
' polos
senx#=tan(e#*ra#)*tan(la.c*ra#)
tanx#=senx#/sqr(1-senx#^2)
ax#=atn(tanx#) ' ax en radianes
if ax#<0 then ax#=ax#+2*pi#

for n=1 to 5

if n=1 then casa$="Casa XI" :h$(n)=h11$
if n=2 then casa$="Casa XII":h$(n)=h12$
if n=3 then casa$="Casa I" :h$(n)=h1$
if n=4 then casa$="Casa II" :h$(n)=h2$
if n=5 then casa$="Casa III":h$(n)=h3$
a=0
t=n*30

a=armc+t 'Ascension Oblicua de la cuspide

if n=1 or n=5 then tanpolo#=sin(ax#/3)/tan(e#*ra#) ' casas XI y III
if n=3 then tanpolo#=tan(la*ra#) ' casa I
if n=2 or n=4 then tanpolo#=sin(2*ax#/3)/tan(e#*ra#) ' casas XII y II

if a= 90 then a= 90.0010
if hs#=12 and a=270 then a=270.0026
if a=270 then a=270.0001
if a>=360 then a=a-360

'Longitud de la cuspide
AA#= (atn(tanpolo#/cos(a*ra#)))/ra#
cus(n)= (atn(tan(a*ra#)*cos(AA#*ra#)/cos((AA#+e#)*ra#)))/ra#

cus(1)=cus(1)-180
cus(2)=cus(2)-180

if cus(n)<0 and hs#>12 then cus(n)=cus(n)+180

do
if cus(n)<0 then cus(n)=cus(n)+360 else exit loop
loop
' print n;"cus ";cus(n)

if hs#=<6 and cus(1)> 90 and cus(1)<180 then cus(1)=cus(1)+180
if hs#=<6 and cus(2)> 90 and cus(2)<180 then cus(2)=cus(2)+180
if hs#>6 and hs#=<12 and cus(1)>180 and cus(1)<270 then cus(1)=cus(1)+180
if hs#>6 and hs#=<12 and cus(2)>180 and cus(2)<270 then cus(2)=cus(2)+180
if hs#>6 and hs#=<12 and cus(4)>270 and cus(4)<360 then cus(4)=cus(4)+180
if hs#>6 and hs#=<12 and cus(5)>270 and cus(5)<360 then cus(5)=cus(5)+180
if hs#>12 and hs#<18 and cus(1)>270 and cus(1)<360 then cus(1)=cus(1)+180
if hs#>12 and hs#<18 and cus(2)>270 and cus(2)<360 then cus(2)=cus(2)+180
if hs#>12 and hs#<18 and cus(4)> 0 and cus(4)< 90 then cus(4)=cus(4)+180
if hs#>12 and hs#<18 and cus(5)> 0 and cus(5)< 90 then cus(5)=cus(5)+180
if hs#>=18 and cus(1)>270 and cus(1)<360 then cus(1)=cus(1)+180
if hs#>=18 and cus(2)>270 and cus(2)<360 then cus(2)=cus(2)+180
if hs#>18 and cus(1)>0 and cus(1)< 90 then cus(1)=cus(1)+180
if hs#>18 and cus(2)>0 and cus(2)< 90 then cus(2)=cus(2)+180
if hs#>=18 and cus(4)>0 and cus(4)< 90 then cus(4)=cus(4)+180
if hs#>=18 and cus(5)>0 and cus(5)< 90 then cus(5)=cus(5)+180
if hs#>18 and cus(4)>90 and cus(4)<180 then cus(4)=cus(4)+180
if hs#>18 and cus(5)>90 and cus(5)<180 then cus(5)=cus(5)+180
' print n;"cus ";cus(n)

cus(n)=cus(n)+180
if cus(n)>=360 then cus(n)=cus(n)-360
if hs#>=18 then cus(3)=cus(3)-180

call trans(cus(n),cus$(n))
call trans(MC,MC$)

'..................
d$= "\ \ ###,.## \ \\ \"
if n=1 then print "latitud: ";la$;" hora sideral: ";hs$
if n=1 then print "----------------------------------------------"
if n=1 then print using d$;"Medio Cielo ";MC;MC$;h10$
print using d$;casa$ ;cus(n);cus$(n);h$(n)
if n=5 then print ".............................................."
next
'next
end
'..................

sub s: while inkey$="":wend: end sub
'..................
sub trans(x,x$)
s=int(x/30)+1
ga=x-(s-1)*30:g=int(ga)
m=cint((ga-g)*60)
if m=60 then m=0:g=g+1
if g=30 then g=0:s=s+1
if s=13 then s=1
x$=str$(g)+"("+str$(s)+")"+str$(m)
end sub

Determinacion del Medio Cielo y Ascendente

2011-09-06T21:40:00.000+02:00

'casa9.bas

for n=1 to 21:read hs$,h11$,h12$,hs1$,hs2$,hs3$:next

' Hora Sid Casa X Casa I CASAXI CASA XII CASA II CASA III
data "00:00:00"," 0( 1)","26( 4)10"," 8.6( 2)","22.0( 3)","12.4( 5)"," 2.4( 5)"
data "00:03:40"," 1( 1)","26( 4)51"," 9.6( 2)","22.9( 3)","13.1( 5)"," 3.4( 6)"
data "02:03:08"," 3( 2)","18( 5)12","12.4( 3)","19.2( 4)"," 6.6( 6)"," 0.7( 7)"
data "04:16:27"," 6( 3)","11( 6)39","13.1( 3)","15.2( 5)"," 3.3( 7)"," 1.4( 8)"
data "05:55:38","29( 3)","29( 6)13"," 5.0( 5)"," 5.1( 6)","23.9( 7)","23.2( 8)"
data "06:00:00"," 0( 4)"," 0( 7)00"," 5.9( 5)"," 6.0( 7)","24.0( 7)","24.1( 8)"
data "06:04:22"," 1( 4)"," 0( 7)47"," 6.8( 5)"," 6.9( 6)","24.9( 7)","25.0( 8)"
data "08:12:54"," 1( 5)","23( 7)31"," 5.4( 6)"," 2.6( 7)","20.1( 8)","23.2( 9)"
data "10:23:35"," 4( 6)","16( 8)30"," 5.4( 7)","28.6( 7)","16.3( 9)","24.4(10)"
data "11:56:20","29( 6)"," 3( 9)09","26.8( 7)","16.9( 8)"," 7.1(10)","20.4(11)"
data "12:00:00"," 0( 7)"," 3( 9)50","27.6( 7)","17.6( 8)"," 8.0(10)","21.4(11)"
data "12:03:40"," 1( 7)"," 4( 9)30","28.3( 7)","18.3( 8)"," 8.9(10)","22.6(11)"
data "14:30:21","10( 8)"," 4(10)55"," 0.7( 9)","17.7( 9)","23.6(10)","10.0( 1)"
data "16:46:16","13( 9)","20(11)09"," 0.7(10)","12.9(10)","19.4( 1)","22.0( 2)"
data "17:55:38","29( 9)","27(12)27","17.2(10)","11.2(11)","15.8( 2)","10.5( 3)"
data "18:00:00"," 0(10)"," 0( 1)00","18.3(10)","12.7(11)","17.3( 2)","11.7( 3)"
data "18:04:22"," 1(10)"," 2( 1)33","19.5(10)","14.2(11)","18.8( 2)","12.8( 3)"
data "20:41:41"," 8(11)","12( 3)18"," 4.4(12)","17.5( 1)"," 1.8( 3)","19.0( 3)"
data "22:49:53","11(12)","12( 4)43","16.4( 1)"," 3.1( 3)","28.7( 4)","16.8( 5)"
data "23:56:20","29(12)","25( 4)30"," 7.4( 2)","21.1( 3)","11.7( 5)"," 1.7( 6)"
data "24:00:00"," 0( 1)","26( 4)10"," 8.6( 2)","22.0( 3)","12.4( 5)"," 2.4( 6)"
la$="51N00"
'..................
cls
e$="23:28:00"
pi#=4*atn(1)
ra#=2*pi#/360
'......................
hs=val(left$(hs$,2))+val(mid$(hs$,4,2))/60+val(right$(hs$,2))/3600
if hs=24 then hs=0
la=val(left$(la$,2))+val(right$(la$,2))/60
e#=val(left$(e$,2))+val(mid$(e$,4,2))/60+val(right$(e$,2))/3600
' ................

armc=hs*360/24 'Ascension Recta del Medio Cielo
MC=atn(tan((armc)*ra#)/cos(e#*ra#))/ra# 'Longitud del Medio Cielo

if MC<0 then MC=MC+360
MC=MC+180
if MC>=360 then MC=MC-360

if hs>=0 and hs=<6 then MC=MC-180
if hs>18 and hs=<24 then MC=MC+180
'................

'a=armc+90 'Ascension Oblicua del Ascendente
a=armc+120

if a= 90 then a= 90.0010
if hs=12 and a=270 then a=270.0026
if a=270 then a=270.0001
if a>=360 then a=a-360

AA= (atn((tan(la*ra#))/cos(a*ra#)))/ra#
H1= (atn(tan(a*ra#)*cos(AA*ra#)/cos((AA+e#)*ra#)))/ra# 'Longitud del Ascendente

if H1<0 and hs>12 then H1=H1+180
if H1<0 then H1=H1+360
H1=H1+180
if H1>=360 then H1=H1-360

if hs>=18 then H1=H1-180

senx=tan(e#*ra#)*tan(la*ra#)
cosx=sqr(1-senx^2)
tanx=senx/cosx
x=atn(tanx) ' x en radianes
tanpolo11=sin(2*x/3)*(1/tan(e#*ra#))

AA= (atn(tanpolo11/cos(a*ra#)))/ra#
H11= (atn(tan(a*ra#)*cos(AA*ra#)/cos((AA+e#)*ra#)))/ra# 'Longitud de la casa 11
if H11<0 then H11=H11+360
print "H11";H11

'..................
call trans(H1,H1$)
call trans(MC,MC$)
print "latitud: ";la$;" hora sideral: ";hs$
d$= "\ \ ###,.## \ \ \ \"
print using d$;"Medio Cielo ";MC;MC$;h10$
print using d$;"Ascendente ";H1;H1$;h1$
end
'..................
sub s: while inkey$="":wend: end sub
'..................
sub trans(x,x$)
s=int(x/30)+1
ga=x-(s-1)*30:g=int(ga)
m=cint((ga-g)*60)
if m=60 then m=0:g=g+1
if g=30 then g=0:s=s+1
if s=13 then s=1
x$=str$(g)+"("+str$(s)+")"+str$(m)
end sub

Determinacion de cuspides. Sistema de Placidus

2011-09-05T18:51:00.000+02:00

'placidus.bas

' Hora Sid Casa X CASAXI CASA XII Casa I CASA II CASA III
data "00:00:00"," 0( 1)"," 8.6( 2)","22.0( 3)","26( 4)10","12.4( 5)"," 2.4( 6)"
data "00:03:40"," 1( 1)"," 9.6( 2)","22.9( 3)","26( 4)51","13.1( 5)"," 3.4( 6)"
data "02:03:08"," 3( 2)","12.4( 3)","19.2( 4)","18( 5)12"," 6.6( 6)"," 0.7( 7)"
data "04:16:27"," 6( 3)","13.4( 4)","15.6( 5)","11( 6)39"," 3.3( 7)"," 1.3( 8)"
data "05:55:38","29( 3)"," 5.0( 5)"," 5.1( 6)","29( 6)13","23.9( 7)","23.2( 8)"
data "06:00:00"," 0( 4)"," 5.9( 5)"," 6.0( 6)"," 0( 7)00","24.0( 7)","24.1( 8)"
data "06:04:22"," 1( 4)"," 6.8( 5)"," 6.9( 6)"," 0( 7)47","24.9( 7)","25.0( 8)"
data "08:12:54"," 1( 5)"," 5.4( 6)"," 2.6( 7)","23( 7)31","20.1( 8)","23.2( 9)"
data "10:23:35"," 4( 6)"," 5.4( 7)","28.6( 7)","16( 8)30","16.3( 9)","24.4(10)"
data "11:56:20","29( 6)","26.8( 7)","16.9( 8)"," 3( 9)09"," 7.1(10)","20.4(11)"
data "12:00:00"," 0( 7)","27.6( 7)","17.6( 8)"," 3( 9)50"," 8.0(10)","21.4(11)"
data "12:03:40"," 1( 7)","28.3( 7)","18.3( 8)"," 4( 9)30"," 8.9(10)","22.6(11)"
data "14:30:21","10( 8)"," 0.7( 9)","17.7( 9)"," 4(10)55","23.6(11)","10.0( 1)"
data "16:46:16","13( 9)"," 0.7(10)","12.9(10)","20(11)09","19.4( 1)","22.0( 2)"
data "17:55:38","29( 9)","17.2(10)","11.2(11)","27(12)27","15.8( 2)","10.5( 3)"
data "18:00:00"," 0(10)","18.3(10)","12.7(11)"," 0( 1)00","17.3( 2)","11.7( 3)"
data "18:04:22"," 1(10)","19.5(10)","14.2(11)"," 2( 1)33","18.8( 2)","12.8( 3)"
data "20:41:41"," 8(11)"," 4.4(12)","17.5( 1)","12( 3)18"," 1.8( 4)","19.0( 4)"
data "22:49:53","11(12)","16.4( 1)"," 3.1( 3)","12( 4)43","28.7( 4)","16.8( 5)"
data "23:56:20","29(12)"," 7.4( 2)","21.1( 3)","25( 4)30","11.7( 5)"," 1.7( 6)"
data "24:00:00"," 0( 1)"," 8.6( 2)","22.0( 3)","26( 4)10","12.4( 5)"," 2.4( 6)"
'..................
for nn=1 to 2:read hs$,h10$,h11$,h12$,h1$,h2$,h3$: next
la$="51N00"
'..................
cls

e$="23:27:15"
pi#=4*atn(1)
ra#=2*pi#/360
'......................
hs#=val(left$(hs$,2))+val(mid$(hs$,4,2))/60+val(right$(hs$,2))/3600
if hs#=24 then hs#=0
la=val(left$(la$,2))+val(right$(la$,2))/60
la.c=la-11/60-29/(60*60) 'correccion de la latitud geografica

e#=val(left$(e$,2))+val(mid$(e$,4,2))/60+val(right$(e$,2))/3600

' ................

'Ascension Recta del Medio Cielo
armc=hs#*360/24

'Longitud del Medio Cielo
MC=atn(tan((armc)*ra#)/cos(e#*ra#))/ra#

if MC<0 then MC=MC+360
MC=MC+180
if MC>=360 then MC=MC-360

if hs#>=0 and hs#=<6 then MC=MC-180
if hs#>18 and hs#=<24 then MC=MC+180
'................

for n=1 to 2
if n=1 then casa$="Casa XI" :h$(n)=h11$
if n=2 then casa$="Casa XII":h$(n)=h12$
' if n=3 then casa$="Casa I" :h$(n)=h1$
' if n=4 then casa$="Casa II" :h$(n)=h2$
' if n=5 then casa$="Casa III":h$(n)=h3$

ra1=armc+30*n
ra2=ra1
f=3/n
' do
for kk=1 to 20
print "ra1 ";ra1
' if abs(ra2-ra3)<0.01 then exit loop
print abs(ra2-ra3)
ra3=ra1
cose=-sin(ra1*ra#)*tan(e#*ra#)*tan(la*ra#)

seno=sqr(1-cose^2)
tang=seno/cose
deno=(atn(tang))/ra#
ra2=armc+deno/f
if ra2<0 then ra2=ra2+360
ra1=ra2
' loop
print "ra1 ";ra1
call s
next kk
print "ra2 ";ra2
cus(n)=(atn(tan(ra2*ra#)/cos(e#*ra#)))/ra#
call trans(cus(n),cus$(n))
call trans(MC,MC$)

'..................
d$= "\ \ ###,.## \ \\ \"
if n=1 then print "latitud: ";la$;" hora sideral: ";hs$
if n=1 then print "----------------------------------------------"
if n=1 then print using d$;"Medio Cielo ";MC;MC$;h10$
print using d$;casa$ ;cus(n);cus$(n);h$(n)
if n=5 then print ".............................................."
next

end
'''..................

sub s: while inkey$="":wend: end sub
'..................
sub trans(x,x$)
s=int(x/30)+1
ga=x-(s-1)*30:g=int(ga)
m=cint((ga-g)*60)
if m=60 then m=0:g=g+1
if g=30 then g=0:s=s+1
if s=13 then s=1
x$=str$(g)+"("+str$(s)+")"+str$(m)
end sub

Programa en Basic para calculo de cuspides mediante interpolacion (Version 2)

2011-09-04T16:39:00.003+02:00

Programa en Basic para calculo de cuspides mediante interpolacion

2011-09-04T16:34:00.000+02:00

cls
d1$=" ## ## ## ## "
d2$=" MC ## ## ## "
pi#=4*atn(1)
a#=2*pi#/360
ce#=(23+27/60)*a# ' constante ecliptica en radianes

input "Horasideral ",hs$
'hs$="163913"
hs=val(left$(hs$,2))+(val(mid$(hs$,3,2)))/60+(val(right$(hs$,2)))/3600
input "latitud ",la$
la=val(left$(la$,2))+(val(mid$(la$,3,2)))/60+(val(right$(la$,2)))/3600
if la>60 then la=60

x=hs*360/24
print "Hora sideral en grados ";x

' calculo del medio cielo
aa=x+2.5*2*sin(x*a#)*(cos(x*a#))
a0=aa
gosub signos
print using d2$;s+1;g;m

'ascendente a 0 grados
'primera correccion
x0=x-2.5*2*sin(x*a#)*(cos(x*a#))
aa=x0
g0=x0+3*30
gosub signos
print " Ascendnte"
print "Latitud Signo Grado Minuto"

print using d1$;0;s+4;g;m

open "casadato.dat" for input as #1
input #1, nd
for n=1 to nd
' g0a=g0p
g2a=g2p
g4a=g4p
g5a=g5p
g6a=g6p
ha =hp 'hacemos igual la hora anterior a la posterior leida antes

input #1, ht$ 'hora sideral de la tabla de datos
ht=val(left$(ht$,2))+(val(mid$(ht$,3,2)))/60+(val(right$(ht$,2)))/3600
input #1, g0$ ' grado 0
' g0=val((left$(g0$,2))-1)*30+(val(mid$(g0$,3,2)))+(val(right$(g0$,2)))/60
input #1, g2$ ' grado 20
g2=(val(left$(g2$,2))-1)*30+(val(mid$(g2$,3,2)))+(val(right$(g2$,2)))/60
input #1, g4$ ' grado 40
g4=(val(left$(g4$,2))-1)*30+(val(mid$(g4$,3,2)))+(val(right$(g4$,2)))/60
input #1, g5$ ' grado 50
g5=(val(left$(g5$,2))-1)*30+(val(mid$(g5$,3,2)))+(val(right$(g5$,2)))/60
input #1, g6$ ' grado 60
g6=(val(left$(g6$,2))-1)*30+(val(mid$(g6$,3,2)))+(val(right$(g6$,2)))/60

' g0p=g0
' if g0p<g0a then g0p=g0p+12*30
g2p=g2
if g2p<g2a then g2p=g2p+12*30
g4p=g4
if g4p<g4a then g4p=g4p+12*30
g5p=g5
if g5p<g5a then g5p=g5p+12*30
g6p=g6
if g6p<g6a then g6p=g6p+12*30

hp=ht ' hacemos igual la hora posterior a la hora leida en ese momento
' si la hora sidral esta comprendida entre dos consecutivas de la tabla
' podemos hacer el calculo de interpolacion con esos dos datos.

if hs>=ha and hs<=hp then
' call inter(g0a,g0p,hp,ha,hs)
' a0=aa
aa=g0
a0=g0
gosub signos
print "Ascendente a 0 grados ";aa$

call inter(g2a,g2p,hp,ha,hs)
a2=aa
print "Ascendente a 20 grados ";aa$

call inter(g4a,g4p,hp,ha,hs)
a4=aa
print "Ascendente a 40 grados ";aa$

call inter(g5a,g5p,hp,ha,hs)
a5=aa
print "Ascendente a 50 grados ";aa$

call inter(g6a,g6p,hp,ha,hs)
a6=aa
print "Ascendente a 60 grados ";aa$

if la<=20 then call inter(a0,a2,20,0 ,la)
if la>20 and la<=40 then call inter(a2,a4,40,20,la)
if la>40 and la<=50 then call inter(a4,a5,50,40,la)
if la>50 then call inter(a5,a6,60,50,la)

al=aa
print "-------------------------------------"
print "ASCENDENTE A";int(la);"GRADOS ";aa$

goto fuera
end if
next n
close #1
fuera:

end

sub inter(ga,gp,p,a,s)
shared gs,aa$,aa
i=p-a
if abs(gp)<=abs(ga) then aa=((ga-gp)*(p-s))/i+gp
if abs(gp)> abs(ga) then aa=((gp-ga)*(s-a))/i+ga
if aa>12*30 then aa=aa-12*30
call s
gosub signos
end sub

signos:
s=int(aa/30) ' signo
g1=(aa-s*30)
g=int(g1) ' grado
m=int((g1-g)*60) ' minuto
aa$=str$(s+1)+str$(g)+str$(m)
return

sub s
while inkey$="":wend
end sub

Ascensión Recta y Longitud de un Planeta

2011-08-30T20:20:00.000+02:00

Si en lugar de utilizar la tabla de efemérides, utilizamos las tablas que nos ofrecen los observatorios, por ejemplo esta página del observatorio astronómico de Madrid, vemos que en lugar de las longitudes, nos viene la "Ascención Recta" . La Ascensión Recta es un ángulo similar a la longitud, pero medido en el Ecuador en lugar de la Eclíptica. Como entre estos dos círculos existe una relación constante, podemos comprobar que ambos datos coinciden.

Ejemplo: 14 de Junio de 1997.

Las tablas de efemérides nos dan una longitud para el Sol a las 0 horas de Greenwich de 22 Gem 59.5 (Esto quiere decir, 22º y 59.5' deGéminis)

Y las tablas del observatorio astronómico de Madrid nos da una Ascensión Recta de 5h 29m 27.4s., comprobaremos que ambos datos son coincidentes.

(Para hacer este ejercicio necesitaremos una calculadora con funciones trigonométricas, si no se dispone de una serian necesario unas tablas trigonométricas).

Pasemos los dos ángulos a una medida homogenea, por ejemplo en grados.

a) En primer lugar pasaremos la longitud a grados:

LO = 22 Gem 59.5 = 60 + 22 + 59.5 / 60 = 82.992 º

Aclaración: Hemos puesto en primer lugar 60º, que es la distancia del Punto 0 de Géminis con el punto 0 de Aries, pues son dos signos de 30º cada uno)

b) Y ahora pasaremos la Ascensión Recta también a grados:

5h 29m 27.4s = 5 + 29/60 + 27.4/3600 = 5 + 0.483 + 0.008 h = 5.491 h

y pasemos las horas a grados:

24 h -------- 360 º

5.491 h ---- AR º

AR = 5.491 x 360 / 24 = 82.365 g = 82g 22 m

c) Ahora calculemos la Ascensión Recta a partir de la longitud eclíptica que nos da las tablas de efemérides.

llamemos:

AR = Ascensión Recta.

LO = Longitud.

e = Inclinación de la eclíptica

En primer lugar calculemos la inclinación de la eclíptica para el 14 de junio de 1997:

La formula mas exacta para este calculo es esta a partir de la inclinación del 1 de enero de 1997:

e = (23.439681 - 0.00000036 x d)

(en esta ecuación d es el numero de días desde el 1 de enero de 1997. Si quisiéramos calcular la inclinación de la eclíptica para una fecha anterior al 1-1-97 pondremos un signo + en lugar de restar)

e = (23.439681 - 0.00000036 * 165)º = 23.439622.

La relación entre la Ascensión Recta y la longitud es la siguiente:

AR = arctan [ tan (LO) x cos (e) ]

(Esto se lee de la forma siguiente: La Ascensión Recta es igual al arco cuya tangente (arctan) es el producto de la tangente de la longitud por el coseno de la inclinación de la eclíptica)

cos(e)= 0.9174797

tan(LO)= 8.135

y el producto de estas dos cantidades es = 7.4637

por lo que el arco cuya tangente es el numero anterior es : 82.369 g.

que equivale a 82g 22m,

que es el mismo valor que nos daba la tabla del observatorio astronómico.

(se produce una pequeña desviación del resultado en el tercer decimal como consecuencia de las operaciones trigonométricas.

Si en lugar de utilizar la tabla de efemérides, utilizamos las tablas que nos ofreces los observatorios. P.ej: Esta página del observatorio astronomico de Madrid.

Vemos que en lugar de las longitudes, nos viene la "AscensiónRecta" . La Ascensión Recta es un ángulo similar a la longitud, pero medido en el Ecuador en lugar de la Eclíptica. Como entre estos dos círculos existe una relación, podemos comprobar que ambos datos coinciden.

Ejemplo: 14 de Junio de 1997. Las efemérides nos dan una longitud para el Sol a las 0 horas de Greenwich de 22 Gem 59.5 (Esto quiere decir, 22º y 59.5' de Géminis)

Y el observatorio astronómico nos da una Ascensión Recta de 5h 29m 27.4s., comprobar que ambos datos son coincidentes.

(Para hacer este ejercicio necesitaremos una calculadora con funciones trigonométricas, si no se dispone de una serian necesario unas tablas trigonométricas).

Pasemos los dos ángulos a una medida homogenea, por ejemplo en grados.

a) En primer lugar pasaremos la longitud a grados:

LO = 22 Gem 59.5 = 60 + 22 + 59.5 / 60 = 82.992 º
Aclaración: Hemos puesto en primer lugar 60º, que es la distancia del Punto 0 de Géminis con el punto 0 de Aries, pues son dos signos de 30º cada uno)

b)Y ahora pasaremos la Ascensión Recta también a grados:

5h 29m 27.4s = 5 + 29/60 + 27.4/3600 = 5 + 0.483 + 0.008 h = 5.491 h

y pasemos las horas a grados:

24 h -------- 360 º
5.491 h ---- AR º

AR = 5.491 x 360 / 24 = 82.365 g = 82g 22 m

c) Ahora calculemos la Ascensión Recta a partir de la longitudeclíptica que nos da las tablas de efemérides.

llamemos:
AR = Ascensión Recta.
LO = Longitud.
e = Inclinación de la eclíptica

En primer lugar calculemos la inclinación de la eclíptica para el 14 de junio de 1997:

La formula mas exacta para este calculo es esta a partir de lainclinación del 1 de enero de 1997:
e = (23.439681 - 0.00000036 x d)
(en esta ecuación d es el numero de días desde el 1 de enero de 1997. Si quisieramos calcular la inclinación de la eclíptica para una fecha anterior al 1-1-97 pondremos un signo + en lugar de restar)

e = (23.439681 - 0.00000036 * 165)º = 23.439622.

La relación entre la Ascensión Recta y la longitud es la siguiente:

AR = arctan [ tan (LO) x cos (e) ]
(Esto se lee de la forma siguiente: La Ascensión Recta es igual al arco cuya tangente (arctan) es el producto de la tangente de la longitud por el coseno de la inclinación de la eclíptica)

cos(e)= 0.9174797
tan(LO)= 8.135

y el producto de estas dos cantidades es = 7.4637

por lo que el arco cuya tangente es el numero anterior es : 82.369 g.
que equivale a 82g 22m,
que es el mismo valor que nos daba la tabla del observatorio astronómico.

(se produce una pequeña desviación del resultado en el tercer decimal como consecuencia de las operaciones trigonométricas.

+++++++++++++++++++++

Calculo de la cúspide del Medio Cielo diréctamente sin ayuda de tablas:

Como complemento de lo anterior explicare la forma de calcular matemáticamente la cúspide de la casa 10, llamada Medio Cielo, astronómicamente no es sino la hora sideral expresada en coordenadas eclípticas, es decir que si tenemos la hora sideral en horas y minutos, bastara con pasarla a grados (24 horas equivalen a 360 grados) y tendremos esa cúspide en coordenada ecuatorial es decir la Ascensión Recta (RA) , si pasamos esa coordenada a eclíptica tendremos la Longitud del Medio Cielo:

En nuestro ejemplo tenemos una Hora Sideral de13h 7 m, y por tanto un RA de:

13h 7.56m = 13 + 7.56/ 60 = 13.126 horas

13.126 x 360 /24 = 196.89 grados

196.89 = 196 g + 0.88x 60 m = 196g 54 m

RA = 196g 54 m

Para hallar su longitud tendremos en cuenta la siguiente relación astronómica:

tan(long)= tan(RA) / cos(O)

donde O es el ángulo de inclinación de la eclíptica : O = 23g 26m 21s
O = 23+26/60+21/360= 23+0.4333+0.0058=23.4393
y como cos(O) = 0.9175
y en este caso RA = 196 g.54m por lo que tan(RA)=0.304
y 0.304 / 0.9175 = 0.331

long = arc.tan(tan(RA) / 0.9175 ) =arc.tan(0.331)=18.31

,por lo que tendremos que la longitud del MC sera 18 grados y 0.31 x 60 =19 m , y la longitud del MC sera:

MC = 18 Li 19

Como puede observarse para realizar el cálculo ha sido necesario introducir la corrección de aceleración que es muy pequeña pero afecta mucho al resultado. Centésimas de minuto en esa corrección alteran varios minutos en el resultado de la longitud del Medio Cielo. Como nadie sabe con esa exactitud su hora de nacimiento, resulta totalmente inapropiado en un calculo manual. Esta corrección tan solo se debe utilizar en un programa informático ya que este lo realiza automáticamente.

Ejemplo 2 de calculo del Ascendente (Hemisferio Norte)

2011-07-30T14:54:00.001+02:00

Ejemplo de Calculo del Ascendente en una ciudad de Hemisferio Norte. Tambien se calcula el Medio Cielo utilizando la curva central para cualquier latitud pero se leen los resultados en las escalas inferior y lateral izquierda.

Para calcular el Ascendente se procede la la forma siguiente:

1.- Vamos a la escala superior y se busca la hora sideral de nacimiento que hemos conseguido antes en unas efemerides o bien con la tabla auxiliar que tenemos en esta web.

2.- Bajamos una linea hasta la curva representativa de la latitud de la ciudad. Si no existe la curva se buscara un punto intermedio entre dos curvas que guarde cierta proporcionalidad con el valor de la latitud.

3.- Seguimos con una linea horizontal hacia la derecha y leemos directamente el valor del ascendente. Tendremos en cuenta que el ascendente se leerá en la escala que va unida con la escala de horas siderales.

En la tabla que tenemos bajo el grafico esta calculado analiticamente lo mismo que se se busco en el grafico.

Para calcular el Medio Cielo se procede la la forma siguiente:

1.- Vamos a la escala inferior y se busca la hora sideral de nacimiento que hemos conseguido antes en unas efemerides o bien con la tabla auxiliar que tenemos en esta web.

2.- Subimos una linea hasta la curva central. (Que es casi recta)

3.- Seguimos con una linea horizontal hacia la izquierda y leemos directamente el valor del Medio Cielo. Tendremos en cuenta que el Medio Cielo se leerá en la escala que va unida con la escala de horas siderales.

Ejemplo de calculo del Ascendente

2011-07-30T14:44:00.000+02:00

Ejemplo de cálculo del Ascendente

Ejemplo de Cálculo del Ascendente en una ciudad del Hemisferio Norte. También se calcula el Medio Cielo utilizando la curva central para cualquier latitud pero se leen los resultados en las escalas inferior y lateral izquierda.

Pulsando sobre el gráfico se agranda para mayor claridad. Puede imprimirse el gráfico para utilizarlo con mas comodidad.

Para calcular tu Ascendente procede en la forma siguiente:

1.- Ve a la escala superior y busca la Hora Sideral de nacimiento que has conseguido antes en unas efemérides o bien con la tabla auxiliar que tenemos en esta web.

2.- Baja una linea recta vertical hasta la curva representativa de la latitud de la ciudad en que naciste. Si no existe la curva busca un punto intermedio entre dos curvas que guarde cierta proporcionalidad con el valor de la latitud.

3.- Seguimos con una linea horizontal hacia la derecha y leemos diréctamente el valor del Ascendente. Tendremos en cuenta que el Ascendente se leerá en la escala que va unida con la escala de Horas Siderales.

En la tabla que tenemos bajo el grafico esta calculado analíticamente lo mismo que se se buscó en el grafico.

Para calcular tu Medio Cielo procede en la forma siguiente:

1.- Ve a la escala inferior y se busca la Hora Sideral de nacimiento que hemos conseguido antes en unas efemerides o bien con la tabla auxiliar que tenemos en esta web.

2.- Subimos una linea hasta la curva central. (Que es casi recta)

Tabla de Casas de Regiomontano (Programa en Basic)

2011-07-30T14:40:00.000+02:00

'regio.bas

'Este programa produce una Tabla de Casas por el sistema de Regiomontano que se graba en un fichero con el nombre Regiomon.doc. Este fichero puede abrirse con cualquier procesador de textos.

'Las instrucciones utilizadas en el codigo son de Basic y pueden ser facilmente traducibles a cualquier lenguaje de programacion.

'Se producen las longitudes de las Casas 10, 11, 12, 1, 2 y 3, las otras no son necesarias pues basta con añadirles 6 signos.

cls
pi=3.14159296535
cdr= pi/180
gecl=23:mecl=27:secl=15
obl=(gecl+mecl/60+secl/3600)*cdr

'print tan(105/cdr)
'input rrr

dim z(12)
dim p(12)
dim v(12)
dim x(12)
dim y(12)
dim l(12)
dim sic(12)
dim sic$(12)
dim gr(12)
dim mi(12)
dim po$(12)
open "regiomon.doc" for output as #1
for si=1 to 12
for grado=0 to 29
g=grado
gr=int(g)
l=(si-1)*30+g                           'l=longitud del medio cielo
ss=si:gosub signos:si$=ss$     'signo del MC en letra
ra=atn(tan(l*cdr)*cos(obl))/cdr 'ascension recta en grados
if ra<0 then ra=ra+360
if ra>360 then ra=ra-360

if l=0   then ra=0.0001
if l=90 then ra=90
if l>90 and l<270 then ra=ra+180
if l=180 then ra=180
if l=270 then ra=270
if l>270 then ra=ra+360
if l>360 then ra=ra-360
if l=360 then ra=360

rag=int(ra):ram=(ra-rag)*60 ' ascension recta en grados y minutos
st=ra*24/360                         ' st=hora sideral

' h.s. en horas,minutos y segundos
sth=int(st)
stmm=(st-sth)*60 :stm=int(stmm)
stss=(stmm-stm)*60:sts=cint(stss)

'for la=22 to 37
la=22
' como ejemplo se crean las Casas para la latitud 22N, basta con sustituir este valor para
' crear la tabla de otra latitud. Tambien puede introducirse un bucle para varias latitudes
' ver la linea anterior donde se crea un bucle de 22 a 37 grados de latitud.

lat=la*cdr

for n=1 to 12 ' bucle de las 12 casas

   z(n)=60+n*30 : if z(n)>360 or z(n)=360 then z(n)=z(n)-360
   p(n)=atn(tan(z(n)*cdr)*cos(lat))/cdr
   sz=sin(z(n))
   cz=cos(z(n))
   if sz>0 and cz<0 then l(n)=l(n)+180        ' cuadrante 2

   if sz<0 and cz<0 then l(n)=l(n)+180        ' cuadrante 3
   if sz<0 and cz>0 then l(n)=l(n)+360        ' cuadrante 4

   if n<7 then then p(n)=p(n)-180
   if p(n)<0 then p(n)=p(n)+360

v(n)=sin(p(n)*cdr) * sin(lat) * sin(obl)

x(n)=cos(ra*cdr+z(n)*cdr) * cos(lat) * cos(obl)- v(n)
y(n)=sin(ra*cdr+z(n)*cdr) * cos(lat)

t=y(n)/(x(n))
l(n)=atn(t)/cdr

if y(n)>0 and x(n)<0 then l(n)=l(n)+180
if y(n)<0 and x(n)<0 then l(n)=l(n)+180
if y(n)<0 and x(n)>0 then l(n)=l(n)+360

if l(n)<0 then l(n)=l(n)+360
if l(n)>360 then l(n)=l(n)-360

sic(n)=int(l(n)/30)+1 'signo de la Casa
ss=sic(n):gosub signos:sic$(n)=ss$
gc=l(n)-30*(sic(n)-1):gr(n)=int(gc)
m=(gc-gr(n))*60:mi(n)=cint(m)

next n

if g=0 then
   print #1, "Tabla de Casas de Regiomontano. Alba de Hermes"
   print #1, " C10 AR   HS"
   print #1, " s g gr mi h m s Lo Casa11 Casa12 Casa 1 Casa 2 Casa 3"
end if
print #1, using "\ \ ## ### ##,.# ## ## ## ## ##\\## ##\\## ##\\## ##\\## ##\\##";_
si$;g;rag;ram;sth;stm;sts;la;_
gr(11);sic$(11);mi(11);gr(12);sic$(12);mi(12);gr(1);_
sic$(1);mi(1);gr(2);sic$(2);mi(2);gr(3);sic$(3);mi(3)
' next la
next grado
next si
close #1
end

signos:
select case ss
case 1: ss$="ar"
case 2: ss$="ta"
case 3: ss$="ge"
case 4: ss$="cn"
case 5: ss$="le"
case 6: ss$="vi"
case 7: ss$="li"
case 8: ss$="es"
case 9: ss$="sa"
case 10:ss$="cp"
case 11:ss$="ac"
case 12:ss$="pi"
end select
return
sub s

while inkey$=""
wend
end sub

efemerides planetarias. programa en basic v2

2011-06-17T19:17:00.001+02:00

' EFEMERIDES PLANETARIAS
cls

dim b$(32)
c%=360
pi#=4*atn(1)

cd#=360/365.2422 ' grados recorridos en un dia medio
ra#=180/pi# ' factor para transformar radianes en grados multiplicando
' grados en radianes dividiendo
rr#=360/pi# ' el doble de ra#

for f=1 to 8
read a$(f),t(f),w(f),p(f),e(f),s(f),i(f),no(f)
next f

principio:
'input "Fecha aaaa,mm,dd", y,m,d
y=1984:m=7:d=12

' calculo del periodo de tiempo en dias desde el 1 de enero de 1980
if m<3 then k=m else k=0
j=int(365.25*(y-k))+int(30.4375*(m+1+12*k))+d-int(y/100)+_
int(int((y/100)/4))+1720996.5

ed%=j-2444238.5+1 ' periodo de dias desde el dia con datos iniciales
print ed%
' POSICION DE LA TIERRA
'ed%=1655
we#=98.83354 ' longitud del periodo
' 102.596403= longitud del perihelio
ee#=0.016718 ' excentricidad de la tierra

' cd#*ed% = grados que recorre en los dias del periodo considerado
' en su orbita, referido al
' 1.00004 = periodo orbital de la tierra en a¤os sidereos
z#=cd#*ed%/1.00004 :gosub rango:ne#=z#

'le sumamos la longitud del punto de partida y le restamos la del perihelio.
' y tenemos la distancia desde el punto a calcular hasta el perihelio
z#=ne#+we#-102.596403 :gosub rango:me#=z#

z#=ne#+rr#*ee#*sin(me#/ra#)+we#:gosub rango:le#=z#
ve#=le#-102.596403 ' anomalia verdadera
re#=(1-ee#^2)/(1+ee#*cos(ve#/ra#)) ' radio vector

print " EFEMERIDES PLANETARIAS"
print
print " Planeta Ascension Recta Declinacion"
d1$= "\ \ ### \\ ## \\ ## \\ ### \\ ## \\ ## \\ ## ## ##"

for n=1 to 8
z#=cd#*ed%/t(n) :gosub rango:np#=z#
z#=np#+w(n)-p(n) :gosub rango:mp#=z#
z#=np#+rr#*e(n)*sin(mp#/ra#)+w(n):gosub posit:lp#=z#
vp#=lp#-p(n)
rp#=s(n)*(1-e(n)^2)/(1+e(n)*cos(vp#/ra#))

pps=sin((lp#-no(n))/ra#)*sin(i(n)/ra#)
ppc=sqr(1-pps^2)
psi=ra#*atn(pps/ppc)
yx=lp#-no(n)
y=sin(yx/ra#)*cos(i(n)/ra#)
x#=cos(yx/ra#)
ct=ra#*atn(y/x#)
q=ct: gosub cuadrante: ct=q
l1=ct+no(n)
r1=rp#*cos(psi/ra#)

if n>2 then goto 430 ' Planetas interiores Mercurio y venus
ll=le#-l1
a=ra#*atn((r1*sin(ll/ra#))/(re#-r1*cos(ll/ra#)))
z#=180+le#+a:gosub posit:lam=z#
goto 470

430 ' planetas exteriores jupiter - pluton
ll=l1-le#
' longitud ecliptica
z#=l1+ra#*atn((re#*sin(ll/ra#)/(r1-re#*cos(ll/ra#)))):gosub posit:lam=z#
470 ' latitud ecliptica
bet=ra#*atn(r1*tan(psi/ra#)*sin((lam-l1)/ra#)/(re#*sin((l1-le#)/ra#)))

' coordenadas ecliptico-ecuatoriales
la=lam/ra#
ec=23.441884/ra#
be=bet/ra#
y=sin(la)*cos(ec)-tan(be)*sin(ec)
x#=cos(la)
rh#=ra#*atn((sin(la)*cos(ec)-tan(be)*sin(ec))/(cos(la)))

' va al cuadrante
q=rh# : gosub cuadrante: rh#=q
' formato en signos, grados y minutos
si=int(rh#/30)+1
grr=rh#-(si-1)*30
gr=int(grr)
mi=(grr-gr)*60
' ascension recta en formato en horas, m. y s.
rx=rh#/15
ry=int(rx)
mx=60*(rx-ry)
my=int(mx)
sx=cint(60*(mx-my))
' declinacion en formato en grados, m. y s.
pps#=sin(be)*cos(ec)+cos(be)*sin(ec)*sin(la)
ppc#=sqr(1-pps#^2)
dc=ra#*atn(pps#/ppc#)
sig=1:if dc<0 then sig=-1
dc1=abs(dc)
dh=int(dc1)
dm=60*(dc1-dh)
dmi=int(dm)
ds=cint(60*(dm-dmi))
dh=sig*dh
print using d1$;a$(n);ry;"h";my;"m";sx;"s";dh;"§";dmi;"m";ds;"s";si;gr;mi
next n

' calculo de longitud y latitud
for n=1 to 8

next

' dibuja los planetas en un atlas estelar
print "t para terminar , s para otra fecha"
a$=input$(1)
cls
if a$="t" then end else goto principio

end

' rango de 0-360 grados
rango:
z#=360*(z#/360-int(z#/360))
return

' numeros positivos
posit:
if z#>360 then z#=z#-360
if z#<0 then z#=z#+360
return

cuadrante:
if x#<0 then q=q+180:return
if y>0 then return
q=q+360
return

sub s
print "para seguir pulsa una tecla"
while inkey$="" :wend
end sub

' (1) (2) (3) (4) (5) (6) (7)
' t we# p ee# s i n
' periodo long. peri. l.perihelio excentr. semie ma. inclinac. nodo asc.
data "Mercu", 0.24085,231.2973 , 77.1442128,.2056306,0.3870986,7.0043579 ,48.0941733
data "Venus", 0.61521,355.73352, 131.2895792,.0067826,0.7233316,3.394435 ,76.4997524
data "Marte", 1.88089,126.30783, 335.6908166,.0933865,1.5236883,1.8498011 ,49.4032001
data "Jupit", 11.86224,146.966365, 14.0095493,.0484658,5.202561 ,1.3041819 ,100.2520175
data "Satur", 29.45771,165.322242, 92.6653974,.0556155,9.554747 ,2.4893741 ,113.4888341
data "Urano", 84.01247,228.070855, 172.7363288,.0463232,19.21814 ,0.7729895 ,73.8768642
data "Neptu",164.79558,260.3578998, 47.8672148,.0090021,30.10957 , 1.7716017,131.5606494
data "Pluto",250.9 ,209.439 ,222.972 ,.25387 ,39.78459 ,17.137 ,109.941
' tierra 1.00004 98.83354 102.596403 .016718 1 0 0
' (1) periodo orbital en a¤os sidereos
' (2) longitud el 1 de enero de 1980
' (3) longitud del perihelio
' (4) excentricidad de la orbita
' (5) valor medio del semieje mayor en relacion con el valor de la tierra
' (6) declinacion de la orbita a la ecliptica
' (7) nodo ascendente

Programa en basic de calculo de efemerides planetarias. (poca aproximacion)

2011-06-04T14:52:00.001+02:00

' EFE.BAS
cls

dim b$(32)
c%=360
pi#=4*atn(1)

cd#=360/365.2422 ' grados recorridos en un dia medio
ra#=180/pi# ' factor para transformar radianes en grados
rr#=360/pi# ' el doble de ra#

for f=1 to 8
read a$(f),t(f),w(f),p(f),e(f),s(f),i(f),no(f)
next f

principio:
'input "Fecha aaaa,mm,dd", y,m,d
y=1984:m=7:d=12
'y=1980:m=1:d=1

' calculo del periodo de tiempo en dias desde el 1 de enero de 1980
if m<3 then k=m else k=0
j=int(365.25*(y-k))+int(30.4375*(m+1+12*k))+d-int(y/100)+_
int(int((y/100)/4))+1720996.5

ed%=j-2444238.5+1 ' periodo de dias desde el dia con datos iniciales
print ed%
' POSICION DE LA TIE
'ed%=1655
we#=98.83354 ' longitud del periodo
' 102.596403= longitud del perihelio
ee#=0.016718 ' excentricidad de la tie

' cd#*ed% = grados que recorre en los dias del periodo considerado
' en su orbita, referido al
' 1.00004 = periodo orbital de la tie en a¤os sidereos
z#=cd#*ed%/1.00004 :gosub rango:ne#=z#

'le sumamos la longitud del punto de partida y le restamos la del perihelio.
' y tenemos la distancia desde el punto a calcular hasta el perihelio
z#=ne#+we#-102.596403 :gosub rango:me#=z#

z#=ne#+rr#*ee#*sin(me#/ra#)+we#:gosub rango:le#=z#
ve#=le#-102.596403 ' anomalia verdadera
re#=(1-ee#^2)/(1+ee#*cos(ve#/ra#)) ' radio vector

print " EFEMERIDES "
print
ca$= " p Ascension Recta Declinacion"
d1$= "\ \ ### \\ ## \\ ## \\ ### \\ ## \\ ## \\ ## ## ##"
ca1$= " p Latitud Longitud"
d2$= "\ \ ### \\ ## \\ ## ## ##"
print ca1$
for n=1 to 8
z#=cd#*ed%/t(n) :gosub rango:np#=z#
z#=np#+w(n)-p(n) :gosub rango:mp#=z#
z#=np#+rr#*e(n)*sin(mp#/ra#)+w(n):gosub posit:lp#=z#
vp#=lp#-p(n)
rp#=s(n)*(1-e(n)^2)/(1+e(n)*cos(vp#/ra#))

pps=sin((lp#-no(n))/ra#)*sin(i(n)/ra#)
ppc=sqr(1-pps^2)
psi=ra#*atn(pps/ppc)
yx=lp#-no(n)
y=sin(yx/ra#)*cos(i(n)/ra#)
x#=cos(yx/ra#)
ct=ra#*atn(y/x#)
q=ct: gosub cuadrante: ct=q
l1=ct+no(n)
r1=rp#*cos(psi/ra#)

' longitud ecliptica
if n>2 then goto 430 ' p. interiores mer y ven
ll=le#-l1
a=ra#*atn((r1*sin(ll/ra#))/(re#-r1*cos(ll/ra#)))
z#=180+le#+a:gosub posit:lam=z#
goto 470
430 ' p. exteriores jup - plu
ll=l1-le#
z#=l1+ra#*atn((re#*sin(ll/ra#)/(r1-re#*cos(ll/ra#)))):gosub posit:lam=z#

470 ' latitud ecliptica
bet=ra#*atn(r1*tan(psi/ra#)*sin((lam-l1)/ra#)/(re#*sin((l1-le#)/ra#)))

' q=lam : gosub cuadrante: lam=q
' formato en signos, grados y minutos
si=int(lam/30)+1
grr=lam-(si-1)*30
gr=int(grr)
mi=int((grr-gr)*60)
' declinacion en formato g.m.s.
sig=1:if bet<0 then sig=-1
dc1=abs(bet)
dh =int(dc1)
dm=int(60*(dc1-dh))
dh=sig*dh
' print ca1$

print using d2$;a$(n);dh;"§";dm;"m";si;gr;mi

goto fin

' coordenadas ecliptico-ecuatoriales
la=lam/ra#
ec=23.441884/ra#
be=bet/ra#
y=sin(la)*cos(ec)-tan(be)*sin(ec)
x#=cos(la)
rh#=ra#*atn((sin(la)*cos(ec)-tan(be)*sin(ec))/(cos(la)))

' va al cuadrante
q=rh# : gosub cuadrante: rh#=q
' formato en signos, grados y minutos
si=int(rh#/30)+1
grr=rh#-(si-1)*30
gr=int(grr)
mi=(grr-gr)*60
' ascension recta en formato h.m.s.
rx=rh#/15
ry=int(rx)
mx=60*(rx-ry)
my=int(mx)
sx=cint(60*(mx-my))
' declinacion en formato g.m.s.
pps#=sin(be)*cos(ec)+cos(be)*sin(ec)*sin(la)
ppc#=sqr(1-pps#^2)
dc=ra#*atn(pps#/ppc#)
sig=1:if dc<0 then sig=-1
dc1=abs(dc)
dh=int(dc1)
dm=60*(dc1-dh)
dmi=int(dm)
ds=cint(60*(dm-dmi))
dh=sig*dh
' print ca$
' print using d1$;a$(n);ry;"h";my;"m";sx;"s";dh;"§";dmi;"m";ds;"s";si;gr;mi
fin:
next n

' calculo de longitud y latitud
for n=1 to 8

next

' dibuja los planetas en un atlas estelar
print "t para terminar , s para otra fecha"
a$=input$(1)
cls
if a$="t" then end else goto principio

end

' rango de 0-360 grados
rango:
z#=360*(z#/360-int(z#/360))
return

' numeros positivos
posit:
if z#>360 then z#=z#-360
if z#<0 then z#=z#+360
return

cuadrante:
if x#<0 then q=q+180:return
if y>0 then return
q=q+360
return

sub s
print "para seguir pulsa una tecla"
while inkey$="" :wend
end sub

' (1) (2) (3) (4) (5) (6) (7)
' t we# p ee# s i n
' periodo long. peri. l.perihelio excentr. semie ma. inclinac. nodo asc.
data "Mercu", 0.24085,231.2973 , 77.1442128,.2056306,0.3870986,7.0043579 ,48.0941733
data "Venus", 0.61521,355.73352, 131.2895792,.0067826,0.7233316,3.394435 ,76.4997524
data "Marte", 1.88089,126.30783, 335.6908166,.0933865,1.5236883,1.8498011 ,49.4032001
data "Jupit", 11.86224,146.966365, 14.0095493,.0484658,5.202561 ,1.3041819 ,100.2520175
data "Satur", 29.45771,165.322242, 92.6653974,.0556155,9.554747 ,2.4893741 ,113.4888341
data "Urano", 84.01247,228.070855, 172.7363288,.0463232,19.21814 ,0.7729895 ,73.8768642
data "Neptu",164.79558,260.3578998, 47.8672148,.0090021,30.10957 , 1.7716017,131.5606494
data "Pluto",250.9 ,209.439 ,222.972 ,.25387 ,39.78459 ,17.137 ,109.941
' tierra 1.00004 98.83354 102.596403 .016718 1 0 0

' (1) periodo orbital en a¤os sidereos
' (2) longitud el 1 de enero de 1980
' (3) longitud del perihelio
' (4) excentricidad de la orbita
' (5) valor medio del semieje mayor en relacion con el valor de la tierra
' (6) declinacion de la orbita a la ecliptica
' (7) nodo ascendente