Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/684/1143 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-18
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/684/1143 BHB candidates in the Milky Way (Xue+, 2008)
*================================================================================
*The Milky Way's circular velocity curve to 60 kpc and an estimate of the dark
*matter halo mass from the kinematics of ~2400 SDSS blue horizontal-branch stars.
* Xue X.X., Rix H.W., Zhao G., Re Fiorentin P., Naab T., Steinmetz M.,
* van den Bosch F.C., Beers T.C., Lee Y.S., Bell E.F., Rockosi C., Yanny B.,
* Newberg H., Wilhelm R., Kang X., Smith M.C., Schneider D.P.
* <Astrophys. J., 684, 1143-1158 (2008)>
* =2008ApJ...684.1143X
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! List of 10224 BHB candidates selected by color cut
integer*4 nr__
parameter (nr__=10224) ! Number of records
character*282 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*14 SpID ! Spectrum name (MJD-Plate-Fiber)
character*19 SDSS ! SDSS-DR6 Name (JHHMMSS.ss+DDMMSS.s) (1)
real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
real*8 GLON ! (deg) Galactic Longitude
real*8 GLAT ! (deg) Galactic Latitude
real*4 gmag ! (mag) Extinction-corrected g band magnitude
real*4 g_r ! (mag) The g-r color index after extinction correction
real*4 u_g ! (mag) The u-g color index after extinction correction
real*4 e_umag ! (mag) Error of u band magnitude
real*4 e_gmag ! (mag) Error of g band magnitude
real*4 e_rmag ! (mag) Error of r band magnitude
real*4 HRV ! (km/s) Heliocentic radial velocity
real*4 e_HRV ! (km/s) Error in HRV
real*4 RVgal ! (km/s) Galactocentric radial velocity (2)
real*4 fm ! ?=-9.99 Flux of H{delta} line relative to the
* continuum at the line core
real*8 D0_2 ! (0.1nm) ?=-9.99 Width of H{delta} line at 20% below the
* local continuum
real*8 cgam ! ?=-9.99 The c of H{gamma} line (G1)
real*8 bgam ! (0.1nm) ?=-9.99 The b of H{gamma} line in Angstroms (G1)
character*3 Type ! Object type: BHB, BS or MS (3)
*Note (1): This sample contains some repeated stars.
*Note (2): Adopting solar motion (10.0,5.2,7.2)km/s and local rotation of 220km/s
*Note (3): Type is classified only by D0.2 & fm method, as;
* BHB = Blue Horizontal Branch;
* BS = Blue Straggler;
* MS = Main Sequence.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! List of 2558 BHB stars selected from SDSS DR6
integer*4 nr__1
parameter (nr__1=2558) ! Number of records
character*280 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*14 SpID_1 ! Spectrum name (MJD-Plate-Fiber)
character*19 SDSS_1 ! SDSS-DR6 Name (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_1 ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 ! (deg) Declination in decimal degrees (J2000)
real*8 GLON_1 ! (deg) Galactic Longitude
real*8 GLAT_1 ! (deg) Galactic Latitude
real*4 gmag_1 ! (mag) Extinction-corrected g band magnitude
real*4 u_g_1 ! (mag) The u-g color index after extinction correction
real*4 g_r_1 ! (mag) The g-r color index after extinction correction
real*4 gMAG_2 ! (mag) Absolute g band magnitude
real*4 D0_2_1 ! (0.1nm) Width of H{delta} line, in Angstroms, at 20%
* below the local continuum
real*4 fm_1 ! Flux of H{delta} line relative to the
* continuum at the line core
real*4 cgam_1 ! The c of H{gamma} line (G1)
real*4 bgam_1 ! (0.1nm) The b of H{gamma} line in Angstroms (G1)
real*4 Dist ! (kpc) Distance to the Sun
real*4 r ! (kpc) Distance to the Galactic center
real*4 x ! (kpc) The x galactocentric coordinate (toward Sun)
real*4 y ! (kpc) The y galactocentric coordinate
* (direction of rotation)
real*4 z ! (kpc) The z galactocentric coordinate (toward NGP)
real*4 HRV_1 ! (km.s-1) Heliocentic radial velocity
real*4 e_HRV_1 ! (km.s-1) Error in HRV
real*4 Vgal ! (km.s-1) Galactocentric radial velocity (1)
integer*4 Teff ! (K) ?=-9999 Effective temperature by Ron Wilhelm
real*4 logg ! ([cm/s2]) ?=-9.99 Surface gravity by Ron Wilhelm
real*4 FeH ! ([-]) ?=-9.99 Metallicity by Ron Wilhelm
character*3 Type_1 ! [BHB] Source type (Blue Horizontal Branch) (2)
*Note (1): Adopting solar motion (10.0,5.2,7.2)km/s and local rotation of 220km/s
*Note (2): BHB stars are classified by combination of D0.2 & fm method and the
* scale width vs shape method.
C=============================================================================
C Loading file 'table1.dat' ! List of 10224 BHB candidates selected by color cut
C Format for file interpretation
1 format(
+ A14,9X,A19,5X,F10.6,5X,F10.6,5X,F10.6,5X,F10.6,3X,F5.2,3X,
+ F5.2,4X,F4.2,4X,F4.2,3X,F5.2,4X,F4.2,9X,F6.1,10X,F5.1,9X,F6.1,
+ 15X,F5.2,8X,F12.2,12X,F8.2,11X,F9.2,4X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,10224
read(1,'(A282)')ar__
read(ar__,1)
+ SpID,SDSS,RAdeg,DEdeg,GLON,GLAT,gmag,g_r,u_g,e_umag,e_gmag,
+ e_rmag,HRV,e_HRV,RVgal,fm,D0_2,cgam,bgam,Type
c ..............Just test output...........
write(6,1)
+ SpID,SDSS,RAdeg,DEdeg,GLON,GLAT,gmag,g_r,u_g,e_umag,e_gmag,
+ e_rmag,HRV,e_HRV,RVgal,fm,D0_2,cgam,bgam,Type
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! List of 2558 BHB stars selected from SDSS DR6
C Format for file interpretation
2 format(
+ A14,9X,A19,10X,F10.6,10X,F10.6,10X,F10.6,10X,F10.6,3X,F5.2,4X,
+ F4.2,3X,F5.2,4X,F4.2,3X,F5.2,4X,F4.2,4X,F4.2,3X,F5.2,4X,F4.1,
+ 4X,F4.1,3X,F5.1,3X,F5.1,3X,F5.1,2X,F6.1,2X,F6.1,2X,F6.1,2X,I6,
+ 3X,F5.2,3X,F5.2,3X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,2558
read(1,'(A280)')ar__1
read(ar__1,2)
+ SpID_1,SDSS_1,RAdeg_1,DEdeg_1,GLON_1,GLAT_1,gmag_1,u_g_1,
+ g_r_1,gMAG_2,D0_2_1,fm_1,cgam_1,bgam_1,Dist,r,x,y,z,HRV_1,
+ e_HRV_1,Vgal,Teff,logg,FeH,Type_1
c ..............Just test output...........
write(6,2)
+ SpID_1,SDSS_1,RAdeg_1,DEdeg_1,GLON_1,GLAT_1,gmag_1,u_g_1,
+ g_r_1,gMAG_2,D0_2_1,fm_1,cgam_1,bgam_1,Dist,r,x,y,z,HRV_1,
+ e_HRV_1,Vgal,Teff,logg,FeH,Type_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end