Newsgroup: it.scienza.astronomia.amatoriale
Titolo: Le meraviglie del Solex e l'eclisse di Luna

Allego in questa mia per gl'interessati al calcolo: Cristiano, Daniele 
Gasparri ed eventuali, i risultati che si possono ottenere con l'utilizzo 
del precisissimo Solex.
In pratica con il solo uso della barra spaziatrice e col settaggio del 
passo temporale (step scaling --> tasto s) e dell'unità di tempo (time unit 
--> tasto u) è possibile spostarsi in unità da 100 fino a 0.01. Semplice 
per chi è nato col DOS, meno semplice per i mouse-dipendenti :-))
Il calcolo si riferisce alla congiunzione Sole<->Luna di 2007 Marzo 03
*********************************************************************
Teoria JPL/NASA DE200  (Solex 8.54)
Geocentric    Angular  Equatorial Coordinates  of Sun      
True Equinox of Date  (Corrected for Aberr. & L-time)
Time is TDT
  2007/03/03  23:01:50  22 57 16.045    - 6 41  5.65     0.99138619
Geocentric    Angular  Equatorial Coordinates  of Moon     
True Equinox of Date  (Corrected for Aberr. & L-time)
Time is TDT
  2007/03/03  23:01:50  10 57 16.043    + 7  0 44.83     0.40217686
###################################################################
Teoria JPL/NASA DE406 (Solex 9.03)
Geocentric    Angular  Equatorial Coordinates  of Sun      
True Equinox of Date  (Corrected for Aberr. & L-time)
Time is TDT
  2007/03/03  23:01:50  22 57 16.045    - 6 41  5.67     0.99138619
Geocentric    Angular  Equatorial Coordinates  of Moon     
True Equinox of Date  (Corrected for Aberr. & L-time)
Time is TDT
  2007/03/03  23:01:50  10 57 16.046    + 7  0 44.81     0.40217686
==============
Un ringraziamento all'Autore (prof. Aldo Vitagliano 
http://main.chemistry.unina.it/~alvitagl/
 )
che ci regala un generatore di effemeridi così accurato.
Eliseo
-- 
"Il n'est pire sourd que celui qui ne veut pas entendre."

 Teoria JPL/NASA DE200  (Solex 8.54)
Vorrei fare una considerazione puramente accademica (così, per
chiacchierare).
Ho notato varie volte che nomini le effemeridi DExxx usando "Teoria".
Secondo me è un termine fuorviante perché normalmente vengono così chiamate
le soluzioni analitiche tipo ELP 2000-82B, VSOPxx, TOP, ecc che introducono
nei loro calcoli termini periodici "aleatori" ovvero teorici (pur basandosi,
in qualche modo, su posizioni planetarie note).
Le effemeridi DExxx (o soluzioni numeriche) sono elenchi di coefficienti di
polinomi interpolatori che consentono di ottenere (quasi istantaneamente) il
valore numerico di un parametro planetario che è stato a suo tempo calcolato
mediante il lento processo d'integrazione numerica (l'integrazione delle
DE414, per esempio, ha richiesto 16 minuti per ogni secolo integrato).
Tenendo conto di un gran numero di corpi perturbatori, le DE/LE non
introducono calcoli arbitrari, per cui gli errori sono dovuti in massima
parte all'incertezza sugli elementi osculatori utilizzati all'inizio
dell'integrazione numerica.
Ciao
Cristiano 

 Cristiano
Ciao


 partenza e una di arrivo
Che vuoi dire? Con le DExxx leggi i coefficienti del polinomio validi per il 
"blocco" temporale che t'interessa e calcoli i corrispondenti valori 
cercati, il tutto molto più velocemente rispetto a quando si usano le 
soluzioni analitiche (composte da migliaia di seni e coseni).
Ciao
Cristiano 

 Cristiano 
