Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/161/295 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-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 ! [1/10] Run
character*10 WDS ! Washington Double Star number
character*10 Design ! Binary discover designation
integer*4 HIP ! Hipparcos Catalogue number
real*8 Date ! (yr) [2016/2020] Fractional Julian Year of orbit
real*4 PAA ! (deg) [-0.6/0.7] Position angle, component A
* ({Delta}{theta}_A_)
real*4 SepA ! (mas) [-3.2/3] Separation, component A ({Delta}{rho}_A_)
real*4 PAB ! (deg) [-0.6/0.5] Position angle, component B
* ({Delta}{theta}_B_)
real*4 SepB ! (mas) [-2.8/3] Separation, component B ({Delta}{rho}_B_)
character*19 Ref ! 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 ! Washington Double Star number
character*11 OName ! Other designations
character*13 Design_1 ! Binary discover designation
character*6 HIP_1 ! Hipparcos Catalogue number
real*8 Date_1 ! (yr) [16/20] Julian year of the observation;
* offset 2000+
real*4 PA ! (deg) [2.2/359] Position angle of the secondary star
* relative to the primary (1)
real*4 Sep ! (arcsec) [0.01/3.24] Separation of the two stars
character*1 l_dmag ! [<] Limit flag on dmag
real*4 dmag ! (mag) [-0.04/7.44] Magnitude difference of the pair in
* the filter used
integer*4 lambda ! (nm) [562/880] Center wavelength of the filter used
integer*4 FWHM ! (nm) [40/50] Full width at half maximum of the filter
* transmission
character*1 Flag ! 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 ! Washington Double Star number
character*12 Design_2 ! Binary discover designation
integer*4 HIP_2 ! Hipparcos Catalogue number
character*7 SpT ! Spectral Type from SIMBAD
real*8 Per ! (yr) [3.72/209] Orbital period
real*4 a ! (arcsec) [0.08/1.2] Semi-major axis
real*8 i ! (deg) [22.9/166] Orbital inclination
real*4 Omega ! (deg) [7.7/341] Longitude of ascending node, {Omega}
real*8 T0 ! (yr) [1999/2188] Time of periastron passage
real*4 e ! [0.1/0.9] Eccentricity
real*4 omega_1 ! (deg) [25.8/386] Argument of periastron, {omega}
real*4 e_Per ! (yr) [0.01/18] Uncertainty in Per
real*4 e_a ! (arcsec) [0.001/0.05] Uncertainty in a
real*4 e_i ! (deg) [0.07/14] Uncertainty in i
real*4 e_Omega ! (deg) [0.3/19] Uncertainty in Omega
real*4 e_T0 ! (yr) [0.02/21] Uncertainty in T0
real*4 e_e ! [0.007/0.2] Uncertainty in e
real*4 e_omega_1 ! (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 ! (d) [52073/59021] Modified Julian Date; MJD-2400000.5
real*8 RVAa ! (km/s) [-82.9/-19.2]? Measured Radial Velocity of the Aa
* Component
real*8 ModRVAa ! (km/s) [-82.3/-19.1]? Model Radial Velocity of the Aa
* Component
real*4 OCAa ! (km/s) [-2.3/2]? Observed minus Calculated Residual for
* the Aa Component
real*4 SigmaAa ! (km/s) [0.37/1.1]? Sigma of the Aa Component
real*8 RVAb ! (km/s) [-82.6/-17.3]? Measured Radial Velocity of the Ab
* Component
real*8 ModRVAb ! (km/s) [-83.1/-17.7]? Model Radial Velocity of the Ab
* Component
real*4 OCAb ! (km/s) [-2.1/2]? Observed minus Calculated Residual for
* the Ab Component
real*4 SigmaAb ! (km/s) [0.42/0.88]? Sigma of the Ab Component
real*8 RVB ! (km/s) [-66.7/-23.9]? Measured Radial Velocity of the B
* Component
real*8 ModRVB ! (km/s) [-66.3/-25.3]? Model Radial Velocity of the B
* Component
real*4 OCB ! (km/s) [-2.7/1.6]? Observed minus Calculated Residual for
* the B Component
real*4 SigmaB ! (km/s) [0.52/1.52]? Sigma of the B Component
integer*4 Set ! [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,WDS,Design,HIP,Date,PAA,SepA,PAB,SepB,Ref
c ..............Just test output...........
write(6,1)Run,WDS,Design,HIP,Date,PAA,SepA,PAB,SepB,Ref
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,OName,Design_1,HIP_1,Date_1,PA,Sep,l_dmag,dmag,lambda,
+ FWHM,Flag
c ..............Just test output...........
write(6,2)
+ WDS_1,OName,Design_1,HIP_1,Date_1,PA,Sep,l_dmag,dmag,lambda,
+ FWHM,Flag
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,Design_2,HIP_2,SpT,Per,a,i,Omega,T0,e,omega_1,e_Per,e_a,
+ e_i,e_Omega,e_T0,e_e,e_omega_1
c ..............Just test output...........
write(6,3)
+ WDS_2,Design_2,HIP_2,SpT,Per,a,i,Omega,T0,e,omega_1,e_Per,e_a,
+ e_i,e_Omega,e_T0,e_e,e_omega_1
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,RVAa,ModRVAa,OCAa,SigmaAa,RVAb,ModRVAb,OCAb,SigmaAb,RVB,
+ ModRVB,OCB,SigmaB,Set
if(ar__3(12:18) .EQ. '') RVAa = rNULL__
if(ar__3(20:26) .EQ. '') ModRVAa = rNULL__
if(ar__3(28:33) .EQ. '') OCAa = rNULL__
if(ar__3(35:39) .EQ. '') SigmaAa = rNULL__
if(ar__3(41:47) .EQ. '') RVAb = rNULL__
if(ar__3(49:55) .EQ. '') ModRVAb = rNULL__
if(ar__3(57:62) .EQ. '') OCAb = rNULL__
if(ar__3(64:68) .EQ. '') SigmaAb = rNULL__
if(ar__3(70:76) .EQ. '') RVB = rNULL__
if(ar__3(78:84) .EQ. '') ModRVB = rNULL__
if(ar__3(86:91) .EQ. '') OCB = rNULL__
if(ar__3(93:97) .EQ. '') SigmaB = rNULL__
c ..............Just test output...........
write(6,4)
+ MJD,RVAa,ModRVAa,OCAa,SigmaAa,RVAb,ModRVAb,OCAb,SigmaAb,RVB,
+ ModRVB,OCB,SigmaB,Set
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end