FORTRAN Generation
(/./ftp/cats/J/ApJ/692/758)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/692/758 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-29
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/692/758               BAL QSOs in SDSS-DR5                 (Gibson+, 2009)
*================================================================================
*A catalog of broad absorption line quasars in Sloan Digital Sky Survey Data
*Release 5.
*    Gibson R.R., Jiang L., Brandt W.N., Hall P.B., Shen Y., Wu J.,
*    Anderson S.F., Schneider D.P., Vanden Berk D., Gallagher S.C., Fan X.,
*    York D.G.
*   <Astrophys. J., 692, 758-777 (2009)>
*   =2009ApJ...692..758G
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! DR5 BAL Catalog

      integer*4 nr__
      parameter (nr__=5035)	! Number of records
      character*265 ar__   	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*18  SDSS       (nr__) ! SDSS-DR5 identification (HHMMSS.ss+DDMMSS.s)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000)
      real*4        z          (nr__) ! SDSS-DR5 redshift (1)
      integer*4     MJD        (nr__) ! (d) ? Modified Julian date of observation (only
*                                   for table 2)
      integer*4     Plate      (nr__) ! ? SDSS plate number (only for table 2)
      integer*4     Fiber      (nr__) ! ? SDSS fiber number (only for table 2)
      real*8        BI_Si      (nr__) ! (km/s) ? SiIV balnicity index
      real*8        BIO_Si     (nr__) ! (km/s) ? SiIV modified balnicity index; see text
      real*4        EW_Si      (nr__) ! (0.1nm) ? SiIV rest-frame absorption equivalent width
      integer*4     vmin_Si    (nr__) ! (km/s) ? SiIV minimum velocity
      integer*4     vmax_Si    (nr__) ! (km/s) ? SiIV maximum velocity
      real*4        fd_Si      (nr__) ! ? SiIV BAL bin fraction below 50% of continuum
      real*8        BI_C       (nr__) ! (km/s) ? CIV balnicity index
      real*8        BIO_C      (nr__) ! (km/s) ? CIV modified balnicity index; see text
      real*4        EW_C       (nr__) ! (0.1nm) ? CIV rest-frame absorption equivalent width
      integer*4     vmin_C     (nr__) ! (km/s) ? CIV minimum velocity
      integer*4     vmax_C     (nr__) ! (km/s) ? CIV maximum velocity
      real*4        fd_C       (nr__) ! ? CIV BAL bin fraction below 50% of continuum
      real*8        BI_Al      (nr__) ! (km/s) ? AlIII balnicity index
      real*8        BIO_Al     (nr__) ! (km/s) ? AlIII modified balnicity index; see text
      real*4        EW_Al      (nr__) ! (0.1nm) ? AlIII rest-frame absorption equivalent width
      integer*4     vmin_Al    (nr__) ! (km/s) ? AlIII minimum velocity
      integer*4     vmax_Al    (nr__) ! (km/s) ? AlIII maximum velocity
      real*4        fd_Al      (nr__) ! ? AlIII BAL bin fraction below 50% of continuum
      real*8        BI_Mg      (nr__) ! (km/s) ? MgII balnicity index
      real*8        BIO_Mg     (nr__) ! (km/s) ? MgII modified balnicity index; see text
      real*4        EW_Mg      (nr__) ! (0.1nm) ? MgII rest-frame absorption equivalent width
      integer*4     vmin_Mg    (nr__) ! (km/s) ? MgII minimum velocity
      integer*4     vmax_Mg    (nr__) ! (km/s) ? MgII maximum velocity
      real*4        fd_Mg      (nr__) ! ? MgII BAL bin fraction below 50% of continuum
      integer*4     f_SiIVEmL  (nr__) ! [0/1] SiIV EmLost flag (2)
      integer*4     f_CIVEmL   (nr__) ! [0/1] CIV EmLost flag (2)
      integer*4     f_AlIIIEmL (nr__) ! [0/1] AlIII EmLost flag (2)
      integer*4     f_MgIIEmL  (nr__) ! [0/1] MgII EmLost flag (2)
      integer*4     f_SiIVBB   (nr__) ! [0/1]? SiIV BALManyBadBins flag (2)
      integer*4     f_CIVBB    (nr__) ! [0/1]? CIV BALManyBadBins flag (2)
      integer*4     f_AlIIIBB  (nr__) ! [0/1]? AlIII BALManyBadBins flag (2)
      integer*4     f_MgIIBB   (nr__) ! [0/1]? MgII BALManyBadBins flag (2)
      integer*4     f_CIV4BWA  (nr__) ! [0/1] CIV BlueWingAbs flag (2)
      real*4        SN1700     (nr__) ! ? Signal-to-noise (3)
      real*4        logF1400   (nr__) ! ([mW/m2/Hz]) ? Log of 1400 monochromatic luminosity (4)
      real*4        logF2500   (nr__) ! ([mW/m2/Hz]) ? Log of 2500 monochromatic luminosity (4)