Su Web ho trovato fra l'altro la Teoria INPOP06 fatta sulle DE414 
dell'Imcce di Paris 
http://www.imcce.fr/Equipes/ASD/inpop/inpop06_preprint.pdf
Le massime differenze si hanno per Marte (al massimo 4/100 di arco-secondo 
in longitudine da -400 anni a +200 anni dal J2000.0) rispetto alla DE405 
(p.5 di 17)
e poi qualcosa di meccanica orbitale in Fortran (occorre ovviamente il 
compilatore) nel sito
http://www.cdeagle.com/html/omfortran.html
io ho prelevato questo demo (fornisce le coordinate rettangolari)
DESCRIPTION :
    The ephemerides SLP96 are available in 8 sequential ASCII files slp96.
    Each one has the same structure and contains Chebychev ephemerides
    covering a century within the interval 1550-2350.
    Fortran programs allow the user to install the files and to compute
    coordinates of the bodies.
    A file of test values is also provided for checking the coordinates
    after the installation.
    Explanations about files SLP96 and sofware package are given in the
    file slp96.doc.


 rispetto al ref J2000 te ne do atto.
Credo che non hai chiaro il metodo attraverso il quale si utilizzano le 
DExxx.
Con le DExxx accedi direttamente ai valori di cui hai bisogno effettuando 
solo le moltiplicazioni e le addizioni che servono per calcolare il dato 
mediante un polinomio di Chebyshev che varia al variare del blocco temporale 
per il quale fai i calcoli.
Il tempo per ottenere il valore cercato è indipendente dall'istante per il 
quale effettui il calcolo.
Ciao
Cristiano 

 
Eh, Eh, ... Per forza! 
Solex non "usa" le DE405, le
produce! 
 :D
Beh, non esattamente ... Ma quasi ...
---------------------------
Posted via 
http://www.AstroForum.it
-- 
Posted via a free Usenet account from 
http://www.teranews.com

 Beh, non esattamente ... Ma quasi ...
Per fortuna l'hai precisato, altrimenti saresti sembrato un alchimista del 
medioevo che trasformava il piombo in l'oro. :-)
Pur essendo convinto della precisione piů che buona del tuo programma 
(verificabile, ad esempio, usando Horizons), francamente credo che ciň che 
dici č molto azzardato.
Ragionando per analogia con le PL404 (che io utilizzo quando non ho la 
copertura delle DE414), le PL406 sono un "condensato" delle DE406 
(consentono, cioč, di ridurre drasticamente la lunghezza dei file che 
contengono i coefficienti per i calcoli).
A loro volta, le DE406 sono una versione volutamente peggiorata delle DE405 
(anche in questo caso, le DE406 sono state create per occupare meno spazio).
Alla luce di questi fatti *oggettivi*, č evidente che la frase (ricomponendo 
il tuo post) "Solex, le DE405 quasi le produce" sembra alquanto azzardata, 
ma proprio tanto! :-)
Io direi semplicemente che Solex ha un grado di precisione piů che adeguato 
per i potenziali utenti di quel tipo di programmi. Perché voler dare a 
credere quello che non č?
Poi, se proprio dovessi darti un suggerimento, ti spronerei a cambiare il 
tipo d'interfaccia utente adattandola alle attuali GUI, molto piů funzionali 
e alle quali siamo ormai abituati da anni. Tu fa' come ti pare. :-)
Cristiano 

 
Grazie del suggerimento. E' una questione che ho più volte
considerato, ma alla mia età non è facile imparare una logica
diversa, e poi convertirci dodicimila righe di codice ... 
Tuttavia ti faccio notare che è solo una questione di pura apparenza.
Al 95% la sola differenza da una "normale" interfaccia
Windows è che i "bottoni" non sono "bottoni" ma
parole di testo. Se ci clicchi sopra col mouse è come avere cliccato
un bottone.


 
Stiamo parlando di due questioni divergenti (?). Io ti dico che dal tempo 
della ILE1954 le posizioni dei corpi celesti vengono pre-calcolate solo in 
alcuni poli, distanziati di giorni od ore a seconda della soluzione 
proposta e poi rappresentate all'interno di un polinomio di Chebyshev che 
dovrà scegliersi fra i tanti possibili. Nessuno si sognerebbe di eseguire 
OGNI VOLTA da capo l'integrazione numerica... Vabbé calcolatori potenti, ma 
con giudizio.
Poi le DE200 e in genere le altre DExxx forniscono solo un dato grezzo 
astrometrico che va affinato con le ultime teorie sulla precessione (per 
es. la IAU2000).
Ecco a che alludo e su cui sembra sussistere la nostra attuale 
incomprensione:
Compressed planetary and lunar ephemerides 
http://www.springerlink.com/content/px231t26j2h68646/
Abstract  A package of FORTRAN software has been developed which provides 
planetary and lunar positions, with respect to the solar system barycenter, 
for all times in the interval 1801–2049; positions agree to 1 
milliarcsecond with those generated by Jet Propulsion Laboratory 
Development Ephemeris 200 (DE200). The system consists of approximately 800 
kilobytes of ephemeris files and 40 kilobytes of programs, totalling 5% of 
the storage required by DE200. After removal of reference orbits, segments 
of DE200 positions were fitted by finite Chebyshev series of degree 40. The 
Chebyshev coefficients were rounded to integer multiples of a suitable unit 
and packed to form the ephemeris files.


 A me sa che ci stiamo perdendo in sofismi terminologici.
