Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/696/580 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-Mar-28 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/696/580 Properties of weak emission-line QSOs (WLQ) (Shemmer+, 2009) *================================================================================ *X-ray insights into the nature of weak emission-line quasars at high redshift. * Shemmer O., Brandt W.N., Anderson S.F., Diamond-Stanic A.M., Fan X., * Richards G.T., Schneider D.P., Strauss M.A. * <Astrophys. J., 696, 580-590 (2009)> * =2009ApJ...696..580S C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Chandra and XMM observation log of WLQs integer*4 nr__ parameter (nr__=11) ! Number of records character*54 ar__ ! Full-size record character*18 SDSS (nr__) ! SDSS identification (HHMMSS.ss+DDMMSS.s) character*1 f_SDSS (nr__) ! [hi] Note about the X-ray data (G1) real*4 z (nr__) ! Redshift real*4 rOX (nr__) ! (arcsec) ? Angular distance between the optical and * X-ray position (2) character*1 f_ObsDate (nr__) ! [d] Merged image (G1) character*10 ObsDate (nr__) ! ("YYYY/MM/DD") Date of observation integer*4 ObsID (nr__) ! ? Chandra observation identification number real*4 ExpT (nr__) ! (ks) Chandra exposure time * (corrected for detector dead time) integer*4 Ref (nr__) ! ? Reference number in which the source were * first identified as WLQ (Weak emission-Line * quasar) or BL Lac candidate (5) *Note (2): The optical positions of the quasars have been obtained from the * reference given in the seventh column, and the X-ray positions have * been obtained with wavdetect. *Note (5): Position reference as follows: * 1 = Collinge et al. 2005, Cat. J/AJ/129/2542 * 2 = this work; * 3 = Schneider et al. 2005, SDSS-DR3 quasar catalog, Cat. VII/243 (catalog * superseded by VII/260); * 4 = Fan et al. 2006AJ....131.1203F; * 5 = Schneider et al. 2007, SDSS-DR5 quasar catalog, Cat. VII/252 (catalog * superseded by VII/260). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! X-ray, optical, and radio properties of WLQs integer*4 nr__1 parameter (nr__1=11) ! Number of records character*131 ar__1 ! Full-size record character*18 SDSS_1 (nr__1) ! SDSS identification (HHMMSS.ss+DDMMSS.s) character*1 f_SDSS_1 (nr__1) ! [hi] Note about the X-ray data (G1) real*4 NH (nr__1) ! (10+20cm-2) Neutral galactic absorption column density * from Dickey & Lockman (1990ARA&A..28..215D) real*4 AB1450 (nr__1) ! (mag) Monochromatic AB magnitude at a rest-frame * wavelength of 1450{AA} character*1 f_AB1450 (nr__1) ! [j] Spectrum from Keck (G1) real*4 BMAG (nr__1) ! (mag) Absolute B-band magnitude real*4 F2500 (nr__1) ! (10uJy) Flux density at rest-frame 2500{AA} * in units of 10^-28^erg/cm^2^/s/Hz (10^-5^Jy) real*4 L2500 (nr__1) ! ([10-7W]) Log monochromatic luminosity at a rest-frame * wavelength of 2500{AA} character*1 l_CRate (nr__1) ! Limit flag on CRate real*4 CRate (nr__1) ! (ct/ks) Observed count rate in the 0.5-2keV band real*4 E_CRate (nr__1) ! (ct/ks) ? Positive error on CRate real*4 e_CRate_1 (nr__1) ! (ct/ks) ? Negative error on CRate character*1 f_CRate (nr__1) ! [k] corrected for exposure map (G1) character*1 l_Fx (nr__1) ! Limit flag on Fx real*4 Fx (nr__1) ! (aW/m2) Galactic absorption-corrected flux in the * 0.5-2keV band (in 10^-15^erg/cm^2^/s) real*4 E_Fx (nr__1) ! (aW/m2) ? Positive error on Fx real*4 e_Fx_1 (nr__1) ! (aW/m2) ? Negative error on Fx character*1 l_F2keV (nr__1) ! Limit flag on F2keV real*4 F2keV (nr__1) ! (nJy) Flux density at rest-frame 2keV in units of * 10^-32^erg/cm^2^/s/Hz (10^-9^Jy) character*1 l_L2keV (nr__1) ! Limit flag on log2keV real*4 L2keV (nr__1) ! ([10-7W]) Log monochromatic luminosity at rest-frame 2keV * (log({nu}L_{nu}_)) character*1 l_logL (nr__1) ! Limit flag on logL real*4 logL (nr__1) ! ([10-7W]) Log luminosity in the rest-frame 2-10keV band character*1 l_aox (nr__1) ! Limit flag on aox real*4 aox (nr__1) ! Optical-to-X-ray power-law slope (4) character*1 l_Daox (nr__1) ! Limit flag on Daox real*4 Daox (nr__1) ! Difference between measured and predicted * {alpha}_ox_ (5) character*1 l_Sigox (nr__1) ! Limit flag on Sigox real*4 Sigox (nr__1) ! Significance of the difference between * measured and predicted {alpha}_ox_ (6) character*1 l_aro (nr__1) ! Limit flag on aro real*4 aro (nr__1) ! Radio-to-optical power-law slope (7) character*1 f_aro (nr__1) ! [l] Radio flux density from NVSS (G1) *Note (4): {alpha}_ox_, defined as: * {alpha}_ox_=log(f_2keV_/f_2500{AA}_)/log({nu}_2keV_/{nu}_2500{AA}_) * where f_2keV_ and f_2500{AA}_ are the flux densities at rest-frame * 2keV and 2500{AA}, respectively. *Note (5): {Delta}{alpha}_ox_ is the difference between measured (aOX) and * predicted {alpha}_ox_, given the UV luminosity of log2500, based on * the observed {alpha}_ox_-L_{nu}_(2500{AA}) relation in AGNs (given as * Equation (3) of Just et al. 2007ApJ...665.1004J) *Note (6): The statistical significance of the difference (DaOX) is also * given in units of {sigma}, where * {sigma}=0.146 for 31<logL_{nu}_(2500{AA})<32, and * {sigma}=0.131 for 32<logL_{nu}_(2500{AA})<33 * (see Table 5 of Steffen et al. 2006AJ....131.2826S; these differences * do not account for flux-density measurement errors). *Note (7): Unless otherwise noted, radio-to-optical flux ratios * ({alpha}_ro_; see Section 3.1 for the conversion of these values to * the radio-loudness parameter R) involve radio flux densities at an * observed-frame frequency of 1.4GHz taken from the FIRST survey * (Becker et al. 1995ApJ...450..559B); upper limits on {alpha}_ro_ are * calculated from the 3{sigma} FIRST detection threshold at the source * position. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! X-ray, optical, and radio properties of the BL Lac sample integer*4 nr__2 parameter (nr__2=279) ! Number of records character*92 ar__2 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*24 Name (nr__2) ! Source name real*8 RAdeg (nr__2) ! (deg) Right Ascension in decimal degrees (J2000) real*8 DEdeg (nr__2) ! (deg) Declination in decimal degrees (J2000) real*4 z_1 (nr__2) ! Redshift real*4 log2keV (nr__2) ! ([10-7W]) Log monochromatic luminosity at 2keV in erg/s integer*4 f_log2keV (nr__2) ! [0/1] 0 = upper limit in log2keV (1) real*4 log5500 (nr__2) ! ([10-7W]) Log monochromatic luminosity at 5500{AA} integer*4 f_log5500 (nr__2) ! [0/1] 0 = upper limit on log5500 (1) real*4 log5GHz (nr__2) ! ([10-7W]) Log monochromatic luminosity at 5GHz in erg/s integer*4 f_log5GHz (nr__2) ! [0/1] 0 = upper limit on log5GHz (1) real*4 aox_1 (nr__2) ! Optical/X-ray spectral index integer*4 f_aox (nr__2) ! [0/1] 0 = upper limit on aox (1) real*4 aro_1 (nr__2) ! Radio/Optical spectral index integer*4 f_aro_1 (nr__2) ! [0/1] 0 = upper limit on aro (1) integer*4 Ref_1 (nr__2) ! Reference for multiwavelength data (2) *Note (1): 0 = upper limit, 1 = measured value. *Note (2): Reference as follows: * 1 = Plotkin et al. (2008, Cat. J/AJ/135/2453); * 2 = Laurent-Muehleisen et al. (1999, Cat. J/ApJ/525/127); * 3 = Collinge et al. (2005, Cat. J/AJ/129/2542); * 4 = Perlman et al. (1998, Cat. J/AJ/115/1253); * 5 = Rector et al. (2000AJ....120.1626R); * 6 = Perlman et al. (1996ApJS..104..251P); * 7 = Landt et al. (2001, Cat. J/MNRAS/323/757). C============================================================================= C Loading file 'table1.dat' ! Chandra and XMM observation log of WLQs C Format for file interpretation 1 format( + A18,1X,A1,1X,F3.1,1X,F3.1,1X,A1,1X,A10,1X,I4,1X,F5.2,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,11 read(1,'(A54)')ar__ read(ar__,1) + SDSS(i__),f_SDSS(i__),z(i__),rOX(i__),f_ObsDate(i__), + ObsDate(i__),ObsID(i__),ExpT(i__),Ref(i__) if(ar__(26:28) .EQ. '') rOX(i__) = rNULL__ if(ar__(43:46) .EQ. '') ObsID(i__) = iNULL__ if(ar__(54:54) .EQ. '') Ref(i__) = iNULL__ c ..............Just test output........... write(6,1) + SDSS(i__),f_SDSS(i__),z(i__),rOX(i__),f_ObsDate(i__), + ObsDate(i__),ObsID(i__),ExpT(i__),Ref(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! X-ray, optical, and radio properties of WLQs C Format for file interpretation 2 format( + A18,1X,A1,1X,F4.2,1X,F4.1,1X,A1,1X,F5.1,1X,F4.1,1X,F4.1,1X,A1, + F5.2,1X,F4.2,1X,F4.2,1X,A1,1X,A1,F4.1,1X,F4.1,1X,F4.1,1X,A1, + F5.2,1X,A1,F4.1,1X,A1,F4.1,1X,A1,F5.2,1X,A1,F5.2,1X,A1,F3.1, + 2X,A1,F5.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,11 read(1,'(A131)')ar__1 read(ar__1,2) + SDSS_1(i__),f_SDSS_1(i__),NH(i__),AB1450(i__),f_AB1450(i__), + BMAG(i__),F2500(i__),L2500(i__),l_CRate(i__),CRate(i__), + E_CRate(i__),e_CRate_1(i__),f_CRate(i__),l_Fx(i__),Fx(i__), + E_Fx(i__),e_Fx_1(i__),l_F2keV(i__),F2keV(i__),l_L2keV(i__), + L2keV(i__),l_logL(i__),logL(i__),l_aox(i__),aox(i__), + l_Daox(i__),Daox(i__),l_Sigox(i__),Sigox(i__),l_aro(i__), + aro(i__),f_aro(i__) if(ar__1(57:60) .EQ. '') E_CRate(i__) = rNULL__ if(ar__1(62:65) .EQ. '') e_CRate_1(i__) = rNULL__ if(ar__1(75:78) .EQ. '') E_Fx(i__) = rNULL__ if(ar__1(80:83) .EQ. '') e_Fx_1(i__) = rNULL__ c ..............Just test output........... write(6,2) + SDSS_1(i__),f_SDSS_1(i__),NH(i__),AB1450(i__),f_AB1450(i__), + BMAG(i__),F2500(i__),L2500(i__),l_CRate(i__),CRate(i__), + E_CRate(i__),e_CRate_1(i__),f_CRate(i__),l_Fx(i__),Fx(i__), + E_Fx(i__),e_Fx_1(i__),l_F2keV(i__),F2keV(i__),l_L2keV(i__), + L2keV(i__),l_logL(i__),logL(i__),l_aox(i__),aox(i__), + l_Daox(i__),Daox(i__),l_Sigox(i__),Sigox(i__),l_aro(i__), + aro(i__),f_aro(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! X-ray, optical, and radio properties of the * BL Lac sample C Format for file interpretation 3 format( + A24,1X,F9.5,1X,F9.5,1X,F5.3,1X,F5.2,1X,I1,1X,F5.2,1X,I1,1X, + F5.2,1X,I1,1X,F5.2,1X,I1,1X,F5.2,1X,I1,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,279 read(1,'(A92)')ar__2 read(ar__2,3) + Name(i__),RAdeg(i__),DEdeg(i__),z_1(i__),log2keV(i__), + f_log2keV(i__),log5500(i__),f_log5500(i__),log5GHz(i__), + f_log5GHz(i__),aox_1(i__),f_aox(i__),aro_1(i__),f_aro_1(i__), + Ref_1(i__) c ..............Just test output........... write(6,3) + Name(i__),RAdeg(i__),DEdeg(i__),z_1(i__),log2keV(i__), + f_log2keV(i__),log5500(i__),f_log5500(i__),log5GHz(i__), + f_log5GHz(i__),aox_1(i__),f_aox(i__),aro_1(i__),f_aro_1(i__), + Ref_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end