Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/158/222 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-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 (nr__) ! WDS identification (HHMMm+DDMM, J2000)
character*2 Comp (nr__) ! Component(s) designation
integer*4 HIP (nr__) ! [14194/111598] Hipparcos number
integer*4 HD (nr__) ! [18795/214169] HD number
character*6 SpType (nr__) ! Spectral type
real*4 Vmag (nr__) ! (mag) [5.63/11.97] V band magnitude
real*4 V_K (nr__) ! (mag) [0.19/2.89]? V-K color index
integer*4 pmRA (nr__) ! (mas/yr) [-279/166] Proper motion in right ascension
* direction (pmRA*cosDE) (1)
integer*4 pmDE (nr__) ! (mas/yr) [-124/172] Proper motion in declination
* direction (1)
real*4 RV (nr__) ! (km/s) [-89.87/36.12] Radial velocity
character*1 u_RV (nr__) ! [:] Uncertainty flag on RV
real*4 Plx (nr__) ! (mas) [6.17/29.67] Parallax (1)
character*1 n_Plx (nr__) ! [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 (nr__1) ! [14194/111598] Hipparcos number
character*2 Comp_1 (nr__1) ! Component(s) designation
integer*4 N (nr__1) ! [1/24] Number of averaged measurements
real*4 a (nr__1) ! [0.014/0.522] Dip amplitude
real*4 sigma (nr__1) ! (km/s) [3.168/24.05] Dispersion {sigma} of a
real*4 asigma (nr__1) ! (km/s) [0.043/2.188] Product a{sigma} proportional
* to the dip area
real*4 Vsini (nr__1) ! (km/s) [0/43] Projected rotation velocity
character*1 u_Vsini (nr__1) ! [:] Uncertainty flag on Vsini
character*2 Li (nr__1) ! 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 (nr__2) ! [14194/111598] Hipparcos number
character*5 Syst (nr__2) ! System designation
real*8 Per (nr__2) ! (d) [2.19189/1136] Period
real*8 e_Per (nr__2) ! (d) [1e-05/5.8] Uncertainty in Per
real*8 T (nr__2) ! (d) [53759.4/58560.7] Julian Date (JD-2400000)
real*4 e_T (nr__2) ! (d) [0.0007/9.11] Uncertainty in T
real*4 e (nr__2) ! [0/0.575] Eccentricity
real*4 e_e (nr__2) ! [0.0005/0.038]? Uncertainty in e
character*1 f_e (nr__2) ! [f] f for fixed value
real*4 omegaA (nr__2) ! (deg) [0/308.1] Longitude of periastron for the
* primary component {omega}_A_
real*4 e_omegaA (nr__2) ! (deg) [0.01/2.1]? Uncertainty in omegaA
character*1 f_omegaA (nr__2) ! [f] f for fixed
real*4 K1 (nr__2) ! (km/s) [10.027/67.484] Radial velocity amplitude of
* the primary component
real*4 e_K1 (nr__2) ! (km/s) [0.0039/0.754] Uncertainty in K1
real*4 K2 (nr__2) ! (km/s) [10.241/71.379]? Radial velocity amplitude of
* the secondary component
real*4 e_K2 (nr__2) ! (km/s) [0.026/1.874]? Uncertainty in K2
real*8 gamma (nr__2) ! (km/s) [-89.867/36.121]? Center-of-mass velocity
* {gamma}
real*4 e_gamma (nr__2) ! (km/s) [0.01/0.341]? Uncertainty in gamma
real*4 rms1 (nr__2) ! (km/s) [0.03/1.67]? Weighted rms residual for the
* primary component
real*4 rms2 (nr__2) ! (km/s) [0.03/0.37]? Weighted rms residual for the
* secondary component
real*4 M1sin3i (nr__2) ! (Msun) [0.27/1.49] Mass of the primary for double-lined
* system, M_1_sin^3^i
character*1 u_M1sin3i (nr__2) ! [:] Uncertainty flag in M1sin3i
character*1 l_M2sin3i (nr__2) ! [>] Limit flag on M2sin3i
real*4 M2sin3i (nr__2) ! (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 (nr__3) ! [14194/111598] Hipparcos number
character*5 Syst_1 (nr__3) ! System designation
real*8 Date (nr__3) ! (d) [50301.354/58702.7466] Julian Date (JD-2400000)
real*8 RV_1 (nr__3) ! (km/s) [-116.04/94.47] Radial velocity
real*4 e_RV (nr__3) ! (km/s) [0.05/20] RV error
real*4 O_C (nr__3) ! (km/s) [-9.04/6.14] Residual to orbit
character*1 Comp_2 (nr__3) ! [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 (nr__4) ! [40523/66438] Hipparcos number
character*1 Comp_3 (nr__4) ! [AB] Component designation
real*8 Date_1 (nr__4) ! (d) [58193.7/58699.6] Julian Date (JD-2400000)
real*4 RV_2 (nr__4) ! (km/s) [-6.4/30.203] Radial velocity
character*1 u_RV_1 (nr__4) ! [:] 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(i__),Comp(i__),HIP(i__),HD(i__),SpType(i__),Vmag(i__),
+ V_K(i__),pmRA(i__),pmDE(i__),RV(i__),u_RV(i__),Plx(i__),
+ n_Plx(i__)
if(ar__(42:45) .EQ. '') V_K(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ WDS(i__),Comp(i__),HIP(i__),HD(i__),SpType(i__),Vmag(i__),
+ V_K(i__),pmRA(i__),pmDE(i__),RV(i__),u_RV(i__),Plx(i__),
+ n_Plx(i__)
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(i__),Comp_1(i__),N(i__),a(i__),sigma(i__),asigma(i__),
+ Vsini(i__),u_Vsini(i__),Li(i__)
c ..............Just test output...........
write(6,2)
+ HIP_1(i__),Comp_1(i__),N(i__),a(i__),sigma(i__),asigma(i__),
+ Vsini(i__),u_Vsini(i__),Li(i__)
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(i__),Syst(i__),Per(i__),e_Per(i__),T(i__),e_T(i__),
+ e(i__),e_e(i__),f_e(i__),omegaA(i__),e_omegaA(i__),
+ f_omegaA(i__),K1(i__),e_K1(i__),K2(i__),e_K2(i__),gamma(i__),
+ e_gamma(i__),rms1(i__),rms2(i__),M1sin3i(i__),u_M1sin3i(i__),
+ l_M2sin3i(i__),M2sin3i(i__)
if(ar__2(58:63) .EQ. '') e_e(i__) = rNULL__
if(ar__2(72:75) .EQ. '') e_omegaA(i__) = rNULL__
if(ar__2(92:97) .EQ. '') K2(i__) = rNULL__
if(ar__2(99:103) .EQ. '') e_K2(i__) = rNULL__
if(ar__2(105:111) .EQ. '') gamma(i__) = rNULL__
if(ar__2(113:117) .EQ. '') e_gamma(i__) = rNULL__
if(ar__2(119:122) .EQ. '') rms1(i__) = rNULL__
if(ar__2(124:127) .EQ. '') rms2(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ HIP_2(i__),Syst(i__),Per(i__),e_Per(i__),T(i__),e_T(i__),
+ e(i__),e_e(i__),f_e(i__),omegaA(i__),e_omegaA(i__),
+ f_omegaA(i__),K1(i__),e_K1(i__),K2(i__),e_K2(i__),gamma(i__),
+ e_gamma(i__),rms1(i__),rms2(i__),M1sin3i(i__),u_M1sin3i(i__),
+ l_M2sin3i(i__),M2sin3i(i__)
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(i__),Syst_1(i__),Date(i__),RV_1(i__),e_RV(i__),O_C(i__),
+ Comp_2(i__)
c ..............Just test output...........
write(6,4)
+ HIP_3(i__),Syst_1(i__),Date(i__),RV_1(i__),e_RV(i__),O_C(i__),
+ Comp_2(i__)
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(i__),Comp_3(i__),Date_1(i__),RV_2(i__),u_RV_1(i__)
c ..............Just test output...........
write(6,5)
+ HIP_4(i__),Comp_3(i__),Date_1(i__),RV_2(i__),u_RV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end