Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/161/295 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-Oct-07
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/161/295 Obs. with Differential Speckle Survey Instrument. X. (Horch+, 2021)
*================================================================================
*Observations with the Differential Speckle Survey Instrument.
*X. Preliminary orbits of K-dwarf binaries and other stars.
* Horch E.P., Broderick K.G., Casetti-Dinescu D.I., Henry T.J., Fekel F.C.,
* Muterspaugh M.W., Willmarth D.W., Winters J.G., van Belle G.T., Clark C.A.,
* Everett M.E.
* <Astron. J., 161, 295 (2021)>
* =2021AJ....161..295H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Orbits and residuals used in the scale determinations
integer*4 nr__
parameter (nr__=21) ! Number of records
character*81 ar__ ! Full-size record
integer*4 Run (nr__) ! [1/10] Run
character*10 WDS (nr__) ! Washington Double Star number
character*10 Design (nr__) ! Binary discover designation
integer*4 HIP (nr__) ! Hipparcos Catalogue number
real*8 Date (nr__) ! (yr) [2016/2020] Fractional Julian Year of orbit
real*4 PAA (nr__) ! (deg) [-0.6/0.7] Position angle, component A
* ({Delta}{theta}_A_)
real*4 SepA (nr__) ! (mas) [-3.2/3] Separation, component A ({Delta}{rho}_A_)
real*4 PAB (nr__) ! (deg) [-0.6/0.5] Position angle, component B
* ({Delta}{theta}_B_)
real*4 SepB (nr__) ! (mas) [-2.8/3] Separation, component B ({Delta}{rho}_B_)
character*19 Ref (nr__) ! Orbit Reference (bibcode)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Binary star Speckle measures
integer*4 nr__1
parameter (nr__1=378) ! Number of records
character*81 ar__1 ! Full-size record
character*10 WDS_1 (nr__1) ! Washington Double Star number
character*11 OName (nr__1) ! Other designations
character*13 Design_1 (nr__1) ! Binary discover designation
character*6 HIP_1 (nr__1) ! Hipparcos Catalogue number
real*8 Date_1 (nr__1) ! (yr) [16/20] Julian year of the observation;
* offset 2000+
real*4 PA (nr__1) ! (deg) [2.2/359] Position angle of the secondary star
* relative to the primary (1)
real*4 Sep (nr__1) ! (arcsec) [0.01/3.24] Separation of the two stars
character*1 l_dmag (nr__1) ! [<] Limit flag on dmag
real*4 dmag (nr__1) ! (mag) [-0.04/7.44] Magnitude difference of the pair in
* the filter used
integer*4 lambda (nr__1) ! (nm) [562/880] Center wavelength of the filter used
integer*4 FWHM (nr__1) ! (nm) [40/50] Full width at half maximum of the filter
* transmission
character*1 Flag (nr__1) ! catalog entry flags (2)
*Note (1): North through East defining the positive sense of PA
*Note (2): Flags as follows:
* a = Photometry for this observation appears as an upper limit because the
* observation may be affected by speckle decorrelation as discussed in
* the text (115 occurrences)
* b = Quadrant inconsistent with other published measures (28 occurrences)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Visual orbital elements for 25 systems
integer*4 nr__2
parameter (nr__2=25) ! Number of records
character*138 ar__2 ! Full-size record
character*10 WDS_2 (nr__2) ! Washington Double Star number
character*12 Design_2 (nr__2) ! Binary discover designation
integer*4 HIP_2 (nr__2) ! Hipparcos Catalogue number
character*7 SpT (nr__2) ! Spectral Type from SIMBAD
real*8 Per (nr__2) ! (yr) [3.72/209] Orbital period
real*4 a (nr__2) ! (arcsec) [0.08/1.2] Semi-major axis
real*8 i (nr__2) ! (deg) [22.9/166] Orbital inclination
real*4 Omega (nr__2) ! (deg) [7.7/341] Longitude of ascending node, {Omega}
real*8 T0 (nr__2) ! (yr) [1999/2188] Time of periastron passage
real*4 e (nr__2) ! [0.1/0.9] Eccentricity
real*4 omega_1 (nr__2) ! (deg) [25.8/386] Argument of periastron, {omega}
real*4 e_Per (nr__2) ! (yr) [0.01/18] Uncertainty in Per
real*4 e_a (nr__2) ! (arcsec) [0.001/0.05] Uncertainty in a
real*4 e_i (nr__2) ! (deg) [0.07/14] Uncertainty in i
real*4 e_Omega (nr__2) ! (deg) [0.3/19] Uncertainty in Omega
real*4 e_T0 (nr__2) ! (yr) [0.02/21] Uncertainty in T0
real*4 e_e (nr__2) ! [0.007/0.2] Uncertainty in e
real*4 e_omega_1 (nr__2) ! (deg) [0.2/40] Uncertainty in omega
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Radial velocity of HD173093
integer*4 nr__3
parameter (nr__3=324) ! Number of records
character*99 ar__3 ! Full-size record
real*8 MJD (nr__3) ! (d) [52073/59021] Modified Julian Date; MJD-2400000.5
real*8 RVAa (nr__3) ! (km/s) [-82.9/-19.2]? Measured Radial Velocity of the Aa
* Component
real*8 ModRVAa (nr__3) ! (km/s) [-82.3/-19.1]? Model Radial Velocity of the Aa
* Component
real*4 OCAa (nr__3) ! (km/s) [-2.3/2]? Observed minus Calculated Residual for
* the Aa Component
real*4 SigmaAa (nr__3) ! (km/s) [0.37/1.1]? Sigma of the Aa Component
real*8 RVAb (nr__3) ! (km/s) [-82.6/-17.3]? Measured Radial Velocity of the Ab
* Component
real*8 ModRVAb (nr__3) ! (km/s) [-83.1/-17.7]? Model Radial Velocity of the Ab
* Component
real*4 OCAb (nr__3) ! (km/s) [-2.1/2]? Observed minus Calculated Residual for
* the Ab Component
real*4 SigmaAb (nr__3) ! (km/s) [0.42/0.88]? Sigma of the Ab Component
real*8 RVB (nr__3) ! (km/s) [-66.7/-23.9]? Measured Radial Velocity of the B
* Component
real*8 ModRVB (nr__3) ! (km/s) [-66.3/-25.3]? Model Radial Velocity of the B
* Component
real*4 OCB (nr__3) ! (km/s) [-2.7/1.6]? Observed minus Calculated Residual for
* the B Component
real*4 SigmaB (nr__3) ! (km/s) [0.52/1.52]? Sigma of the B Component
integer*4 Set (nr__3) ! [1/3] Observational Set (1)
*Note (1): Sets are drawn from the following sources. Blank entries in
* Set 1 indicate that the velocity for a given component was not used in
* the fit due to extensive line blending. Set 1 = KPNO velocities
* obtained by D. Willmarth; Set 2 = KPNO velocities obtained by F.
* Fekel; Set 3 = Fairborn Observatory.
C=============================================================================
C Loading file 'table2.dat' ! Orbits and residuals used in the scale determinations
C Format for file interpretation
1 format(
+ I2,1X,A10,1X,A10,1X,I6,1X,F9.4,1X,F4.1,1X,F4.1,1X,F4.1,1X,
+ F4.1,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,21
read(1,'(A81)')ar__
read(ar__,1)
+ Run(i__),WDS(i__),Design(i__),HIP(i__),Date(i__),PAA(i__),
+ SepA(i__),PAB(i__),SepB(i__),Ref(i__)
c ..............Just test output...........
write(6,1)
+ Run(i__),WDS(i__),Design(i__),HIP(i__),Date(i__),PAA(i__),
+ SepA(i__),PAB(i__),SepB(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Binary star Speckle measures
C Format for file interpretation
2 format(
+ A10,1X,A11,1X,A13,2X,A6,1X,F7.4,1X,F5.1,1X,F6.4,1X,A1,F5.2,1X,
+ I3,1X,I2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,378
read(1,'(A81)')ar__1
read(ar__1,2)
+ WDS_1(i__),OName(i__),Design_1(i__),HIP_1(i__),Date_1(i__),
+ PA(i__),Sep(i__),l_dmag(i__),dmag(i__),lambda(i__),FWHM(i__),
+ Flag(i__)
c ..............Just test output...........
write(6,2)
+ WDS_1(i__),OName(i__),Design_1(i__),HIP_1(i__),Date_1(i__),
+ PA(i__),Sep(i__),l_dmag(i__),dmag(i__),lambda(i__),FWHM(i__),
+ Flag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Visual orbital elements for 25 systems
C Format for file interpretation
3 format(
+ A10,1X,A12,1X,I6,1X,A7,1X,F7.3,1X,F6.4,1X,F7.3,1X,F6.2,1X,
+ F8.3,1X,F6.4,1X,F6.2,1X,F6.3,1X,F6.4,1X,F6.3,1X,F5.2,1X,F6.3,
+ 1X,F6.4,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,25
read(1,'(A138)')ar__2
read(ar__2,3)
+ WDS_2(i__),Design_2(i__),HIP_2(i__),SpT(i__),Per(i__),a(i__),
+ i(i__),Omega(i__),T0(i__),e(i__),omega_1(i__),e_Per(i__),
+ e_a(i__),e_i(i__),e_Omega(i__),e_T0(i__),e_e(i__),
+ e_omega_1(i__)
c ..............Just test output...........
write(6,3)
+ WDS_2(i__),Design_2(i__),HIP_2(i__),SpT(i__),Per(i__),a(i__),
+ i(i__),Omega(i__),T0(i__),e(i__),omega_1(i__),e_Per(i__),
+ e_a(i__),e_i(i__),e_Omega(i__),e_T0(i__),e_e(i__),
+ e_omega_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Radial velocity of HD173093
C Format for file interpretation
4 format(
+ F10.4,1X,F7.3,1X,F7.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F7.3,1X,F6.3,
+ 1X,F5.3,1X,F7.3,1X,F7.3,1X,F6.3,1X,F5.3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,324
read(1,'(A99)')ar__3
read(ar__3,4)
+ MJD(i__),RVAa(i__),ModRVAa(i__),OCAa(i__),SigmaAa(i__),
+ RVAb(i__),ModRVAb(i__),OCAb(i__),SigmaAb(i__),RVB(i__),
+ ModRVB(i__),OCB(i__),SigmaB(i__),Set(i__)
if(ar__3(12:18) .EQ. '') RVAa(i__) = rNULL__
if(ar__3(20:26) .EQ. '') ModRVAa(i__) = rNULL__
if(ar__3(28:33) .EQ. '') OCAa(i__) = rNULL__
if(ar__3(35:39) .EQ. '') SigmaAa(i__) = rNULL__
if(ar__3(41:47) .EQ. '') RVAb(i__) = rNULL__
if(ar__3(49:55) .EQ. '') ModRVAb(i__) = rNULL__
if(ar__3(57:62) .EQ. '') OCAb(i__) = rNULL__
if(ar__3(64:68) .EQ. '') SigmaAb(i__) = rNULL__
if(ar__3(70:76) .EQ. '') RVB(i__) = rNULL__
if(ar__3(78:84) .EQ. '') ModRVB(i__) = rNULL__
if(ar__3(86:91) .EQ. '') OCB(i__) = rNULL__
if(ar__3(93:97) .EQ. '') SigmaB(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ MJD(i__),RVAa(i__),ModRVAa(i__),OCAa(i__),SigmaAa(i__),
+ RVAb(i__),ModRVAb(i__),OCAb(i__),SigmaAb(i__),RVB(i__),
+ ModRVB(i__),OCB(i__),SigmaB(i__),Set(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end