Conversion of standardized ReadMe file for
file /./ftp/cats/J/other/A_ARV/18.67 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 2024-Apr-17
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/other/A+ARV/18.67 Accurate masses and radii of normal stars (Torres+, 2010)
*================================================================================
*Accurate masses and radii of normal stars: modern results and applications.
* Torres G., Andersen J., Gimenez A.
* <Astron. Astrophys. Rev., 18, 67-126 (2010)>
* =2010A&ARv..18...67T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observed parameters and derived quantities for
95 binary systems, with mean errors
integer*4 nr__
parameter (nr__=190) ! Number of records
character*121 ar__ ! Full-size record
integer*4 Seq (nr__) ! [1/95] System sequential number
character*15 Name (nr__) ! Star name of the component
real*4 Per (nr__) ! (d) Period of the component
character*1 u_Per (nr__) ! Uncertainty flag on Per
character*1 x_Per (nr__) ! [y] y if period in years
character*1 Comp (nr__) ! [AB] Component
character*8 SpType (nr__) ! MK spectral type of the component
real*8 Mass (nr__) ! (solMass) Mass of the component
real*8 e_Mass (nr__) ! (solMass) rms uncertainty on Mass
real*8 Rad (nr__) ! (solRad) Radius of the component
real*4 e_Rad (nr__) ! (solRad) rms uncertainty on Rad
integer*4 Teff (nr__) ! (K) Effective temperature of the component
integer*4 e_Teff (nr__) ! (K) rms uncertainty on Teff
real*4 logg (nr__) ! ([cm/s2]) Surface gravity of the component
real*8 e_logg (nr__) ! ([cm/s2]) rms uncertainty on logg
real*4 logL (nr__) ! ([solLum]) Luminosity of the component
real*4 e_logL (nr__) ! ([solLum]) rms uncertainty on logL
real*4 VMAG (nr__) ! (mag) Absolute V magnitude of the component
real*4 e_VMAG (nr__) ! (mag) rms uncertainty on on VMAG
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Other parameters and references for the 95
binary systems
integer*4 nr__1
parameter (nr__1=95) ! Number of records
character*190 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! [1/95] System sequential number
character*14 Name_1 (nr__1) ! Name of the system
character*1 n_Name (nr__1) ! [*] Pre-main-sequence phase (1)
real*8 Dist (nr__1) ! (pc) Heliocentric distance
real*8 e_Dist (nr__1) ! (pc) rms uncertainty on Dist
real*4 E_B_V (nr__1) ! (mag) ? Interstellar redenning
real*4 e_E_B_V (nr__1) ! (mag) ? rms uncertainty on E(B-V)
real*4 vAsini (nr__1) ! (km/s) ? Rotation velocity of A component
real*4 e_vAsini (nr__1) ! (km/s) ? rms uncertainty on vAsini
real*4 vBsini (nr__1) ! (km/s) ? Rotation velocity of B component
real*4 e_vBsini (nr__1) ! (km/s) ? rms uncertainty on vBsini
real*4 v_Fe_H_ (nr__1) ! ([Sun]) ? Metallicity
real*4 e__Fe_H_ (nr__1) ! ([Sun]) ? rms uncertainty on [Fe/H]
character*1 n__Fe_H_ (nr__1) ! [*] Adopted value of [Fe/H] (2)
real*4 logAge (nr__1) ! ([yr]) Age of the system
character*1 u_logAge (nr__1) ! Uncertainty flag on Age
character*19 Ref1 (nr__1) ! Reference (BibCode)
character*19 Ref2 (nr__1) ! Reference (BibCode)
character*19 Ref3 (nr__1) ! Reference (BibCode)
character*19 Ref4 (nr__1) ! Reference (BibCode)
character*19 Ref5 (nr__1) ! Reference (BibCode)
*Note (1): An asterisk after the name indicates one or both components are
* in the pre-main sequence phase.
*Note (2): An asterisk after the metallicity indicates a value adopted from
* the binary's parent population (cluster or galaxy).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Systems with eccentric orbits, and systems
with apsidal motion
integer*4 nr__2
parameter (nr__2=44) ! Number of records
character*115 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! [1/95] System sequential number
character*14 Name_2 (nr__2) ! Name of the system
real*8 Ecc (nr__2) ! Eccentricity
real*8 e_Ecc (nr__2) ! rms uncertainty on Ecc
real*8 dw_dt (nr__2) ! ? Aspidal motion rate in deg/cycle unit
real*8 e_dw_dt (nr__2) ! ? rms uncertainty in dw/dt in deg/cycle unit
real*4 logK2 (nr__2) ! ? Internal structure constant (1)
real*4 e_logK2 (nr__2) ! ? rms uncertainty on logK2
integer*4 Flag (nr__2) ! [0/3]? Peculiarities (2)
character*19 Ref (nr__2) ! Reference (BibCode)
character*32 n_Ref (nr__2) ! Complete reference if having not a bibcode
*Note (1): indicative of the degree of central density concentration
* of the component stars.
*Note (2): Flag as follows:
* 0 = Apsidal motion not detected
* 1 = Rotational axes misaligned
* 2 = Apsidal motion dominated by General Relativity component
* 3 = Significance of apsidal motion too low for reliable
* determination of logk2
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Parameters for 23 astrometric binary systems
integer*4 nr__3
parameter (nr__3=23) ! Number of records
character*336 ar__3 ! Full-size record
integer*4 Seq_3 (nr__3) ! [1/23] Sequential number
character*14 Name1 (nr__3) ! Name of the primary
real*4 Vmag_1 (nr__3) ! (mag) V magnitude of the system
real*8 Plx (nr__3) ! (mas) Orbital parallax
real*4 e_Plx (nr__3) ! (mas) rms uncertainty on Plx
real*4 Dist_1 (nr__3) ! (pc) Heliocentric distance of the system
real*4 e_Dist_1 (nr__3) ! (pc) rms uncertainty on Dist
character*2 cp1 (nr__3) ! [A Aa] Component for primary
character*5 SpT1 (nr__3) ! MK spectral type of primary
real*4 Mass1 (nr__3) ! (solMass) Mass of the primary
real*4 e_Mass1 (nr__3) ! (solMass) rms uncertainty on Mass1
real*4 logL1 (nr__3) ! ([solLum]) Luminosity of the primary
character*1 u_logL1 (nr__3) ! Uncertainty flag on logL1
real*4 e_logL1 (nr__3) ! ([solLum]) ? rms uncertainty on logL1
integer*4 Teff1 (nr__3) ! (K) ? Effective temperature of the primary
integer*4 e_Teff1 (nr__3) ! (K) ? rms uncertainty on Teff
real*4 Rad1 (nr__3) ! (solRad) ? Radius of the primary
real*4 e_Rad1 (nr__3) ! (solRad) ? rms uncertainty on Rad1
character*19 Ref_1 (nr__3) ! Reference (Bibcode)
character*19 Ref2_1 (nr__3) ! Reference (Bibcode)
character*11 Name2 (nr__3) ! Name of the secondary
real*8 Per_1 (nr__3) ! (d) Period of the system
character*1 x_Per_1 (nr__3) ! [y] y if period in years
real*8 a (nr__3) ! (mas) ? Semi-major axis of the system
real*4 e_a (nr__3) ! (mas) ? rms uncertainty on a
real*4 v_Fe_H__1 (nr__3) ! ? Metallicity of the system
real*4 e__Fe_H__1 (nr__3) ! ? rms uncertainty on [Fe/H]
character*2 cp2 (nr__3) ! [Ab B Ba] Component for secondary
character*5 SpT2 (nr__3) ! MK spectral type of the secondary
real*4 Mass2 (nr__3) ! (solMass) Mass of the secondary
real*4 e_Mass2 (nr__3) ! (solMass) rms uncertainty on Mass2
real*4 logL2 (nr__3) ! ([solLum]) ? Luminosity of the secondary
character*1 u_logL2 (nr__3) ! Uncertainty flag on logL2
real*4 e_logL2 (nr__3) ! ([solLum]) ? rms uncertainty on the logL2
integer*4 Teff2 (nr__3) ! (K) ? Effective temperature of the secondary
integer*4 e_Teff2 (nr__3) ! (K) ? rms uncertainty on Teff2
real*4 Rad2 (nr__3) ! (solRad) ? Radius of the secondary
real*4 e_Rad2 (nr__3) ! (solRad) ? rms uncertainty of the secondary
character*11 Name3 (nr__3) ! Name of the the tertiary
real*4 Per3 (nr__3) ! (d) ? Period of the tertiary
character*1 x_Per3 (nr__3) ! [y] y if period in years
real*4 a3 (nr__3) ! (mas) ? Semi-major axis of the tertiary
real*4 e_a3 (nr__3) ! (mas) ? rms uncertainty on a3
character*2 cp3 (nr__3) ! [B Bb] Component for tertiary
character*5 SpT3 (nr__3) ! MK spectral type of the tertiary
real*4 Mass3 (nr__3) ! (solMass) ? Mass of the tertiary
real*4 e_Mass3 (nr__3) ! (solMass) ? rms uncertainty on Mass3
real*4 logL3 (nr__3) ! ([solLum]) ? Luminosity of the tertiary
character*1 u_logL3 (nr__3) ! Uncertainty flag on logL3
real*4 e_logL3 (nr__3) ! ([solLum]) ? rms uncertainty on logL3
integer*4 Teff3 (nr__3) ! (K) ? Effective temperature of the tertiary
integer*4 e_Teff3 (nr__3) ! (K) ? rms uncertainty on the tertiary
real*4 Rad3 (nr__3) ! (solRad) ? Radius of the tertiary
real*4 e_Rad3 (nr__3) ! (solRad) ? rms uncertainty on Rad3
C=============================================================================
C Loading file 'table1.dat' ! Observed parameters and derived quantities for
* 95 binary systems, with mean errors
C Format for file interpretation
1 format(
+ I2,1X,A15,1X,F6.2,A1,A1,1X,A1,1X,A8,2X,F8.5,1X,F7.5,2X,F7.4,
+ 1X,F6.4,1X,I5,1X,I4,1X,F6.4,1X,F7.5,F6.3,1X,F5.3,1X,F5.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,190
read(1,'(A121)')ar__
read(ar__,1)
+ Seq(i__),Name(i__),Per(i__),u_Per(i__),x_Per(i__),Comp(i__),
+ SpType(i__),Mass(i__),e_Mass(i__),Rad(i__),e_Rad(i__),
+ Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),logL(i__),
+ e_logL(i__),VMAG(i__),e_VMAG(i__)
c ..............Just test output...........
write(6,1)
+ Seq(i__),Name(i__),Per(i__),u_Per(i__),x_Per(i__),Comp(i__),
+ SpType(i__),Mass(i__),e_Mass(i__),Rad(i__),e_Rad(i__),
+ Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),logL(i__),
+ e_logL(i__),VMAG(i__),e_VMAG(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Other parameters and references for the 95
* binary systems
C Format for file interpretation
2 format(
+ I2,1X,A14,A1,F9.3,1X,F8.3,2X,F5.3,1X,F5.3,1X,F5.1,1X,F4.1,2X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F4.2,A1,1X,F4.1,A1,1X,A19,1X,A19,1X,
+ A19,1X,A19,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,95
read(1,'(A190)')ar__1
read(ar__1,2)
+ Seq_1(i__),Name_1(i__),n_Name(i__),Dist(i__),e_Dist(i__),
+ E_B_V(i__),e_E_B_V(i__),vAsini(i__),e_vAsini(i__),vBsini(i__),
+ e_vBsini(i__),v_Fe_H_(i__),e__Fe_H_(i__),n__Fe_H_(i__),
+ logAge(i__),u_logAge(i__),Ref1(i__),Ref2(i__),Ref3(i__),
+ Ref4(i__),Ref5(i__)
if(ar__1(39:43) .EQ. '') E_B_V(i__) = rNULL__
if(ar__1(45:49) .EQ. '') e_E_B_V(i__) = rNULL__
if(ar__1(51:55) .EQ. '') vAsini(i__) = rNULL__
if(ar__1(57:60) .EQ. '') e_vAsini(i__) = rNULL__
if(ar__1(63:67) .EQ. '') vBsini(i__) = rNULL__
if(ar__1(69:72) .EQ. '') e_vBsini(i__) = rNULL__
if(ar__1(74:78) .EQ. '') v_Fe_H_(i__) = rNULL__
if(ar__1(80:83) .EQ. '') e__Fe_H_(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),Name_1(i__),n_Name(i__),Dist(i__),e_Dist(i__),
+ E_B_V(i__),e_E_B_V(i__),vAsini(i__),e_vAsini(i__),vBsini(i__),
+ e_vBsini(i__),v_Fe_H_(i__),e__Fe_H_(i__),n__Fe_H_(i__),
+ logAge(i__),u_logAge(i__),Ref1(i__),Ref2(i__),Ref3(i__),
+ Ref4(i__),Ref5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Systems with eccentric orbits, and systems
* with apsidal motion
C Format for file interpretation
3 format(
+ I2,1X,A14,F7.5,1X,F7.5,2X,F7.5,1X,F7.5,1X,F5.2,1X,F4.2,1X,I1,
+ 1X,A19,1X,A32)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,44
read(1,'(A115)')ar__2
read(ar__2,3)
+ Seq_2(i__),Name_2(i__),Ecc(i__),e_Ecc(i__),dw_dt(i__),
+ e_dw_dt(i__),logK2(i__),e_logK2(i__),Flag(i__),Ref(i__),
+ n_Ref(i__)
if(ar__2(35:41) .EQ. '') dw_dt(i__) = rNULL__
if(ar__2(43:49) .EQ. '') e_dw_dt(i__) = rNULL__
if(ar__2(51:55) .EQ. '') logK2(i__) = rNULL__
if(ar__2(57:60) .EQ. '') e_logK2(i__) = rNULL__
if(ar__2(62:62) .EQ. '') Flag(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),Name_2(i__),Ecc(i__),e_Ecc(i__),dw_dt(i__),
+ e_dw_dt(i__),logK2(i__),e_logK2(i__),Flag(i__),Ref(i__),
+ n_Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Parameters for 23 astrometric binary systems
C Format for file interpretation
4 format(
+ I2,1X,A14,2X,F5.2,2X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A2,1X,A5,
+ 1X,F6.4,1X,F6.4,1X,F6.3,A1,F5.3,2X,I4,1X,I3,1X,F5.2,1X,F4.2,
+ 1X,A19,1X,A19,1X,A11,F8.3,A1,1X,F7.3,1X,F5.3,1X,F5.2,1X,F4.2,
+ 1X,A2,1X,A5,1X,F6.4,1X,F6.4,1X,F6.3,A1,F5.3,1X,I4,1X,I3,1X,
+ F4.2,1X,F4.2,1X,A11,1X,F6.2,A1,1X,F5.1,1X,F3.1,1X,A2,1X,A5,1X,
+ F6.4,1X,F6.4,2X,F6.3,A1,F5.3,1X,I4,1X,I3,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,23
read(1,'(A336)')ar__3
read(ar__3,4)
+ Seq_3(i__),Name1(i__),Vmag_1(i__),Plx(i__),e_Plx(i__),
+ Dist_1(i__),e_Dist_1(i__),cp1(i__),SpT1(i__),Mass1(i__),
+ e_Mass1(i__),logL1(i__),u_logL1(i__),e_logL1(i__),Teff1(i__),
+ e_Teff1(i__),Rad1(i__),e_Rad1(i__),Ref_1(i__),Ref2_1(i__),
+ Name2(i__),Per_1(i__),x_Per_1(i__),a(i__),e_a(i__),
+ v_Fe_H__1(i__),e__Fe_H__1(i__),cp2(i__),SpT2(i__),Mass2(i__),
+ e_Mass2(i__),logL2(i__),u_logL2(i__),e_logL2(i__),Teff2(i__),
+ e_Teff2(i__),Rad2(i__),e_Rad2(i__),Name3(i__),Per3(i__),
+ x_Per3(i__),a3(i__),e_a3(i__),cp3(i__),SpT3(i__),Mass3(i__),
+ e_Mass3(i__),logL3(i__),u_logL3(i__),e_logL3(i__),Teff3(i__),
+ e_Teff3(i__),Rad3(i__),e_Rad3(i__)
if(ar__3(84:88) .EQ. '') e_logL1(i__) = rNULL__
if(ar__3(91:94) .EQ. '') Teff1(i__) = iNULL__
if(ar__3(96:98) .EQ. '') e_Teff1(i__) = iNULL__
if(ar__3(100:104) .EQ. '') Rad1(i__) = rNULL__
if(ar__3(106:109) .EQ. '') e_Rad1(i__) = rNULL__
if(ar__3(172:178) .EQ. '') a(i__) = rNULL__
if(ar__3(180:184) .EQ. '') e_a(i__) = rNULL__
if(ar__3(186:190) .EQ. '') v_Fe_H__1(i__) = rNULL__
if(ar__3(192:195) .EQ. '') e__Fe_H__1(i__) = rNULL__
if(ar__3(220:225) .EQ. '') logL2(i__) = rNULL__
if(ar__3(227:231) .EQ. '') e_logL2(i__) = rNULL__
if(ar__3(233:236) .EQ. '') Teff2(i__) = iNULL__
if(ar__3(238:240) .EQ. '') e_Teff2(i__) = iNULL__
if(ar__3(242:245) .EQ. '') Rad2(i__) = rNULL__
if(ar__3(247:250) .EQ. '') e_Rad2(i__) = rNULL__
if(ar__3(264:269) .EQ. '') Per3(i__) = rNULL__
if(ar__3(272:276) .EQ. '') a3(i__) = rNULL__
if(ar__3(278:280) .EQ. '') e_a3(i__) = rNULL__
if(ar__3(291:296) .EQ. '') Mass3(i__) = rNULL__
if(ar__3(298:303) .EQ. '') e_Mass3(i__) = rNULL__
if(ar__3(306:311) .EQ. '') logL3(i__) = rNULL__
if(ar__3(313:317) .EQ. '') e_logL3(i__) = rNULL__
if(ar__3(319:322) .EQ. '') Teff3(i__) = iNULL__
if(ar__3(324:326) .EQ. '') e_Teff3(i__) = iNULL__
if(ar__3(328:331) .EQ. '') Rad3(i__) = rNULL__
if(ar__3(333:336) .EQ. '') e_Rad3(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Seq_3(i__),Name1(i__),Vmag_1(i__),Plx(i__),e_Plx(i__),
+ Dist_1(i__),e_Dist_1(i__),cp1(i__),SpT1(i__),Mass1(i__),
+ e_Mass1(i__),logL1(i__),u_logL1(i__),e_logL1(i__),Teff1(i__),
+ e_Teff1(i__),Rad1(i__),e_Rad1(i__),Ref_1(i__),Ref2_1(i__),
+ Name2(i__),Per_1(i__),x_Per_1(i__),a(i__),e_a(i__),
+ v_Fe_H__1(i__),e__Fe_H__1(i__),cp2(i__),SpT2(i__),Mass2(i__),
+ e_Mass2(i__),logL2(i__),u_logL2(i__),e_logL2(i__),Teff2(i__),
+ e_Teff2(i__),Rad2(i__),e_Rad2(i__),Name3(i__),Per3(i__),
+ x_Per3(i__),a3(i__),e_a3(i__),cp3(i__),SpT3(i__),Mass3(i__),
+ e_Mass3(i__),logL3(i__),u_logL3(i__),e_logL3(i__),Teff3(i__),
+ e_Teff3(i__),Rad3(i__),e_Rad3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end