Infatti. E leggendo tutto il resto, mi sembra che la nostra chiacchierata 
stia prendendo una brutta piega; tanto più brutta se pensi che stiamo 
parlando di "aria fritta". :-)
Inoltre, mi sembra che non riusciamo proprio a capirci (è probabile che non 
sappia spiegarmi), per cui è inutile continuare.
Ciao e amici come prima
Cristiano 

 oppure freeware, fanno altrettanto ?
Qui mi fai pensare che tu credi che io stia mettendo in dubbio la serietà 
tua e del tuo programma del tuo programma, ma io non l'ho mai fatto, tengo a 
precisarlo.
Cristiano 

 Cristiano 
Ma almeno sei riuscito a trovare il motivo per cui ottieni scarti cosě 
elevati rispetto a Espenak, relativi all'eclisse di marzo? Guarda che mi 
sono accorto solo oggi della tua risposta a me indirizzata del 8/3 e la mia 
la trovi in questo ng.
Torna pure quando ti pare ma ricordati di portare con te gli esempi 
numerici ben dettagliati, cosě ci capiamo prima e meglio.
Ciao.


 
Se ti interessa davvero, ti posso rendere disponibili il codice
sorgente della versione base per DOS (il nucleo storico) che serve
per capire il principio e il metodo di funzionamento. 
L'integratore è un Bulirsh-Stoer polinomiale (modificato) di 16°
ordine. Una Ferrari truccata ... I corpi in considerazione sono
normalmente tutti i pianeti, la Luna, Cerere, Pallade e Vesta. 
Si può avere anche un Sistema Solare "ridotto" senza Plutone
e senza asteroidi, oppure anche senza la Luna.
Il principio generale è che l'integratore si va a leggere dalla
libreria le condizioni iniziali più vicine all'epoca richiesta, e poi
parte da lì e raggiunge l'epoca richiesta e gli eventuali step
successivi. La libreria a sua volta è stata creata con una versione
più precisa dello stesso integratore (che lavora a 19 decimali
anzichè a 16), a partire da un unico set di condizioni iniziali
ottimizzato con i minimi quadrati, usando l'intero arco delle
DE405/406 (o DE409 nella opzione corrispondente) come se fossero
delle "osservazioni virtuali".
 L'errore di integrazione è completamente trascurabile all'interno del
