Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/190/1 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-19
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