Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/160/169 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-16
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/160/169 WOCS. LXXXII. Orbital parameters & RVs in NGC 7789 (Nine+, 2020)
*================================================================================
*WIYN open cluster study.
*LXXXII. Radial-velocity measurements and spectroscopic binary orbits in the
*open cluster NGC 7789.
* Nine A.C., Milliman K.E., Mathieu R.D., Geller A.M., Leiner E.M.,
* Platais I., Tofflemire B.M.
* <Astron. J., 160, 169 (2020)>
* =2020AJ....160..169N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Radial-velocity measurements
integer*4 nr__
parameter (nr__=8275) ! Number of records
character*62 ar__ ! Full-size record
integer*4 ID ! [1001/72035] Identifier
real*8 HJD ! (d) [7565/58643] Heliocentric Julian Date, HJD-2400000
real*4 RV1 ! (km/s) [-162.5/101.5] Primary radial velocity
real*8 CH1 ! [0.4/5200000] Primary Correlation Height
real*4 O_C1 ! (km/s) [-8.6/9.8]? Primary residual
real*4 RV2 ! (km/s) [-171.7/71.1]? Secondary radial velocity
real*4 CH2 ! [0.14/0.9]? Secondary Correlation Height
real*4 O_C2 ! (km/s) [-8.6/9.8]? Secondary residual
real*4 Phase ! [0.002/0.998]? Orbital phase
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Radial-velocity summary table
integer*4 nr__1
parameter (nr__1=1206) ! Number of records
character*116 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID_1 ! [1001/107032] Identifier
integer*4 RAh ! (h) [23] Hour of Right Ascension (J2000)
integer*4 RAm ! (min) [55/59] Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! [+] Sign of declination (J2000)
integer*4 DEd ! (deg) [56/57] Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 Gmag ! (mag) [7.43/14.99]? Gaia DR2 G band magnitude
real*4 GBP_GRP ! (mag) [0.28/4.36]? Gaia DR2 G_BP_-G_RP_ color
integer*4 N ! [1/49]? Number of observations
real*4 RVel ! (km/s) [-162.3/54.4]? Mean radial velocity
real*4 e_RVel ! (km/s) [0/45]? Standard uncertainty in <RVel>
real*4 e_i ! [0/112.5]? Eccentricity to inclination ratio
integer*4 PPM ! (%) [0/99]? Proper-motion membership probability
integer*4 PRV ! (%) [0/95]? Radial velocity membership probability
character*7 Class ! Membership classification
real*8 gamma ! (km/s) [-56.8/-51.2]? Center-of-mass velocity
real*4 e_gamma ! (km/s) [0.05/0.7]? Standard uncertainty in gamma
character*26 Com ! Additional comment(s)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Orbital parameters for NGC 7789 single-lined
binaries
integer*4 nr__2
parameter (nr__2=60) ! Number of records
character*138 ar__2 ! Full-size record
integer*4 ID_2 ! [1006/47010] Identifier
real*8 Porb ! (d) [1.88/4190] Orbital period
real*8 e_Porb ! (d) [1e-05/230] The 1{sigma} uncertainty in Porb
real*4 Cycle ! [0.9/1823] Number of orbital cycles observed
real*4 gamma_1 ! (km/s) [-56.21/-51.4] Center-of-mass radial velocity
real*4 e_gamma_1 ! (km/s) [0.05/0.5] Standard uncertainty in gamma
real*4 K ! (km/s) [2.5/58.4] Orbital amplitude
real*4 e_K ! (km/s) [0.07/5] The 1{sigma} uncertainty in K
real*4 e ! [0.006/0.9] Eccentricity
real*4 e_e ! [0.003/0.16] The 1{sigma} uncertainty in e
real*4 omega ! [6.0/348] Longitude of periastron
real*4 e_omega ! [0.7/120] The 1{sigma} uncertainty in omega
real*8 T0 ! (d) [54180/58390] Heliocentric Julian Date of
* periastron passage JD-2400000
real*4 e_T0 ! (d) [0.04/160] The 1{sigma} uncertainty in T0
real*8 asini ! (10+6km) [0.364/253] Projected semi-major axis
real*4 e_asini ! (10+6km) [0.008/60] The 1{sigma} uncertainty in asini
real*4 fm ! (Msun) [0.00011/0.223] The mass function
real*4 e_fm ! (solMass) [2e-05/0.1] The 1{sigma} uncertainty in fm
real*4 sigma ! (km/s) [0.14/1.99]? RMS residual velocity (1)
integer*4 Num ! [12/44]? Number of observations (2)
*Note (1): For the primary from the orbital solution.
*Note (2): Used in the orbital solutions after excluding measurements where
* the two velocities could not be resolved.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Orbital parameters for NGC 7789 double-dined
binaries
integer*4 nr__3
parameter (nr__3=42) ! Number of records
character*141 ar__3 ! Full-size record
integer*4 ID_3 ! [7003/40008] Identifier
real*8 Porb_1 ! (d) [2.36601/191.64]? Orbital period
real*8 e_Porb_1 ! (d) [6e-06/0.23]? The 1{sigma} uncertainty in Porb
real*4 Cycle_1 ! [8.6/1753]? Number of orbital cycles observed
real*4 gamma_2 ! (km/s) [-56.5/-51.2]? Center-of-mass radial velocity
real*4 e_gamma_2 ! (km/s) [0.08/0.7]? Standard uncertainty in gamma
real*4 K_1 ! (km/s) [26.7/122.2] Orbital amplitude
real*4 e_K_1 ! (km/s) [0.1/2.2] The 1{sigma} uncertainty in K
real*4 e_1 ! [0.007/0.573]? Eccentricity
real*4 e_e_1 ! [0.002/0.04]? The 1{sigma} uncertainty in e
real*4 omega_1 ! [24/330]? Longitude of periastron
real*4 e_omega_1 ! [0.5/50]? The 1{sigma} uncertainty in omega
real*8 T0_1 ! (d) [54550/57123]? Heliocentric Julian Date of
* periastron passage JD-2400000
real*4 e_T0_1 ! (d) [0.018/2]? The 1{sigma} uncertainty in T0
real*4 asini_1 ! (10+6km) [2.199/65.5] Projected semi-major axis
real*4 e_asini_1 ! (10+6km) [0.014/1.7] The 1{sigma} uncertainty in asini
real*4 msin3i ! (solMass) [0.193/1.89] Projected mass
real*4 e_msin3i ! (solMass) [0.005/0.114] The 1{sigma} uncertainty in msin3i
real*4 q ! [0.685/1.02]? Mass ratio
real*4 e_q ! [0.006/0.05]? The 1{sigma} uncertainty in q
real*4 sigma_1 ! (km/s) [0.25/7.68]? RMS residual velocity (1)
integer*4 Num_1 ! [12/35]? Number of radial velocity measurements
* (2)
*Note (1): For the primary from the orbital solution.
*Note (2): Used in the orbital solutions after excluding measurements where
* the two velocities could not be resolved.
C=============================================================================
C Loading file 'table1.dat' ! Radial-velocity measurements
C Format for file interpretation
1 format(
+ I5,1X,F10.4,1X,F6.1,1X,F10.2,1X,F4.1,1X,F6.1,1X,F4.2,1X,F4.1,
+ 1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,8275
read(1,'(A62)')ar__
read(ar__,1)ID,HJD,RV1,CH1,O_C1,RV2,CH2,O_C2,Phase
if(ar__(36:39) .EQ. '') O_C1 = rNULL__
if(ar__(41:46) .EQ. '') RV2 = rNULL__
if(ar__(48:51) .EQ. '') CH2 = rNULL__
if(ar__(53:56) .EQ. '') O_C2 = rNULL__
if(ar__(58:62) .EQ. '') Phase = rNULL__
c ..............Just test output...........
write(6,1)ID,HJD,RV1,CH1,O_C1,RV2,CH2,O_C2,Phase
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Radial-velocity summary table
C Format for file interpretation
2 format(
+ I6,1X,I2,1X,I2,1X,F5.2,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,F4.2,1X,
+ I2,1X,F6.1,1X,F4.1,1X,F5.1,1X,I2,1X,I2,1X,A7,1X,F7.3,1X,F5.3,
+ 1X,A26)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,1206
read(1,'(A116)')ar__1
read(ar__1,2)
+ ID_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Gmag,GBP_GRP,N,RVel,e_RVel,
+ e_i,PPM,PRV,Class,gamma,e_gamma,Com
if(ar__1(31:35) .EQ. '') Gmag = rNULL__
if(ar__1(37:40) .EQ. '') GBP_GRP = rNULL__
if(ar__1(42:43) .EQ. '') N = iNULL__
if(ar__1(45:50) .EQ. '') RVel = rNULL__
if(ar__1(52:55) .EQ. '') e_RVel = rNULL__
if(ar__1(57:61) .EQ. '') e_i = rNULL__
if(ar__1(63:64) .EQ. '') PPM = iNULL__
if(ar__1(66:67) .EQ. '') PRV = iNULL__
if(ar__1(77:83) .EQ. '') gamma = rNULL__
if(ar__1(85:89) .EQ. '') e_gamma = rNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,2)
+ ID_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Gmag,GBP_GRP,N,RVel,e_RVel,
+ e_i,PPM,PRV,Class,gamma,e_gamma,Com
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Orbital parameters for NGC 7789 single-lined
* binaries
C Format for file interpretation
3 format(
+ I6,1X,F11.6,1X,F10.6,1X,F6.1,1X,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F5.3,1X,F5.3,1X,F5.1,1X,F5.1,1X,F8.2,1X,F6.2,1X,F7.3,1X,
+ F6.3,1X,E7.2,1X,E6.1,1X,F4.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,60
read(1,'(A138)')ar__2
read(ar__2,3)
+ ID_2,Porb,e_Porb,Cycle,gamma_1,e_gamma_1,K,e_K,e,e_e,omega,
+ e_omega,T0,e_T0,asini,e_asini,fm,e_fm,sigma,Num
if(ar__2(131:134) .EQ. '') sigma = rNULL__
if(ar__2(136:138) .EQ. '') Num = iNULL__
c ..............Just test output...........
write(6,3)
+ ID_2,Porb,e_Porb,Cycle,gamma_1,e_gamma_1,K,e_K,e,e_e,omega,
+ e_omega,T0,e_T0,asini,e_asini,fm,e_fm,sigma,Num
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Orbital parameters for NGC 7789 double-dined
* binaries
C Format for file interpretation
4 format(
+ I5,1X,F10.6,1X,F8.6,1X,F6.1,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,
+ 1X,F5.3,1X,F5.3,1X,F5.1,1X,F4.1,1X,F9.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F4.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,42
read(1,'(A141)')ar__3
read(ar__3,4)
+ ID_3,Porb_1,e_Porb_1,Cycle_1,gamma_2,e_gamma_2,K_1,e_K_1,e_1,
+ e_e_1,omega_1,e_omega_1,T0_1,e_T0_1,asini_1,e_asini_1,msin3i,
+ e_msin3i,q,e_q,sigma_1,Num_1
if(ar__3(7:16) .EQ. '') Porb_1 = rNULL__
if(ar__3(18:25) .EQ. '') e_Porb_1 = rNULL__
if(ar__3(27:32) .EQ. '') Cycle_1 = rNULL__
if(ar__3(34:39) .EQ. '') gamma_2 = rNULL__
if(ar__3(41:44) .EQ. '') e_gamma_2 = rNULL__
if(ar__3(58:62) .EQ. '') e_1 = rNULL__
if(ar__3(64:68) .EQ. '') e_e_1 = rNULL__
if(ar__3(70:74) .EQ. '') omega_1 = rNULL__
if(ar__3(76:79) .EQ. '') e_omega_1 = rNULL__
if(ar__3(81:89) .EQ. '') T0_1 = rNULL__
if(ar__3(91:95) .EQ. '') e_T0_1 = rNULL__
if(ar__3(122:126) .EQ. '') q = rNULL__
if(ar__3(128:132) .EQ. '') e_q = rNULL__
if(ar__3(134:137) .EQ. '') sigma_1 = rNULL__
if(ar__3(139:141) .EQ. '') Num_1 = iNULL__
c ..............Just test output...........
write(6,4)
+ ID_3,Porb_1,e_Porb_1,Cycle_1,gamma_2,e_gamma_2,K_1,e_K_1,e_1,
+ e_e_1,omega_1,e_omega_1,T0_1,e_T0_1,asini_1,e_asini_1,msin3i,
+ e_msin3i,q,e_q,sigma_1,Num_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end