Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/133/313 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/AJ/133/313 AGN from RASS and SDSS DR5 (Anderson+, 2007) *================================================================================ *A large, uniform sample of x-ray-emitting active galactic nuclei from the *ROSAT All Sky and Sloan Digital Sky Surveys: the Data Release 5 sample. * Anderson S.F., Margon B., Voges W., Plotkin R.M., Syphers D., Haggard D., * Collinge M.J., Meyer J., Strauss M.A., Agueros M.A., Hall P.B., Homer L., * Ivezic Z., Richards G.T., Richmond M.W., Schneider D.P., Stinson G., * Vanden Berk D.E., York D.G. * <Astron. J., 133, 313-329 (2007)> * =2007AJ....133..313A C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Observed Parameters of Broad-Line RASS/SDSS AGN integer*4 nr__ parameter (nr__=6224) ! Number of records character*108 ar__ ! Full-size record character*15 RASS (nr__) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS (nr__) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 umag (nr__) ! (mag) SDSS u band PSF magnitude real*4 gmag (nr__) ! (mag) SDSS g band PSF magnitude real*4 rmag (nr__) ! (mag) SDSS r band PSF magnitude real*4 imag (nr__) ! (mag) SDSS i band PSF magnitude real*4 zmag (nr__) ! (mag) SDSS z band PSF magnitude integer*4 St (nr__) ! [3/6] SDSS imaging morphology parameter * (6=stellar/unresolved, 3=extended/resolved) real*4 z (nr__) ! ? Redshift real*4 CRate (nr__) ! (ct/s) X-ray count rate (3) integer*4 Exp (nr__) ! (s) The X-ray exposure time real*4 HR1 (nr__) ! The X-ray hardness ratio 1 (4) real*4 HR2 (nr__) ! The X-ray hardness ratio 2 (4) integer*4 L (nr__) ! The X-ray source detection likelihood real*4 FX (nr__) ! (10-16W/m2) Absorbed X-ray flux (5) *Note (3): In the 0.1-2.4keV broadband, corrected for vignetting. *Note (4): The usual ROSAT PSPC HR1 hardness ratio defined as * HR1=(B-A)/(B+A), where * A = count rate in pulse height amplitude (PHA) range 11-41, and * B = count rate in PHA range 52-201. * The usual ROSAT PSPC HR2 hardness ratio defined as HR2=(D-C)/(D+C), where * C = count rate in PHA range 52-90, and * D = count rate in PHA range 91-201. *Note (5): In the 0.1-2.4keV band and in units of 10^-13^erg/s/cm^2^. * Assuming a power-law X-ray spectrum with energy index {alpha}_x_=1.5, * typical of low redshift quasars in the RASS PSPC bandpass (e.g., see * Schartel et al. 1996MNRAS.283.1015S). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Derived Parameters of Broad-Line RASS/SDSS AGN integer*4 nr__1 parameter (nr__1=6224) ! Number of records character*116 ar__1 ! Full-size record character*15 RASS_1 (nr__1) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS_1 (nr__1) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 g0mag (nr__1) ! (mag) SDSS g band PSF magnitude (2) real*4 z_1 (nr__1) ! ? Redshift real*4 FXcor (nr__1) ! (10-16W/m2) Corrected X-ray flux (3) real*4 logLX (nr__1) ! ([10-7W]) ? Log of broadband (0.1-2.4kev) X-ray * luminosity in units of erg/s (4) real*4 LogL2500 (nr__1) ! ([10-7W/Hz]) ? Log of 2500{AA} monochromatic optical * luminosity in units of erg/s/Hz (4) real*4 logL2keV (nr__1) ! ([10-7W/Hz]) ? Log of 2keV monochromatic X-ray * luminosity in units of erg/s/Hz (4) real*4 alpha (nr__1) ! Slope of a hypothetical power law * from 2500{AA} to 2keV character*39 Comm (nr__1) ! Additional comment (5) *Note (2): Extinction corrected according to the reddening maps of Schlegel, * Finkbeiner, & Davis (1998ApJ...500..525S). *Note (3): Corrected for absorption within the Galaxy, with X-ray absorbing * columns estimated using the N_H_ column density measures of the Stark * et al. (1992) 21cm maps. In units of 10^-13^erg/s/cm^2^. *Note (4): We adopt values of H_0_=70 km/s/Mpc, {Omega}_M_=0.3, and * {Omega}_{Lambda}_=0.7. In converting from corrected broadband (X-ray * 0.1-2.4kev, and optical g-band) fluxes to luminosities, we again * assume an X-ray power-law spectrum with energy index {alpha}_x_=1.5, * and an optical power-law with energy index {alpha}_o_=0.5. *Note (5): For example noting selected objects that are radio sources, * listing an alternate name (shortened, if needed, to the first 20 * characters) taken from NED for about 10% of the quasars, and noting * the ~1.5% of the cases where two of our cataloged AGN fall within the * same RASS error circle. (The latter are denoted by `ambigID' in the * comments). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Observed Parameters of RASS/SDSS AGN Having Narrower Permitted Emission integer*4 nr__2 parameter (nr__2=515) ! Number of records character*108 ar__2 ! Full-size record character*15 RASS_2 (nr__2) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS_2 (nr__2) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 umag_1 (nr__2) ! (mag) SDSS u band PSF magnitude real*4 gmag_1 (nr__2) ! (mag) SDSS g band PSF magnitude real*4 rmag_1 (nr__2) ! (mag) SDSS r band PSF magnitude real*4 imag_1 (nr__2) ! (mag) SDSS i band PSF magnitude real*4 zmag_1 (nr__2) ! (mag) SDSS z band PSF magnitude integer*4 St_1 (nr__2) ! [3/6] SDSS imaging morphology parameter * (6=stellar/unresolved, 3=extended/resolved) real*4 z_2 (nr__2) ! ? Redshift real*4 CRate_1 (nr__2) ! (ct/s) X-ray count rate (3) integer*4 Exp_1 (nr__2) ! (s) The X-ray exposure time real*4 HR1_1 (nr__2) ! The X-ray hardness ratio 1 (4) real*4 HR2_1 (nr__2) ! The X-ray hardness ratio 2 (4) integer*4 L_1 (nr__2) ! The X-ray source detection likelihood real*4 FX_1 (nr__2) ! (10-16W/m2) Absorbed X-ray flux (5) *Note (3): In the 0.1-2.4keV broadband, corrected for vignetting. *Note (4): The usual ROSAT PSPC HR1 hardness ratio defined as * HR1=(B-A)/(B+A), where * A = count rate in pulse height amplitude (PHA) range 11-41, and * B = count rate in PHA range 52-201. * The usual ROSAT PSPC HR2 hardness ratio defined as HR2=(D-C)/(D+C), where * C = count rate in PHA range 52-90, and * D = count rate in PHA range 91-201. *Note (5): In the 0.1-2.4keV band and in units of 10^-13^erg/s/cm^2^. * Assuming a power-law X-ray spectrum with energy index {alpha}_x_=1.5, * typical of low redshift quasars in the RASS PSPC bandpass (e.g., see * Schartel et al. 1996MNRAS.283.1015S). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Derived Parameters of RASS/SDSS AGN Having Narrower Permitted Emission integer*4 nr__3 parameter (nr__3=515) ! Number of records character*116 ar__3 ! Full-size record character*15 RASS_3 (nr__3) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS_3 (nr__3) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 g0mag_1 (nr__3) ! (mag) SDSS g band PSF magnitude (2) real*4 z_3 (nr__3) ! ? Redshift real*4 FXcor_1 (nr__3) ! (10-16W/m2) Corrected X-ray flux (3) real*4 logLX_1 (nr__3) ! ([10-7W]) ? Log of broadband (0.1-2.4kev) X-ray * luminosity in units of erg/s (4) real*4 LogL2500_1 (nr__3) ! ([10-7W/Hz]) ? Log of 2500{AA} monochromatic optical * luminosity in units of erg/s/Hz (4) real*4 logL2keV_1 (nr__3) ! ([10-7W/Hz]) ? Log of 2keV monochromatic X-ray * luminosity in units of erg/s/Hz (4) real*4 alpha_1 (nr__3) ! Slope of a hypothetical power law * from 2500{AA} to 2keV character*39 Comm_1 (nr__3) ! Additional comment (5) *Note (2): Extinction corrected according to the reddening maps of Schlegel, * Finkbeiner, & Davis (1998ApJ...500..525S). *Note (3): Corrected for absorption within the Galaxy, with X-ray absorbing * columns estimated using the N_H_ column density measures of the Stark * et al. (1992) 21cm maps. In units of 10^-13^erg/s/cm^2^. *Note (4): We adopt values of H_0_=70 km/s/Mpc, {Omega}_M_=0.3, and * {Omega}_{Lambda}_=0.7. In converting from corrected broadband (X-ray * 0.1-2.4kev, and optical g-band) fluxes to luminosities, we again * assume an X-ray power-law spectrum with energy index {alpha}_x_=1.5, * and an optical power-law with energy index {alpha}_o_=0.5. *Note (5): For example noting selected objects that are radio sources, * listing an alternate name (shortened, if needed, to the first 20 * characters) taken from NED for about 10% of the quasars, and noting * the ~1.5% of the cases where two of our cataloged AGN fall within the * same RASS error circle. (The latter are denoted by `ambigID' in the * comments). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Observed Parameters of RASS/SDSS BL Lac Candidates integer*4 nr__4 parameter (nr__4=266) ! Number of records character*108 ar__4 ! Full-size record character*15 RASS_4 (nr__4) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS_4 (nr__4) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 umag_2 (nr__4) ! (mag) SDSS u band PSF magnitude real*4 gmag_2 (nr__4) ! (mag) SDSS g band PSF magnitude real*4 rmag_2 (nr__4) ! (mag) SDSS r band PSF magnitude real*4 imag_2 (nr__4) ! (mag) SDSS i band PSF magnitude real*4 zmag_2 (nr__4) ! (mag) SDSS z band PSF magnitude integer*4 St_2 (nr__4) ! [3/6] SDSS imaging morphology parameter * (6=stellar/unresolved, 3=extended/resolved) real*4 z_4 (nr__4) ! ? Redshift real*4 CRate_2 (nr__4) ! (ct/s) X-ray count rate (3) integer*4 Exp_2 (nr__4) ! (s) The X-ray exposure time real*4 HR1_2 (nr__4) ! The X-ray hardness ratio 1 (4) real*4 HR2_2 (nr__4) ! The X-ray hardness ratio 2 (4) integer*4 L_2 (nr__4) ! The X-ray source detection likelihood real*4 FX_2 (nr__4) ! (10-16W/m2) Absorbed X-ray flux (5) *Note (3): In the 0.1-2.4keV broadband, corrected for vignetting. *Note (4): The usual ROSAT PSPC HR1 hardness ratio defined as * HR1=(B-A)/(B+A), where * A = count rate in pulse height amplitude (PHA) range 11-41, and * B = count rate in PHA range 52-201. * The usual ROSAT PSPC HR2 hardness ratio defined as HR2=(D-C)/(D+C), where * C = count rate in PHA range 52-90, and * D = count rate in PHA range 91-201. *Note (5): In the 0.1-2.4keV band and in units of 10^-13^erg/s/cm^2^. * Assuming a power-law X-ray spectrum with energy index {alpha}_x_=1.5, * typical of low redshift quasars in the RASS PSPC bandpass (e.g., see * Schartel et al. 1996MNRAS.283.1015S). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Derived Parameters of RASS/SDSS BL Lac Candidates integer*4 nr__5 parameter (nr__5=266) ! Number of records character*116 ar__5 ! Full-size record character*15 RASS_5 (nr__5) ! RASS X-ray source name, HHMMSS.s+DDMMSS (G1) character*18 SDSS_5 (nr__5) ! SDSS optical counterpart name, * HHMMSS.ss+DDMMSS.s (G1) real*4 g0mag_2 (nr__5) ! (mag) SDSS g band PSF magnitude (2) real*4 z_5 (nr__5) ! ? Redshift real*4 FXcor_2 (nr__5) ! (10-16W/m2) Corrected X-ray flux (3) real*4 logLX_2 (nr__5) ! ([10-7W]) ? Log of broadband (0.1-2.4kev) X-ray * luminosity in units of erg/s (4) real*4 LogL2500_2 (nr__5) ! ([10-7W/Hz]) ? Log of 2500{AA} monochromatic optical * luminosity in units of erg/s/Hz (4) real*4 logL2keV_2 (nr__5) ! ([10-7W/Hz]) ? Log of 2keV monochromatic X-ray * luminosity in units of erg/s/Hz (4) real*4 alpha_2 (nr__5) ! Slope of a hypothetical power law * from 2500{AA} to 2keV character*39 Comm_2 (nr__5) ! Additional comment (5) *Note (2): Extinction corrected according to the reddening maps of Schlegel, * Finkbeiner, & Davis (1998ApJ...500..525S). *Note (3): Corrected for absorption within the Galaxy, with X-ray absorbing * columns estimated using the N_H_ column density measures of the Stark * et al. (1992) 21cm maps. In units of 10^-13^erg/s/cm^2^. *Note (4): We adopt values of H_0_=70 km/s/Mpc, {Omega}_M_=0.3, and * {Omega}_{Lambda}_=0.7. In converting from corrected broadband (X-ray * 0.1-2.4kev, and optical g-band) fluxes to luminosities, we again * assume an X-ray power-law spectrum with energy index {alpha}_x_=1.5, * and an optical power-law with energy index {alpha}_o_=0.5. *Note (5): For example noting selected objects that are radio sources, * listing an alternate name (shortened, if needed, to the first 20 * characters) taken from NED for about 10% of the quasars, and noting * the ~1.5% of the cases where two of our cataloged AGN fall within the * same RASS error circle. (The latter are denoted by `ambigID' in the * comments). C============================================================================= C Loading file 'table1.dat' ! Observed Parameters of Broad-Line RASS/SDSS AGN C Format for file interpretation 1 format( + A15,1X,A18,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I1,1X, + F5.3,1X,F6.4,1X,I4,1X,F5.2,1X,F5.2,1X,I4,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,6224 read(1,'(A108)')ar__ read(ar__,1) + RASS(i__),SDSS(i__),umag(i__),gmag(i__),rmag(i__),imag(i__), + zmag(i__),St(i__),z(i__),CRate(i__),Exp(i__),HR1(i__), + HR2(i__),L(i__),FX(i__) if(ar__(68:72) .EQ. '') z(i__) = rNULL__ c ..............Just test output........... write(6,1) + RASS(i__),SDSS(i__),umag(i__),gmag(i__),rmag(i__),imag(i__), + zmag(i__),St(i__),z(i__),CRate(i__),Exp(i__),HR1(i__), + HR2(i__),L(i__),FX(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Derived Parameters of Broad-Line RASS/SDSS AGN C Format for file interpretation 2 format( + A15,1X,A18,1X,F5.2,1X,F5.3,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X, + F4.2,1X,A39) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,6224 read(1,'(A116)')ar__1 read(ar__1,2) + RASS_1(i__),SDSS_1(i__),g0mag(i__),z_1(i__),FXcor(i__), + logLX(i__),LogL2500(i__),logL2keV(i__),alpha(i__),Comm(i__) if(ar__1(42:46) .EQ. '') z_1(i__) = rNULL__ if(ar__1(55:59) .EQ. '') logLX(i__) = rNULL__ if(ar__1(61:65) .EQ. '') LogL2500(i__) = rNULL__ if(ar__1(67:71) .EQ. '') logL2keV(i__) = rNULL__ c ..............Just test output........... write(6,2) + RASS_1(i__),SDSS_1(i__),g0mag(i__),z_1(i__),FXcor(i__), + logLX(i__),LogL2500(i__),logL2keV(i__),alpha(i__),Comm(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Observed Parameters of RASS/SDSS AGN Having * Narrower Permitted Emission C Format for file interpretation 3 format( + A15,1X,A18,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I1,1X, + F5.3,1X,F6.4,1X,I4,1X,F5.2,1X,F5.2,1X,I4,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,515 read(1,'(A108)')ar__2 read(ar__2,3) + RASS_2(i__),SDSS_2(i__),umag_1(i__),gmag_1(i__),rmag_1(i__), + imag_1(i__),zmag_1(i__),St_1(i__),z_2(i__),CRate_1(i__), + Exp_1(i__),HR1_1(i__),HR2_1(i__),L_1(i__),FX_1(i__) if(ar__2(68:72) .EQ. '') z_2(i__) = rNULL__ c ..............Just test output........... write(6,3) + RASS_2(i__),SDSS_2(i__),umag_1(i__),gmag_1(i__),rmag_1(i__), + imag_1(i__),zmag_1(i__),St_1(i__),z_2(i__),CRate_1(i__), + Exp_1(i__),HR1_1(i__),HR2_1(i__),L_1(i__),FX_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Derived Parameters of RASS/SDSS AGN Having * Narrower Permitted Emission C Format for file interpretation 4 format( + A15,1X,A18,1X,F5.2,1X,F5.3,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X, + F4.2,1X,A39) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,515 read(1,'(A116)')ar__3 read(ar__3,4) + RASS_3(i__),SDSS_3(i__),g0mag_1(i__),z_3(i__),FXcor_1(i__), + logLX_1(i__),LogL2500_1(i__),logL2keV_1(i__),alpha_1(i__), + Comm_1(i__) if(ar__3(42:46) .EQ. '') z_3(i__) = rNULL__ if(ar__3(55:59) .EQ. '') logLX_1(i__) = rNULL__ if(ar__3(61:65) .EQ. '') LogL2500_1(i__) = rNULL__ if(ar__3(67:71) .EQ. '') logL2keV_1(i__) = rNULL__ c ..............Just test output........... write(6,4) + RASS_3(i__),SDSS_3(i__),g0mag_1(i__),z_3(i__),FXcor_1(i__), + logLX_1(i__),LogL2500_1(i__),logL2keV_1(i__),alpha_1(i__), + Comm_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Observed Parameters of RASS/SDSS BL Lac Candidates C Format for file interpretation 5 format( + A15,1X,A18,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I1,1X, + F5.3,1X,F6.4,1X,I4,1X,F5.2,1X,F5.2,1X,I4,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,266 read(1,'(A108)')ar__4 read(ar__4,5) + RASS_4(i__),SDSS_4(i__),umag_2(i__),gmag_2(i__),rmag_2(i__), + imag_2(i__),zmag_2(i__),St_2(i__),z_4(i__),CRate_2(i__), + Exp_2(i__),HR1_2(i__),HR2_2(i__),L_2(i__),FX_2(i__) if(ar__4(68:72) .EQ. '') z_4(i__) = rNULL__ c ..............Just test output........... write(6,5) + RASS_4(i__),SDSS_4(i__),umag_2(i__),gmag_2(i__),rmag_2(i__), + imag_2(i__),zmag_2(i__),St_2(i__),z_4(i__),CRate_2(i__), + Exp_2(i__),HR1_2(i__),HR2_2(i__),L_2(i__),FX_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Derived Parameters of RASS/SDSS BL Lac Candidates C Format for file interpretation 6 format( + A15,1X,A18,1X,F5.2,1X,F5.3,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X, + F4.2,1X,A39) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,266 read(1,'(A116)')ar__5 read(ar__5,6) + RASS_5(i__),SDSS_5(i__),g0mag_2(i__),z_5(i__),FXcor_2(i__), + logLX_2(i__),LogL2500_2(i__),logL2keV_2(i__),alpha_2(i__), + Comm_2(i__) if(ar__5(42:46) .EQ. '') z_5(i__) = rNULL__ if(ar__5(55:59) .EQ. '') logLX_2(i__) = rNULL__ if(ar__5(61:65) .EQ. '') LogL2500_2(i__) = rNULL__ if(ar__5(67:71) .EQ. '') logL2keV_2(i__) = rNULL__ c ..............Just test output........... write(6,6) + RASS_5(i__),SDSS_5(i__),g0mag_2(i__),z_5(i__),FXcor_2(i__), + logLX_2(i__),LogL2500_2(i__),logL2keV_2(i__),alpha_2(i__), + Comm_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end