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.5, on 2013-May-24
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__
parame++++Only one proper motion (pm) found for: table[13].dat
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJS/180/67 Photometric selection of quasars from SDSS. II. (Richards+, 2009)
*================================================================================
*Efficient photometric selection of quasars from the Sloan Digital Sky Survey.
*II. ~1,000,000 quasars from Data Release 6.
* Richards G.T., Myers A.D., Gray A.G., Riegel R.N., Nichol R.C.,
* Brunner R.J., Szalay A.S., Schneider D.P., Anderson S.F.
* <Astrophys. J. Suppl. Ser., 180, 67-83 (2009)>
* =2009ApJS..180...67R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Nonparametric Bayes classifier kernel density
estimate (NBC KDE) quasar candidate catalog
integer*4 nr__
parameter (nr__=1015082) ! Number of records
character*271 ar__ ! Full-size record
integer*4 Seq (nr__) ! Unique catalog number
character*18 SDSS (nr__) ! SDSS name (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
character*19 objID (nr__) ! SDSS object identification
real*8 zph (nr__) ! Photometric redshift; see Weinstein et al.,
* 2004ApJS..155..243W
real*4 zph_L (nr__) ! Lower limit of zphot
real*4 zph_U (nr__) ! Upper limit of zphot
real*4 zpprob (nr__) ! Photometric redshift range probability
real*8 umag (nr__) ! (mag) u band PSF ubercalibrated asinh magnitude (1)
real*4 gmag (nr__) ! (mag) g band PSF ubercalibrated asinh magnitude (1)
real*4 rmag (nr__) ! (mag) r band PSF ubercalibrated asinh magnitude (1)
real*4 imag (nr__) ! (mag) i band PSF ubercalibrated asinh magnitude (1)
real*4 zmag (nr__) ! (mag) z band PSF ubercalibrated asinh magnitude (1)
real*4 e_umag (nr__) ! (mag) Error in umag
real*4 e_gmag (nr__) ! (mag) Error in gmag
real*4 e_rmag (nr__) ! (mag) Error in rmag
real*4 e_imag (nr__) ! (mag) Error in imag
real*4 e_zmag (nr__) ! (mag) Error in zmag
real*8 E_B_V (nr__) ! (mag) The (B-V) extinction (2)
real*8 c (nr__) ! (mag) Concentration for star/galaxy separation (3)
real*8 F20cm (nr__) ! (mJy) ?=-1.0 the 20cm (1.5GHz) flux density
real*8 FX (nr__) ! (ct/s) ?=-9.0 RASS full-band count rate (0.1-2keV)
real*8 pm (nr__) ! (mas/yr) ?=-1.0 Proper motion
integer*4 mov (nr__) ! [-1/1] Flag to indicate possible moving
* object (4)
integer*4 fz (nr__) ! [0/1] Full z range flag, 95% star prior
integer*4 lz (nr__) ! [0/1] Low z range flag (z<=2.2), 98% star prior
integer*4 mz (nr__) ! [0/1] Mid z range flag (2.2<z<3.5), 98% star
* prior
integer*4 hz (nr__) ! [0/1] High z range flag (z>=3.5), 98% star
* prior
integer*4 uv (nr__) ! [0/1] UV-excess flag, 88% star prior;
* see Paper I, Cat. <J/ApJS/155/257>
real*8 log_q (nr__) ! ([-]) Kernel density estimate (KDE) quasar probability
real*8 log_s (nr__) ! ([-]) Kernel density estimate (KDE) star probability
integer*4 good (nr__) ! [-6/6] Quality flag (5)
character*16 Class (nr__) ! Previous catalog object classification (6)
real*4 z (nr__) ! ?=-1.0 Previous catalog object redshift
*Note (1): Corrected for Galactic extinction.
*Note (2): A_u_, A_g_, A_r_, A_i_, and A_z_=(5.155, 3.793, 2.751, 2.086, and
* 1.479)*E(B-V).
*Note (3): Equal to PSFMag_i-modelMag_i.
*Note (4): Equal to 1 if moving.
*Note (5): Where 6 is most robust and -6 is least robust.
* Objects with good<0 are in table3 (rejected quasar candidates)
*Note (6): Previous catalog classification from sources:
* 2SLAQ = SDSS-2dF LRG and QSO Survey (2SLAQ) Early Data Release quasar catalog
* (Croom et al., 2009MNRAS.392...19C)
* DR5 = DR5 quasar catalog (Schneider et al., 2007, Cat. <VII/252>)
* 2QZ = 2QZ quasar catalog (Croom et al., 2004, Cat. <VII/241>)
* DR6 = SDSS-DR6 spectroscopic database (Adelman-McCarthy et al.,
* 2008ApJS..175..297A)
* U = Unknown
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Rejected quasar candidates
integer*4 nr__1
parameter (nr__1=157075) ! Number of records
character*271 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! Unique catalog number
character*18 SDSS_1 (nr__1) ! SDSS name (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
character*19 objID_1 (nr__1) ! SDSS object identification
real*8 zph_1 (nr__1) ! Photometric redshift; see Weinstein et al.,
* 2004ApJS..155..243W
real*4 zph_L_1 (nr__1) ! Lower limit of zphot
real*4 zph_U_1 (nr__1) ! Upper limit of zphot
real*4 zpprob_1 (nr__1) ! Photometric redshift range probability
real*8 umag_1 (nr__1) ! (mag) u band PSF ubercalibrated asinh magnitude (1)
real*4 gmag_1 (nr__1) ! (mag) g band PSF ubercalibrated asinh magnitude (1)
real*4 rmag_1 (nr__1) ! (mag) r band PSF ubercalibrated asinh magnitude (1)
real*4 imag_1 (nr__1) ! (mag) i band PSF ubercalibrated asinh magnitude (1)
real*4 zmag_1 (nr__1) ! (mag) z band PSF ubercalibrated asinh magnitude (1)
real*4 e_umag_1 (nr__1) ! (mag) Error in umag
real*4 e_gmag_1 (nr__1) ! (mag) Error in gmag
real*4 e_rmag_1 (nr__1) ! (mag) Error in rmag
real*4 e_imag_1 (nr__1) ! (mag) Error in imag
real*4 e_zmag_1 (nr__1) ! (mag) Error in zmag
real*8 E_B_V_1 (nr__1) ! (mag) The (B-V) extinction (2)
real*8 c_1 (nr__1) ! (mag) Concentration for star/galaxy separation (3)
real*8 F20cm_1 (nr__1) ! (mJy) ?=-1.0 the 20cm (1.5GHz) flux density
real*8 FX_1 (nr__1) ! (ct/s) ?=-9.0 RASS full-band count rate (0.1-2keV)
real*8 pm_1 (nr__1) ! (mas/yr) ?=-1.0 Proper motion
integer*4 mov_1 (nr__1) ! [-1/1] Flag to indicate possible moving
* object (4)
integer*4 fz_1 (nr__1) ! [0/1] Full z range flag, 95% star prior
integer*4 lz_1 (nr__1) ! [0/1] Low z range flag (z<=2.2), 98% star prior
integer*4 mz_1 (nr__1) ! [0/1] Mid z range flag (2.2<z<3.5), 98% star
* prior
integer*4 hz_1 (nr__1) ! [0/1] High z range flag (z>=3.5), 98% star
* prior
integer*4 uv_1 (nr__1) ! [0/1] UV-excess flag, 88% star prior;
* see Paper I, Cat. <J/ApJS/155/257>
real*8 log_q_1 (nr__1) ! ([-]) Kernel density estimate (KDE) quasar probability
real*8 log_s_1 (nr__1) ! ([-]) Kernel density estimate (KDE) star probability
integer*4 good_1 (nr__1) ! [-6/6] Quality flag (5)
character*16 Class_1 (nr__1) ! Previous catalog object classification (6)
real*4 z_1 (nr__1) ! ?=-1.0 Previous catalog object redshift
*Note (1): Corrected for Galactic extinction.
*Note (2): A_u_, A_g_, A_r_, A_i_, and A_z_=(5.155, 3.793, 2.751, 2.086, and
* 1.479)*E(B-V).
*Note (3): Equal to PSFMag_i-modelMag_i.
*Note (4): Equal to 1 if moving.
*Note (5): Where 6 is most robust and -6 is least robust.
* Objects with good<0 are in table3 (rejected quasar candidates)
*Note (6): Previous catalog classification from sources:
* 2SLAQ = SDSS-2dF LRG and QSO Survey (2SLAQ) Early Data Release quasar catalog
* (Croom et al., 2009MNRAS.392...19C)
* DR5 = DR5 quasar catalog (Schneider et al., 2007, Cat. <VII/252>)
* 2QZ = 2QZ quasar catalog (Croom et al., 2004, Cat. <VII/241>)
* DR6 = SDSS-DR6 spectroscopic database (Adelman-McCarthy et al.,
* 2008ApJS..175..297A)
* U = Unknown
C=============================================================================
C Loading file 'table1.dat' ! Nonparametric Bayes classifier kernel density
* estimate (NBC KDE) quasar candidate catalog
C Format for file interpretation
1 format(
+ I7,1X,A18,1X,F12.7,1X,F11.7,2X,A19,1X,F7.3,1X,F6.3,1X,F6.3,1X,
+ F6.3,1X,F7.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F7.3,1X,F7.3,1X,F8.2,1X,F7.4,1X,
+ F8.2,1X,I2,2X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.3,1X,F8.3,2X,
+ I2,1X,A16,1X,F6.3)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,1015082
read(1,'(A271)')ar__
read(ar__,1)
+ Seq(i__),SDSS(i__),RAdeg(i__),DEdeg(i__),objID(i__),zph(i__),
+ zph_L(i__),zph_U(i__),zpprob(i__),umag(i__),gmag(i__),
+ rmag(i__),imag(i__),zmag(i__),e_umag(i__),e_gmag(i__),
+ e_rmag(i__),e_imag(i__),e_zmag(i__),E_B_V(i__),c(i__),
+ F20cm(i__),FX(i__),pm(i__),mov(i__),fz(i__),lz(i__),mz(i__),
+ hz(i__),uv(i__),log_q(i__),log_s(i__),good(i__),Class(i__),
+ z(i__)
c ..............Just test output...........
write(6,1)
+ Seq(i__),SDSS(i__),RAdeg(i__),DEdeg(i__),objID(i__),zph(i__),
+ zph_L(i__),zph_U(i__),zpprob(i__),umag(i__),gmag(i__),
+ rmag(i__),imag(i__),zmag(i__),e_umag(i__),e_gmag(i__),
+ e_rmag(i__),e_imag(i__),e_zmag(i__),E_B_V(i__),c(i__),
+ F20cm(i__),FX(i__),pm(i__),mov(i__),fz(i__),lz(i__),mz(i__),
+ hz(i__),uv(i__),log_q(i__),log_s(i__),good(i__),Class(i__),
+ z(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Rejected quasar candidates
C Format for file interpretation
2 format(
+ I7,1X,A18,1X,F12.7,1X,F11.7,2X,A19,1X,F7.3,1X,F6.3,1X,F6.3,1X,
+ F6.3,1X,F7.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F7.3,1X,F7.3,1X,F8.2,1X,F7.4,1X,
+ F8.2,1X,I2,2X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,F9.3,1X,F8.3,2X,
+ I2,1X,A16,1X,F6.3)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,157075
read(1,'(A271)')ar__1
read(ar__1,2)
+ Seq_1(i__),SDSS_1(i__),RAdeg_1(i__),DEdeg_1(i__),objID_1(i__),
+ zph_1(i__),zph_L_1(i__),zph_U_1(i__),zpprob_1(i__),
+ umag_1(i__),gmag_1(i__),rmag_1(i__),imag_1(i__),zmag_1(i__),
+ e_umag_1(i__),e_gmag_1(i__),e_rmag_1(i__),e_imag_1(i__),
+ e_zmag_1(i__),E_B_V_1(i__),c_1(i__),F20cm_1(i__),FX_1(i__),
+ pm_1(i__),mov_1(i__),fz_1(i__),lz_1(i__),mz_1(i__),hz_1(i__),
+ uv_1(i__),log_q_1(i__),log_s_1(i__),good_1(i__),Class_1(i__),
+ z_1(i__)
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),SDSS_1(i__),RAdeg_1(i__),DEdeg_1(i__),objID_1(i__),
+ zph_1(i__),zph_L_1(i__),zph_U_1(i__),zpprob_1(i__),
+ umag_1(i__),gmag_1(i__),rmag_1(i__),imag_1(i__),zmag_1(i__),
+ e_umag_1(i__),e_gmag_1(i__),e_rmag_1(i__),e_imag_1(i__),
+ e_zmag_1(i__),E_B_V_1(i__),c_1(i__),F20cm_1(i__),FX_1(i__),
+ pm_1(i__),mov_1(i__),fz_1(i__),lz_1(i__),mz_1(i__),hz_1(i__),
+ uv_1(i__),log_q_1(i__),log_s_1(i__),good_1(i__),Class_1(i__),
+ z_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end