Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/324/366 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 2026-Jun-08
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/324/366 Theory of motion & ephemerides of Hyperion (Duriez+ 1997)
*================================================================================
*Theory of motion and ephemerides of Hyperion
* Duriez L., Vienne A.
* <Astron. Astrophys. 324, 366 (1997)>
* =1997A&A...324..366D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3' ! Fundamental arguments of the theory
integer*4 nr__
parameter (nr__=7) ! Number of records
character*74 ar__ ! Full-size record
integer*4 Num ! Number of the argument
character*5 Arg ! Argument, see note (1)
real*8 Freq ! (rad/d) Frequency
real*8 Phas ! (rad) Phase
real*8 Per ! (d) ? Period
*Note (1): psi: Synodic argument between Titan and Hyperion
* tau: argument of the libration
* pi7: longitude of the proper pericentre of Hyperion
* pi6: longitude of the proper pericentre of Titan
* Om7: longitude of the proper node of Hyperion
* Om6: longitude of the proper node of Titan
* Om0: longitude of the node of the invariable plane
* Each argument is (Freq * t + Phas) where: t = Julian Date - 2451545.0
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4' ! Series for element p of Hyperion
integer*4 nr__1
parameter (nr__1=105) ! Number of records
character*70 ar__1 ! Full-size record
integer*4 Num_1 ! Number of the term
real*8 Ampl ! (rad) Amplitude
integer*4 N1 ! Argument (1)
integer*4 N2 ! Argument (1)
integer*4 N3 ! Argument (1)
integer*4 N4 ! Argument (1)
integer*4 N5 ! Argument (1)
integer*4 N6 ! Argument (1)
integer*4 N7 ! Argument (1)
real*8 Per_1 ! (d) ? Period
*Note (1): The argument of each term has to be computed as:
* N1*psi + N2*tau + N3*pi7 + N4*pi6 + N5*Om7 + N6*Om6 + N7*Om0
* where psi, tau, pi7, pi6, Om7, Om6, Om0 are given in table3.
* This is the argument of a cosine in table4, of a sine in table5
* and of a complex exponential in table6 and table7
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5' ! Series for element q of Hyperion
integer*4 nr__2
parameter (nr__2=214) ! Number of records
character*70 ar__2 ! Full-size record
integer*4 Num_2 ! Number of the term
real*8 Ampl_1 ! (rad) Amplitude
integer*4 N1_1 ! Argument (1)
integer*4 N2_1 ! Argument (1)
integer*4 N3_1 ! Argument (1)
integer*4 N4_1 ! Argument (1)
integer*4 N5_1 ! Argument (1)
integer*4 N6_1 ! Argument (1)
integer*4 N7_1 ! Argument (1)
real*8 Per_2 ! (d) ? Period
*Note (1): The argument of each term has to be computed as:
* N1*psi + N2*tau + N3*pi7 + N4*pi6 + N5*Om7 + N6*Om6 + N7*Om0
* where psi, tau, pi7, pi6, Om7, Om6, Om0 are given in table3.
* This is the argument of a cosine in table4, of a sine in table5
* and of a complex exponential in table6 and table7
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6' ! Series for element z of Hyperion
integer*4 nr__3
parameter (nr__3=179) ! Number of records
character*70 ar__3 ! Full-size record
integer*4 Num_3 ! Number of the term
real*8 Ampl_2 ! (rad) Amplitude
integer*4 N1_2 ! Argument (1)
integer*4 N2_2 ! Argument (1)
integer*4 N3_2 ! Argument (1)
integer*4 N4_2 ! Argument (1)
integer*4 N5_2 ! Argument (1)
integer*4 N6_2 ! Argument (1)
integer*4 N7_2 ! Argument (1)
real*8 Per_3 ! (d) ? Period
*Note (1): The argument of each term has to be computed as:
* N1*psi + N2*tau + N3*pi7 + N4*pi6 + N5*Om7 + N6*Om6 + N7*Om0
* where psi, tau, pi7, pi6, Om7, Om6, Om0 are given in table3.
* This is the argument of a cosine in table4, of a sine in table5
* and of a complex exponential in table6 and table7
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7' ! Series for element zeta of Hyperion
integer*4 nr__4
parameter (nr__4=52) ! Number of records
character*70 ar__4 ! Full-size record
integer*4 Num_4 ! Number of the term
real*8 Ampl_3 ! (rad) Amplitude
integer*4 N1_3 ! Argument (1)
integer*4 N2_3 ! Argument (1)
integer*4 N3_3 ! Argument (1)
integer*4 N4_3 ! Argument (1)
integer*4 N5_3 ! Argument (1)
integer*4 N6_3 ! Argument (1)
integer*4 N7_3 ! Argument (1)
real*8 Per_4 ! (d) ? Period
*Note (1): The argument of each term has to be computed as:
* N1*psi + N2*tau + N3*pi7 + N4*pi6 + N5*Om7 + N6*Om6 + N7*Om0
* where psi, tau, pi7, pi6, Om7, Om6, Om0 are given in table3.
* This is the argument of a cosine in table4, of a sine in table5
* and of a complex exponential in table6 and table7
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8' ! Solar and short period perturbations of Hyperion
integer*4 nr__5
parameter (nr__5=46) ! Number of records
character*92 ar__5 ! Full-size record
character*3 Ser ! Element of Hyperion (1)
real*8 Ampl_4 ! (rad) Amplitude
real*8 Freq_1 ! (rad/d) Frequency
real*8 Phas_1 ! (rad) Phase
real*8 Per_5 ! (d) ? Period
*Note (1): p7 : perturbation of element p of Hyperion (series in cosine)
* q7 : perturbation of element q of Hyperion (series in sine)
* z7 : perturbation of element z of Hyperion (series in complex
* exponential)
* zt7: perturbation of zeta of Hyperion (series in complex exponential)
* The argument of each term is: (Freq * t + Phas),
* where: t = Julian Date - 2451545.0
C=============================================================================
C Loading file 'table3' ! Fundamental arguments of the theory
C Format for file interpretation
1 format(I3,A5,D25.15,D25.15,F16.6)
C Effective file loading
open(unit=1,status='old',file=
+'table3')
write(6,*) '....Loading file: table3'
do i__=1,7
read(1,'(A74)')ar__
read(ar__,1)Num,Arg,Freq,Phas,Per
if(ar__(59:74) .EQ. '') Per = rNULL__
c ..............Just test output...........
write(6,1)Num,Arg,Freq,Phas,Per
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4' ! Series for element p of Hyperion
C Format for file interpretation
2 format(I3,D25.15,I4,I4,I4,I4,I4,I4,I4,F14.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4')
write(6,*) '....Loading file: table4'
do i__=1,105
read(1,'(A70)')ar__1
read(ar__1,2)Num_1,Ampl,N1,N2,N3,N4,N5,N6,N7,Per_1
if(ar__1(57:70) .EQ. '') Per_1 = rNULL__
c ..............Just test output...........
write(6,2)Num_1,Ampl,N1,N2,N3,N4,N5,N6,N7,Per_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5' ! Series for element q of Hyperion
C Format for file interpretation
3 format(I3,D25.15,I4,I4,I4,I4,I4,I4,I4,F14.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5')
write(6,*) '....Loading file: table5'
do i__=1,214
read(1,'(A70)')ar__2
read(ar__2,3)
+ Num_2,Ampl_1,N1_1,N2_1,N3_1,N4_1,N5_1,N6_1,N7_1,Per_2
if(ar__2(57:70) .EQ. '') Per_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ Num_2,Ampl_1,N1_1,N2_1,N3_1,N4_1,N5_1,N6_1,N7_1,Per_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6' ! Series for element z of Hyperion
C Format for file interpretation
4 format(I3,D25.15,I4,I4,I4,I4,I4,I4,I4,F14.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6')
write(6,*) '....Loading file: table6'
do i__=1,179
read(1,'(A70)')ar__3
read(ar__3,4)
+ Num_3,Ampl_2,N1_2,N2_2,N3_2,N4_2,N5_2,N6_2,N7_2,Per_3
if(ar__3(57:70) .EQ. '') Per_3 = rNULL__
c ..............Just test output...........
write(6,4)
+ Num_3,Ampl_2,N1_2,N2_2,N3_2,N4_2,N5_2,N6_2,N7_2,Per_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7' ! Series for element zeta of Hyperion
C Format for file interpretation
5 format(I3,D25.15,I4,I4,I4,I4,I4,I4,I4,F14.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7')
write(6,*) '....Loading file: table7'
do i__=1,52
read(1,'(A70)')ar__4
read(ar__4,5)
+ Num_4,Ampl_3,N1_3,N2_3,N3_3,N4_3,N5_3,N6_3,N7_3,Per_4
if(ar__4(57:70) .EQ. '') Per_4 = rNULL__
c ..............Just test output...........
write(6,5)
+ Num_4,Ampl_3,N1_3,N2_3,N3_3,N4_3,N5_3,N6_3,N7_3,Per_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8' ! Solar and short period perturbations of Hyperion
C Format for file interpretation
6 format(A3,D25.15,D25.15,D25.15,F14.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8')
write(6,*) '....Loading file: table8'
do i__=1,46
read(1,'(A92)')ar__5
read(ar__5,6)Ser,Ampl_4,Freq_1,Phas_1,Per_5
if(ar__5(79:92) .EQ. '') Per_5 = rNULL__
c ..............Just test output...........
write(6,6)Ser,Ampl_4,Freq_1,Phas_1,Per_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end