Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/400/1145 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 2024-Apr-20
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/400/1145 Celestial Intermediate Pole & Ephemeris Origin (Capitaine+, 2003)
*================================================================================
*Expressions for the Celestial Intermediate Pole and Celestial Ephemeris Origin
*consistent with the IAU 2000A precession-nutation model.
* Capitaine N., Chapront J., Lambert S., Wallace P.T.
* <Astron. Astrophys. 400, 1145 (2003)>
* =2003A&A...400.1145C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea.dat' ! Non-polynomial terms in the development of the
coordinate X(t) of the Celestial Intermediate Pole
in the ICRS compatible with IAU 2000A
precession-nutation model (unit mas)
integer*4 nr__
parameter (nr__=1600) ! Number of records
character*105 ar__ ! Full-size record
integer*4 N (nr__) ! Number of records of order j
integer*4 j (nr__) ! Order of the Poisson term (power of t)
integer*4 i (nr__) ! Term number i
real*8 asji (nr__) ! (uarcsec) a_{s,j})_i coefficient, in micro-arcsec
real*8 acji (nr__) ! (uarcsec) a_{c,j})_i coefficient, in micro-arcsec
integer*4 nl (nr__) ! Mean anomaly of the Moon coefficient
integer*4 nl_ (nr__) ! Mean anomaly of the Sun coefficient
integer*4 nF (nr__) ! L - Omega (L: Mean longitude of the Moon)
* coefficient
integer*4 nD (nr__) ! Mean elongation from the Moon to the Sun
* coefficient
integer*4 nOm (nr__) ! Mean longitude of the ascending node of the
* Moon (Omega) coefficient
integer*4 nLMe (nr__) ! Mean longitude of Mercure coefficient
integer*4 nLV (nr__) ! Mean longitude of Venus coefficient
integer*4 nLE (nr__) ! Mean longitude of the Earth coefficient
integer*4 nLMa (nr__) ! Mean longitude of Mars coefficient
integer*4 nLJ (nr__) ! Mean longitude of Jupiter coefficient
integer*4 nLS (nr__) ! Mean longitude of Saturn coefficient
integer*4 nLU (nr__) ! Mean longitude of Uranus coefficient
integer*4 nLN (nr__) ! Mean longitude of Neptune coefficient
integer*4 npa (nr__) ! Accumulated general precession in longitude
* coefficient
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb.dat' ! Non-polynomial terms in the development of the
coordinate Y(t) of the Celestial Intermediate Pole
in the ICRS compatible with IAU 2000A
precession-nutation model (unit mas)
integer*4 nr__1
parameter (nr__1=1275) ! Number of records
character*105 ar__1 ! Full-size record
integer*4 N_1 (nr__1) ! Number of records of order j
integer*4 j_1 (nr__1) ! Order of the Poisson term (power of t)
integer*4 i_1 (nr__1) ! Term number i
real*8 bsji (nr__1) ! (uarcsec) b_{s,j})_i coefficient, in micro-arcsec
real*8 bcji (nr__1) ! (uarcsec) b_{c,j})_i coefficient, in micro-arcsec
integer*4 nl_1 (nr__1) ! Mean anomaly of the Moon coefficient
integer*4 nl__1 (nr__1) ! Mean anomaly of the Sun coefficient
integer*4 nF_1 (nr__1) ! L - Omega (L: Mean longitude of the Moon)
* coefficient
integer*4 nD_1 (nr__1) ! Mean elongation from the Moon to the Sun
* coefficient
integer*4 nOm_1 (nr__1) ! Mean longitude of the ascending node of the
* Moon (Omega) coefficient
integer*4 nLMe_1 (nr__1) ! Mean longitude of Mercure coefficient
integer*4 nLV_1 (nr__1) ! Mean longitude of Venus coefficient
integer*4 nLE_1 (nr__1) ! Mean longitude of the Earth coefficient
integer*4 nLMa_1 (nr__1) ! Mean longitude of Mars coefficient
integer*4 nLJ_1 (nr__1) ! Mean longitude of Jupiter coefficient
integer*4 nLS_1 (nr__1) ! Mean longitude of Saturn coefficient
integer*4 nLU_1 (nr__1) ! Mean longitude of Uranus coefficient
integer*4 nLN_1 (nr__1) ! Mean longitude of Neptune coefficient
integer*4 npa_1 (nr__1) ! Accumulated general precession in longitude
* coefficient
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec.dat' ! Non-polynomial terms in the development of the
quantity s(t)+XY/2 compatible with IAU 2000A
precession-nutation model (unit mas)
integer*4 nr__2
parameter (nr__2=66) ! Number of records
character*105 ar__2 ! Full-size record
integer*4 N_2 (nr__2) ! Number of records of order j
integer*4 j_2 (nr__2) ! Order of the Poisson term (power of t)
integer*4 i_2 (nr__2) ! Term number i
real*8 Csji (nr__2) ! (uarcsec) C_{s,j})_i coefficient, in micro-arcsec
real*8 Ccji (nr__2) ! (uarcsec) C_{c,j})_i coefficient, in micro-arcsec
integer*4 nl_2 (nr__2) ! Mean anomaly of the Moon coefficient
integer*4 nl__2 (nr__2) ! Mean anomaly of the Sun coefficient
integer*4 nF_2 (nr__2) ! L - Omega (L: Mean longitude of the Moon)
* coefficient
integer*4 nD_2 (nr__2) ! Mean elongation from the Moon to the Sun
* coefficient
integer*4 nOm_2 (nr__2) ! Mean longitude of the ascending node of the
* Moon (Omega) coefficient
integer*4 nLMe_2 (nr__2) ! Mean longitude of Mercure coefficient
integer*4 nLV_2 (nr__2) ! Mean longitude of Venus coefficient
integer*4 nLE_2 (nr__2) ! Mean longitude of the Earth coefficient
integer*4 nLMa_2 (nr__2) ! Mean longitude of Mars coefficient
integer*4 nLJ_2 (nr__2) ! Mean longitude of Jupiter coefficient
integer*4 nLS_2 (nr__2) ! Mean longitude of Saturn coefficient
integer*4 nLU_2 (nr__2) ! Mean longitude of Uranus coefficient
integer*4 nLN_2 (nr__2) ! Mean longitude of Neptune coefficient
integer*4 npa_2 (nr__2) ! Accumulated general precession in longitude
* coefficient
C=============================================================================
C Loading file 'tablea.dat' ! Non-polynomial terms in the development of the
* coordinate X(t) of the Celestial Intermediate Pole
* in the ICRS compatible with IAU 2000A
* precession-nutation model (unit mas)
C Format for file interpretation
1 format(
+ I4,1X,I1,1X,I5,1X,F11.2,1X,F10.2,3X,I2,3X,I2,3X,I2,3X,I2,3X,
+ I2,3X,I2,2X,I3,2X,I3,2X,I3,3X,I2,3X,I2,3X,I2,3X,I2,3X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea.dat')
write(6,*) '....Loading file: tablea.dat'
do i__=1,1600
read(1,'(A105)')ar__
read(ar__,1)
+ N(i__),j(i__),i(i__),asji(i__),acji(i__),nl(i__),nl_(i__),
+ nF(i__),nD(i__),nOm(i__),nLMe(i__),nLV(i__),nLE(i__),
+ nLMa(i__),nLJ(i__),nLS(i__),nLU(i__),nLN(i__),npa(i__)
c ..............Just test output...........
write(6,1)
+ N(i__),j(i__),i(i__),asji(i__),acji(i__),nl(i__),nl_(i__),
+ nF(i__),nD(i__),nOm(i__),nLMe(i__),nLV(i__),nLE(i__),
+ nLMa(i__),nLJ(i__),nLS(i__),nLU(i__),nLN(i__),npa(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb.dat' ! Non-polynomial terms in the development of the
* coordinate Y(t) of the Celestial Intermediate Pole
* in the ICRS compatible with IAU 2000A
* precession-nutation model (unit mas)
C Format for file interpretation
2 format(
+ I4,1X,I1,1X,I5,1X,F11.2,1X,F10.2,3X,I2,3X,I2,3X,I2,3X,I2,3X,
+ I2,3X,I2,2X,I3,2X,I3,2X,I3,3X,I2,3X,I2,3X,I2,3X,I2,3X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb.dat')
write(6,*) '....Loading file: tableb.dat'
do i__=1,1275
read(1,'(A105)')ar__1
read(ar__1,2)
+ N_1(i__),j_1(i__),i_1(i__),bsji(i__),bcji(i__),nl_1(i__),
+ nl__1(i__),nF_1(i__),nD_1(i__),nOm_1(i__),nLMe_1(i__),
+ nLV_1(i__),nLE_1(i__),nLMa_1(i__),nLJ_1(i__),nLS_1(i__),
+ nLU_1(i__),nLN_1(i__),npa_1(i__)
c ..............Just test output...........
write(6,2)
+ N_1(i__),j_1(i__),i_1(i__),bsji(i__),bcji(i__),nl_1(i__),
+ nl__1(i__),nF_1(i__),nD_1(i__),nOm_1(i__),nLMe_1(i__),
+ nLV_1(i__),nLE_1(i__),nLMa_1(i__),nLJ_1(i__),nLS_1(i__),
+ nLU_1(i__),nLN_1(i__),npa_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec.dat' ! Non-polynomial terms in the development of the
* quantity s(t)+XY/2 compatible with IAU 2000A
* precession-nutation model (unit mas)
C Format for file interpretation
3 format(
+ I4,1X,I1,1X,I5,1X,F11.2,1X,F10.2,3X,I2,3X,I2,3X,I2,3X,I2,3X,
+ I2,3X,I2,2X,I3,2X,I3,2X,I3,3X,I2,3X,I2,3X,I2,3X,I2,3X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec.dat')
write(6,*) '....Loading file: tablec.dat'
do i__=1,66
read(1,'(A105)')ar__2
read(ar__2,3)
+ N_2(i__),j_2(i__),i_2(i__),Csji(i__),Ccji(i__),nl_2(i__),
+ nl__2(i__),nF_2(i__),nD_2(i__),nOm_2(i__),nLMe_2(i__),
+ nLV_2(i__),nLE_2(i__),nLMa_2(i__),nLJ_2(i__),nLS_2(i__),
+ nLU_2(i__),nLN_2(i__),npa_2(i__)
c ..............Just test output...........
write(6,3)
+ N_2(i__),j_2(i__),i_2(i__),Csji(i__),Ccji(i__),nl_2(i__),
+ nl__2(i__),nF_2(i__),nD_2(i__),nOm_2(i__),nLMe_2(i__),
+ nLV_2(i__),nLE_2(i__),nLMa_2(i__),nLJ_2(i__),nLS_2(i__),
+ nLU_2(i__),nLN_2(i__),npa_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end