*Note (1): The redshifts for our QSOs are taken from the DR5 QSO catalog
*     (Cat. VII/252, superseded by VII/260), with  the exception of two
*     sources, J100424.88+122922.2 and J153029.05+553247.9, which we
*     consider to be misclassified. We assign these sources redshifts of
*     4.66 and 1.73, respectively.
*Note (2): There are cases, such as strongly absorbed emission-line
*     profiles, for which we were unable to accurately determine the
*     underlying emission. We define a set of flags that we use to indicate
*     when the fit model is particularly uncertain due to such effects.
*  EmLost = The profile of the emission line at (nearly) zero velocity has been
*     strongly absorbed and may not have been accurately reconstructed.
*  BlueWingAbs = Additional "blue wing" emission occurs on the blue side
*     of the C emission line, and absorption is apparent between the blue
*     wing and the line core. Since we fit only the line core and cannot
*     accurately reconstruct the "blue wing" emission, we could potentially
*     miss BAL absorption in these cases.
*  BALManyBadBins = Some putative BALs consist of a large number of
*     spectral bins that have been flagged as "bad" in the SDSS pipeline.
*  See text for further details.
*Note (3): Median of the flux divided by the noise for all spectral bins in
*     the 1650-1750{AA} range. See text for further details.
*Note (4): Only for table 1. Corrected for Galactic reddening.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! DR5 BAL catalog duplicate spectra

      integer*4 nr__1
      parameter (nr__1=569)	! Number of records
      character*265 ar__1  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*18  SDSS_1     (nr__1) ! SDSS-DR5 identification (HHMMSS.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)
      real*4        z_1        (nr__1) ! SDSS-DR5 redshift (1)
      integer*4     MJD_1      (nr__1) ! (d) ? Modified Julian date of observation (only
*                                   for table 2)
      integer*4     Plate_1    (nr__1) ! ? SDSS plate number (only for table 2)
      integer*4     Fiber_1    (nr__1) ! ? SDSS fiber number (only for table 2)
      real*8        BI_Si_1    (nr__1) ! (km/s) ? SiIV balnicity index
      real*8        BIO_Si_1   (nr__1) ! (km/s) ? SiIV modified balnicity index; see text
      real*4        EW_Si_1    (nr__1) ! (0.1nm) ? SiIV rest-frame absorption equivalent width
      integer*4     vmin_Si_1  (nr__1) ! (km/s) ? SiIV minimum velocity
      integer*4     vmax_Si_1  (nr__1) ! (km/s) ? SiIV maximum velocity
      real*4        fd_Si_1    (nr__1) ! ? SiIV BAL bin fraction below 50% of continuum
      real*8        BI_C_1     (nr__1) ! (km/s) ? CIV balnicity index
      real*8        BIO_C_1    (nr__1) ! (km/s) ? CIV modified balnicity index; see text
      real*4        EW_C_1     (nr__1) ! (0.1nm) ? CIV rest-frame absorption equivalent width
      integer*4     vmin_C_1   (nr__1) ! (km/s) ? CIV minimum velocity
      integer*4     vmax_C_1   (nr__1) ! (km/s) ? CIV maximum velocity
      real*4        fd_C_1     (nr__1) ! ? CIV BAL bin fraction below 50% of continuum
      real*8        BI_Al_1    (nr__1) ! (km/s) ? AlIII balnicity index
      real*8        BIO_Al_1   (nr__1) ! (km/s) ? AlIII modified balnicity index; see text
      real*4        EW_Al_1    (nr__1) ! (0.1nm) ? AlIII rest-frame absorption equivalent width
      integer*4     vmin_Al_1  (nr__1) ! (km/s) ? AlIII minimum velocity
      integer*4     vmax_Al_1  (nr__1) ! (km/s) ? AlIII maximum velocity
      real*4        fd_Al_1    (nr__1) ! ? AlIII BAL bin fraction below 50% of continuum
      real*8        BI_Mg_1    (nr__1) ! (km/s) ? MgII balnicity index
      real*8        BIO_Mg_1   (nr__1) ! (km/s) ? MgII modified balnicity index; see text
      real*4        EW_Mg_1    (nr__1) ! (0.1nm) ? MgII rest-frame absorption equivalent width
      integer*4     vmin_Mg_1  (nr__1) ! (km/s) ? MgII minimum velocity
      integer*4     vmax_Mg_1  (nr__1) ! (km/s) ? MgII maximum velocity
      real*4        fd_Mg_1    (nr__1) ! ? MgII BAL bin fraction below 50% of continuum
      integer*4     f_SiIVEmL_1(nr__1) ! [0/1] SiIV EmLost flag (2)
      integer*4     f_CIVEmL_1 (nr__1) ! [0/1] CIV EmLost flag (2)
      integer*4     f_AlIIIEmL_1(nr__1) ! [0/1] AlIII EmLost flag (2)
      integer*4     f_MgIIEmL_1(nr__1) ! [0/1] MgII EmLost flag (2)
      integer*4     f_SiIVBB_1 (nr__1) ! [0/1]? SiIV BALManyBadBins flag (2)
      integer*4     f_CIVBB_1  (nr__1) ! [0/1]? CIV BALManyBadBins flag (2)
      integer*4     f_AlIIIBB_1(nr__1) ! [0/1]? AlIII BALManyBadBins flag (2)
      integer*4     f_MgIIBB_1 (nr__1) ! [0/1]? MgII BALManyBadBins flag (2)
      integer*4     f_CIV4BWA_1(nr__1) ! [0/1] CIV BlueWingAbs flag (2)
      real*4        SN1700_1   (nr__1) ! ? Signal-to-noise (3)
      real*4        logF1400_1 (nr__1) ! ([mW/m2/Hz]) ? Log of 1400 monochromatic luminosity (4)
      real*4        logF2500_1 (nr__1) ! ([mW/m2/Hz]) ? Log of 2500 monochromatic luminosity (4)
