Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/160/169 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-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 (nr__) ! [1001/72035] Identifier
real*8 HJD (nr__) ! (d) [7565/58643] Heliocentric Julian Date, HJD-2400000
real*4 RV1 (nr__) ! (km/s) [-162.5/101.5] Primary radial velocity
real*8 CH1 (nr__) ! [0.4/5200000] Primary Correlation Height
real*4 O_C1 (nr__) ! (km/s) [-8.6/9.8]? Primary residual
real*4 RV2 (nr__) ! (km/s) [-171.7/71.1]? Secondary radial velocity
real*4 CH2 (nr__) ! [0.14/0.9]? Secondary Correlation Height
real*4 O_C2 (nr__) ! (km/s) [-8.6/9.8]? Secondary residual
real*4 Phase (nr__) ! [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 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID_1 (nr__1) ! [1001/107032] Identifier
integer*4 RAh (nr__1) ! (h) [23] Hour of Right Ascension (J2000)
integer*4 RAm (nr__1) ! (min) [55/59] Minute of Right Ascension (J2000)
real*4 RAs (nr__1) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__1) ! [+] Sign of declination (J2000)
integer*4 DEd (nr__1) ! (deg) [56/57] Degree of Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*4 Gmag (nr__1) ! (mag) [7.43/14.99]? Gaia DR2 G band magnitude
real*4 GBP_GRP (nr__1) ! (mag) [0.28/4.36]? Gaia DR2 G_BP_-G_RP_ color
integer*4 N (nr__1) ! [1/49]? Number of observations
real*4 RVel (nr__1) ! (km/s) [-162.3/54.4]? Mean radial velocity
real*4 e_RVel (nr__1) ! (km/s) [0/45]? Standard uncertainty in <RVel>
real*4 e_i (nr__1) ! [0/112.5]? Eccentricity to inclination ratio
integer*4 PPM (nr__1) ! (%) [0/99]? Proper-motion membership probability
integer*4 PRV (nr__1) ! (%) [0/95]? Radial velocity membership probability
character*7 Class (nr__1) ! Membership classification
real*8 gamma (nr__1) ! (km/s) [-56.8/-51.2]? Center-of-mass velocity
real*4 e_gamma (nr__1) ! (km/s) [0.05/0.7]? Standard uncertainty in gamma
character*26 Com (nr__1) ! 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 (nr__2) ! [1006/47010] Identifier
real*8 Porb (nr__2) ! (d) [1.88/4190] Orbital period
real*8 e_Porb (nr__2) ! (d) [1e-05/230] The 1{sigma} uncertainty in Porb
real*4 Cycle (nr__2) ! [0.9/1823] Number of orbital cycles observed
real*4 gamma_1 (nr__2) ! (km/s) [-56.21/-51.4] Center-of-mass radial velocity
real*4 e_gamma_1 (nr__2) ! (km/s) [0.05/0.5] Standard uncertainty in gamma
real*4 K (nr__2) ! (km/s) [2.5/58.4] Orbital amplitude
real*4 e_K (nr__2) ! (km/s) [0.07/5] The 1{sigma} uncertainty in K
real*4 e (nr__2) ! [0.006/0.9] Eccentricity
real*4 e_e (nr__2) ! [0.003/0.16] The 1{sigma} uncertainty in e
real*4 omega (nr__2) ! [6.0/348] Longitude of periastron
real*4 e_omega (nr__2) ! [0.7/120] The 1{sigma} uncertainty in omega
real*8 T0 (nr__2) ! (d) [54180/58390] Heliocentric Julian Date of
* periastron passage JD-2400000
real*4 e_T0 (nr__2) ! (d) [0.04/160] The 1{sigma} uncertainty in T0
real*8 asini (nr__2) ! (10+6km) [0.364/253] Projected semi-major axis
real*4 e_asini (nr__2) ! (10+6km) [0.008/60] The 1{sigma} uncertainty in asini
real*4 fm (nr__2) ! (Msun) [0.00011/0.223] The mass function
real*4 e_fm (nr__2) ! (solMass) [2e-05/0.1] The 1{sigma} uncertainty in fm
real*4 sigma (nr__2) ! (km/s) [0.14/1.99]? RMS residual velocity (1)
integer*4 Num (nr__2) ! [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 (nr__3) ! [7003/40008] Identifier
real*8 Porb_1 (nr__3) ! (d) [2.36601/191.64]? Orbital period
real*8 e_Porb_1 (nr__3) ! (d) [6e-06/0.23]? The 1{sigma} uncertainty in Porb
real*4 Cycle_1 (nr__3) ! [8.6/1753]? Number of orbital cycles observed
real*4 gamma_2 (nr__3) ! (km/s) [-56.5/-51.2]? Center-of-mass radial velocity
real*4 e_gamma_2 (nr__3) ! (km/s) [0.08/0.7]? Standard uncertainty in gamma
real*4 K_1 (nr__3) ! (km/s) [26.7/122.2] Orbital amplitude
real*4 e_K_1 (nr__3) ! (km/s) [0.1/2.2] The 1{sigma} uncertainty in K
real*4 e_1 (nr__3) ! [0.007/0.573]? Eccentricity
real*4 e_e_1 (nr__3) ! [0.002/0.04]? The 1{sigma} uncertainty in e
real*4 omega_1 (nr__3) ! [24/330]? Longitude of periastron
real*4 e_omega_1 (nr__3) ! [0.5/50]? The 1{sigma} uncertainty in omega
real*8 T0_1 (nr__3) ! (d) [54550/57123]? Heliocentric Julian Date of
* periastron passage JD-2400000
real*4 e_T0_1 (nr__3) ! (d) [0.018/2]? The 1{sigma} uncertainty in T0
real*4 asini_1 (nr__3) ! (10+6km) [2.199/65.5] Projected semi-major axis
real*4 e_asini_1 (nr__3) ! (10+6km) [0.014/1.7] The 1{sigma} uncertainty in asini
real*4 msin3i (nr__3) ! (solMass) [0.193/1.89] Projected mass
real*4 e_msin3i (nr__3) ! (solMass) [0.005/0.114] The 1{sigma} uncertainty in msin3i
real*4 q (nr__3) ! [0.685/1.02]? Mass ratio
real*4 e_q (nr__3) ! [0.006/0.05]? The 1{sigma} uncertainty in q
real*4 sigma_1 (nr__3) ! (km/s) [0.25/7.68]? RMS residual velocity (1)
integer*4 Num_1 (nr__3) ! [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(i__),HJD(i__),RV1(i__),CH1(i__),O_C1(i__),RV2(i__),
+ CH2(i__),O_C2(i__),Phase(i__)
if(ar__(36:39) .EQ. '') O_C1(i__) = rNULL__
if(ar__(41:46) .EQ. '') RV2(i__) = rNULL__
if(ar__(48:51) .EQ. '') CH2(i__) = rNULL__
if(ar__(53:56) .EQ. '') O_C2(i__) = rNULL__
if(ar__(58:62) .EQ. '') Phase(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),HJD(i__),RV1(i__),CH1(i__),O_C1(i__),RV2(i__),
+ CH2(i__),O_C2(i__),Phase(i__)
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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Gmag(i__),GBP_GRP(i__),N(i__),RVel(i__),
+ e_RVel(i__),e_i(i__),PPM(i__),PRV(i__),Class(i__),gamma(i__),
+ e_gamma(i__),Com(i__)
if(ar__1(31:35) .EQ. '') Gmag(i__) = rNULL__
if(ar__1(37:40) .EQ. '') GBP_GRP(i__) = rNULL__
if(ar__1(42:43) .EQ. '') N(i__) = iNULL__
if(ar__1(45:50) .EQ. '') RVel(i__) = rNULL__
if(ar__1(52:55) .EQ. '') e_RVel(i__) = rNULL__
if(ar__1(57:61) .EQ. '') e_i(i__) = rNULL__
if(ar__1(63:64) .EQ. '') PPM(i__) = iNULL__
if(ar__1(66:67) .EQ. '') PRV(i__) = iNULL__
if(ar__1(77:83) .EQ. '') gamma(i__) = rNULL__
if(ar__1(85:89) .EQ. '') e_gamma(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ ID_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Gmag(i__),GBP_GRP(i__),N(i__),RVel(i__),
+ e_RVel(i__),e_i(i__),PPM(i__),PRV(i__),Class(i__),gamma(i__),
+ e_gamma(i__),Com(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),Porb(i__),e_Porb(i__),Cycle(i__),gamma_1(i__),
+ e_gamma_1(i__),K(i__),e_K(i__),e(i__),e_e(i__),omega(i__),
+ e_omega(i__),T0(i__),e_T0(i__),asini(i__),e_asini(i__),
+ fm(i__),e_fm(i__),sigma(i__),Num(i__)
if(ar__2(131:134) .EQ. '') sigma(i__) = rNULL__
if(ar__2(136:138) .EQ. '') Num(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ ID_2(i__),Porb(i__),e_Porb(i__),Cycle(i__),gamma_1(i__),
+ e_gamma_1(i__),K(i__),e_K(i__),e(i__),e_e(i__),omega(i__),
+ e_omega(i__),T0(i__),e_T0(i__),asini(i__),e_asini(i__),
+ fm(i__),e_fm(i__),sigma(i__),Num(i__)
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(i__),Porb_1(i__),e_Porb_1(i__),Cycle_1(i__),gamma_2(i__),
+ e_gamma_2(i__),K_1(i__),e_K_1(i__),e_1(i__),e_e_1(i__),
+ omega_1(i__),e_omega_1(i__),T0_1(i__),e_T0_1(i__),
+ asini_1(i__),e_asini_1(i__),msin3i(i__),e_msin3i(i__),q(i__),
+ e_q(i__),sigma_1(i__),Num_1(i__)
if(ar__3(7:16) .EQ. '') Porb_1(i__) = rNULL__
if(ar__3(18:25) .EQ. '') e_Porb_1(i__) = rNULL__
if(ar__3(27:32) .EQ. '') Cycle_1(i__) = rNULL__
if(ar__3(34:39) .EQ. '') gamma_2(i__) = rNULL__
if(ar__3(41:44) .EQ. '') e_gamma_2(i__) = rNULL__
if(ar__3(58:62) .EQ. '') e_1(i__) = rNULL__
if(ar__3(64:68) .EQ. '') e_e_1(i__) = rNULL__
if(ar__3(70:74) .EQ. '') omega_1(i__) = rNULL__
if(ar__3(76:79) .EQ. '') e_omega_1(i__) = rNULL__
if(ar__3(81:89) .EQ. '') T0_1(i__) = rNULL__
if(ar__3(91:95) .EQ. '') e_T0_1(i__) = rNULL__
if(ar__3(122:126) .EQ. '') q(i__) = rNULL__
if(ar__3(128:132) .EQ. '') e_q(i__) = rNULL__
if(ar__3(134:137) .EQ. '') sigma_1(i__) = rNULL__
if(ar__3(139:141) .EQ. '') Num_1(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ ID_3(i__),Porb_1(i__),e_Porb_1(i__),Cycle_1(i__),gamma_2(i__),
+ e_gamma_2(i__),K_1(i__),e_K_1(i__),e_1(i__),e_e_1(i__),
+ omega_1(i__),e_omega_1(i__),T0_1(i__),e_T0_1(i__),
+ asini_1(i__),e_asini_1(i__),msin3i(i__),e_msin3i(i__),q(i__),
+ e_q(i__),sigma_1(i__),Num_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end