Conversion of standardized ReadMe file for
catalog 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 2024-Apr-24
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/ApJS/190/1 A survey of stellar families (Raghavan+, 2010)
*================================================================================
*A survey of stellar families: multiplicity of solar-type stars.
* Raghavan D., McAlister H.A., Henry T.J., Latham D.W., Marcy G.W.,
* Mason B.D., Gies D.R., White R.J., Ten Brummelaar T.A.
* <Astrophys. J. Suppl. Ser., 190, 1-42 (2010)>
* =2010ApJS..190....1R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Volume-limited sample of 454 solar-type stars
integer*4 nr__
parameter (nr__=454) ! Number of records
character*89 ar__ ! 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 RAh ! (h) ? Hour of Right Ascension (J2000)
integer*4 RAm ! (min) ? Minute of Right Ascension (J2000)
real*4 RAs ! (s) ? Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination (J2000)
integer*4 DEd ! (deg) ? Degree of Declination (J2000)
integer*4 DEm ! (arcmin) ? Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) ? Arcsecond of Declination (J2000)
integer*4 HIP ! Hipparcos name (Cat. I/239) (1)
integer*4 HD ! ? HD name (Cat. III/135)
real*4 Vmag ! (mag) ? Hipparcos V band magnitude
real*4 B_V ! (mag) Hipparcos (B-V) color
real*4 plx ! (mas) ? Hipparcos parallax
integer*4 r_plx ! [12]? Reference for plx
real*4 e_plx ! (mas) ? The 1{sigma} uncertainty in plx
real*4 vLplx ! (mas) ? van Leeuwen 2007 (Cat. I/311) parallax
real*4 e_vLplx ! (mas) ? The 1{sigma} uncertainty in vLplx
character*8 SpT1 ! MK Spectral type
integer*4 r_SpT1 ! ? Spectral type reference. See refs.dat file
*Note (1): 0 for Sun with no coordinates.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! CPM (Common Proper Motion) companions identified
integer*4 nr__1
parameter (nr__1=88) ! Number of records
character*90 ar__1 ! Full-size record
character*10 Type ! CPM companion type: Associated or Refuted
character*10 WDS ! The WDS identifier (Cat. B/wds)
character*7 Disc ! WDS discoverer designation
character*4 m_Disc ! WDS pair identifier
character*10 PName ! Primary name (HD or HIP)
character*1 f_PName ! [a] New companion discovered by this effort.
integer*4 rho ! (arcsec) CPM candidate companion separation
integer*4 theta ! (deg) CPM candidate companion position angle
real*8 Epoch ! CPM candidate companion epoch of observation
character*15 Name ! CPM candidate companion name
integer*4 Code ! Reason code for association/refutation (3)
integer*4 HIP_1 ! Hipparcos name of the system (added by CDS)
*Note (3): Reason code as follows:
* 1 = Photometric distance to the CPM candidate matches the Hipparcos distance
* to the primary (see Table 5).
* 2 = Published parallax and proper motion for the CPM candidate matches the
* corresponding primary's values from Hipparcos
* 3 = Spectroscopic distance to the CPM candidate matches the Hipparcos
* distance to the primary.
* 4 = Known companion with a published orbit (each of these were seen as
* comoving diffraction spikes).
* 5 = Published evidence, see Section 4.3.
* 6 = Companionship implied by proximity to the primary and a matching,
* large proper motion.
* 7 = Measurements of the pair in the WDS confirm orbital motion.
* 8 = While the proper motion is similar enough to enable selection as a
* possible CPM companion, the numerical value for the candidate in
* catalogs (Hog et al., 1998, Cat. I/259; Salim & Gould, 2003,
* Cat. J/ApJ/582/1011; Zacharias et al., 2004, Cat. I/289;
* Lepine & Shara, 2005, Cat. I/298; Hipparcos, Cat. I/239) are
* significantly different from the primary's Hipparcos values, ruling out
* a physical association.
* 9 = The candidate companion is a non-stellar artifact such as a plate
* defect.
* 10 = Photometric distance to the CPM candidate is significantly different
* than the primary's Hipparcos distance (see Table 5).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Spectral type, proper motion, and photometry of
CPM (Common Proper Motion) candidates
integer*4 nr__2
parameter (nr__2=38) ! Number of records
character*109 ar__2 ! Full-size record
character*3 Type_1 ! CPM companion type (1)
integer*4 HD_1 ! Primary name
integer*4 rho_1 ! (arcsec) ? Separation
character*8 SpType ! CPM candidate MK spectral type
integer*4 r_SpType ! ? SpType reference. See refs.dat file.
integer*4 pmRA ! (mas/yr) ? CPM candidate proper motion along RA
integer*4 pmDE ! (mas/yr) ? CPM candidate proper motion along DE
integer*4 r_pmRA ! ? Proper motion reference. See refs.dat file.
real*4 Vmag_1 ! (mag) CCD V-band magnitude of the CPM candidate
integer*4 r_Vmag ! Vmag reference. See refs.dat file.
real*4 Rmag ! (mag) ? CCD R-band magnitude of the CPM candidate
integer*4 r_Rmag ! ? Rmag reference. See refs.dat file.
real*4 Imag ! (mag) ? CCD I-band magnitude of the CPM candidate
integer*4 r_Imag ! ? Imag reference. See refs.dat file.
integer*4 No ! ? Observation number
real*4 Jmag ! (mag) CPM candidate 2MASS J-band magnitude
real*4 Hmag ! (mag) CPM candidate 2MASS H-band magnitude
real*4 Ksmag ! (mag) CPM candidate 2MASS Ks-band magnitude
real*4 Dist ! (pc) CPM candidate distance
character*1 f_Dist ! [abc] Flag on Dist (2)
real*4 e_Dist ! (pc) Dist error
integer*4 HIP_2 ! Hipparcos name of the system (added by CDS)
*Note (1): Type as follows:
* PAC = Physically associated CPM companions
* RCC = Refuted CPM candidates
*Note (2): Flag as follows:
* a = While this distance is too low compared to the primary's Hipparcos
* distance of 24.6pc, the companion is a roughly equal-brightness binary.
* Adjusting the Hipparcos and 2MASS magnitudes accordingly changes the
* distance estimate to 24.0+/-3.7, a much better match with the distance
* to the primary.
* b = See Section 4.3 for a discussion of these photometric distance estimates
* and the status of these companions.
* c = Even though the companion is HIP 110719, its Hipparcos astrometry has
* large errors, necessitating the photometric distance check.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Optical WDS entries
integer*4 nr__3
parameter (nr__3=298) ! Number of records
character*62 ar__3 ! Full-size record
character*10 WDS_1 ! The WDS identification (Cat. B/wds)
character*7 Disc_1 ! WDS discoverer designation
character*4 m_WDS ! WDS pair identifier
character*1 f_WDS ! [a] Flag for HD 200560 (1)
character*10 PName_1 ! Primary name (HD or HIP)
integer*4 Nobs ! Number of observations
integer*4 theta_1 ! (deg) Companion position angle
real*4 rho_2 ! (arcsec) Companion separation
integer*4 Epoch_1 ! Companion epoch of observation
integer*4 HIP_3 ! Hipparcos name of the system (added by CDS)
*Note (1): a = While the AB pair has an orbital solution in ORB6 and is likely
* physically bound, it is listed here to identify that component B is only
* optically associated with the sample star, which is component C. As noted in
* the next line, the blinking of archival images also helps identify AC as an
* optical pair. Thus, the physical pair AB is not associated with sample star C.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Accelerating proper motion solutions
integer*4 nr__4
parameter (nr__4=91) ! Number of records
character*89 ar__4 ! Full-size record
integer*4 HD_2 ! ? HD name
character*1 f_HD ! [cde] Flag on HD (1)
integer*4 HIP_4 ! Hipparcos name of primary
character*1 H59 ! [CGOX] Primary identification of companions in
* the Hipparcos catalog (2)
real*8 pmRA_1 ! (mas/yr) Hipparcos proper motion along RA (3)
real*4 e_pmRA ! (mas/yr) Uncertainty in pmRA
real*8 pmDE_1 ! (mas/yr) Hipparcos proper motion along DE
real*4 e_pmDE ! (mas/yr) Uncertainty in pmDE
real*4 pmRA2 ! (mas/yr) ? Tycho-2 proper motion along RA (3)
real*4 e_pmRA2 ! (mas/yr) ? Uncertainty in pmRA2
real*8 pmDE2 ! (mas/yr) ? Tycho-2 proper motion along DE
real*4 e_pmDE2 ! (mas/yr) ? Uncertainty in pmDE2
real*4 sigma ! ? Significance of difference between Hipparcos
* Tycho-2 proper motions (4)
character*1 MK05 ! [Y] if proper-motion binary in Makarov & Kaplan
* 2005, Cat. J/AJ/129/2420
character*1 F07 ! [Y] if proper-motion binary in Frankowski et
* al., 2007, Cat. J/A+A/464/377)
character*3 St ! Companion status (5)
integer*4 Note ! Companion reason note (6)
integer*4 HIPc ! ? HIP number of companion (added at CDS)
*Note (1): Flag as follows:
* c = Wide companion, 319" away from HD 7693.
* d = Wide companion, 185" away from HD 53705.
* e = Wide companion, 756" away from HIP 36357.
*Note (2): Multiplicity flag as follows:
* C = Double stars resolved by Hipparcos as separated components that could be
* modeled as single stars, usually with an assumed common parallax.
* G = Accelerating proper motion or higher-order terms.
* O = Orbital solutions
* X = Stochastic solutions (probably astrometric binaries with short period)
* Further details are given in the Double and Multiple Systems Annex of
* Hipparcos. See also section 3.2 for further details.
*Note (3): {mu}_{alpha}_*cos[{delta}].
*Note (4): Computed as the root sum squared of the difference in each axis
* divided by the corresponding larger error.
*Note (5): Status as follows:
* YES = physically associated companion;
* MAY = unconfirmed candidate retained for further investigations.
*Note (6): Reason note as follows:
* 1 = Visual and/or spectroscopic binary with an orbital solution.
* 2 = Nearby companion is likely responsible for the proper motion
* acceleration (see Section 4.3).
* 3 = Hipparcos G flag and the {chi}^2^ test in Frankowski et al. (2007,
* Cat. J/A+A/464/377) suggest an unseen companion, but because the
* Hipparcos and Tycho-2 proper motions differ by less than 3{sigma}, this
* is retained as a candidate for further investigations.
* 4 = See individual system notes in Section 4.3.
* 5 = Hipparcos "G" flag and a greater than 3{sigma} difference in proper
* motion indicate an unseen companion.
* 6 = Nearby companion with evidence of orbital motion based on WDS
* measurements is likely responsible for the proper motion acceleration.
* 7 = Radial velocity variations indicate a spectroscopic binary, but not
* enough observations exist to derive an orbit.
* 8 = Greater than 3{sigma} difference in proper motion is the only evidence
* of a companion. The companion is considered physical if it also passed
* the {chi}^2^ test in Frankowski et al. (2007, Cat. J/A+A/464/377),
* otherwise is retained as a candidate.
* 9 = Nearby companion with matching photometric distance is likely
* responsible for the proper motion acceleration (see Table 10).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Visual orbit solutions
integer*4 nr__5
parameter (nr__5=98) ! Number of records
character*96 ar__5 ! Full-size record
character*19 Type_2 ! Visual orbit type
character*10 WDS_2 ! The WDS identification
integer*4 HD_3 ! HD name
integer*4 HIP_5 ! Hipparcos name
character*14 Disc_2 ! WDS discoverer designation
integer*4 Nobs_1 ! ? Number of observations
integer*4 DelT ! (yr) ? Observational time span
real*8 Per ! ? Orbital period
character*1 x_Per ! Units for Per (y = years, d = days)
integer*4 G ! Orbit grade from ORB6, (1=definitive to
* 5=indeterminate) (1)
character*8 r_G ! Orbit reference from ORB6, see refs.dat file (1)
integer*4 SB ! [1,2]? single-lined (1) or double-lined (2) SB
character*3 Status ! Companion status (3)
*Note (1): ORB6 = Sixth Catalog of Orbits of Visual Binary Stars, available on:
* http://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/wds/orb6
*Note (3): Status as follows:
* YES = physically associated companion;
* NO = unrelated field star;
* MAY = unconfirmed candidate retained for further investigations.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Survey stars and their stellar, brown dwarf, and
planetary companions
integer*4 nr__6
parameter (nr__6=799) ! Number of records
character*105 ar__6 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh_1 ! (h) ? Sample star Hour of Right Ascension (J2000) (1)
integer*4 RAm_1 ! (min) ? Sample star Minute of Right Ascension (J2000)
real*4 RAs_1 ! (s) ? Sample star Second of Right Ascension (J2000)
character*1 DE__1 ! Sign of sample star Declination (J2000) (1)
integer*4 DEd_1 ! (deg) ? Sample star Degree of Declination (J2000)
integer*4 DEm_1 ! (arcmin) ? Sample star Arcminute of Declination (J2000)
real*4 DEs_1 ! (arcsec) ? Sample star Arcsecond of Declination (J2000)
integer*4 HD_4 ! ? HD name
character*15 OName ! Other name (1)
character*1 N ! [N] System has specific note in Section 4.3.
* See notes.dat file
character*7 Comp ! Companion identification
character*1 f_Comp ! [abc] Flag on Comp (2)
real*8 Per_1 ! ? Orbital period of the pair
character*1 x_Per_1 ! Units for Per (y = years, d = days, h = hours)
real*8 Asep ! (arcsec) ? Angular separation between the components
real*8 Lsep ! (AU) ? Projected linear separation (3)
character*1 St_1 ! [YM] Companionship status:
* confirmed (Y) or candidate (M)
character*1 VB ! [OPU] Visual binary orbit quality code (5)
character*1 SB_1 ! [12V] Spectroscopic binary code (6)
character*1 CP ! [M-T] Common proper motion code (7)
character*1 OT ! [ELM] Other unresolved companion code (8)
character*1 CH ! [SV] Companions seen by CHARA LBOI (9)
integer*4 HIP_6 ! Hipparcos name of the system (added by CDS) (1)
*Note (1): Sun, Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, and Neptune
* have "0" for HIP number and no coordinates.
*Note (2): Flag as follows:
* a = The sample star is not the system's primary, which is identified as
* component A below.
* b = The brightest component of the system is HD 25998, but is designated as
* component E in the WDS. Component A is the wide CPM companion,
* HD 25893, which is about 2 magnitudes fainter and itself a visual
* binary. We have retained the component designations of the WDS, so the
* fainter visual pair is AB and the wide CPM companion is E. WDS
* components C and D are optical, and E itself might have a close
* companion, as evidenced by its accelerating proper motion (see Table 8.)
* c = WDS lists these entries for HD 200595, a bright binary 153" away from
* the sample star HD 200560, but one that is not physically associated
* with it. HD 200560 is itself is a close CPM pair and listed in the WDS
* as CD. We have retained the WDS designations, which makes C and D the
* only physically associated components of this system.
*Note (3): If Asep is empty, Lsep contains the semi-major axis from orbital
* solutions; it corresponds to Asep otherwise.
*Note (5): Code as follows:
* O = robust orbit of grade 1, 2, 3, or 8;
* P = preliminary orbit of grade 4 or 5;
* U = unresolved photocentric-motion orbit.
*Note (6): Code as follows:
* 1 = single-lined;
* 2 = double-lined;
* V = radial-velocity variations indicating a companion, but no
* orbital solution.
*Note (7): Code as follows:
* M = close pair with matching proper motions;
* O = pair with evidence of orbital motion;
* P = companion with matching proper motion and photometric distance;
* R = close pair with published evidence of companionship;
* S = companion with matching proper motion and spectral type
* identification that is consistent with the primary's distance;
* T = pair with independently-measured matching proper motions and
* trigonometric parallaxes.
*Note (8): Code as follows:
* E = eclipsing binary;
* L = companion indicated by an overluminous star;
* M = companion implied by proper motion acceleration.
*Note (9): Code as follows:
* S = SFP (separated fringe packet) (Raghavan et al. 2010a, in prep.);
* V = visibility-modulation binaries (Raghavan et al. 2009ApJ...690..394R,
* 2010b, in prep.).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table17.dat' ! Physical parameters of the sample stars
integer*4 nr__7
parameter (nr__7=454) ! Number of records
character*55 ar__7 ! Full-size record
integer*4 HD_5 ! ? HD name
integer*4 HIP_7 ! Hipparcos name (0 for Sun)
character*2 Comp_1 ! Component identification
character*8 SpT2 ! MK spectral type
character*2 r_SpT2 ! Reference for SpType. See refs.dat file
real*4 Mass ! (solMass) ? Stellar mass
character*2 r_Mass ! Reference for Mass. See refs.dat file
real*4 v_Fe_H_ ! ([Sun]) ? Metallicity
character*2 r__Fe_H_ ! Reference for [Fe/H]. See refs.dat file
real*4 logRHK ! ([---]) ? Chromospheric emission measure log(R'_HK_)
character*2 r_logRHK ! Reference for logRHK. See refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table18.dat' ! Spectral types and masses of the companions
integer*4 nr__8
parameter (nr__8=213) ! Number of records
character*64 ar__8 ! Full-size record
character*13 CompName ! Companion name (HD or HIP NNNNNN AA)
character*15 OName_1 ! Alternative name
character*14 SpType_1 ! MK spectral type
character*2 r_SpType_1 ! Reference for SpType. See refs.dat file
real*4 Mass_1 ! (solMass) ? Companion mass
character*1 f_Mass ! [a] Flag on Mass (1)
character*2 r_Mass_1 ! Reference for Mass. See refs.dat file
integer*4 HIP_8 ! Hipparcos name of the system (added by CDS)
*Note (1):
* a = Estimated using mass-sum of Aa+Ab from DM91 (Duquennoy & Mayor,
* 1991A&A...248..485D) and mass of Aa from Nordstrom et al. (2004,
* Cat V/117).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References for table 1, 5, 11, 17 and 18
integer*4 nr__9
parameter (nr__9=135) ! Number of records
character*199 ar__9 ! Full-size record
character*8 Ref ! Reference identification
character*19 BibCode ! Bibcode
character*21 Aut ! Author's name
character*148 Com ! Comment
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Notes on individual systems (from section 4.3)
integer*4 nr__10
parameter (nr__10=905) ! Number of records
character*84 ar__10 ! Full-size record
integer*4 HIP_9 ! Hipparcos name of the system (added at CDS)
integer*4 HD_6 ! HD number of the star
character*70 Text ! Text of note
C=============================================================================
C Loading file 'table1.dat' ! Volume-limited sample of 454 solar-type stars
C Format for file interpretation
1 format(
+ I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,I6,1X,I6,1X,F5.2,
+ 1X,F5.3,1X,F6.2,1X,I2,1X,F4.2,1X,F6.2,1X,F4.2,1X,A8,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,454
read(1,'(A89)')ar__
read(ar__,1)
+ RAh,RAm,RAs,DE_,DEd,DEm,DEs,HIP,HD,Vmag,B_V,plx,r_plx,e_plx,
+ vLplx,e_vLplx,SpT1,r_SpT1
if(ar__(1:2) .EQ. '') RAh = iNULL__
if(ar__(4:5) .EQ. '') RAm = iNULL__
if(ar__(7:11) .EQ. '') RAs = rNULL__
if(ar__(14:15) .EQ. '') DEd = iNULL__
if(ar__(17:18) .EQ. '') DEm = iNULL__
if(ar__(20:23) .EQ. '') DEs = rNULL__
if(ar__(32:37) .EQ. '') HD = iNULL__
if(ar__(39:43) .EQ. '') Vmag = rNULL__
if(ar__(51:56) .EQ. '') plx = rNULL__
if(ar__(58:59) .EQ. '') r_plx = iNULL__
if(ar__(61:64) .EQ. '') e_plx = rNULL__
if(ar__(66:71) .EQ. '') vLplx = rNULL__
if(ar__(73:76) .EQ. '') e_vLplx = rNULL__
if(ar__(88:89) .EQ. '') r_SpT1 = iNULL__
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,1)
+ RAh,RAm,RAs,DE_,DEd,DEm,DEs,HIP,HD,Vmag,B_V,plx,r_plx,e_plx,
+ vLplx,e_vLplx,SpT1,r_SpT1
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! CPM (Common Proper Motion) companions identified
C Format for file interpretation
2 format(
+ A10,1X,A10,1X,A7,1X,A4,1X,A10,1X,A1,1X,I4,1X,I3,1X,F7.2,1X,
+ A15,1X,I2,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,88
read(1,'(A90)')ar__1
read(ar__1,2)
+ Type,WDS,Disc,m_Disc,PName,f_PName,rho,theta,Epoch,Name,Code,
+ HIP_1
c ..............Just test output...........
write(6,2)
+ Type,WDS,Disc,m_Disc,PName,f_PName,rho,theta,Epoch,Name,Code,
+ HIP_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Spectral type, proper motion, and photometry of
* CPM (Common Proper Motion) candidates
C Format for file interpretation
3 format(
+ A3,4X,I6,1X,I3,1X,A8,1X,I2,1X,I4,1X,I4,1X,I2,1X,F5.2,1X,I2,1X,
+ F5.2,1X,I2,1X,F5.2,1X,I2,1X,I1,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.1,1X,A1,1X,F4.1,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,38
read(1,'(A109)')ar__2
read(ar__2,3)
+ Type_1,HD_1,rho_1,SpType,r_SpType,pmRA,pmDE,r_pmRA,Vmag_1,
+ r_Vmag,Rmag,r_Rmag,Imag,r_Imag,No,Jmag,Hmag,Ksmag,Dist,f_Dist,
+ e_Dist,HIP_2
if(ar__2(15:17) .EQ. '') rho_1 = iNULL__
if(ar__2(28:29) .EQ. '') r_SpType = iNULL__
if(ar__2(31:34) .EQ. '') pmRA = iNULL__
if(ar__2(36:39) .EQ. '') pmDE = iNULL__
if(ar__2(41:42) .EQ. '') r_pmRA = iNULL__
if(ar__2(53:57) .EQ. '') Rmag = rNULL__
if(ar__2(59:60) .EQ. '') r_Rmag = iNULL__
if(ar__2(62:66) .EQ. '') Imag = rNULL__
if(ar__2(68:69) .EQ. '') r_Imag = iNULL__
if(ar__2(71:71) .EQ. '') No = iNULL__
c ..............Just test output...........
write(6,3)
+ Type_1,HD_1,rho_1,SpType,r_SpType,pmRA,pmDE,r_pmRA,Vmag_1,
+ r_Vmag,Rmag,r_Rmag,Imag,r_Imag,No,Jmag,Hmag,Ksmag,Dist,f_Dist,
+ e_Dist,HIP_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Optical WDS entries
C Format for file interpretation
4 format(
+ A10,1X,A7,1X,A4,1X,A1,1X,A10,1X,I3,1X,I3,1X,F5.1,1X,I4,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,298
read(1,'(A62)')ar__3
read(ar__3,4)
+ WDS_1,Disc_1,m_WDS,f_WDS,PName_1,Nobs,theta_1,rho_2,Epoch_1,
+ HIP_3
c ..............Just test output...........
write(6,4)
+ WDS_1,Disc_1,m_WDS,f_WDS,PName_1,Nobs,theta_1,rho_2,Epoch_1,
+ HIP_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Accelerating proper motion solutions
C Format for file interpretation
5 format(
+ I6,1X,A1,1X,I6,1X,A1,1X,F7.2,1X,F4.2,1X,F8.2,1X,F4.2,1X,F6.1,
+ 1X,F3.1,1X,F7.1,1X,F3.1,1X,F4.1,1X,A1,1X,A1,1X,A3,1X,I1,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,91
read(1,'(A89)')ar__4
read(ar__4,5)
+ HD_2,f_HD,HIP_4,H59,pmRA_1,e_pmRA,pmDE_1,e_pmDE,pmRA2,e_pmRA2,
+ pmDE2,e_pmDE2,sigma,MK05,F07,St,Note,HIPc
if(ar__4(1:6) .EQ. '') HD_2 = iNULL__
if(ar__4(46:51) .EQ. '') pmRA2 = rNULL__
if(ar__4(53:55) .EQ. '') e_pmRA2 = rNULL__
if(ar__4(57:63) .EQ. '') pmDE2 = rNULL__
if(ar__4(65:67) .EQ. '') e_pmDE2 = rNULL__
if(ar__4(69:72) .EQ. '') sigma = rNULL__
if(ar__4(84:89) .EQ. '') HIPc = iNULL__
c ..............Just test output...........
write(6,5)
+ HD_2,f_HD,HIP_4,H59,pmRA_1,e_pmRA,pmDE_1,e_pmDE,pmRA2,e_pmRA2,
+ pmDE2,e_pmDE2,sigma,MK05,F07,St,Note,HIPc
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Visual orbit solutions
C Format for file interpretation
6 format(
+ A19,1X,A10,1X,I6,1X,I6,1X,A14,1X,I4,1X,I3,1X,F8.2,1X,A1,1X,I1,
+ 1X,A8,1X,I1,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,98
read(1,'(A96)')ar__5
read(ar__5,6)
+ Type_2,WDS_2,HD_3,HIP_5,Disc_2,Nobs_1,DelT,Per,x_Per,G,r_G,SB,
+ Status
if(ar__5(61:64) .EQ. '') Nobs_1 = iNULL__
if(ar__5(66:68) .EQ. '') DelT = iNULL__
if(ar__5(70:77) .EQ. '') Per = rNULL__
if(ar__5(92:92) .EQ. '') SB = iNULL__
c ..............Just test output...........
write(6,6)
+ Type_2,WDS_2,HD_3,HIP_5,Disc_2,Nobs_1,DelT,Per,x_Per,G,r_G,SB,
+ Status
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Survey stars and their stellar, brown dwarf, and
* planetary companions
C Format for file interpretation
7 format(
+ I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,I6,1X,A15,1X,A1,1X,
+ A7,1X,A1,1X,F8.2,1X,A1,1X,F7.2,1X,F8.2,1X,A1,1X,A1,1X,A1,1X,
+ A1,1X,A1,1X,A1,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,799
read(1,'(A105)')ar__6
read(ar__6,7)
+ RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,HD_4,OName,N,Comp,
+ f_Comp,Per_1,x_Per_1,Asep,Lsep,St_1,VB,SB_1,CP,OT,CH,HIP_6
if(ar__6(1:2) .EQ. '') RAh_1 = iNULL__
if(ar__6(4:5) .EQ. '') RAm_1 = iNULL__
if(ar__6(7:11) .EQ. '') RAs_1 = rNULL__
if(ar__6(14:15) .EQ. '') DEd_1 = iNULL__
if(ar__6(17:18) .EQ. '') DEm_1 = iNULL__
if(ar__6(20:23) .EQ. '') DEs_1 = rNULL__
if(ar__6(25:30) .EQ. '') HD_4 = iNULL__
if(ar__6(60:67) .EQ. '') Per_1 = rNULL__
if(ar__6(71:77) .EQ. '') Asep = rNULL__
if(ar__6(79:86) .EQ. '') Lsep = rNULL__
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,7)
+ RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,HD_4,OName,N,Comp,
+ f_Comp,Per_1,x_Per_1,Asep,Lsep,St_1,VB,SB_1,CP,OT,CH,HIP_6
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table17.dat' ! Physical parameters of the sample stars
C Format for file interpretation
8 format(
+ I6,1X,I6,1X,A2,1X,A8,1X,A2,1X,F4.2,1X,A2,1X,F5.2,1X,A2,1X,
+ F6.3,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table17.dat')
write(6,*) '....Loading file: table17.dat'
do i__=1,454
read(1,'(A55)')ar__7
read(ar__7,8)
+ HD_5,HIP_7,Comp_1,SpT2,r_SpT2,Mass,r_Mass,v_Fe_H_,r__Fe_H_,
+ logRHK,r_logRHK
if(ar__7(1:6) .EQ. '') HD_5 = iNULL__
if(ar__7(30:33) .EQ. '') Mass = rNULL__
if(ar__7(38:42) .EQ. '') v_Fe_H_ = rNULL__
if(ar__7(47:52) .EQ. '') logRHK = rNULL__
c ..............Just test output...........
write(6,8)
+ HD_5,HIP_7,Comp_1,SpT2,r_SpT2,Mass,r_Mass,v_Fe_H_,r__Fe_H_,
+ logRHK,r_logRHK
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table18.dat' ! Spectral types and masses of the companions
C Format for file interpretation
9 format(A13,1X,A15,1X,A14,1X,A2,1X,F4.2,1X,A1,1X,A2,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table18.dat')
write(6,*) '....Loading file: table18.dat'
do i__=1,213
read(1,'(A64)')ar__8
read(ar__8,9)
+ CompName,OName_1,SpType_1,r_SpType_1,Mass_1,f_Mass,r_Mass_1,
+ HIP_8
if(ar__8(49:52) .EQ. '') Mass_1 = rNULL__
c ..............Just test output...........
write(6,9)
+ CompName,OName_1,SpType_1,r_SpType_1,Mass_1,f_Mass,r_Mass_1,
+ HIP_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References for table 1, 5, 11, 17 and 18
C Format for file interpretation
10 format(A8,1X,A19,1X,A21,1X,A148)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,135
read(1,'(A199)')ar__9
read(ar__9,10)Ref,BibCode,Aut,Com
c ..............Just test output...........
write(6,10)Ref,BibCode,Aut,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Notes on individual systems (from section 4.3)
C Format for file interpretation
11 format(I6,1X,I6,1X,A70)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,905
read(1,'(A84)')ar__10
read(ar__10,11)HIP_9,HD_6,Text
c ..............Just test output...........
write(6,11)HIP_9,HD_6,Text
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end