Astronomia Teorica e Meccanica Celeste

Since summer 2010 this website, which involves Astrodynamics and Celestial Mechanics, has a limited section for English Speaking Users (ESU), where new enhancements of SOFA routines are shown. As usually done up to now, much worked examples will be presented with their complete codes (in Fortran90 and C) alongwith the related outputs in text format.
To download SOFA routines/functions please take the last version directly from its Official Site. And now, get a glance to: ESU Main Page

(13 Dicembre 2016). Operazione di restyling del codice asteroidale aster16.c pubblicato sotto. Per la lettura delle stringhe del database non viene usato il laborioso comando fseek ma fscanf che si avvale di un file di deposito esterno per copiarvi (in una riga) i dati del record intercettato. Ciò permette di modificare rapidamente il codice C in caso di variazione del template dei campi (larghezza e numero). L'algoritmo è stato provato sul file AsterSelect.c ed applicato successivamente allo script finale AsterDIC16.c, entrambi reperibili qui, insieme allo stesso dbase di 20,000 pianetini.

(28 Settembre 2016). Una interessante variante del software asteroidale è questa. I due Database degli elementi orbitali, richiamati dai codici AstGeocP.c e AstNEAp.c, contengono rispettivamente 20,000 asteroidi del catalogo MPC e 14,690 NEA, acronimo di Near Earth Asteroids; essi vengono scandagliati scegliendo come parametro di screening la distanza geocentrica; per una migliore consultazione, i risultati ottenuti sono copiati su file testo. Nel caso dei NEA, si percepisce subito quali siano i più vicini alla Terra alla data dell'effemeride, come nell'esempio seguente:


Distanza Geocentrica limite assunta= 0.02500
Data Effemeride= 2016  9 24

	            	  2016 RM20       
	Distanza GEOCENTRICA d=   0.01991001   (UA)
	Distanza Eliocentr.  r=   1.01175656   (UA)
	Long.Eclitt.Eliocen. l=    1.07918  (gradi) 
	Lati.Eclitt.Eliocen. b=   -1.01969  (gradi) 
	  Ascensione Retta  AR=   36.12525  (gradi) =  2H 24M 30S
	      Declinazione  DE=  -58.36825  (gradi) 
	                   psi= -114.98831  (gradi) Elongaz. Ovest
	   Angolo di Fase beta=   63.98968  (gradi) 
	           Magnitudine=   20.0          

	            	  2016 RD34       
	Distanza GEOCENTRICA d=   0.01327780   (UA)
	Distanza Eliocentr.  r=   1.00928176   (UA)
	Long.Eclitt.Eliocen. l=    1.82229  (gradi) 
	Lati.Eclitt.Eliocen. b=   -0.01367  (gradi) 
	  Ascensione Retta  AR=   62.36386  (gradi) =  4H  9M 27S
	      Declinazione  DE=   19.96857  (gradi) 
	                   psi= -116.99746  (gradi) Elongaz. Ovest
	   Angolo di Fase beta=   62.33090  (gradi) 
	           Magnitudine=   20.1          

	   =============================
	    N.ro corpi selezionati=   2
	   =============================


	
	
Altro esempio, per gentile concessione dell'autore di AstGeoc.f90 [Fortran 90],
l'astrofilo piemontese Aldo Nicola.

        Date and Time UTC   
    ==============================
  Input : Year , Month, Day (yyyy,mm,dd) 
