FORTRAN Generation
(/./ftp/cats/J/AJ/133/313)

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