Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/158/222 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-May-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/AJ/158/222 Orbits of subsystems in multiple stars. VI. (Tokovinin, 2019)
*================================================================================
*Spectroscopic orbits of subsystems in multiple stars. VI.
* Tokovinin A.
* <Astron. J., 158, 222 (2019)>
* =2019AJ....158..222T (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Basic parameters of observed multiple systems
integer*4 nr__
parameter (nr__=20) ! Number of records
character*70 ar__ ! Full-size record
character*10 WDS ! WDS identification (HHMMm+DDMM, J2000)
character*2 Comp ! Component(s) designation
integer*4 HIP ! [14194/111598] Hipparcos number
integer*4 HD ! [18795/214169] HD number
character*6 SpType ! Spectral type
real*4 Vmag ! (mag) [5.63/11.97] V band magnitude
real*4 V_K ! (mag) [0.19/2.89]? V-K color index
integer*4 pmRA ! (mas/yr) [-279/166] Proper motion in right ascension
* direction (pmRA*cosDE) (1)
integer*4 pmDE ! (mas/yr) [-124/172] Proper motion in declination
* direction (1)
real*4 RV ! (km/s) [-89.87/36.12] Radial velocity
character*1 u_RV ! [:] Uncertainty flag on RV
real*4 Plx ! (mas) [6.17/29.67] Parallax (1)
character*1 n_Plx ! [b] Note on Plx (2)
*Note (1): Proper motions and parallaxes are taken from the Gaia DR2 (Gaia
* Collaboration et al. 2018, Cat. I/345), where available.
*Note (2): Note as follows:
* b = Hipparcos parallax (van Leeuwen 2007, Cat. I/311).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Cross-correlation function (CCF) parameters
integer*4 nr__1
parameter (nr__1=29) ! Number of records
character*40 ar__1 ! Full-size record
integer*4 HIP_1 ! [14194/111598] Hipparcos number
character*2 Comp_1 ! Component(s) designation
integer*4 N ! [1/24] Number of averaged measurements
real*4 a ! [0.014/0.522] Dip amplitude
real*4 sigma ! (km/s) [3.168/24.05] Dispersion {sigma} of a
real*4 asigma ! (km/s) [0.043/2.188] Product a{sigma} proportional
* to the dip area
real*4 Vsini ! (km/s) [0/43] Projected rotation velocity
character*1 u_Vsini ! [:] Uncertainty flag on Vsini
character*2 Li ! Presence or absence of the lithium 6708 {AA}
* line in individual components
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Spectroscopic orbits
integer*4 nr__2
parameter (nr__2=13) ! Number of records
character*139 ar__2 ! Full-size record
integer*4 HIP_2 ! [14194/111598] Hipparcos number
character*5 Syst ! System designation
real*8 Per ! (d) [2.19189/1136] Period
real*8 e_Per ! (d) [1e-05/5.8] Uncertainty in Per
real*8 T ! (d) [53759.4/58560.7] Julian Date (JD-2400000)
real*4 e_T ! (d) [0.0007/9.11] Uncertainty in T
real*4 e ! [0/0.575] Eccentricity
real*4 e_e ! [0.0005/0.038]? Uncertainty in e
character*1 f_e ! [f] f for fixed value
real*4 omegaA ! (deg) [0/308.1] Longitude of periastron for the
* primary component {omega}_A_
real*4 e_omegaA ! (deg) [0.01/2.1]? Uncertainty in omegaA
character*1 f_omegaA ! [f] f for fixed
real*4 K1 ! (km/s) [10.027/67.484] Radial velocity amplitude of
* the primary component
real*4 e_K1 ! (km/s) [0.0039/0.754] Uncertainty in K1
real*4 K2 ! (km/s) [10.241/71.379]? Radial velocity amplitude of
* the secondary component
real*4 e_K2 ! (km/s) [0.026/1.874]? Uncertainty in K2
real*8 gamma ! (km/s) [-89.867/36.121]? Center-of-mass velocity
* {gamma}
real*4 e_gamma ! (km/s) [0.01/0.341]? Uncertainty in gamma
real*4 rms1 ! (km/s) [0.03/1.67]? Weighted rms residual for the
* primary component
real*4 rms2 ! (km/s) [0.03/0.37]? Weighted rms residual for the
* secondary component
real*4 M1sin3i ! (Msun) [0.27/1.49] Mass of the primary for double-lined
* system, M_1_sin^3^i
character*1 u_M1sin3i ! [:] Uncertainty flag in M1sin3i
character*1 l_M2sin3i ! [>] Limit flag on M2sin3i
real*4 M2sin3i ! (Msun) [0.27/1.28] Mass of the secondary M_2_sin^3^i
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Radial velocities and residuals
integer*4 nr__3
parameter (nr__3=337) ! Number of records
character*45 ar__3 ! Full-size record
integer*4 HIP_3 ! [14194/111598] Hipparcos number
character*5 Syst_1 ! System designation
real*8 Date ! (d) [50301.354/58702.7466] Julian Date (JD-2400000)
real*8 RV_1 ! (km/s) [-116.04/94.47] Radial velocity
real*4 e_RV ! (km/s) [0.05/20] RV error
real*4 O_C ! (km/s) [-9.04/6.14] Residual to orbit
character*1 Comp_2 ! [abc] Component label (a=primary, b=secondary,
* c=tertiary)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Radial velocities of other components
integer*4 nr__4
parameter (nr__4=21) ! Number of records
character*26 ar__4 ! Full-size record
integer*4 HIP_4 ! [40523/66438] Hipparcos number
character*1 Comp_3 ! [AB] Component designation
real*8 Date_1 ! (d) [58193.7/58699.6] Julian Date (JD-2400000)
real*4 RV_2 ! (km/s) [-6.4/30.203] Radial velocity
character*1 u_RV_1 ! [:] Uncertainty flag on RV
C=============================================================================
C Loading file 'table1.dat' ! Basic parameters of observed multiple systems
C Format for file interpretation
1 format(
+ A10,1X,A2,1X,I6,1X,I6,1X,A6,1X,F5.2,1X,F4.2,1X,I4,1X,I4,1X,
+ F6.2,A1,1X,F5.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,20
read(1,'(A70)')ar__
read(ar__,1)
+ WDS,Comp,HIP,HD,SpType,Vmag,V_K,pmRA,pmDE,RV,u_RV,Plx,n_Plx
if(ar__(42:45) .EQ. '') V_K = rNULL__
c ..............Just test output...........
write(6,1)
+ WDS,Comp,HIP,HD,SpType,Vmag,V_K,pmRA,pmDE,RV,u_RV,Plx,n_Plx
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Cross-correlation function (CCF) parameters
C Format for file interpretation
2 format(I6,1X,A2,1X,I2,1X,F5.3,1X,F6.3,1X,F5.3,1X,F4.1,A1,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,29
read(1,'(A40)')ar__1
read(ar__1,2)HIP_1,Comp_1,N,a,sigma,asigma,Vsini,u_Vsini,Li
c ..............Just test output...........
write(6,2)HIP_1,Comp_1,N,a,sigma,asigma,Vsini,u_Vsini,Li
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Spectroscopic orbits
C Format for file interpretation
3 format(
+ I6,1X,A5,1X,F10.5,1X,F7.5,1X,F10.4,1X,F6.4,1X,F6.4,1X,F6.4,A1,
+ 1X,F5.1,1X,F4.2,A1,1X,F6.3,1X,F6.4,1X,F6.3,1X,F5.3,1X,F7.3,1X,
+ F5.3,1X,F4.2,1X,F4.2,1X,F4.2,A1,1X,A1,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,13
read(1,'(A139)')ar__2
read(ar__2,3)
+ HIP_2,Syst,Per,e_Per,T,e_T,e,e_e,f_e,omegaA,e_omegaA,f_omegaA,
+ K1,e_K1,K2,e_K2,gamma,e_gamma,rms1,rms2,M1sin3i,u_M1sin3i,
+ l_M2sin3i,M2sin3i
if(ar__2(58:63) .EQ. '') e_e = rNULL__
if(ar__2(72:75) .EQ. '') e_omegaA = rNULL__
if(ar__2(92:97) .EQ. '') K2 = rNULL__
if(ar__2(99:103) .EQ. '') e_K2 = rNULL__
if(ar__2(105:111) .EQ. '') gamma = rNULL__
if(ar__2(113:117) .EQ. '') e_gamma = rNULL__
if(ar__2(119:122) .EQ. '') rms1 = rNULL__
if(ar__2(124:127) .EQ. '') rms2 = rNULL__
c ..............Just test output...........
write(6,3)
+ HIP_2,Syst,Per,e_Per,T,e_T,e,e_e,f_e,omegaA,e_omegaA,f_omegaA,
+ K1,e_K1,K2,e_K2,gamma,e_gamma,rms1,rms2,M1sin3i,u_M1sin3i,
+ l_M2sin3i,M2sin3i
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Radial velocities and residuals
C Format for file interpretation
4 format(I6,1X,A5,1X,F10.4,1X,F7.2,1X,F5.2,1X,F5.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,337
read(1,'(A45)')ar__3
read(ar__3,4)HIP_3,Syst_1,Date,RV_1,e_RV,O_C,Comp_2
c ..............Just test output...........
write(6,4)HIP_3,Syst_1,Date,RV_1,e_RV,O_C,Comp_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Radial velocities of other components
C Format for file interpretation
5 format(I5,1X,A1,1X,F10.4,1X,F6.3,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,21
read(1,'(A26)')ar__4
read(ar__4,5)HIP_4,Comp_3,Date_1,RV_2,u_RV_1
c ..............Just test output...........
write(6,5)HIP_4,Comp_3,Date_1,RV_2,u_RV_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end