Conversion of standardized ReadMe file for
catalog 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/MNRAS/352/285 u(bj)r & BV photometry of blue HB stars (Clewley+, 2004)
*================================================================================
*Distant field blue horizontal branch stars and the mass of the Galaxy -
*II. Photometry and spectroscopy of UKST candidates 16<B<19.5, 11<R<52kpc
* Clewley L., Warren S.J., Hewett P.C., Norris J.E., Evans N.W.
* <Mon. Not. R. Astron. Soc., 352, 285-298 (2004)>
* =2004MNRAS.352..285C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Photometric catalogue of A-type stars selected
from six UKST fields
integer*4 nr__
parameter (nr__=461) ! Number of records
character*64 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*9 v_CWH2004_ (nr__) ! Name
integer*4 RAh (nr__) ! (h) Right ascension (J2000.0)
integer*4 RAm (nr__) ! (min) Right ascension (J2000.0)
real*4 RAs (nr__) ! (s) Right ascension (J2000.0)
character*1 DE_ (nr__) ! Declination sign (J2000.0)
integer*4 DEd (nr__) ! (deg) Declination (J2000.0)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000.0)
real*4 DEs (nr__) ! (arcsec) Declination (J2000.0)
real*4 bjmag (nr__) ! (mag) bj magnitude (1)
real*4 u_bj (nr__) ! (mag) u-bj colour index (1)
real*4 bj_or (nr__) ! (mag) bj-or colour index (1)
real*4 e_u_bj (nr__) ! (mag) rms uncertainty on u-bj
real*4 e_bj_or (nr__) ! (mag) rms uncertainty on bj-or
*Note (1): UKST filters:
* ---------------------------------------------------------
* Name Wavelength (nm) at transmittances of
* 1/10 and 9/10 maximum values
* ---------------------------------------------------------
* bj : GG395 365 415
* u : UG1 315,390 330,375
* or : OG590 580 600
* ---------------------------------------------------------
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Photometric data for the blue horizontal-branch
star candidates, corrected for extinction
integer*4 nr__1
parameter (nr__1=280) ! Number of records
character*68 ar__1 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*9 v_CWH2004__1(nr__1) ! Name (1)
character*1 m__CWH2004_(nr__1) ! [qv] Multiplicity index on [CWH2004]
integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000.0)
integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000.0)
real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000.0)
character*1 DE__1 (nr__1) ! Declination sign (J2000.0)
integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000.0)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000.0)
real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000.0)
real*4 Vmag (nr__1) ! (mag) V magnitude
real*4 v_B_V_0 (nr__1) ! (mag) Dereddend B-V colour index
real*4 e_Vmag (nr__1) ! (mag) rms uncertainty on Vmag
real*4 e__B_V_0 (nr__1) ! (mag) rms uncertainty on (B-V)0
real*4 E_B_V (nr__1) ! (mag) Reddening in B-V
integer*4 N (nr__1) ! Number of photometric nights (2)
*Note (1): Suspected variables are labelled 'V'.
* Stars that were subsequently found to be quasars are labelled 'q'.
* Stars in the complete sample have names prefixed 'C', otherwise the
* stars are prefixed 'I'.
*Note (2): The dates and total number of photometric nights 'N' for the
* photometric survey in each of the six fields presented in Table 1.
* Observations were made using the 1.0-m Jacobus Kapteyn Telescope
* (JKT), the 2.5-m Isaac Newton Telescope (INT) and the 2.3-m Telescope
* at Siding Spring Observatories (SSO).
* ----------------------------------------------------------------------
* Telescope Date N Fields
* -----------------------------------------------------------------------
* SSO 1998 Sep 13-16 4 SGP
* JKT 1998 Oct 15, 17-20 5 SA94, MTF
* SSO 1999 Mar 12-18 7 F854, F789, SGP
* JKT 1999 April 4, 7-9 4 F854, F789
* SSO 1999 Oct 11-17 7 SGP
* JKT 1999 Oct 30-31, Sep 1, 3-5 6 SA94, MTF
* SSO 1999 Nov 17-20 4 F358
* SSO 2000 Feb 24 4 F854, F789
* JKT 2000 Mar 23, 26, 28 3 F854, F789
* INT 2001 Mar 29-30, Apr 1 3 SA94, MTF
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Spectroscopic data for the horizontal-branch
star candidates
integer*4 nr__2
parameter (nr__2=142) ! Number of records
character*88 ar__2 ! Full-size record
character*9 v_CWH2004__2(nr__2) ! Name
character*1 n__CWH2004_(nr__2) ! [Nbv] Multiplicity index on [CWH2004]
real*4 D0_15_gd (nr__2) ! (0.1nm) D0.15 parameter (1)
real*4 e_D0_15_gd (nr__2) ! (0.1nm) rms uncertainty on D0.15
real*4 b_gd (nr__2) ! (0.1nm) b parameter (1)
real*4 c_gd (nr__2) ! (0.1nm) c parameter (1)
real*4 A (nr__2) ! (0.1nm) Semi-major axes of the error ellipse
* in the b-c plane (2)
real*4 B (nr__2) ! (0.1nm) Semi-minor axes of the error ellipse
* in the b-c plane (2)
real*4 theta (nr__2) ! (rad) Orientation of the semi-major axis measured
* anti-clockwise from the b-axis (2)
real*4 v_Fe_H_ (nr__2) ! ([Sun]) ? Measured metallicity
real*4 e__Fe_H_ (nr__2) ! ([Sun]) ? rms uncertainty on [Fe/H]
real*8 HRV (nr__2) ! (km/s) Heliocentric radial velocity
real*4 e_HRV (nr__2) ! (km/s) rms uncertainty on HRV
real*4 Dist (nr__2) ! (kpc) Heliocentric distance
real*4 e_Dist (nr__2) ! (kpc) rms uncertainty on Dist
character*5 Class (nr__2) ! Classification as BHB (blue horizontal-branch)
* or A/BS (field blue straggler)
*Note (1): Two methods were developed in paper I (Clemens et al.,
* 2002MNRAS.337...87C) to separate the two populations of A stars. Both
* methods employ a Sersic function fit to the H{gamma} and H{delta}
* absorption lines.
* -- The first method, the 'D0.15-colour' method, plots the average of
* the width of the two Balmer lines against (B-V)0 colour.
* -- The second method, called the scale-width-shape method, plots two
* parameters of the Sersic fit, the scale-width b, against the
* power-law exponent c.
*Note (2): The errors on the parameters b(gd) and c(gd) are provided with
* the A, B and theta parameters of the error ellipse.
C=============================================================================
C Loading file 'table4.dat' ! Photometric catalogue of A-type stars selected
* from six UKST fields
C Format for file interpretation
1 format(
+ A9,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,461
read(1,'(A64)')ar__
read(ar__,1)
+ v_CWH2004_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),bjmag(i__),u_bj(i__),bj_or(i__),e_u_bj(i__),
+ e_bj_or(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ v_CWH2004_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),bjmag(i__),u_bj(i__),bj_or(i__),e_u_bj(i__),
+ e_bj_or(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Photometric data for the blue horizontal-branch
* star candidates, corrected for extinction
C Format for file interpretation
2 format(
+ A9,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,280
read(1,'(A68)')ar__1
read(ar__1,2)
+ v_CWH2004__1(i__),m__CWH2004_(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ Vmag(i__),v_B_V_0(i__),e_Vmag(i__),e__B_V_0(i__),E_B_V(i__),
+ N(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ v_CWH2004__1(i__),m__CWH2004_(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ Vmag(i__),v_B_V_0(i__),e_Vmag(i__),e__B_V_0(i__),E_B_V(i__),
+ N(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Spectroscopic data for the horizontal-branch
* star candidates
C Format for file interpretation
3 format(
+ A9,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.3,1X,F5.3,
+ 1X,F4.1,1X,F3.1,1X,F7.2,1X,F4.1,1X,F4.1,1X,F3.1,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,142
read(1,'(A88)')ar__2
read(ar__2,3)
+ v_CWH2004__2(i__),n__CWH2004_(i__),D0_15_gd(i__),
+ e_D0_15_gd(i__),b_gd(i__),c_gd(i__),A(i__),B(i__),theta(i__),
+ v_Fe_H_(i__),e__Fe_H_(i__),HRV(i__),e_HRV(i__),Dist(i__),
+ e_Dist(i__),Class(i__)
if(ar__2(52:55) .EQ. '') v_Fe_H_(i__) = rNULL__
if(ar__2(57:59) .EQ. '') e__Fe_H_(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ v_CWH2004__2(i__),n__CWH2004_(i__),D0_15_gd(i__),
+ e_D0_15_gd(i__),b_gd(i__),c_gd(i__),A(i__),B(i__),theta(i__),
+ v_Fe_H_(i__),e__Fe_H_(i__),HRV(i__),e_HRV(i__),Dist(i__),
+ e_Dist(i__),Class(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end