Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/687/1264 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-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/ApJ/687/1264 Age estimation for solar-type dwarfs (Mamajek+, 2008)
*================================================================================
*Improved age estimation for solar-type dwarfs using activity-rotation
*diagnostics.
* Mamajek E.E., Hillenbrand L.A.
* <Astrophys. J., 687, 1264-1293 (2008)>
* =2008ApJ...687.1264M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! logR'_HK_ for color-separated solar-type
dwarf binaries
integer*4 nr__
parameter (nr__=21) ! Number of records
character*58 ar__ ! Full-size record
character*10 StarA (nr__) ! Name of primary
character*10 StarB (nr__) ! Name of secondary
real*4 B_V_A (nr__) ! (mag) B-V colour index for A component
real*4 B_V_B (nr__) ! (mag) B-V colour index for B component
real*4 logR_HKA (nr__) ! R'HK index for A component
real*4 logR_HKB (nr__) ! R'HK index for B component
character*13 Ref (nr__) ! References, in refs.dat file (the '*' indicates
* values converted to logR'_HK_ by the author
* following Noyes et al. (1984ApJ...279..763N))
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! logR'_HK_ for near-identical solar-type
dwarf binaries
integer*4 nr__1
parameter (nr__1=14) ! Number of records
character*58 ar__1 ! Full-size record
character*10 StarA_1 (nr__1) ! Name of primary
character*10 StarB_1 (nr__1) ! Name of secondary
real*4 B_V_A_1 (nr__1) ! (mag) B-V colour index for A component
real*4 B_V_B_1 (nr__1) ! (mag) B-V colour index for B component
real*4 logR_HKA_1 (nr__1) ! R'HK index for A component
real*4 logR_HKB_1 (nr__1) ! R'HK index for B component
character*13 Ref_1 (nr__1) ! References, in refs.dat file (the '*' indicates
* values converted to logR'_HK_ by the author
* following Noyes et al. (1984ApJ...279..763N))
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Field binaries with rotation periods
integer*4 nr__2
parameter (nr__2=5) ! Number of records
character*57 ar__2 ! Full-size record
character*10 StarA_2 (nr__2) ! Name of primary
character*10 StarB_2 (nr__2) ! Name of secondary
real*4 B_V_A_2 (nr__2) ! (mag) B-V colour index for A component (1)
real*4 B_V_B_2 (nr__2) ! (mag) B-V colour index for B component
real*4 Per_A (nr__2) ! (d) Period of A component
real*4 Per_B (nr__2) ! (d) Period of B component
character*13 Ref_2 (nr__2) ! References, in refs.dat file
*Note (1): The period for HD 128620 (alpha Cen A) is a mean (25.6days) from
* values given by E. Guinan (2008, priv. com.; 22+/-3days) and Hallam et
* al. (1991ApJ...372..610H; 28.8+/-2.5days) and is consistent within the
* constraints from vsini and p-mode rotational splitting (Fletcher et
* al. 2006MNRAS.371..935F; Bazot et al. 2007A&A...470..295B).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Members of stellar aggregates with logR'_HK_
measurements
integer*4 nr__3
parameter (nr__3=274) ! Number of records
character*108 ar__3 ! Full-size record
character*17 Name (nr__3) ! Name
character*23 SName (nr__3) ! Simbad Name
character*10 OName (nr__3) ! HD name
real*4 B_V (nr__3) ! (mag) ? B-V colour index
integer*4 r_B_V (nr__3) ! ? Reference for B-V, in refs.dat file
real*4 E_B_V (nr__3) ! (mag) Interstellar reddening
character*6 r_E_B_V (nr__3) ! Reference for E(B-V), in refs.dat file
real*4 v_B_V_0 (nr__3) ! (mag) Intrinsic B-V colour index
character*6 r__B_V_0 (nr__3) ! Reference for (B-V)0, in refs.dat file
real*4 logR_HK (nr__3) ! R'HK index
integer*4 o_logR_HK (nr__3) ! ? Number of measurements
integer*4 r_logR_HK (nr__3) ! Reference for logR'HK, in refs.dat file
character*9 Group (nr__3) ! Group name (1)
*Note (1): the groups are (from table6 of paper):
* -----------------------------------------------------------------------
* Group Age log RHK 68% N* Usual name
* Name Myr median
* -----------------------------------------------------------------------
* US 5 -4.05+/-0.03 0.13 9 Upper Sco association
* beta Pic 12 -4.03+/-0.13 0.23 6 {beta} Pic association
* UCL 16 -4.04+/-0.01 0.07 10 Upper Cen-Lup and Lower Cen-Cru
* Tuc 30 -4.16+/-0.13 0.16 8 Tuc-Hor association
* alpha Per 85 -4.16+/-0.08 0.27 13 {alpha} Per cluster
* Pleiades 130 -4.33+/-0.04 0.24 56 Pleiades cluster
* UMa 500 -4.48+/-0.03 0.09 17 UMa Group
* Hyades 625 -4.47+/-0.01 0.09 87 Hyades cluster
* M67 4000 -4.84+/-0.01 0.11 76 M67 open cluster
* -----------------------------------------------------------------------
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Old solar-type dwarfs from Valenti et al. (2005,
Cat. J/ApJS/159/141) with age uncertainties of <20%
integer*4 nr__4
parameter (nr__4=23) ! Number of records
character*32 ar__4 ! Full-size record
integer*4 HD (nr__4) ! HD number (Cat. III/135)
real*4 B_V_1 (nr__4) ! (mag) B-V color from Perryman et al.
* (HIP, 1997, Cat. I/239)
real*4 tau (nr__4) ! (Gyr) Isochronal age (1)
real*4 logR_HK_1 (nr__4) ! R'_HK_ Chromospheric activity
integer*4 r_logR_HK_1(nr__4) ! Reference for logR'HK, in refs.dat file
real*4 DVMAG (nr__4) ! (mag) Difference between stellar absolute magnitude
* and that for an MS star of the same B-V color
*Note (1): from Valenti et al. (2005, Cat. J/ApJS/159/141>), uncertainties <20%.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Activity-gyro ages for solar-type binaries
integer*4 nr__5
parameter (nr__5=17) ! Number of records
character*47 ar__5 ! Full-size record
character*10 StarA_3 (nr__5) ! Name of primary
character*10 StarB_3 (nr__5) ! Name of secondary
real*4 logtA (nr__5) ! ([yr]) Activity-gyro age for component A
real*4 logtB (nr__5) ! ([yr]) Activity-gyro age for component B
real*4 logt (nr__5) ! ([yr]) Mean gyro age for the system
real*4 e_logt (nr__5) ! ([yr]) rms uncertainty on logt
character*1 n_StarA (nr__5) ! [a] Note for HD 137763 (1)
*Note (1): a: HD 137763 is a pathological case discussed in footnote 10.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Activity ages for the solar-type dwarfs within 16pc
integer*4 nr__6
parameter (nr__6=103) ! Number of records
character*110 ar__6 ! Full-size record
integer*4 HD_1 (nr__6) ! HD number (Cat. III/135)
integer*4 HIP (nr__6) ! HIP number (Cat. I/311)
character*7 GJ (nr__6) ! Nearby Star Catalog designation (Cat. V/70)
character*11 OName_1 (nr__6) ! Other name
real*4 plx (nr__6) ! (mas) Parallax
real*4 e_plx (nr__6) ! (mas) rms uncertainty on plx
real*4 B_V_2 (nr__6) ! (mag) B-V colour index
integer*4 r_B_V_1 (nr__6) ! Reference for B-V, in refs.dat file
real*4 logR_HK_2 (nr__6) ! R'_HK_ index
integer*4 r_logR_HK_2(nr__6) ! Reference for logR'HK, in refs.dat file
real*4 Vmag (nr__6) ! (mag) V magnitude
real*4 VMAG_1 (nr__6) ! (mag) Absolute V magnitude
real*4 DVMAG_1 (nr__6) ! (mag) Difference between stellar absolute magnitude
* and that for an MS star of the same B-V color
character*15 SpType (nr__6) ! MK spectral type
integer*4 r_SpType (nr__6) ! Reference for SpType, in refs.dat file
character*1 l_tau1 (nr__6) ! Limit flag on tau1
real*4 tau1 (nr__6) ! (Gyr) Age from using the revised activity-age
* relation
character*1 l_tau2 (nr__6) ! Limit flag on tau2
real*4 tau2 (nr__6) ! (Gyr) Age inferred from converting the chromospheric
* activity levels to a rotation period via the
* Rossby number and then converting the
* rotation period to an age using the revised
* gyro relation
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References and notes
integer*4 nr__7
parameter (nr__7=64) ! Number of records
character*185 ar__7 ! Full-size record
integer*4 Ref_3 (nr__7) ! Reference code
character*19 BibCode (nr__7) ! BibCode
character*21 Aut (nr__7) ! Author's name
character*139 Com (nr__7) ! Comments
C=============================================================================
C Loading file 'table2.dat' ! logR'_HK_ for color-separated solar-type
* dwarf binaries
C Format for file interpretation
1 format(A10,1X,A10,2X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,21
read(1,'(A58)')ar__
read(ar__,1)
+ StarA(i__),StarB(i__),B_V_A(i__),B_V_B(i__),logR_HKA(i__),
+ logR_HKB(i__),Ref(i__)
c ..............Just test output...........
write(6,1)
+ StarA(i__),StarB(i__),B_V_A(i__),B_V_B(i__),logR_HKA(i__),
+ logR_HKB(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! logR'_HK_ for near-identical solar-type
* dwarf binaries
C Format for file interpretation
2 format(A10,1X,A10,2X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,14
read(1,'(A58)')ar__1
read(ar__1,2)
+ StarA_1(i__),StarB_1(i__),B_V_A_1(i__),B_V_B_1(i__),
+ logR_HKA_1(i__),logR_HKB_1(i__),Ref_1(i__)
c ..............Just test output...........
write(6,2)
+ StarA_1(i__),StarB_1(i__),B_V_A_1(i__),B_V_B_1(i__),
+ logR_HKA_1(i__),logR_HKB_1(i__),Ref_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Field binaries with rotation periods
C Format for file interpretation
3 format(A10,1X,A10,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,5
read(1,'(A57)')ar__2
read(ar__2,3)
+ StarA_2(i__),StarB_2(i__),B_V_A_2(i__),B_V_B_2(i__),
+ Per_A(i__),Per_B(i__),Ref_2(i__)
c ..............Just test output...........
write(6,3)
+ StarA_2(i__),StarB_2(i__),B_V_A_2(i__),B_V_B_2(i__),
+ Per_A(i__),Per_B(i__),Ref_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Members of stellar aggregates with logR'_HK_
* measurements
C Format for file interpretation
4 format(
+ A17,A23,1X,A10,1X,F4.2,1X,I2,1X,F4.2,1X,A6,2X,F4.2,1X,A6,2X,
+ F5.2,1X,I3,1X,I2,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,274
read(1,'(A108)')ar__3
read(ar__3,4)
+ Name(i__),SName(i__),OName(i__),B_V(i__),r_B_V(i__),
+ E_B_V(i__),r_E_B_V(i__),v_B_V_0(i__),r__B_V_0(i__),
+ logR_HK(i__),o_logR_HK(i__),r_logR_HK(i__),Group(i__)
if(ar__3(53:56) .EQ. '') B_V(i__) = rNULL__
if(ar__3(58:59) .EQ. '') r_B_V(i__) = iNULL__
if(ar__3(93:95) .EQ. '') o_logR_HK(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ Name(i__),SName(i__),OName(i__),B_V(i__),r_B_V(i__),
+ E_B_V(i__),r_E_B_V(i__),v_B_V_0(i__),r__B_V_0(i__),
+ logR_HK(i__),o_logR_HK(i__),r_logR_HK(i__),Group(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Old solar-type dwarfs from Valenti et al. (2005,
* Cat. J/ApJS/159/141) with age uncertainties of <20%
C Format for file interpretation
5 format(I6,1X,F5.3,1X,F4.1,1X,F5.2,1X,I2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,23
read(1,'(A32)')ar__4
read(ar__4,5)
+ HD(i__),B_V_1(i__),tau(i__),logR_HK_1(i__),r_logR_HK_1(i__),
+ DVMAG(i__)
c ..............Just test output...........
write(6,5)
+ HD(i__),B_V_1(i__),tau(i__),logR_HK_1(i__),r_logR_HK_1(i__),
+ DVMAG(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Activity-gyro ages for solar-type binaries
C Format for file interpretation
6 format(A10,3X,A10,3X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,17
read(1,'(A47)')ar__5
read(ar__5,6)
+ StarA_3(i__),StarB_3(i__),logtA(i__),logtB(i__),logt(i__),
+ e_logt(i__),n_StarA(i__)
c ..............Just test output...........
write(6,6)
+ StarA_3(i__),StarB_3(i__),logtA(i__),logtB(i__),logt(i__),
+ e_logt(i__),n_StarA(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Activity ages for the solar-type dwarfs within 16pc
C Format for file interpretation
7 format(
+ I6,1X,I6,1X,A7,1X,A11,1X,F6.2,1X,F4.2,1X,F5.3,1X,I1,1X,F6.3,
+ 1X,I2,1X,F5.2,1X,F4.2,1X,F5.2,1X,A15,1X,I2,1X,A1,F3.1,1X,A1,
+ F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,103
read(1,'(A110)')ar__6
read(ar__6,7)
+ HD_1(i__),HIP(i__),GJ(i__),OName_1(i__),plx(i__),e_plx(i__),
+ B_V_2(i__),r_B_V_1(i__),logR_HK_2(i__),r_logR_HK_2(i__),
+ Vmag(i__),VMAG_1(i__),DVMAG_1(i__),SpType(i__),r_SpType(i__),
+ l_tau1(i__),tau1(i__),l_tau2(i__),tau2(i__)
c ..............Just test output...........
write(6,7)
+ HD_1(i__),HIP(i__),GJ(i__),OName_1(i__),plx(i__),e_plx(i__),
+ B_V_2(i__),r_B_V_1(i__),logR_HK_2(i__),r_logR_HK_2(i__),
+ Vmag(i__),VMAG_1(i__),DVMAG_1(i__),SpType(i__),r_SpType(i__),
+ l_tau1(i__),tau1(i__),l_tau2(i__),tau2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References and notes
C Format for file interpretation
8 format(I2,1X,A19,1X,A21,2X,A139)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,64
read(1,'(A185)')ar__7
read(ar__7,8)Ref_3(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,8)Ref_3(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end