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