FORTRAN Generation
(/./ftp/cats/J/A_A/355/398)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/355/398 into FORTRAN code for reading data files line by line.

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 2024-Mar-28
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           ! [1/269]+ Number i
      integer*4     l           ! Mean anomaly of the Moon coefficient
      integer*4     l_          ! Mean anomaly of the Sun coefficient
      integer*4     F           ! L - Omega (L: Mean longitude of the Moon)
*                                 coefficient
      integer*4     D           ! Mean elongation from the Moon to the Sun
*                                 coefficient
      integer*4     Omega       ! Mean longitude of the ascending node of the Moon
*                                 coefficient
      real*8        Per         ! (d) Period of the nutation
      integer*4     as0i        ! (uas) (a_s,0_)_i_ coefficient (1)
      integer*4     as1i        ! (uas/hyr) (a_s,1_)_i_ coefficient (2)
      integer*4     as2i        ! (uas/hyr2) (a_s,2_)_i_ coefficient (2)
      integer*4     ac0i        ! (uas) (a_c,0_)_i_ coefficient
      integer*4     ac1i        ! (uas/hyr) (a_c,1_)_i_ coefficient
      integer*4     ac2i        ! (uas/hyr2) (a_c,2_)_i_ coefficient
      integer*4     bc0i        ! (uas) (b_c,0_)_i_ coefficient (1)
      integer*4     bc1i        ! (uas/hyr) (b_c,1_)_i_ coefficient (2)
      integer*4     bc2i        ! (uas/hyr2) (b_c,2_)_i_ coefficient (2)
      integer*4     bs0i        ! (uas) (b_s,0_)_i_ coefficient
      integer*4     bs1i        ! (uas/hyr) (b_s,1_)_i_ coefficient
      integer*4     bs2i        ! (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         ! [1/112]+ Number i
      integer*4     lVe         ! Mean longitude of Venus coefficient
      integer*4     lE          ! Mean longitude of the Earth coefficient
      integer*4     lMa         ! Mean longitude of Mars coefficient
      integer*4     lJ          ! Mean longitude of Jupiter coefficient
      integer*4     lSa         ! Mean longitude of Saturn coefficient
      integer*4     pa          ! General precession in longitude coefficient
      integer*4     D_1         ! Mean elongation from the Moon to the Sun
*                                 coefficient
      integer*4     F_1         ! L - Omega (L: Mean longitude of the Moon)
*                                 coefficient
      integer*4     l_1         ! Mean anomaly of the Moon coefficient
      integer*4     Omega_1     ! Mean longitude of the ascending node of the Moon
*                                 coefficient
      real*8        Per_1       ! (d) Period of the nutation
      integer*4     as0i_1      ! (uas) (a_s,0_)_i_ coefficient
      integer*4     as1i_1      ! (uas/hyr) (a_s,1_)_i_ coefficient
      integer*4     as2i_1      ! (uas/hyr2) (a_s,2_)_i_ coefficient
      integer*4     ac0i_1      ! (uas) (a_c,0_)_i_ coefficient
      integer*4     ac1i_1      ! (uas/hyr) (a_c,1_)_i_ coefficient
      integer*4     ac2i_1      ! (uas/hyr2) (a_c,2_)_i_ coefficient
      integer*4     bc0i_1      ! (uas) (b_c,0_)_i_ coefficient
      integer*4     bc1i_1      ! (uas/hyr) (b_c,1_)_i_ coefficient
      integer*4     bc2i_1      ! (uas/hyr2) (b_c,2_)_i_ coefficient
      integer*4     bs0i_1      ! (uas) (b_s,0_)_i_ coefficient
      integer*4     bs1i_1      ! (uas/hyr) (b_s,1_)_i_ coefficient
      integer*4     bs2i_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,l,l_,F,D,Omega,Per,as0i,as1i,as2i,ac0i,ac1i,ac2i,bc0i,bc1i,
     +  bc2i,bs0i,bs1i,bs2i
c    ..............Just test output...........
        write(6,1)
     +  i,l,l_,F,D,Omega,Per,as0i,as1i,as2i,ac0i,ac1i,ac2i,bc0i,bc1i,
     +  bc2i,bs0i,bs1i,bs2i
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,lVe,lE,lMa,lJ,lSa,pa,D_1,F_1,l_1,Omega_1,Per_1,as0i_1,
     +  as1i_1,as2i_1,ac0i_1,ac1i_1,ac2i_1,bc0i_1,bc1i_1,bc2i_1,
     +  bs0i_1,bs1i_1,bs2i_1
c    ..............Just test output...........
        write(6,2)
     +  i_1,lVe,lE,lMa,lJ,lSa,pa,D_1,F_1,l_1,Omega_1,Per_1,as0i_1,
     +  as1i_1,as2i_1,ac0i_1,ac1i_1,ac2i_1,bc0i_1,bc1i_1,bc2i_1,
     +  bs0i_1,bs1i_1,bs2i_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end