Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/154/62 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Apr-13
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/AJ/154/62 Orbital parameters of Kuiper Belt objects (Volk+, 2017)
*================================================================================
*The curiously warped mean plane of the Kuiper Belt.
* Volk K., Malhotra R.
* <Astron. J., 154, 62-62 (2017)>
* =2017AJ....154...62V (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! List objects used to calculate mean planes
integer*4 nr__
parameter (nr__=906) ! Number of records
character*108 ar__ ! Full-size record
character*7 MPC (nr__) ! Minor Planet Center designation (1)
real*4 a (nr__) ! (AU) [36.3/149.85] Best-fit barycentric semi-major
* axis (2)
real*4 e_a (nr__) ! (AU) [0.0005/4.14] The 1-{sigma} uncertainty in a
* (da) (3)
real*4 e (nr__) ! [0.0018/0.76] Best-fit eccentricity (2)
real*4 i (nr__) ! (rad) [0.0018/0.88] Best-fit inclination (2)
real*4 Omega (nr__) ! (rad) [0.0022/6.28] Best-fit longitude of ascending
* node ({Omega}) (2)
real*4 omega_1 (nr__) ! (rad) [0.0014/6.27] Best-fit argument of perihelion
* ({omega}) (2)
real*4 lambda (nr__) ! (rad) [-3.14/3.14] Ecliptic longitude ({lambda}) (2)
real*4 beta (nr__) ! (rad) [-0.844/0.805] Ecliptic latitude ({beta}) (2)
real*4 Dist (nr__) ! (AU) [31.544/96.255] Heliocentric distance (r_h_)
real*8 Epoch (nr__) ! (d) Epoch of the orbit fit and sky position
* (in Julian Date)
*Note (1): The MPC designations are given in their packed format (see
* http://www.minorplanetcenter.net/iau/info/PackedDes.html).
*Note (2): The orbital elements are the best-fit barycentric elements from a
* Bernstein & Khushalani (2000AJ....120.3323B) orbit fit to the astrometry
* available for each object from the Minor Planet Center (MPC). All elements
* are referenced to the J2000 ecliptic coordinate system.
*Note (3): Taken from the orbit-fit covariance matrix.
C=============================================================================
C Loading file 'table1.dat' ! List objects used to calculate mean planes
C Format for file interpretation
1 format(
+ A7,1X,E11.5,1X,E9.3,1X,E11.5,1X,E11.5,1X,E11.5,1X,E11.5,1X,
+ F6.3,1X,F6.3,1X,F6.3,1X,F9.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,906
read(1,'(A108)')ar__
read(ar__,1)
+ MPC(i__),a(i__),e_a(i__),e(i__),i(i__),Omega(i__),
+ omega_1(i__),lambda(i__),beta(i__),Dist(i__),Epoch(i__)
c ..............Just test output...........
write(6,1)
+ MPC(i__),a(i__),e_a(i__),e(i__),i(i__),Omega(i__),
+ omega_1(i__),lambda(i__),beta(i__),Dist(i__),Epoch(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end