FORTRAN Generation
(/./ftp/cats/J/ApJ/696/580)

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