*Note (1): The redshifts for our QSOs are taken from the DR5 QSO catalog
*     (Cat. VII/252, superseded by VII/260), with  the exception of two
*     sources, J100424.88+122922.2 and J153029.05+553247.9, which we
*     consider to be misclassified. We assign these sources redshifts of
*     4.66 and 1.73, respectively.
*Note (2): There are cases, such as strongly absorbed emission-line
*     profiles, for which we were unable to accurately determine the
*     underlying emission. We define a set of flags that we use to indicate
*     when the fit model is particularly uncertain due to such effects.
*  EmLost = The profile of the emission line at (nearly) zero velocity has been
*     strongly absorbed and may not have been accurately reconstructed.
*  BlueWingAbs = Additional "blue wing" emission occurs on the blue side
*     of the C emission line, and absorption is apparent between the blue
*     wing and the line core. Since we fit only the line core and cannot
*     accurately reconstruct the "blue wing" emission, we could potentially
*     miss BAL absorption in these cases.
*  BALManyBadBins = Some putative BALs consist of a large number of
*     spectral bins that have been flagged as "bad" in the SDSS pipeline.
*  See text for further details.
*Note (3): Median of the flux divided by the noise for all spectral bins in
*     the 1650-1750{AA} range. See text for further details.
*Note (4): Only for table 1. Corrected for Galactic reddening.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! DR5 non-BALs with CIV BlueWingAbs flag set

      integer*4 nr__2
      parameter (nr__2=118)	! Number of records
      character*51 ar__2  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*18  SDSS_2     (nr__2) ! SDSS-DR5 name (HHMMSS.ss+DDMMSS.s)
      real*8        RAdeg_2    (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination in decimal degrees (J2000)
      real*4        z_2        (nr__2) ! Redshift
      integer*4     SiIV       (nr__2) ! [0/1]? SiIV contamination? (only for table 4)
      integer*4     AlIII      (nr__2) ! [0/1]? AlIII contamination? (only for table 4)
      integer*4     MgII       (nr__2) ! [0/1]? MgII contamination? (only for table 4)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4.dat'	! BALs with absorption from other ions

      integer*4 nr__3
      parameter (nr__3=92)	! Number of records
      character*51 ar__3  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*18  SDSS_3     (nr__3) ! SDSS-DR5 name (HHMMSS.ss+DDMMSS.s)
      real*8        RAdeg_3    (nr__3) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination in decimal degrees (J2000)
      real*4        z_3        (nr__3) ! Redshift
      integer*4     SiIV_1     (nr__3) ! [0/1]? SiIV contamination? (only for table 4)
      integer*4     AlIII_1    (nr__3) ! [0/1]? AlIII contamination? (only for table 4)
      integer*4     MgII_1     (nr__3) ! [0/1]? MgII contamination? (only for table 4)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table5.dat'	! X-ray source information

      integer*4 nr__4
      parameter (nr__4=73)	! Number of records
      character*90 ar__4  	! Full-size record

      character*18  SDSS_4     (nr__4) ! SDSS source (HHMMSS.ss+DDMMSS.s)
      character*10  ObsID      (nr__4) ! Chandra or XMM-Newton observation
*                                   identification number (1)
      real*4        Exp        (nr__4) ! (ks) X-ray exposure time (2)
      character*7   Sct        (nr__4) ! (ct) 0.5-2keV counts source/BG
      character*7   Hct        (nr__4) ! (ct) 2-8keV counts source/BG
      integer*4     Det        (nr__4) ! [0/1] X-ray detected? (0 for no detection)
      character*1   l_logR_    (nr__4) ! Limit flag on logR*
      real*4        logR_      (nr__4) ! ([-]) Radio loudness parameter (3)
      character*1   l_logLx    (nr__4) ! Limit flag on log(L2keV)
      real*4        logLx      (nr__4) ! ([10-7W/m2/Hz]) The 2keV monochromatic luminosity in erg/s/Hz
      character*1   l_aox      (nr__4) ! Limit flag on aox
      real*4        aox        (nr__4) ! Optical-to-X-ray spectral index (4)
      real*4        e_aox      (nr__4) ! ? aox uncertainty (5)
      character*1   l_Daox     (nr__4) ! Limit flag on Daox
      real*4        Daox       (nr__4) ! Difference in X-ray luminosity (6)
      real*4        e_Daox     (nr__4) ! ? Delaox uncertainty (5)
*Note (1): Ten-digit numbers correspond to XMM-Newton observations, while
*     shorter numbers correspond to Chandra observations.
*Note (2): The effective X-ray exposure reported by the CIAO or SAS
*     toolchains for the extraction of the source region.
*Note (3): log(R*)=log(L_{nu}_(5GHz)/L_2500_). See section 4.6 for further
*     details.
*Note (4): {alpha}_ox_ is usually expressed in terms of the logarithm of the
*     ratio between the monochromatic luminosities L_{nu}_ at 2keV and
*     2500{AA} with the parameter:
*     {alpha}_ox_=0.3838log(L_{nu}_(2keV)/L_{nu}_(2500{AA})). For typical
*     QSOs, {alpha}_ox_ is a function of L_2500_, the monochromatic
*     luminosity at 2500{AA}. Just et al. (2007ApJ...665.1004J) found the
*     equation (3):
*     {alpha}_ox_(L_2500_)=(-0.140+/-0.007)log(L_2500_)+(2.705+/-0.212).
*Note (5): The error is assumed to be dominated by the error on X-ray luminosity.
*Note (6): Equation (3) enables quantitative measurement of how the X-ray
*     luminosity of a given QSO differs from that expected for a typical
*     QSO of the same UV luminosity. This difference is parameterized by:
*     {Delta}{alpha}_ox_={alpha}_ox_-{alpha}_ox_(L_2500_), where {alpha}_ox_
*     is the value observed for a particular source observed to have a UV
*     luminosity of L_2500_, and {alpha}_ox_(L_2500_) is determined for a
*     typical QSO from Equation (3).

C=============================================================================

C  Loading file 'table1.dat'	! DR5 BAL Catalog

C  Format for file interpretation

    1 format(
     +  A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,I5,1X,I4,1X,I3,1X,F7.1,1X,
     +  F7.1,1X,F5.1,1X,I6,1X,I6,1X,F4.2,1X,F7.1,1X,F7.1,1X,F6.1,1X,
     +  I6,1X,I6,1X,F4.2,1X,F7.1,1X,F7.1,1X,F6.1,1X,I6,1X,I6,1X,F4.2,
     +  1X,F7.1,1X,F7.1,1X,F6.1,1X,I6,1X,I6,1X,F4.2,1X,I1,1X,I1,1X,I1,
     +  1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,F5.1,1X,F6.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,5035
        read(1,'(A265)')ar__
        read(ar__,1)
     +  SDSS(i__),RAdeg(i__),DEdeg(i__),z(i__),MJD(i__),Plate(i__),
     +  Fiber(i__),BI_Si(i__),BIO_Si(i__),EW_Si(i__),vmin_Si(i__),
     +  vmax_Si(i__),fd_Si(i__),BI_C(i__),BIO_C(i__),EW_C(i__),
     +  vmin_C(i__),vmax_C(i__),fd_C(i__),BI_Al(i__),BIO_Al(i__),
     +  EW_Al(i__),vmin_Al(i__),vmax_Al(i__),fd_Al(i__),BI_Mg(i__),
     +  BIO_Mg(i__),EW_Mg(i__),vmin_Mg(i__),vmax_Mg(i__),fd_Mg(i__),
     +  f_SiIVEmL(i__),f_CIVEmL(i__),f_AlIIIEmL(i__),f_MgIIEmL(i__),
     +  f_SiIVBB(i__),f_CIVBB(i__),f_AlIIIBB(i__),f_MgIIBB(i__),
     +  f_CIV4BWA(i__),SN1700(i__),logF1400(i__),logF2500(i__)
        if(ar__(47:51) .EQ. '') MJD(i__) = iNULL__
        if(ar__(53:56) .EQ. '') Plate(i__) = iNULL__
        if(ar__(58:60) .EQ. '') Fiber(i__) = iNULL__
        if(ar__(62:68) .EQ. '') BI_Si(i__) = rNULL__
        if(ar__(70:76) .EQ. '') BIO_Si(i__) = rNULL__
        if(ar__(78:82) .EQ. '') EW_Si(i__) = rNULL__
        if(ar__(84:89) .EQ. '') vmin_Si(i__) = iNULL__
        if(ar__(91:96) .EQ. '') vmax_Si(i__) = iNULL__
        if(ar__(98:101) .EQ. '') fd_Si(i__) = rNULL__
        if(ar__(103:109) .EQ. '') BI_C(i__) = rNULL__
        if(ar__(111:117) .EQ. '') BIO_C(i__) = rNULL__
        if(ar__(119:124) .EQ. '') EW_C(i__) = rNULL__
        if(ar__(126:131) .EQ. '') vmin_C(i__) = iNULL__
        if(ar__(133:138) .EQ. '') vmax_C(i__) = iNULL__
        if(ar__(140:143) .EQ. '') fd_C(i__) = rNULL__
        if(ar__(145:151) .EQ. '') BI_Al(i__) = rNULL__
        if(ar__(153:159) .EQ. '') BIO_Al(i__) = rNULL__
        if(ar__(161:166) .EQ. '') EW_Al(i__) = rNULL__
        if(ar__(168:173) .EQ. '') vmin_Al(i__) = iNULL__
        if(ar__(175:180) .EQ. '') vmax_Al(i__) = iNULL__
        if(ar__(182:185) .EQ. '') fd_Al(i__) = rNULL__
        if(ar__(187:193) .EQ. '') BI_Mg(i__) = rNULL__
        if(ar__(195:201) .EQ. '') BIO_Mg(i__) = rNULL__
        if(ar__(203:208) .EQ. '') EW_Mg(i__) = rNULL__
        if(ar__(210:215) .EQ. '') vmin_Mg(i__) = iNULL__
        if(ar__(217:222) .EQ. '') vmax_Mg(i__) = iNULL__
        if(ar__(224:227) .EQ. '') fd_Mg(i__) = rNULL__
        if(ar__(237:237) .EQ. '') f_SiIVBB(i__) = iNULL__
        if(ar__(239:239) .EQ. '') f_CIVBB(i__) = iNULL__
        if(ar__(241:241) .EQ. '') f_AlIIIBB(i__) = iNULL__
        if(ar__(243:243) .EQ. '') f_MgIIBB(i__) = iNULL__
        if(ar__(247:251) .EQ. '') SN1700(i__) = rNULL__
        if(ar__(253:258) .EQ. '') logF1400(i__) = rNULL__
        if(ar__(260:265) .EQ. '') logF2500(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  SDSS(i__),RAdeg(i__),DEdeg(i__),z(i__),MJD(i__),Plate(i__),
     +  Fiber(i__),BI_Si(i__),BIO_Si(i__),EW_Si(i__),vmin_Si(i__),
     +  vmax_Si(i__),fd_Si(i__),BI_C(i__),BIO_C(i__),EW_C(i__),
     +  vmin_C(i__),vmax_C(i__),fd_C(i__),BI_Al(i__),BIO_Al(i__),
     +  EW_Al(i__),vmin_Al(i__),vmax_Al(i__),fd_Al(i__),BI_Mg(i__),
     +  BIO_Mg(i__),EW_Mg(i__),vmin_Mg(i__),vmax_Mg(i__),fd_Mg(i__),
     +  f_SiIVEmL(i__),f_CIVEmL(i__),f_AlIIIEmL(i__),f_MgIIEmL(i__),
     +  f_SiIVBB(i__),f_CIVBB(i__),f_AlIIIBB(i__),f_MgIIBB(i__),
     +  f_CIV4BWA(i__),SN1700(i__),logF1400(i__),logF2500(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! DR5 BAL catalog duplicate spectra

C  Format for file interpretation

    2 format(
     +  A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,I5,1X,I4,1X,I3,1X,F7.1,1X,
     +  F7.1,1X,F5.1,1X,I6,1X,I6,1X,F4.2,1X,F7.1,1X,F7.1,1X,F6.1,1X,
     +  I6,1X,I6,1X,F4.2,1X,F7.1,1X,F7.1,1X,F6.1,1X,I6,1X,I6,1X,F4.2,
     +  1X,F7.1,1X,F7.1,1X,F6.1,1X,I6,1X,I6,1X,F4.2,1X,I1,1X,I1,1X,I1,
     +  1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,F5.1,1X,F6.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,569
        read(1,'(A265)')ar__1
        read(ar__1,2)
     +  SDSS_1(i__),RAdeg_1(i__),DEdeg_1(i__),z_1(i__),MJD_1(i__),
     +  Plate_1(i__),Fiber_1(i__),BI_Si_1(i__),BIO_Si_1(i__),
     +  EW_Si_1(i__),vmin_Si_1(i__),vmax_Si_1(i__),fd_Si_1(i__),
     +  BI_C_1(i__),BIO_C_1(i__),EW_C_1(i__),vmin_C_1(i__),
     +  vmax_C_1(i__),fd_C_1(i__),BI_Al_1(i__),BIO_Al_1(i__),
     +  EW_Al_1(i__),vmin_Al_1(i__),vmax_Al_1(i__),fd_Al_1(i__),
     +  BI_Mg_1(i__),BIO_Mg_1(i__),EW_Mg_1(i__),vmin_Mg_1(i__),
     +  vmax_Mg_1(i__),fd_Mg_1(i__),f_SiIVEmL_1(i__),f_CIVEmL_1(i__),
     +  f_AlIIIEmL_1(i__),f_MgIIEmL_1(i__),f_SiIVBB_1(i__),
     +  f_CIVBB_1(i__),f_AlIIIBB_1(i__),f_MgIIBB_1(i__),
     +  f_CIV4BWA_1(i__),SN1700_1(i__),logF1400_1(i__),logF2500_1(i__)
        if(ar__1(47:51) .EQ. '') MJD_1(i__) = iNULL__
        if(ar__1(53:56) .EQ. '') Plate_1(i__) = iNULL__
        if(ar__1(58:60) .EQ. '') Fiber_1(i__) = iNULL__
        if(ar__1(62:68) .EQ. '') BI_Si_1(i__) = rNULL__
        if(ar__1(70:76) .EQ. '') BIO_Si_1(i__) = rNULL__
        if(ar__1(78:82) .EQ. '') EW_Si_1(i__) = rNULL__
        if(ar__1(84:89) .EQ. '') vmin_Si_1(i__) = iNULL__
        if(ar__1(91:96) .EQ. '') vmax_Si_1(i__) = iNULL__
        if(ar__1(98:101) .EQ. '') fd_Si_1(i__) = rNULL__
        if(ar__1(103:109) .EQ. '') BI_C_1(i__) = rNULL__
        if(ar__1(111:117) .EQ. '') BIO_C_1(i__) = rNULL__
        if(ar__1(119:124) .EQ. '') EW_C_1(i__) = rNULL__
        if(ar__1(126:131) .EQ. '') vmin_C_1(i__) = iNULL__
        if(ar__1(133:138) .EQ. '') vmax_C_1(i__) = iNULL__
        if(ar__1(140:143) .EQ. '') fd_C_1(i__) = rNULL__
        if(ar__1(145:151) .EQ. '') BI_Al_1(i__) = rNULL__
        if(ar__1(153:159) .EQ. '') BIO_Al_1(i__) = rNULL__
        if(ar__1(161:166) .EQ. '') EW_Al_1(i__) = rNULL__
        if(ar__1(168:173) .EQ. '') vmin_Al_1(i__) = iNULL__
        if(ar__1(175:180) .EQ. '') vmax_Al_1(i__) = iNULL__
        if(ar__1(182:185) .EQ. '') fd_Al_1(i__) = rNULL__
        if(ar__1(187:193) .EQ. '') BI_Mg_1(i__) = rNULL__
        if(ar__1(195:201) .EQ. '') BIO_Mg_1(i__) = rNULL__
        if(ar__1(203:208) .EQ. '') EW_Mg_1(i__) = rNULL__
        if(ar__1(210:215) .EQ. '') vmin_Mg_1(i__) = iNULL__
        if(ar__1(217:222) .EQ. '') vmax_Mg_1(i__) = iNULL__
        if(ar__1(224:227) .EQ. '') fd_Mg_1(i__) = rNULL__
        if(ar__1(237:237) .EQ. '') f_SiIVBB_1(i__) = iNULL__
        if(ar__1(239:239) .EQ. '') f_CIVBB_1(i__) = iNULL__
        if(ar__1(241:241) .EQ. '') f_AlIIIBB_1(i__) = iNULL__
        if(ar__1(243:243) .EQ. '') f_MgIIBB_1(i__) = iNULL__
        if(ar__1(247:251) .EQ. '') SN1700_1(i__) = rNULL__
        if(ar__1(253:258) .EQ. '') logF1400_1(i__) = rNULL__
        if(ar__1(260:265) .EQ. '') logF2500_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  SDSS_1(i__),RAdeg_1(i__),DEdeg_1(i__),z_1(i__),MJD_1(i__),
     +  Plate_1(i__),Fiber_1(i__),BI_Si_1(i__),BIO_Si_1(i__),
     +  EW_Si_1(i__),vmin_Si_1(i__),vmax_Si_1(i__),fd_Si_1(i__),
     +  BI_C_1(i__),BIO_C_1(i__),EW_C_1(i__),vmin_C_1(i__),
     +  vmax_C_1(i__),fd_C_1(i__),BI_Al_1(i__),BIO_Al_1(i__),
     +  EW_Al_1(i__),vmin_Al_1(i__),vmax_Al_1(i__),fd_Al_1(i__),
     +  BI_Mg_1(i__),BIO_Mg_1(i__),EW_Mg_1(i__),vmin_Mg_1(i__),
     +  vmax_Mg_1(i__),fd_Mg_1(i__),f_SiIVEmL_1(i__),f_CIVEmL_1(i__),
     +  f_AlIIIEmL_1(i__),f_MgIIEmL_1(i__),f_SiIVBB_1(i__),
     +  f_CIVBB_1(i__),f_AlIIIBB_1(i__),f_MgIIBB_1(i__),
     +  f_CIV4BWA_1(i__),SN1700_1(i__),logF1400_1(i__),logF2500_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! DR5 non-BALs with CIV BlueWingAbs flag set

C  Format for file interpretation

    3 format(A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,I1,1X,I1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,118
        read(1,'(A51)')ar__2
        read(ar__2,3)
     +  SDSS_2(i__),RAdeg_2(i__),DEdeg_2(i__),z_2(i__),SiIV(i__),
     +  AlIII(i__),MgII(i__)
        if(ar__2(47:47) .EQ. '') SiIV(i__) = iNULL__
        if(ar__2(49:49) .EQ. '') AlIII(i__) = iNULL__
        if(ar__2(51:51) .EQ. '') MgII(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  SDSS_2(i__),RAdeg_2(i__),DEdeg_2(i__),z_2(i__),SiIV(i__),
     +  AlIII(i__),MgII(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table4.dat'	! BALs with absorption from other ions

C  Format for file interpretation

    4 format(A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,I1,1X,I1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,92
        read(1,'(A51)')ar__3
        read(ar__3,4)
     +  SDSS_3(i__),RAdeg_3(i__),DEdeg_3(i__),z_3(i__),SiIV_1(i__),
     +  AlIII_1(i__),MgII_1(i__)
        if(ar__3(47:47) .EQ. '') SiIV_1(i__) = iNULL__
        if(ar__3(49:49) .EQ. '') AlIII_1(i__) = iNULL__
        if(ar__3(51:51) .EQ. '') MgII_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  SDSS_3(i__),RAdeg_3(i__),DEdeg_3(i__),z_3(i__),SiIV_1(i__),
     +  AlIII_1(i__),MgII_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table5.dat'	! X-ray source information

C  Format for file interpretation

    5 format(
     +  A18,1X,A10,1X,F4.1,1X,A7,1X,A7,1X,I1,1X,A1,F5.2,1X,A1,F5.2,1X,
     +  A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,73
        read(1,'(A90)')ar__4
        read(ar__4,5)
     +  SDSS_4(i__),ObsID(i__),Exp(i__),Sct(i__),Hct(i__),Det(i__),
     +  l_logR_(i__),logR_(i__),l_logLx(i__),logLx(i__),l_aox(i__),
     +  aox(i__),e_aox(i__),l_Daox(i__),Daox(i__),e_Daox(i__)
        if(ar__4(75:78) .EQ. '') e_aox(i__) = rNULL__
        if(ar__4(87:90) .EQ. '') e_Daox(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  SDSS_4(i__),ObsID(i__),Exp(i__),Sct(i__),Hct(i__),Det(i__),
     +  l_logR_(i__),logR_(i__),l_logLx(i__),logLx(i__),l_aox(i__),
     +  aox(i__),e_aox(i__),l_Daox(i__),Daox(i__),e_Daox(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end