range di epoche coperto dalle librerie, nel senso che se ad esempio,
nei dati di libreria, la Luna è "sbagliata" di 2 metri,
l'integratore può al massimo, dopo un migliaio di orbite, aggiungerci
di suo un paio di centimetri.
Il punto chiave del tutto è che,  se ci si accontenta della precisione
massima delle osservazioni ottiche (diciamo un centesimo di secondo
d'arco), non è necessario un modello sofisticato come quello del JPL
per riprodurre le orbite, anche su scale temporali di migliaia di
anni.


 numerici ben dettagliati, cosě ci capiamo prima e meglio.
Queste sono le coordinate (in metri) che ottengo per la posizione della 
Terra il 2007-03-03 23:00 CT:
301946736,05316  601219239,233835  244381967,170802
ed č il valore "grezzo" dato dalle DE414 (equatore medio terrestre ed 
equinozio ICRF/J2000.0) e questi sono i valori dati da Horizons (che usa le 
DE405):
3.019467360532904E+05  6.012192392337621E+05  2.443819671707689E+05
č possibile notare l'ottimo accordo tra i due vettori.
A questo punto, sono in grado di apportare tutte le correzioni per 
calcolarmi le coordinate apparenti, ma in fase di rotazione del vettore 
dovrei commettere qualche sciocchezza.
Cristiano 

 successivi. [...]
Avendo a disposizione quella "sventola" d'integratore, mi viene pensato che 
se il programma l'avessi scritto io, avrei messo a disposizione dell'utente 
un'opzione per crearsi il suo sistema planetario, cioè Solex che si va a 
leggere un file in semplice formato testo con i nomi, i vettori di stato 
(posizione e velocità) e il GM di tutti i corpi del sistema.
Questo permetterebbe di usare, per esempio, una versione ridotta delle DExxx 
con le posizioni di tutti i pianeti prese ogni 50 o 100 anni, oppure uno 
potrebbe mettere una sonda per simulare un fly-by (sarebbe molto 
interessante per applicazioni didattiche), io, per esempio, l'ho fatto 
utilizzando il modesto RK del 4° ordine.
Se non fosse per il fatto che il JPL ha usato un integratore del 16° ordine 
in quadrupla precisione, sarei portato a dire che il tuo integratore è 
migliore di quello del JPL.  Tanto vale sfruttarlo! :-)
L'unica cosa che io veramente mi sforzerei di migliorare è la dimensione 
della finestra del programma; con le attuali risoluzioni video dei moderni 
PC, la finestra di Solex appare veramente piccolina; non che sia un 
problema, ma già che ci sei...
Nell'altro messaggio lasci intendere che non vuoi cimentarti nel lungo 
adattamento che sarebbe necessario; se vuoi, io mi offro volontario.
Cristiano 

 
Ma questo si può fare benissimo! Basta andare nel menù
"proprietà" della finestra (mentre si è in modalità testo)
e scegliersi un font più grande di quello in uso. La finestra si
allarga di conseguenza. Inoltre nella modalità grafica esiste
comunque l'opzione di riempire tutto lo schermo con la finestra
grafica.


 Eh, Eh, come si vede che non leggi i manuali ...
E pensare che sono un assiduo lettore di manuali! :-) Ora ho visto 
l'opzione.
Eventualmente il codice sorgente dell'integratore è disponibile o non hai 
piacere a renderlo pubblico? Ho scaricato "NUMERICAL INTEGRATION FOR..." 
sperando di trovare la descrizione dell'integratore, ma...
Cristiano 

 sperando di trovare la descrizione dell'integratore, ma...
Che ne dici se intanto facciamo conoscere anche agli altri Lettori il
lavoro in discorso, sei d'accordo?  ;-)
"Numerical integration for the real time production of fundamental
ephemerides over a wide time span"
by Aldo Vitagliano, Dipartimento di Chimica, Università di Napoli
"Federico II", Via Mezzocannone 4, 80134 Napoli.
Received 21 May 1996; accepted: 4 December 1996
Celestial Mechanics and Dynamical Astronomy *66* 293-308, 1997.
http://chemistry.unina.it/~alvitagl/solex/numint.pdf
Ciao.
Eliseo

 
La primissima versione DOS del programma, il "nocciolo", con
il codice sorgente e la sua libreria basata sulle DE200, è ora
scaricabile dalla pagina web di Solex (in fondoallapagina):
http://chemistry.unina.it/~alvitagl/solex/
