Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/422/527 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/A+A/422/527 Metal-poor star uvby-beta photometry. X. (Schuster+, 2004)
*================================================================================
*uvby-beta photometry of high-velocity and metal-poor stars.
*X. Stars of very low metal abundance: observations, reddenings, metallicities,
*classifications, distances, and relative ages.
* Schuster W.J., Beers T.C., Michel R., Nissen P.E., Garcia G.
* <Astron. Astrophys. 422, 527 (2004)>
* =2004A&A...422..527S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! uvby-{beta} data from SPM
integer*4 nr__
parameter (nr__=245) ! Number of records
character*78 ar__ ! Full-size record
character*9 Name ! Star designation
real*4 Vmag ! (mag) ? Johnson V magnitude
real*4 b_y ! (mag) Stroemgren colour index b-y
real*4 m1 ! (mag) Stroemgren metallicity index m_1_
real*4 c1 ! (mag) Stroemgren Balmer discontinuity index c_1_
integer*4 o_b_y ! Number of independent observations in uvby
real*4 Hbeta ! (mag) ? H{beta} line index
integer*4 o_Hbeta ! ? Number of independent observations in Hbeta
character*30 Notes ! Notes (1)
*Note (1): Notes:
* V,V? = variable(?)
* RG = red subgiant-giant
* fnt(r) = faint(er)
* st = star
* os = outside
* dph = diaphragm
* D,D? = double(?)
* wo = without
* excl = exclude
* obs = observation(s)
* contm = contaminated
* ocntr = offcenter
* js = just
* S = south
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! uvby data from La Silla
integer*4 nr__1
parameter (nr__1=177) ! Number of records
character*78 ar__1 ! Full-size record
character*9 Name_1 ! Star designation
real*4 Vmag_1 ! (mag) ? Johnson V magnitude
real*4 b_y_1 ! (mag) Stroemgren colour index b-y
real*4 m1_1 ! (mag) Stroemgren metallicity index m_1_
real*4 c1_1 ! (mag) Stroemgren Balmer discontinuity index c_1_
integer*4 o_b_y_1 ! Number of independent observations in uvby
real*4 Hbeta_1 ! (mag) ? H{beta} line index
integer*4 o_Hbeta_1 ! ? Number of independent observations in Hbeta
character*30 Notes_1 ! Notes (1)
*Note (1): Notes:
* V,V? = variable(?)
* RG = red subgiant-giant
* fnt(r) = faint(er)
* st = star
* os = outside
* dph = diaphragm
* D,D? = double(?)
* wo = without
* excl = exclude
* obs = observation(s)
* contm = contaminated
* ocntr = offcenter
* js = just
* S = south
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Coordinates, BV photometry, reddening excesses,
(B-V)_0_'s, and Notes
integer*4 nr__2
parameter (nr__2=497) ! Number of records
character*83 ar__2 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000.0
real*8 DEdeg ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*9 Name_2 ! Star name
integer*4 RAh ! (h) Right ascension (J2000.0)
integer*4 RAm ! (min) Right ascension (J2000.0)
real*4 RAs ! (s) Right ascension (J2000.0)
character*1 DE_ ! Declination sign (J2000.0)
integer*4 DEd ! (deg) Declination (J2000.0)
integer*4 DEm ! (arcmin) Declination (J2000.0)
integer*4 DEs ! (arcsec) Declination (J2000.0)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 Vmag_2 ! (mag) ? V magnitude
real*4 B_V ! (mag) ? B-V colour index
real*4 E_B_V_S ! (mag) Interstellar reddening from Schlegel et al.,
* 1998ApJ...500..525S
real*4 E_B_V ! (mag) Adopted interstellar reddening in (B-V)
* from Eq. 1
real*4 v_B_V_0 ! (mag) ? De-reddened Johnson color index B-V
real*4 v_B_V_0b ! (mag) ? Estimated de-reddened Johnson color index
* B-V, derived from (b-y)_0, Fig. 4
character*1 CH ! [*] *: star is likely carbon enhanced,
* from Fig. 5
character*1 OUT ! [*] *: star is considered a photometric
* outlier, from Figs. 1, 4
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Cross-Identifications of VMP Stars
integer*4 nr__3
parameter (nr__3=53) ! Number of records
character*31 ar__3 ! Full-size record
character*9 Name_3 ! Star name
character*9 Name2 ! Second alternative name of star
character*9 Name3 ! Third alternative name of star
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Dereddened photometry, Reddening Excesses,
and Photometric Classifications
integer*4 nr__4
parameter (nr__4=508) ! Number of records
character*61 ar__4 ! Full-size record
character*9 Name_4 ! Star name
real*4 V0 ! (mag) ? Visual magnitude corrected for interstellar
* absorption (1)
real*4 v_b_y_0 ! (mag) De-reddened Stroemgren color index b-y (1)
real*4 m0 ! (mag) De-reddened Stroemgren metallicity index m1 (1)
real*4 c0 ! (mag) De-reddened Stroemgren Balmer discontinuity
* index c1 (1)
real*4 E_b_y ! (mag) ? Stroemgren interstellar reddening estimate
* from intrinsic-color calibration of
* Schuster & Nissen (1989, Cat. <II/180>)
real*4 E_B_V_1 ! (mag) ? Adopted interstellar reddening in (B-V),
* from eq. 1
character*11 Class ! Photometric classification of star, Fig.6 (G1)
*Note (1): V, b-y, m1 and c1 have been corrected for the effect of
* interstellar reddening if E(b-y)>=0.015 or E(B-V)>=0.02.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Line-strength indices from the HK survey
integer*4 nr__5
parameter (nr__5=497) ! Number of records
character*31 ar__5 ! Full-size record
character*9 Name_5 ! Star name
real*4 KP ! (0.1nm) ? Ca II K line index
real*4 HP2 ! (0.1nm) ? Balmer H{delta} line index
real*4 HG2 ! (0.1nm) ? Balmer H{gamma} line index
real*4 GP ! (0.1nm) ? CH G-band line index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! [Fe/H] values from different methods,
plus adopted average
integer*4 nr__6
parameter (nr__6=497) ! Number of records
character*41 ar__6 ! Full-size record
character*9 Name_6 ! Star name
real*4 v_Fe_H_K1 ! ([Sun]) ? Metallicity based on KP, (B-V)0b
* (Beers et al., 1999, Cat. <J/AJ/117/981>)
real*4 v_Fe_H_A1 ! ([Sun]) ? Metallicity based on KP, (B-V)0b, and an
* Artificial Neural Network
real*4 v_Fe_H_K2 ! ([Sun]) ? Metallicity based on KP, (B-V)0
* (Beers et al., 1999, Cat. <J/AJ/117/981>)
real*4 v_Fe_H_A2 ! ([Sun]) ? Metallicity based on KP, (B-V)0, and an
* Artificial Neural Network
real*4 v_Fe_H_ ! ([Sun]) ? Final estimate of metallicity
integer*4 Notes_2 ! [1/3]? Individual notes (1)
*Note (1): Individual notes:
* 1: Late-type star, CaII-K emission
* 2: Composite spectrum (Early + Late)
* 3: CD-38:245
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8a.dat' ! Metal-rich BS (Am) stars - observed quantities
integer*4 nr__7
parameter (nr__7=10) ! Number of records
character*68 ar__7 ! Full-size record
character*9 Name_7 ! Star name
real*4 VJmag ! (mag) Johnson V magnitude
real*4 b_y_2 ! (mag) Stroemgren colour index b-y
real*4 m1_2 ! (mag) Stroemgren metallicity index m_1_
real*4 c1_2 ! (mag) Stroemgren Balmer discontinuity index c_1_
real*4 Hbeta_2 ! (mag) H{beta} line index
real*4 v_c1_ ! (mag) Reddening-free Balmer discontinuity index
real*4 v_m1_ ! (mag) Reddening-free metallicity index
real*4 Vmag_3 ! (mag) ? Johnson broad-band V magnitude
real*4 B_V_1 ! (mag) ? Johnson broad-band B-V colour index
real*4 U_B ! (mag) ? Johnson broad-band U-B colour index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8b.dat' ! Metal-rich BS (Am) stars - derived quantities
integer*4 nr__8
parameter (nr__8=10) ! Number of records
character*59 ar__8 ! Full-size record
character*9 Name_8 ! Star name
real*4 E_B_V_1 ! (mag) Adopted interstellar reddening in (B-V)
* from Eq. 1
real*4 v_B_V_0_1 ! (mag) De-reddened Johnson color index B-V
character*1 u__B_V_0 ! [)] Uncertainty flag on (B-V)0
real*4 v_U_B_0 ! (mag) ? De-reddened Johnson color index U-B
real*4 V0_1 ! (mag) V magnitude corrected for interstellar
* absorption
real*4 v_b_y_0_1 ! (mag) De-reddened Stroemgren color index b-y
real*4 m0_1 ! (mag) De-reddened Stroemgren metallicity index m1
real*4 c0_1 ! (mag) De-reddened Stroemgren Balmer discontinuity
* index c1
character*2 Class_1 ! Photometric classification of star (G1)
real*4 v_Fe_H__1 ! ([Sun]) ? Final estimate of metallicity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Photometric absolute magnitudes and distances
of the VMP stars
integer*4 nr__9
parameter (nr__9=497) ! Number of records
character*69 ar__9 ! Full-size record
character*9 Name_9 ! Star name
character*11 Class_2 ! Photometric classification of star (G1)
real*4 VMAG1 ! (mag) ? Estimated absolute magnitude, assuming primary
* classification (MS for TO stars, as listed
* for others) and UBV photometry
real*4 VMAG2 ! (mag) ? Estimated absolute magnitude, assuming secondary
* classification (SG or RG for TO stars) and UBV
* photometry
integer*4 Dist1 ! (pc) ? Estimated distance, assuming primary
* classification (MS for TO stars, as listed for
* others) and UBV photometry
integer*4 Dist2 ! (pc) ? Estimated distance, assuming secondary
* classification (SG or RG for TO stars) and UBV
* photometry
integer*4 Dist ! (pc) ? Estimated distance based on Stroemgren
* photometry and indices
character*14 Notes_3 ! Comments on individual stars
C=============================================================================
C Loading file 'table1.dat' ! uvby-{beta} data from SPM
C Format for file interpretation
1 format(
+ A9,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,I1,1X,F5.3,1X,I1,1X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,245
read(1,'(A78)')ar__
read(ar__,1)Name,Vmag,b_y,m1,c1,o_b_y,Hbeta,o_Hbeta,Notes
if(ar__(11:16) .EQ. '') Vmag = rNULL__
if(ar__(41:45) .EQ. '') Hbeta = rNULL__
if(ar__(47:47) .EQ. '') o_Hbeta = iNULL__
c ..............Just test output...........
write(6,1)Name,Vmag,b_y,m1,c1,o_b_y,Hbeta,o_Hbeta,Notes
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! uvby data from La Silla
C Format for file interpretation
2 format(
+ A9,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,I1,1X,F5.3,1X,I1,1X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,177
read(1,'(A78)')ar__1
read(ar__1,2)
+ Name_1,Vmag_1,b_y_1,m1_1,c1_1,o_b_y_1,Hbeta_1,o_Hbeta_1,
+ Notes_1
if(ar__1(11:16) .EQ. '') Vmag_1 = rNULL__
if(ar__1(41:45) .EQ. '') Hbeta_1 = rNULL__
if(ar__1(47:47) .EQ. '') o_Hbeta_1 = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,Vmag_1,b_y_1,m1_1,c1_1,o_b_y_1,Hbeta_1,o_Hbeta_1,
+ Notes_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Coordinates, BV photometry, reddening excesses,
* (B-V)_0_'s, and Notes
C Format for file interpretation
3 format(
+ A9,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F5.1,1X,F5.1,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,A1,2X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,497
read(1,'(A83)')ar__2
read(ar__2,3)
+ Name_2,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,Vmag_2,B_V,
+ E_B_V_S,E_B_V,v_B_V_0,v_B_V_0b,CH,OUT
if(ar__2(44:48) .EQ. '') Vmag_2 = rNULL__
if(ar__2(50:54) .EQ. '') B_V = rNULL__
if(ar__2(68:72) .EQ. '') v_B_V_0 = rNULL__
if(ar__2(74:78) .EQ. '') v_B_V_0b = rNULL__
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,3)
+ Name_2,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,Vmag_2,B_V,
+ E_B_V_S,E_B_V,v_B_V_0,v_B_V_0b,CH,OUT
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Cross-Identifications of VMP Stars
C Format for file interpretation
4 format(A9,2X,A9,2X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,53
read(1,'(A31)')ar__3
read(ar__3,4)Name_3,Name2,Name3
c ..............Just test output...........
write(6,4)Name_3,Name2,Name3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Dereddened photometry, Reddening Excesses,
* and Photometric Classifications
C Format for file interpretation
5 format(
+ A9,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F4.2,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,508
read(1,'(A61)')ar__4
read(ar__4,5)Name_4,V0,v_b_y_0,m0,c0,E_b_y,E_B_V_1,Class
if(ar__4(11:16) .EQ. '') V0 = rNULL__
if(ar__4(39:44) .EQ. '') E_b_y = rNULL__
if(ar__4(46:49) .EQ. '') E_B_V_1 = rNULL__
c ..............Just test output...........
write(6,5)Name_4,V0,v_b_y_0,m0,c0,E_b_y,E_B_V_1,Class
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Line-strength indices from the HK survey
C Format for file interpretation
6 format(A9,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,497
read(1,'(A31)')ar__5
read(ar__5,6)Name_5,KP,HP2,HG2,GP
if(ar__5(11:14) .EQ. '') KP = rNULL__
if(ar__5(16:20) .EQ. '') HP2 = rNULL__
if(ar__5(22:26) .EQ. '') HG2 = rNULL__
if(ar__5(28:31) .EQ. '') GP = rNULL__
c ..............Just test output...........
write(6,6)Name_5,KP,HP2,HG2,GP
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! [Fe/H] values from different methods,
* plus adopted average
C Format for file interpretation
7 format(A9,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,497
read(1,'(A41)')ar__6
read(ar__6,7)
+ Name_6,v_Fe_H_K1,v_Fe_H_A1,v_Fe_H_K2,v_Fe_H_A2,v_Fe_H_,Notes_2
if(ar__6(11:15) .EQ. '') v_Fe_H_K1 = rNULL__
if(ar__6(17:21) .EQ. '') v_Fe_H_A1 = rNULL__
if(ar__6(23:27) .EQ. '') v_Fe_H_K2 = rNULL__
if(ar__6(29:33) .EQ. '') v_Fe_H_A2 = rNULL__
if(ar__6(35:39) .EQ. '') v_Fe_H_ = rNULL__
if(ar__6(41:41) .EQ. '') Notes_2 = iNULL__
c ..............Just test output...........
write(6,7)
+ Name_6,v_Fe_H_K1,v_Fe_H_A1,v_Fe_H_K2,v_Fe_H_A2,v_Fe_H_,Notes_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8a.dat' ! Metal-rich BS (Am) stars - observed quantities
C Format for file interpretation
8 format(
+ A9,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8a.dat')
write(6,*) '....Loading file: table8a.dat'
do i__=1,10
read(1,'(A68)')ar__7
read(ar__7,8)
+ Name_7,VJmag,b_y_2,m1_2,c1_2,Hbeta_2,v_c1_,v_m1_,Vmag_3,B_V_1,
+ U_B
if(ar__7(54:58) .EQ. '') Vmag_3 = rNULL__
if(ar__7(60:63) .EQ. '') B_V_1 = rNULL__
if(ar__7(65:68) .EQ. '') U_B = rNULL__
c ..............Just test output...........
write(6,8)
+ Name_7,VJmag,b_y_2,m1_2,c1_2,Hbeta_2,v_c1_,v_m1_,Vmag_3,B_V_1,
+ U_B
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8b.dat' ! Metal-rich BS (Am) stars - derived quantities
C Format for file interpretation
9 format(
+ A9,1X,F4.2,1X,F4.2,A1,1X,F4.2,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,A2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8b.dat')
write(6,*) '....Loading file: table8b.dat'
do i__=1,10
read(1,'(A59)')ar__8
read(ar__8,9)
+ Name_8,E_B_V_1,v_B_V_0_1,u__B_V_0,v_U_B_0,V0_1,v_b_y_0_1,m0_1,
+ c0_1,Class_1,v_Fe_H__1
if(ar__8(22:25) .EQ. '') v_U_B_0 = rNULL__
if(ar__8(55:59) .EQ. '') v_Fe_H__1 = rNULL__
c ..............Just test output...........
write(6,9)
+ Name_8,E_B_V_1,v_B_V_0_1,u__B_V_0,v_U_B_0,V0_1,v_b_y_0_1,m0_1,
+ c0_1,Class_1,v_Fe_H__1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Photometric absolute magnitudes and distances
* of the VMP stars
C Format for file interpretation
10 format(A9,1X,A11,2X,F6.3,1X,F6.3,1X,I5,1X,I5,1X,I5,1X,A14)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,497
read(1,'(A69)')ar__9
read(ar__9,10)
+ Name_9,Class_2,VMAG1,VMAG2,Dist1,Dist2,Dist,Notes_3
if(ar__9(24:29) .EQ. '') VMAG1 = rNULL__
if(ar__9(31:36) .EQ. '') VMAG2 = rNULL__
if(ar__9(38:42) .EQ. '') Dist1 = iNULL__
if(ar__9(44:48) .EQ. '') Dist2 = iNULL__
if(ar__9(50:54) .EQ. '') Dist = iNULL__
c ..............Just test output...........
write(6,10)
+ Name_9,Class_2,VMAG1,VMAG2,Dist1,Dist2,Dist,Notes_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end