2016,10,1

  Input : Maximum Geocentric distance 1.00 [UA]
  
                (2016, revised by G. Matarazzo)
       Time of Ephemeris  ID   2457662.5007891473        Time of epoch  JD    2457600.5000000000     
 --------------------------------------------------------------------------------------------------
                             R.A.2000       Decl.2000   Delta        r      Phase   Mag  Elong.    
  num      body             hh mm ss.ss    dd pp ss.s    A.U.       A.U.     deg     V    deg      
 --------------------------------------------------------------------------------------------------
 00433   Eros               21  0 12.87  -  2 21 50.50  0.846252  1.657753   28.7  12.75  127.4 E     
 01863   Antinous            0  4 18.17  + 16 11 57.34  0.919574  1.905601    7.5  17.30  165.6 E     
 01865   Cerberus            5 55 50.21  + 19 54 23.11  0.685827  1.298868   49.6  18.44   99.0 W     
 01943   Anteros             1 37 11.49  + 25 41 24.22  0.783340  1.733142   15.6  17.27  152.2 W     
 02100   Ra-Shalom           3  7 42.15  - 25 51 57.79  0.162216  1.116143   41.8  14.03  131.9 W     
 03103   Eger                8 12 41.29  - 19 25 44.57  0.520085  0.909234   84.2  16.70   64.6 W     
 03200   Phaethon           13 58  8.90  + 82 51 22.89  0.403303  1.056688   71.0  15.24   86.6 E     
 03352   McAuliffe           3  1 54.05  +  6 45 50.91  0.818840  1.721235   21.0  17.58  141.9 W     
 03753   Cruithne            6 50 56.68  - 13 46 43.52  0.644850  1.133080   61.5  17.11   84.0 W     
 04179   Toutatis           18 52 11.08  - 23 21 37.69  0.633991  1.221312   54.9  16.86   93.9 E     
 04341   Poseidon           23 50 16.60  - 31 14 20.35  0.674437  1.598263   21.4  17.21  144.3 E     
 04581   Asclepius          17  8 41.34  - 19  0 22.81  0.727503  1.015213   67.8  22.43   69.9 E     
 04953   1990 MU             6  5 16.61  -  3 36 31.94  0.773186  1.325591   48.7  15.98   95.8 W     
 05381   Sekhmet             9 35 54.04  - 51 52 52.45  0.649299  0.885928   79.8  18.19   60.6 W     
 05653   Camarillo           4 18 22.87  + 27  1 12.91  0.934055  1.679978   30.9  18.40  120.4 W     
 05786   Talos              23 16 17.22  - 18  4 25.08  0.956216  1.899822   14.3  19.20  152.1 E     
 05836   1993 MF             3 10  7.85  + 30 24 38.95  0.545095  1.431187   30.6  15.58  133.3 W     
 05863   Tara               18  7 28.72  - 12 24 16.33  0.785987  1.206008   55.6  17.71   84.0 E     
 06239   Minos              16  7 49.78  - 23 12 17.01  0.946415  0.921460   64.8  20.50   56.4 E     
 06569   Ondaatje            3 18 54.45  - 34 37 49.91  0.794182  1.596334   30.9  18.24  125.1 W     
 07088   Ishtar              0 51 33.27  - 16 12 50.45  0.708463  1.684057   11.8  17.80  159.9 W     
 07341   1991 VK             0 18 39.83  + 19 14 29.00  0.497297  1.485093   10.9  16.72  163.7 E     
 07888   1993 UC             6 29 48.04  - 12 57  4.86  0.814011  1.281518   51.4  17.09   89.2 W     
 10150   1994 PN            16 20 16.84  - 56  9 18.88  0.904571  1.095014   59.2  17.50   69.9 E     
 10636   1998 QK56          23 10 39.21  - 10 15 51.21  0.539845  1.510947   15.5  18.01  156.2 E     
 11500   Tomaiyowit          5 21 31.17  + 33  9  0.18  0.516005  1.245886   50.6  19.22  105.9 W     
 15745   Yuliya              3  6 15.24  - 11 39 45.48  0.878535  1.757120   22.3  19.31  138.3 W     
 16834   1997 WU22          19 57 10.94  + 14 32 16.72  0.682906  1.407616   41.3  17.13  112.0 E     
  END OF PROGRAM 
	   	

(Settembre 2016).Con il software asteroidale pubblicato qui appresso, voglio rendere omaggio all'amico, collega ed astrofilo Sergio Foglia ben conosciuto, in Italia e all'estero, per la sua poliedrica attività di scopritore di pianetini e cultore di programmi astronomici di rara eleganza e precisione.
Nel lontano 1996 ho avuto il privilegio di ricevere una personal release del codice MPEPH.c per il calcolo delle effemeridi degli asteroidi tramite lettura del Dbase degli elementi orbitali pubblicato dal Minor Planet Center (MPC). Per avere un'idea più dettagliata basta leggere questo suo esaustivo paper sull'argomento.
Va da sé, però, che nel corso degli anni il MPC ha cambiato lo standard del database, aumentando considerevolmente il numero dei campi (fissi) da cui estrarre i dati e le modalità di lettura di alcuni: stringhe invece di numeri. Ho quindi rivisto la parte del codice da cambiare e con il permesso dell'autore procedo alla pubblicazione della sua opera intellettuale, scaricabile da qui; il file zippato comprende il più recente dbase limitato ai primi 20mila asteroidi del catalogo MPC e il nuovo listato aster16.c, per la cui compilazione è sufficiente usare (gcc) di licenza GNU, gratuito e di facile reperibilità sia per sistemi Windows che Linux.
Grazie di cuore, Sergio.



        +-------------------------------------------+
        |         MINOR  PLANET  EPHEMERIS          |
        |            M P E P H   v. 2.1             |
        | 1996, S.Foglia, Serafino Zani Observatory |
        |     (2016, revisited by G. Matarazzo)     |
        +-------------------------------------------+

        Minor Planet catalogue number [1 to 20,000]: 1620

         Anomalia Media M=  322.32014 (gradi)
         Argomento Peri w=  276.86945 (gradi)
         Longitud. Nodo O=  337.21625 (gradi)
           Inclinazione i=   13.33747 (gradi)
          Eccentricita' e=   0.3353926
             Moto Medio n=   0.70916069 (gradi/g)
         Semiasse Magg. a=   1.2453926  (UA)
         Numero Catalogo :    (1620)
          Nome Pianetino :  Geographos
              Parametr