Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/355/398 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-Jun-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/A+A/355/398 Celestial Ephemeris Origin definition (Capitaine+, 2000)
*================================================================================
*Definition of the Celestial Ephemeris Origin and of UT1 in the
*International Celestial Reference Frame.
* Capitaine N., Guinot B., McCarthy D.D.
* <Astron. Astrophys. 355, 398 (2000)>
* =2000A&A...355..398C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1a.dat' ! Numerical Development of the coordinates X(t),
Y(t) of the Celestial Ephemeris Pole in the ICRS
consistent with the IERS 1996 Conventions and
the IERS Annual report for 1997 (unit {mu}as)
integer*4 nr__
parameter (nr__=269) ! Number of records
character*97 ar__ ! Full-size record
integer*4 i (nr__) ! [1/269]+ Number i
integer*4 l (nr__) ! Mean anomaly of the Moon coefficient
integer*4 l_ (nr__) ! Mean anomaly of the Sun coefficient
integer*4 F (nr__) ! L - Omega (L: Mean longitude of the Moon)
* coefficient
integer*4 D (nr__) ! Mean elongation from the Moon to the Sun
* coefficient
integer*4 Omega (nr__) ! Mean longitude of the ascending node of the Moon
* coefficient
real*8 Per (nr__) ! (d) Period of the nutation
integer*4 as0i (nr__) ! (uas) (a_s,0_)_i_ coefficient (1)
integer*4 as1i (nr__) ! (uas/hyr) (a_s,1_)_i_ coefficient (2)
integer*4 as2i (nr__) ! (uas/hyr2) (a_s,2_)_i_ coefficient (2)
integer*4 ac0i (nr__) ! (uas) (a_c,0_)_i_ coefficient
integer*4 ac1i (nr__) ! (uas/hyr) (a_c,1_)_i_ coefficient
integer*4 ac2i (nr__) ! (uas/hyr2) (a_c,2_)_i_ coefficient
integer*4 bc0i (nr__) ! (uas) (b_c,0_)_i_ coefficient (1)
integer*4 bc1i (nr__) ! (uas/hyr) (b_c,1_)_i_ coefficient (2)
integer*4 bc2i (nr__) ! (uas/hyr2) (b_c,2_)_i_ coefficient (2)
integer*4 bs0i (nr__) ! (uas) (b_s,0_)_i_ coefficient
integer*4 bs1i (nr__) ! (uas/hyr) (b_s,1_)_i_ coefficient
integer*4 bs2i (nr__) ! (uas/hyr2) (b_s,2_)_i_ coefficient
*Note (1): The amplitudes (a_s,0_)_i_, (b_c,0_)_i_ are equal to the amplitudes
* A_i_.sin{epsilon}_0_ and B_i_ of the IERS 1996 series for nutation
* in longitude x sin{epsilon}_0_ and obliquity, except for 20 terms in
* each coordinate X and Y in which appears a contribution from
* crossed-nutation effect.
*Note (2): The amplitudes (a_s,j_)_i_, (b_c,j_)_i_, for j=1,2, are due to crossed
* terms between precession and nutation i, of the form t.sin or
* t.cos for j=1 and of the form t^2^.sin or t^2^.cos for j=2.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1b.dat' ! Numerical Development of the planetary nutations
in the coordinates X(t), Y(t) of the CEP in the
ICRS consistent with the IERS 1996 Conventions
and the IERS Annual report for 1997 (unit {mu}as)
integer*4 nr__1
parameter (nr__1=112) ! Number of records
character*82 ar__1 ! Full-size record
integer*4 i_1 (nr__1) ! [1/112]+ Number i
integer*4 lVe (nr__1) ! Mean longitude of Venus coefficient
integer*4 lE (nr__1) ! Mean longitude of the Earth coefficient
integer*4 lMa (nr__1) ! Mean longitude of Mars coefficient
integer*4 lJ (nr__1) ! Mean longitude of Jupiter coefficient
integer*4 lSa (nr__1) ! Mean longitude of Saturn coefficient
integer*4 pa (nr__1) ! General precession in longitude coefficient
integer*4 D_1 (nr__1) ! Mean elongation from the Moon to the Sun
* coefficient
integer*4 F_1 (nr__1) ! L - Omega (L: Mean longitude of the Moon)
* coefficient
integer*4 l_1 (nr__1) ! Mean anomaly of the Moon coefficient
integer*4 Omega_1 (nr__1) ! Mean longitude of the ascending node of the Moon
* coefficient
real*8 Per_1 (nr__1) ! (d) Period of the nutation
integer*4 as0i_1 (nr__1) ! (uas) (a_s,0_)_i_ coefficient
integer*4 as1i_1 (nr__1) ! (uas/hyr) (a_s,1_)_i_ coefficient
integer*4 as2i_1 (nr__1) ! (uas/hyr2) (a_s,2_)_i_ coefficient
integer*4 ac0i_1 (nr__1) ! (uas) (a_c,0_)_i_ coefficient
integer*4 ac1i_1 (nr__1) ! (uas/hyr) (a_c,1_)_i_ coefficient
integer*4 ac2i_1 (nr__1) ! (uas/hyr2) (a_c,2_)_i_ coefficient
integer*4 bc0i_1 (nr__1) ! (uas) (b_c,0_)_i_ coefficient
integer*4 bc1i_1 (nr__1) ! (uas/hyr) (b_c,1_)_i_ coefficient
integer*4 bc2i_1 (nr__1) ! (uas/hyr2) (b_c,2_)_i_ coefficient
integer*4 bs0i_1 (nr__1) ! (uas) (b_s,0_)_i_ coefficient
integer*4 bs1i_1 (nr__1) ! (uas/hyr) (b_s,1_)_i_ coefficient
integer*4 bs2i_1 (nr__1) ! (uas/hyr2) (b_s,2_)_i_ coefficient
C=============================================================================
C Loading file 'table1a.dat' ! Numerical Development of the coordinates X(t),
* Y(t) of the Celestial Ephemeris Pole in the ICRS
* consistent with the IERS 1996 Conventions and
* the IERS Annual report for 1997 (unit {mu}as)
C Format for file interpretation
1 format(
+ I3,1X,I2,1X,I2,1X,I2,1X,I2,1X,I1,1X,F9.2,1X,I8,1X,I5,1X,I4,1X,
+ I4,1X,I6,1X,I2,1X,I7,1X,I4,1X,I5,1X,I4,1X,I6,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table1a.dat')
write(6,*) '....Loading file: table1a.dat'
do i__=1,269
read(1,'(A97)')ar__
read(ar__,1)
+ i(i__),l(i__),l_(i__),F(i__),D(i__),Omega(i__),Per(i__),
+ as0i(i__),as1i(i__),as2i(i__),ac0i(i__),ac1i(i__),ac2i(i__),
+ bc0i(i__),bc1i(i__),bc2i(i__),bs0i(i__),bs1i(i__),bs2i(i__)
c ..............Just test output...........
write(6,1)
+ i(i__),l(i__),l_(i__),F(i__),D(i__),Omega(i__),Per(i__),
+ as0i(i__),as1i(i__),as2i(i__),ac0i(i__),ac1i(i__),ac2i(i__),
+ bc0i(i__),bc1i(i__),bc2i(i__),bs0i(i__),bs1i(i__),bs2i(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1b.dat' ! Numerical Development of the planetary nutations
* in the coordinates X(t), Y(t) of the CEP in the
* ICRS consistent with the IERS 1996 Conventions
* and the IERS Annual report for 1997 (unit {mu}as)
C Format for file interpretation
2 format(
+ I3,1X,I3,1X,I3,1X,I3,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,1X,
+ I2,1X,F9.2,1X,I3,1X,I2,1X,I1,1X,I3,1X,I2,1X,I1,1X,I3,1X,I2,1X,
+ I1,1X,I3,1X,I2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table1b.dat')
write(6,*) '....Loading file: table1b.dat'
do i__=1,112
read(1,'(A82)')ar__1
read(ar__1,2)
+ i_1(i__),lVe(i__),lE(i__),lMa(i__),lJ(i__),lSa(i__),pa(i__),
+ D_1(i__),F_1(i__),l_1(i__),Omega_1(i__),Per_1(i__),
+ as0i_1(i__),as1i_1(i__),as2i_1(i__),ac0i_1(i__),ac1i_1(i__),
+ ac2i_1(i__),bc0i_1(i__),bc1i_1(i__),bc2i_1(i__),bs0i_1(i__),
+ bs1i_1(i__),bs2i_1(i__)
c ..............Just test output...........
write(6,2)
+ i_1(i__),lVe(i__),lE(i__),lMa(i__),lJ(i__),lSa(i__),pa(i__),
+ D_1(i__),F_1(i__),l_1(i__),Omega_1(i__),Per_1(i__),
+ as0i_1(i__),as1i_1(i__),as2i_1(i__),ac0i_1(i__),ac1i_1(i__),
+ ac2i_1(i__),bc0i_1(i__),bc1i_1(i__),bc2i_1(i__),bs0i_1(i__),
+ bs1i_1(i__),bs2i_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end