FORTRAN Generation
(/./ftp/cats/J/ApJS/179/142)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/179/142 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 2026-Aug-11
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/ApJS/179/142     Catalog of X-ray sources in NGC 3379     (Brassington+, 2008)
*================================================================================
*Deep Chandra monitoring observations of NGC 3379: catalog of source properties.
*    Brassington N.J., Fabbiano G., Kim D.-W., Zezas A., Zepf S., Kundu A.,
*    Angelini L., Davies R.L., Gallagher J., Kalogera V., Fragos T., King A.R.,
*    Pellegrini S., Trinchieri G.
*   <Astrophys. J. Suppl. Ser., 179, 142-165 (2008)>
*   =2008ApJS..179..142B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table4.dat'	! Master Source List

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   n__BFK2008_(nr__) ! [<>] Indicates detection with SNR < or > 3
      integer*4     v_BFK2008_ (nr__) ! Master identification number
      character*1   f__BFK2008_(nr__) ! [b] source within optical limits (1)
      character*16  CXOU       (nr__) ! IAU identification (JHHMMSS.s+DDMMSS) (2)
      integer*4     RAh        (nr__) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs        (nr__) ! (s) Second of Right Ascension (J2000)
      character*1   DE_        (nr__) ! Sign of the Declination (J2000)
      integer*4     DEd        (nr__) ! (deg) Degree of Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Arcminute of Declination (J2000)
      real*4        DEs        (nr__) ! (arcsec) Arcsecond of Declination (J2000)
      real*4        Rad        (nr__) ! (arcsec) Source extraction radius
      real*4        e_pos      (nr__) ! (arcsec) ? Position uncertainty; see Eq.1
      real*4        SNR        (nr__) ! ? Signal-to-Noise ratio
      character*2   l_logLX    (nr__) ! [<= ] Limit flag on logLX
      real*4        logLX      (nr__) ! ([10-7W]) Log of 0.3-8.0keV luminosity (3)
      character*2   l_3logLX   (nr__) ! Limit flag on logLX
      real*4        v3logLX    (nr__) ! ([10-7W]) ? The 3{sigma} uncertainty in logLX (3)
      character*3   LT         (nr__) ! Long-term source variability (4)
      character*1   BB         (nr__) ! Short-term variability from Bayesian block
*                                    analysis (5)
      character*1   K_S        (nr__) ! Short-term variability from
*                                    Kolmogorov-Smirnov test (5)
      character*6   Opt        (nr__) ! Optical association (6)
      real*4        DG         (nr__) ! (arcsec) Distance from center of galaxy
      character*1   n_DG       (nr__) ! [b] source within optical limits (1)
      character*7   Note       (nr__) ! Additional information (7)
*Note (1): b = source is within the D_25_ ellipse of the galaxy
*     (duplicated column)
*Note (2): Following the convention "CXOU Jhhmmss.s+/-ddmmss".
*          [BFK2008] JHHMMSS.s+DDMMSS in Simbad.
*Note (3): Assuming D=10.6Mpc. For sources only detected in a single observation,
*          1{sigma} upper limits from the co-added observation are quoted (with
*          3{sigma} upper limits from the detected observation following.
*Note (4): Long-term variability as follows:
*    N = non-variable source;
*    V = variable source;
*   TC = transient candidate;
*  PTC = possible transient.
*Note (5): Short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (6): Optical association as follows:
*   GC = optical source is a confirmed globular cluster;
*   BG = optical source is a background object;
* corr = match has been defined as correlationed;
* exmt = match between 0.6" and 3" in separation;
* none = source inside the field of view of the HST observation, but has no
*        optical counterpart;
*    X = source within the HST FOV but also within 5" of the nucleus.
*Note (7): Note as follows:
*      X = source detected in a single observation only,
*     O1 = single overlapping source;
*     O2 = complicated overlapping sources;
*   BKG? = possible background object;
*double? = possible double source.

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

C  Declarations for 'table5.dat'	! Source counts, hardness ratios and color-color
                             values: Coadded Observation

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

      integer*4     v_BFK2008__1(nr__1) ! Master identification number
      character*2   l_srcB     (nr__1) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB       (nr__1) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB     (nr__1) ! (ct) ? Uncertainty in srcB
      real*4        srcS1      (nr__1) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1    (nr__1) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2      (nr__1) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2    (nr__1) ! (ct) ? Uncertainty in srcS2
      real*4        srcS       (nr__1) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS     (nr__1) ! (ct) ? Uncertainty in srcS
      real*4        srcH       (nr__1) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH     (nr__1) ! (ct) ? Uncertainty in srcH
      real*4        srcSc      (nr__1) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc    (nr__1) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc      (nr__1) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc    (nr__1) ! (ct) ? Uncertainty in srcHc
      real*4        HR         (nr__1) ! ? Hardness ratio (G1)
      real*4        e_HR       (nr__1) ! ? Lower limit uncertainty in HR
      real*4        E_HR_1     (nr__1) ! ? Upper limit uncertainty in HR
      real*4        C21        (nr__1) ! ? X-ray color C21 (log(S1/S2) (G1)
      real*4        e_C21      (nr__1) ! ? Lower limit uncertainty in C21
      real*4        E_C21_1    (nr__1) ! ? Upper limit uncertainty in C21
      real*4        C32        (nr__1) ! ? X-ray color C32 (log(S2/H) (G1)
      real*4        e_C32      (nr__1) ! ? Lower limit uncertainty in C32
      real*4        E_C32_1    (nr__1) ! ? Upper limit uncertainty in C32
      character*2   l_logLX_1  (nr__1) ! [<= ] Upper 68% confidence level flag on 
*                                         logLX
      real*4        logLX_1    (nr__1) ! ([10-7W]) Log of 0.3-8.0keV luminosity (2)
*Note (2): Assuming D=10.6Mpc.

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

C  Declarations for 'table6.dat'	! Source counts, hardness ratios, color-color
                             values and variability: Observation 1

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

      integer*4     v_BFK2008__2(nr__2) ! Master identification number
      character*2   l_srcB_1   (nr__2) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB_1     (nr__2) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB_1   (nr__2) ! (ct) ? Uncertainty in srcB
      real*4        srcS1_1    (nr__2) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1_1  (nr__2) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2_1    (nr__2) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2_1  (nr__2) ! (ct) ? Uncertainty in srcS2
      real*4        srcS_1     (nr__2) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS_1   (nr__2) ! (ct) ? Uncertainty in srcS
      real*4        srcH_1     (nr__2) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH_1   (nr__2) ! (ct) ? Uncertainty in srcH
      real*4        srcSc_1    (nr__2) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc_1  (nr__2) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc_1    (nr__2) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc_1  (nr__2) ! (ct) ? Uncertainty in srcHc
      real*4        HR_1       (nr__2) ! ? Hardness ratio (G1)
      real*4        e_HR_2     (nr__2) ! ? Lower limit uncertainty in HR
      real*4        E_HR_3     (nr__2) ! ? Upper limit uncertainty in HR
      real*4        C21_1      (nr__2) ! ? X-ray color C21 = log(S1/S2) (G1)
      real*4        e_C21_2    (nr__2) ! ? Lower limit uncertainty in C21
      real*4        E_C21_3    (nr__2) ! ? Upper limit uncertainty in C21
      real*4        C32_1      (nr__2) ! ? X-ray color C32 = log(S2/H) (G1)
      real*4        e_C32_2    (nr__2) ! ? Lower limit uncertainty in C32
      real*4        E_C32_3    (nr__2) ! ? Upper limit uncertainty in C32
      character*1   BB_1       (nr__2) ! Short-term variability from Bayesian block
*                                   analysis (2)
      character*1   K_S_1      (nr__2) ! Short-term variability from
*                                   Kolmogorov-Smirnov test (2)
      real*4        Sig        (nr__2) ! ? Significance of change in L_X_ (3)
      character*2   l_logLX_2  (nr__2) ! [<=] Upper 68% confidence level flag on logLX
      real*4        logLX_2    (nr__2) ! ([10-7W]) Log of 0.3-8.0keV luminosity (4)
*Note (2): short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (3): Between the previous observation and the current observation
*          respectively (Equ 2).
*Note (4): Assuming D=10.6 Mpc.

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

C  Declarations for 'table7.dat'	! Source counts, hardness ratios, color-color
                             values and variability: Observation 2

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

      integer*4     v_BFK2008__3(nr__3) ! Master identification number
      character*2   l_srcB_2   (nr__3) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB_2     (nr__3) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB_2   (nr__3) ! (ct) ? Uncertainty in srcB
      real*4        srcS1_2    (nr__3) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1_2  (nr__3) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2_2    (nr__3) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2_2  (nr__3) ! (ct) ? Uncertainty in srcS2
      real*4        srcS_2     (nr__3) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS_2   (nr__3) ! (ct) ? Uncertainty in srcS
      real*4        srcH_2     (nr__3) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH_2   (nr__3) ! (ct) ? Uncertainty in srcH
      real*4        srcSc_2    (nr__3) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc_2  (nr__3) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc_2    (nr__3) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc_2  (nr__3) ! (ct) ? Uncertainty in srcHc
      real*4        HR_2       (nr__3) ! ? Hardness ratio (G1)
      real*4        e_HR_4     (nr__3) ! ? Lower limit uncertainty in HR
      real*4        E_HR_5     (nr__3) ! ? Upper limit uncertainty in HR
      real*4        C21_2      (nr__3) ! ? X-ray color C21 = log(S1/S2) (G1)
      real*4        e_C21_4    (nr__3) ! ? Lower limit uncertainty in C21
      real*4        E_C21_5    (nr__3) ! ? Upper limit uncertainty in C21
      real*4        C32_2      (nr__3) ! ? X-ray color C32 = log(S2/H) (G1)
      real*4        e_C32_4    (nr__3) ! ? Lower limit uncertainty in C32
      real*4        E_C32_5    (nr__3) ! ? Upper limit uncertainty in C32
      character*1   BB_2       (nr__3) ! Short-term variability from Bayesian block
*                                   analysis (2)
      character*1   K_S_2      (nr__3) ! Short-term variability from
*                                   Kolmogorov-Smirnov test (2)
      real*4        Sig_1      (nr__3) ! ? Significance of change in L_X_ (3)
      character*2   l_logLX_3  (nr__3) ! [<=] Upper 68% confidence level flag on logLX
      real*4        logLX_3    (nr__3) ! ([10-7W]) Log of 0.3-8.0keV luminosity (4)
*Note (2): short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (3): Between the previous observation and the current observation
*          respectively (Equ 2).
*Note (4): Assuming D=10.6 Mpc.

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

C  Declarations for 'table8.dat'	! Source counts, hardness ratios, color-color
                             values and variability: Observation 3

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

      integer*4     v_BFK2008__4(nr__4) ! Master identification number
      character*2   l_srcB_3   (nr__4) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB_3     (nr__4) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB_3   (nr__4) ! (ct) ? Uncertainty in srcB
      real*4        srcS1_3    (nr__4) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1_3  (nr__4) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2_3    (nr__4) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2_3  (nr__4) ! (ct) ? Uncertainty in srcS2
      real*4        srcS_3     (nr__4) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS_3   (nr__4) ! (ct) ? Uncertainty in srcS
      real*4        srcH_3     (nr__4) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH_3   (nr__4) ! (ct) ? Uncertainty in srcH
      real*4        srcSc_3    (nr__4) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc_3  (nr__4) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc_3    (nr__4) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc_3  (nr__4) ! (ct) ? Uncertainty in srcHc
      real*4        HR_3       (nr__4) ! ? Hardness ratio (G1)
      real*4        e_HR_6     (nr__4) ! ? Lower limit uncertainty in HR
      real*4        E_HR_7     (nr__4) ! ? Upper limit uncertainty in HR
      real*4        C21_3      (nr__4) ! ? X-ray color C21 = log(S1/S2) (G1)
      real*4        e_C21_6    (nr__4) ! ? Lower limit uncertainty in C21
      real*4        E_C21_7    (nr__4) ! ? Upper limit uncertainty in C21
      real*4        C32_3      (nr__4) ! ? X-ray color C32 = log(S2/H) (G1)
      real*4        e_C32_6    (nr__4) ! ? Lower limit uncertainty in C32
      real*4        E_C32_7    (nr__4) ! ? Upper limit uncertainty in C32
      character*1   BB_3       (nr__4) ! Short-term variability from Bayesian block
*                                   analysis (2)
      character*1   K_S_3      (nr__4) ! Short-term variability from
*                                   Kolmogorov-Smirnov test (2)
      real*4        Sig_2      (nr__4) ! ? Significance of change in L_X_ (3)
      character*2   l_logLX_4  (nr__4) ! [<=] Upper 68% confidence level flag on logLX
      real*4        logLX_4    (nr__4) ! ([10-7W]) Log of 0.3-8.0keV luminosity (4)
*Note (2): short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (3): Between the previous observation and the current observation
*          respectively (Equ 2).
*Note (4): Assuming D=10.6 Mpc.

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

C  Declarations for 'table9.dat'	! Source counts, hardness ratios, color-color
                             values and variability: Observation 4

      integer*4 nr__5
      parameter (nr__5=132)	! Number of records
      character*150 ar__5  	! Full-size record

      integer*4     v_BFK2008__5(nr__5) ! Master identification number
      character*2   l_srcB_4   (nr__5) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB_4     (nr__5) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB_4   (nr__5) ! (ct) ? Uncertainty in srcB
      real*4        srcS1_4    (nr__5) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1_4  (nr__5) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2_4    (nr__5) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2_4  (nr__5) ! (ct) ? Uncertainty in srcS2
      real*4        srcS_4     (nr__5) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS_4   (nr__5) ! (ct) ? Uncertainty in srcS
      real*4        srcH_4     (nr__5) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH_4   (nr__5) ! (ct) ? Uncertainty in srcH
      real*4        srcSc_4    (nr__5) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc_4  (nr__5) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc_4    (nr__5) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc_4  (nr__5) ! (ct) ? Uncertainty in srcHc
      real*4        HR_4       (nr__5) ! ? Hardness ratio (G1)
      real*4        e_HR_8     (nr__5) ! ? Lower limit uncertainty in HR
      real*4        E_HR_9     (nr__5) ! ? Upper limit uncertainty in HR
      real*4        C21_4      (nr__5) ! ? X-ray color C21 = log(S1/S2) (G1)
      real*4        e_C21_8    (nr__5) ! ? Lower limit uncertainty in C21
      real*4        E_C21_9    (nr__5) ! ? Upper limit uncertainty in C21
      real*4        C32_4      (nr__5) ! ? X-ray color C32 = log(S2/H) (G1)
      real*4        e_C32_8    (nr__5) ! ? Lower limit uncertainty in C32
      real*4        E_C32_9    (nr__5) ! ? Upper limit uncertainty in C32
      character*1   BB_4       (nr__5) ! Short-term variability from Bayesian block
*                                   analysis (2)
      character*1   K_S_4      (nr__5) ! Short-term variability from
*                                   Kolmogorov-Smirnov test (2)
      real*4        Sig_3      (nr__5) ! ? Significance of change in L_X_ (3)
      character*2   l_logLX_5  (nr__5) ! [<=] Upper 68% confidence level flag on logLX
      real*4        logLX_5    (nr__5) ! ([10-7W]) Log of 0.3-8.0keV luminosity (4)
*Note (2): short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (3): Between the previous observation and the current observation
*          respectively (Equ 2).
*Note (4): Assuming D=10.6 Mpc.

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

C  Declarations for 'table10.dat'	! Source counts, hardness ratios, color-color
                             values and variability: Observation 5

      integer*4 nr__6
      parameter (nr__6=132)	! Number of records
      character*150 ar__6  	! Full-size record

      integer*4     v_BFK2008__6(nr__6) ! Master identification number
      character*2   l_srcB_5   (nr__6) ! [<= ] Upper 68% confidence level flag on srcB
      real*4        srcB_5     (nr__6) ! (ct) Net source B-band counts (G1)
      real*4        e_srcB_5   (nr__6) ! (ct) ? Uncertainty in srcB
      real*4        srcS1_5    (nr__6) ! (ct) ? Net source S1-band counts (G1)
      real*4        e_srcS1_5  (nr__6) ! (ct) ? Uncertainty in srcS1
      real*4        srcS2_5    (nr__6) ! (ct) ? Net source S2-band counts (G1)
      real*4        e_srcS2_5  (nr__6) ! (ct) ? Uncertainty in srcS2
      real*4        srcS_5     (nr__6) ! (ct) ? Net source S-band counts (G1)
      real*4        e_srcS_5   (nr__6) ! (ct) ? Uncertainty in srcS
      real*4        srcH_5     (nr__6) ! (ct) ? Net source H-band counts (G1)
      real*4        e_srcH_5   (nr__6) ! (ct) ? Uncertainty in srcH
      real*4        srcSc_5    (nr__6) ! (ct) ? Net source Sc-band counts (G1)
      real*4        e_srcSc_5  (nr__6) ! (ct) ? Uncertainty in srcSc
      real*4        srcHc_5    (nr__6) ! (ct) ? Net source Hc-band counts (G1)
      real*4        e_srcHc_5  (nr__6) ! (ct) ? Uncertainty in srcHc
      real*4        HR_5       (nr__6) ! ? Hardness ratio (G1)
      real*4        e_HR_10    (nr__6) ! ? Lower limit uncertainty in HR
      real*4        E_HR_11    (nr__6) ! ? Upper limit uncertainty in HR
      real*4        C21_5      (nr__6) ! ? X-ray color C21 = log(S1/S2) (G1)
      real*4        e_C21_10   (nr__6) ! ? Lower limit uncertainty in C21
      real*4        E_C21_11   (nr__6) ! ? Upper limit uncertainty in C21
      real*4        C32_5      (nr__6) ! ? X-ray color C32 = log(S2/H) (G1)
      real*4        e_C32_10   (nr__6) ! ? Lower limit uncertainty in C32
      real*4        E_C32_11   (nr__6) ! ? Upper limit uncertainty in C32
      character*1   BB_5       (nr__6) ! Short-term variability from Bayesian block
*                                   analysis (2)
      character*1   K_S_5      (nr__6) ! Short-term variability from
*                                   Kolmogorov-Smirnov test (2)
      real*4        Sig_4      (nr__6) ! ? Significance of change in L_X_ (3)
      character*2   l_logLX_6  (nr__6) ! [<=] Upper 68% confidence level flag on logLX
      real*4        logLX_6    (nr__6) ! ([10-7W]) Log of 0.3-8.0keV luminosity (4)
*Note (2): short-term variability as follows:
*    N = non-variable in all observations;
*    V = variable in at least one observation;
*    P = possible variability in at least one observation.
*Note (3): Between the previous observation and the current observation
*          respectively (Equ 2).
*Note (4): Assuming D=10.6 Mpc.

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

C  Declarations for 'table15.dat'	! Raw source and background counts from the
                             co-added observation

      integer*4 nr__7
      parameter (nr__7=132)	! Number of records
      character*80 ar__7  	! Full-size record

      integer*4     v_BFK2008__7(nr__7) ! Master identification number
      real*4        netB       (nr__7) ! (ct) ? Net B counts (G1)
      integer*4     rawB       (nr__7) ! (ct) ? Source B-band counts (G1)
      real*4        bkgB       (nr__7) ! (ct) ? Background B-band counts (G1)
      integer*4     rawS1      (nr__7) ! (ct) ? Source S1-band counts (G1)
      real*4        bkgS1      (nr__7) ! (ct) ? Background S1-band counts (G1)
      integer*4     rawS2      (nr__7) ! (ct) ? Source S2-band counts (G1)
      real*4        bkgS2      (nr__7) ! (ct) ? Background S2-band counts (G1)
      integer*4     rawH       (nr__7) ! (ct) ? Source H-band counts (G1)
      real*4        bkgH       (nr__7) ! (ct) ? Background H-band counts (G1)
      integer*4     rawBc      (nr__7) ! (ct) ? Source Bc-band counts (G1)
      real*4        bkgBc      (nr__7) ! (ct) ? Background Bc-band counts (G1)
      integer*4     rawSc      (nr__7) ! (ct) ? Source Sc-band counts (G1)
      real*4        bkgSc      (nr__7) ! (ct) ? Background Sc-band counts (G1)
      integer*4     rawHc      (nr__7) ! (ct) ? Source Hc-band counts (G1)
      real*4        bkgHc      (nr__7) ! (ct) ? Background Hc-band counts (G1)

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

C  Loading file 'table4.dat'	! Master Source List

C  Format for file interpretation

    1 format(
     +  A1,1X,I3,1X,A1,1X,A16,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,
     +  F4.1,1X,F4.2,1X,F4.2,1X,F4.1,1X,A2,F4.1,1X,A2,F4.1,1X,A3,1X,
     +  A1,1X,A1,1X,A6,1X,F5.1,1X,A1,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,132
        read(1,'(A107)')ar__
        read(ar__,1)
     +  n__BFK2008_(i__),v_BFK2008_(i__),f__BFK2008_(i__),CXOU(i__),
     +  RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),Rad(i__),e_pos(i__),SNR(i__),l_logLX(i__),logLX(i__),
     +  l_3logLX(i__),v3logLX(i__),LT(i__),BB(i__),K_S(i__),Opt(i__),
     +  DG(i__),n_DG(i__),Note(i__)
        if(ar__(54:57) .EQ. '') e_pos(i__) = rNULL__
        if(ar__(59:62) .EQ. '') SNR(i__) = rNULL__
        if(ar__(73:76) .EQ. '') v3logLX(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  n__BFK2008_(i__),v_BFK2008_(i__),f__BFK2008_(i__),CXOU(i__),
     +  RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),Rad(i__),e_pos(i__),SNR(i__),l_logLX(i__),logLX(i__),
     +  l_3logLX(i__),v3logLX(i__),LT(i__),BB(i__),K_S(i__),Opt(i__),
     +  DG(i__),n_DG(i__),Note(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Source counts, hardness ratios and color-color
*                             values: Coadded Observation

C  Format for file interpretation

    2 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,132
        read(1,'(A144)')ar__1
        read(ar__1,2)
     +  v_BFK2008__1(i__),l_srcB(i__),srcB(i__),e_srcB(i__),
     +  srcS1(i__),e_srcS1(i__),srcS2(i__),e_srcS2(i__),srcS(i__),
     +  e_srcS(i__),srcH(i__),e_srcH(i__),srcSc(i__),e_srcSc(i__),
     +  srcHc(i__),e_srcHc(i__),HR(i__),e_HR(i__),E_HR_1(i__),
     +  C21(i__),e_C21(i__),E_C21_1(i__),C32(i__),e_C32(i__),
     +  E_C32_1(i__),l_logLX_1(i__),logLX_1(i__)
        if(ar__1(14:17) .EQ. '') e_srcB(i__) = rNULL__
        if(ar__1(19:24) .EQ. '') srcS1(i__) = rNULL__
        if(ar__1(26:29) .EQ. '') e_srcS1(i__) = rNULL__
        if(ar__1(31:36) .EQ. '') srcS2(i__) = rNULL__
        if(ar__1(38:41) .EQ. '') e_srcS2(i__) = rNULL__
        if(ar__1(43:48) .EQ. '') srcS(i__) = rNULL__
        if(ar__1(50:53) .EQ. '') e_srcS(i__) = rNULL__
        if(ar__1(55:60) .EQ. '') srcH(i__) = rNULL__
        if(ar__1(62:65) .EQ. '') e_srcH(i__) = rNULL__
        if(ar__1(67:72) .EQ. '') srcSc(i__) = rNULL__
        if(ar__1(74:77) .EQ. '') e_srcSc(i__) = rNULL__
        if(ar__1(79:84) .EQ. '') srcHc(i__) = rNULL__
        if(ar__1(86:89) .EQ. '') e_srcHc(i__) = rNULL__
        if(ar__1(91:95) .EQ. '') HR(i__) = rNULL__
        if(ar__1(97:100) .EQ. '') e_HR(i__) = rNULL__
        if(ar__1(102:105) .EQ. '') E_HR_1(i__) = rNULL__
        if(ar__1(107:111) .EQ. '') C21(i__) = rNULL__
        if(ar__1(113:116) .EQ. '') e_C21(i__) = rNULL__
        if(ar__1(118:121) .EQ. '') E_C21_1(i__) = rNULL__
        if(ar__1(123:127) .EQ. '') C32(i__) = rNULL__
        if(ar__1(129:132) .EQ. '') e_C32(i__) = rNULL__
        if(ar__1(134:137) .EQ. '') E_C32_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v_BFK2008__1(i__),l_srcB(i__),srcB(i__),e_srcB(i__),
     +  srcS1(i__),e_srcS1(i__),srcS2(i__),e_srcS2(i__),srcS(i__),
     +  e_srcS(i__),srcH(i__),e_srcH(i__),srcSc(i__),e_srcSc(i__),
     +  srcHc(i__),e_srcHc(i__),HR(i__),e_HR(i__),E_HR_1(i__),
     +  C21(i__),e_C21(i__),E_C21_1(i__),C32(i__),e_C32(i__),
     +  E_C32_1(i__),l_logLX_1(i__),logLX_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Source counts, hardness ratios, color-color
*                             values and variability: Observation 1

C  Format for file interpretation

    3 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A1,1X,A1,1X,F4.1,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,132
        read(1,'(A150)')ar__2
        read(ar__2,3)
     +  v_BFK2008__2(i__),l_srcB_1(i__),srcB_1(i__),e_srcB_1(i__),
     +  srcS1_1(i__),e_srcS1_1(i__),srcS2_1(i__),e_srcS2_1(i__),
     +  srcS_1(i__),e_srcS_1(i__),srcH_1(i__),e_srcH_1(i__),
     +  srcSc_1(i__),e_srcSc_1(i__),srcHc_1(i__),e_srcHc_1(i__),
     +  HR_1(i__),e_HR_2(i__),E_HR_3(i__),C21_1(i__),e_C21_2(i__),
     +  E_C21_3(i__),C32_1(i__),e_C32_2(i__),E_C32_3(i__),BB_1(i__),
     +  K_S_1(i__),Sig(i__),l_logLX_2(i__),logLX_2(i__)
        if(ar__2(14:17) .EQ. '') e_srcB_1(i__) = rNULL__
        if(ar__2(19:23) .EQ. '') srcS1_1(i__) = rNULL__
        if(ar__2(25:28) .EQ. '') e_srcS1_1(i__) = rNULL__
        if(ar__2(30:35) .EQ. '') srcS2_1(i__) = rNULL__
        if(ar__2(37:40) .EQ. '') e_srcS2_1(i__) = rNULL__
        if(ar__2(42:47) .EQ. '') srcS_1(i__) = rNULL__
        if(ar__2(49:52) .EQ. '') e_srcS_1(i__) = rNULL__
        if(ar__2(54:58) .EQ. '') srcH_1(i__) = rNULL__
        if(ar__2(60:63) .EQ. '') e_srcH_1(i__) = rNULL__
        if(ar__2(65:70) .EQ. '') srcSc_1(i__) = rNULL__
        if(ar__2(72:75) .EQ. '') e_srcSc_1(i__) = rNULL__
        if(ar__2(77:81) .EQ. '') srcHc_1(i__) = rNULL__
        if(ar__2(83:86) .EQ. '') e_srcHc_1(i__) = rNULL__
        if(ar__2(88:92) .EQ. '') HR_1(i__) = rNULL__
        if(ar__2(94:97) .EQ. '') e_HR_2(i__) = rNULL__
        if(ar__2(99:102) .EQ. '') E_HR_3(i__) = rNULL__
        if(ar__2(104:108) .EQ. '') C21_1(i__) = rNULL__
        if(ar__2(110:113) .EQ. '') e_C21_2(i__) = rNULL__
        if(ar__2(115:118) .EQ. '') E_C21_3(i__) = rNULL__
        if(ar__2(120:124) .EQ. '') C32_1(i__) = rNULL__
        if(ar__2(126:129) .EQ. '') e_C32_2(i__) = rNULL__
        if(ar__2(131:134) .EQ. '') E_C32_3(i__) = rNULL__
        if(ar__2(140:143) .EQ. '') Sig(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  v_BFK2008__2(i__),l_srcB_1(i__),srcB_1(i__),e_srcB_1(i__),
     +  srcS1_1(i__),e_srcS1_1(i__),srcS2_1(i__),e_srcS2_1(i__),
     +  srcS_1(i__),e_srcS_1(i__),srcH_1(i__),e_srcH_1(i__),
     +  srcSc_1(i__),e_srcSc_1(i__),srcHc_1(i__),e_srcHc_1(i__),
     +  HR_1(i__),e_HR_2(i__),E_HR_3(i__),C21_1(i__),e_C21_2(i__),
     +  E_C21_3(i__),C32_1(i__),e_C32_2(i__),E_C32_3(i__),BB_1(i__),
     +  K_S_1(i__),Sig(i__),l_logLX_2(i__),logLX_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Source counts, hardness ratios, color-color
*                             values and variability: Observation 2

C  Format for file interpretation

    4 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A1,1X,A1,1X,F4.1,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,132
        read(1,'(A150)')ar__3
        read(ar__3,4)
     +  v_BFK2008__3(i__),l_srcB_2(i__),srcB_2(i__),e_srcB_2(i__),
     +  srcS1_2(i__),e_srcS1_2(i__),srcS2_2(i__),e_srcS2_2(i__),
     +  srcS_2(i__),e_srcS_2(i__),srcH_2(i__),e_srcH_2(i__),
     +  srcSc_2(i__),e_srcSc_2(i__),srcHc_2(i__),e_srcHc_2(i__),
     +  HR_2(i__),e_HR_4(i__),E_HR_5(i__),C21_2(i__),e_C21_4(i__),
     +  E_C21_5(i__),C32_2(i__),e_C32_4(i__),E_C32_5(i__),BB_2(i__),
     +  K_S_2(i__),Sig_1(i__),l_logLX_3(i__),logLX_3(i__)
        if(ar__3(14:17) .EQ. '') e_srcB_2(i__) = rNULL__
        if(ar__3(19:23) .EQ. '') srcS1_2(i__) = rNULL__
        if(ar__3(25:28) .EQ. '') e_srcS1_2(i__) = rNULL__
        if(ar__3(30:35) .EQ. '') srcS2_2(i__) = rNULL__
        if(ar__3(37:40) .EQ. '') e_srcS2_2(i__) = rNULL__
        if(ar__3(42:47) .EQ. '') srcS_2(i__) = rNULL__
        if(ar__3(49:52) .EQ. '') e_srcS_2(i__) = rNULL__
        if(ar__3(54:58) .EQ. '') srcH_2(i__) = rNULL__
        if(ar__3(60:63) .EQ. '') e_srcH_2(i__) = rNULL__
        if(ar__3(65:70) .EQ. '') srcSc_2(i__) = rNULL__
        if(ar__3(72:75) .EQ. '') e_srcSc_2(i__) = rNULL__
        if(ar__3(77:81) .EQ. '') srcHc_2(i__) = rNULL__
        if(ar__3(83:86) .EQ. '') e_srcHc_2(i__) = rNULL__
        if(ar__3(88:92) .EQ. '') HR_2(i__) = rNULL__
        if(ar__3(94:97) .EQ. '') e_HR_4(i__) = rNULL__
        if(ar__3(99:102) .EQ. '') E_HR_5(i__) = rNULL__
        if(ar__3(104:108) .EQ. '') C21_2(i__) = rNULL__
        if(ar__3(110:113) .EQ. '') e_C21_4(i__) = rNULL__
        if(ar__3(115:118) .EQ. '') E_C21_5(i__) = rNULL__
        if(ar__3(120:124) .EQ. '') C32_2(i__) = rNULL__
        if(ar__3(126:129) .EQ. '') e_C32_4(i__) = rNULL__
        if(ar__3(131:134) .EQ. '') E_C32_5(i__) = rNULL__
        if(ar__3(140:143) .EQ. '') Sig_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v_BFK2008__3(i__),l_srcB_2(i__),srcB_2(i__),e_srcB_2(i__),
     +  srcS1_2(i__),e_srcS1_2(i__),srcS2_2(i__),e_srcS2_2(i__),
     +  srcS_2(i__),e_srcS_2(i__),srcH_2(i__),e_srcH_2(i__),
     +  srcSc_2(i__),e_srcSc_2(i__),srcHc_2(i__),e_srcHc_2(i__),
     +  HR_2(i__),e_HR_4(i__),E_HR_5(i__),C21_2(i__),e_C21_4(i__),
     +  E_C21_5(i__),C32_2(i__),e_C32_4(i__),E_C32_5(i__),BB_2(i__),
     +  K_S_2(i__),Sig_1(i__),l_logLX_3(i__),logLX_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Source counts, hardness ratios, color-color
*                             values and variability: Observation 3

C  Format for file interpretation

    5 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A1,1X,A1,1X,F4.1,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,132
        read(1,'(A150)')ar__4
        read(ar__4,5)
     +  v_BFK2008__4(i__),l_srcB_3(i__),srcB_3(i__),e_srcB_3(i__),
     +  srcS1_3(i__),e_srcS1_3(i__),srcS2_3(i__),e_srcS2_3(i__),
     +  srcS_3(i__),e_srcS_3(i__),srcH_3(i__),e_srcH_3(i__),
     +  srcSc_3(i__),e_srcSc_3(i__),srcHc_3(i__),e_srcHc_3(i__),
     +  HR_3(i__),e_HR_6(i__),E_HR_7(i__),C21_3(i__),e_C21_6(i__),
     +  E_C21_7(i__),C32_3(i__),e_C32_6(i__),E_C32_7(i__),BB_3(i__),
     +  K_S_3(i__),Sig_2(i__),l_logLX_4(i__),logLX_4(i__)
        if(ar__4(14:17) .EQ. '') e_srcB_3(i__) = rNULL__
        if(ar__4(19:23) .EQ. '') srcS1_3(i__) = rNULL__
        if(ar__4(25:28) .EQ. '') e_srcS1_3(i__) = rNULL__
        if(ar__4(30:35) .EQ. '') srcS2_3(i__) = rNULL__
        if(ar__4(37:40) .EQ. '') e_srcS2_3(i__) = rNULL__
        if(ar__4(42:47) .EQ. '') srcS_3(i__) = rNULL__
        if(ar__4(49:52) .EQ. '') e_srcS_3(i__) = rNULL__
        if(ar__4(54:58) .EQ. '') srcH_3(i__) = rNULL__
        if(ar__4(60:63) .EQ. '') e_srcH_3(i__) = rNULL__
        if(ar__4(65:70) .EQ. '') srcSc_3(i__) = rNULL__
        if(ar__4(72:75) .EQ. '') e_srcSc_3(i__) = rNULL__
        if(ar__4(77:81) .EQ. '') srcHc_3(i__) = rNULL__
        if(ar__4(83:86) .EQ. '') e_srcHc_3(i__) = rNULL__
        if(ar__4(88:92) .EQ. '') HR_3(i__) = rNULL__
        if(ar__4(94:97) .EQ. '') e_HR_6(i__) = rNULL__
        if(ar__4(99:102) .EQ. '') E_HR_7(i__) = rNULL__
        if(ar__4(104:108) .EQ. '') C21_3(i__) = rNULL__
        if(ar__4(110:113) .EQ. '') e_C21_6(i__) = rNULL__
        if(ar__4(115:118) .EQ. '') E_C21_7(i__) = rNULL__
        if(ar__4(120:124) .EQ. '') C32_3(i__) = rNULL__
        if(ar__4(126:129) .EQ. '') e_C32_6(i__) = rNULL__
        if(ar__4(131:134) .EQ. '') E_C32_7(i__) = rNULL__
        if(ar__4(140:143) .EQ. '') Sig_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v_BFK2008__4(i__),l_srcB_3(i__),srcB_3(i__),e_srcB_3(i__),
     +  srcS1_3(i__),e_srcS1_3(i__),srcS2_3(i__),e_srcS2_3(i__),
     +  srcS_3(i__),e_srcS_3(i__),srcH_3(i__),e_srcH_3(i__),
     +  srcSc_3(i__),e_srcSc_3(i__),srcHc_3(i__),e_srcHc_3(i__),
     +  HR_3(i__),e_HR_6(i__),E_HR_7(i__),C21_3(i__),e_C21_6(i__),
     +  E_C21_7(i__),C32_3(i__),e_C32_6(i__),E_C32_7(i__),BB_3(i__),
     +  K_S_3(i__),Sig_2(i__),l_logLX_4(i__),logLX_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Source counts, hardness ratios, color-color
*                             values and variability: Observation 4

C  Format for file interpretation

    6 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A1,1X,A1,1X,F4.1,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,132
        read(1,'(A150)')ar__5
        read(ar__5,6)
     +  v_BFK2008__5(i__),l_srcB_4(i__),srcB_4(i__),e_srcB_4(i__),
     +  srcS1_4(i__),e_srcS1_4(i__),srcS2_4(i__),e_srcS2_4(i__),
     +  srcS_4(i__),e_srcS_4(i__),srcH_4(i__),e_srcH_4(i__),
     +  srcSc_4(i__),e_srcSc_4(i__),srcHc_4(i__),e_srcHc_4(i__),
     +  HR_4(i__),e_HR_8(i__),E_HR_9(i__),C21_4(i__),e_C21_8(i__),
     +  E_C21_9(i__),C32_4(i__),e_C32_8(i__),E_C32_9(i__),BB_4(i__),
     +  K_S_4(i__),Sig_3(i__),l_logLX_5(i__),logLX_5(i__)
        if(ar__5(14:17) .EQ. '') e_srcB_4(i__) = rNULL__
        if(ar__5(19:23) .EQ. '') srcS1_4(i__) = rNULL__
        if(ar__5(25:28) .EQ. '') e_srcS1_4(i__) = rNULL__
        if(ar__5(30:35) .EQ. '') srcS2_4(i__) = rNULL__
        if(ar__5(37:40) .EQ. '') e_srcS2_4(i__) = rNULL__
        if(ar__5(42:47) .EQ. '') srcS_4(i__) = rNULL__
        if(ar__5(49:52) .EQ. '') e_srcS_4(i__) = rNULL__
        if(ar__5(54:58) .EQ. '') srcH_4(i__) = rNULL__
        if(ar__5(60:63) .EQ. '') e_srcH_4(i__) = rNULL__
        if(ar__5(65:70) .EQ. '') srcSc_4(i__) = rNULL__
        if(ar__5(72:75) .EQ. '') e_srcSc_4(i__) = rNULL__
        if(ar__5(77:81) .EQ. '') srcHc_4(i__) = rNULL__
        if(ar__5(83:86) .EQ. '') e_srcHc_4(i__) = rNULL__
        if(ar__5(88:92) .EQ. '') HR_4(i__) = rNULL__
        if(ar__5(94:97) .EQ. '') e_HR_8(i__) = rNULL__
        if(ar__5(99:102) .EQ. '') E_HR_9(i__) = rNULL__
        if(ar__5(104:108) .EQ. '') C21_4(i__) = rNULL__
        if(ar__5(110:113) .EQ. '') e_C21_8(i__) = rNULL__
        if(ar__5(115:118) .EQ. '') E_C21_9(i__) = rNULL__
        if(ar__5(120:124) .EQ. '') C32_4(i__) = rNULL__
        if(ar__5(126:129) .EQ. '') e_C32_8(i__) = rNULL__
        if(ar__5(131:134) .EQ. '') E_C32_9(i__) = rNULL__
        if(ar__5(140:143) .EQ. '') Sig_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  v_BFK2008__5(i__),l_srcB_4(i__),srcB_4(i__),e_srcB_4(i__),
     +  srcS1_4(i__),e_srcS1_4(i__),srcS2_4(i__),e_srcS2_4(i__),
     +  srcS_4(i__),e_srcS_4(i__),srcH_4(i__),e_srcH_4(i__),
     +  srcSc_4(i__),e_srcSc_4(i__),srcHc_4(i__),e_srcHc_4(i__),
     +  HR_4(i__),e_HR_8(i__),E_HR_9(i__),C21_4(i__),e_C21_8(i__),
     +  E_C21_9(i__),C32_4(i__),e_C32_8(i__),E_C32_9(i__),BB_4(i__),
     +  K_S_4(i__),Sig_3(i__),l_logLX_5(i__),logLX_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Source counts, hardness ratios, color-color
*                             values and variability: Observation 5

C  Format for file interpretation

    7 format(
     +  I3,1X,A2,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,
     +  1X,F4.1,1X,F5.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F4.2,1X,A1,1X,A1,1X,F4.1,1X,A2,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,132
        read(1,'(A150)')ar__6
        read(ar__6,7)
     +  v_BFK2008__6(i__),l_srcB_5(i__),srcB_5(i__),e_srcB_5(i__),
     +  srcS1_5(i__),e_srcS1_5(i__),srcS2_5(i__),e_srcS2_5(i__),
     +  srcS_5(i__),e_srcS_5(i__),srcH_5(i__),e_srcH_5(i__),
     +  srcSc_5(i__),e_srcSc_5(i__),srcHc_5(i__),e_srcHc_5(i__),
     +  HR_5(i__),e_HR_10(i__),E_HR_11(i__),C21_5(i__),e_C21_10(i__),
     +  E_C21_11(i__),C32_5(i__),e_C32_10(i__),E_C32_11(i__),
     +  BB_5(i__),K_S_5(i__),Sig_4(i__),l_logLX_6(i__),logLX_6(i__)
        if(ar__6(14:17) .EQ. '') e_srcB_5(i__) = rNULL__
        if(ar__6(19:23) .EQ. '') srcS1_5(i__) = rNULL__
        if(ar__6(25:28) .EQ. '') e_srcS1_5(i__) = rNULL__
        if(ar__6(30:35) .EQ. '') srcS2_5(i__) = rNULL__
        if(ar__6(37:40) .EQ. '') e_srcS2_5(i__) = rNULL__
        if(ar__6(42:47) .EQ. '') srcS_5(i__) = rNULL__
        if(ar__6(49:52) .EQ. '') e_srcS_5(i__) = rNULL__
        if(ar__6(54:58) .EQ. '') srcH_5(i__) = rNULL__
        if(ar__6(60:63) .EQ. '') e_srcH_5(i__) = rNULL__
        if(ar__6(65:70) .EQ. '') srcSc_5(i__) = rNULL__
        if(ar__6(72:75) .EQ. '') e_srcSc_5(i__) = rNULL__
        if(ar__6(77:81) .EQ. '') srcHc_5(i__) = rNULL__
        if(ar__6(83:86) .EQ. '') e_srcHc_5(i__) = rNULL__
        if(ar__6(88:92) .EQ. '') HR_5(i__) = rNULL__
        if(ar__6(94:97) .EQ. '') e_HR_10(i__) = rNULL__
        if(ar__6(99:102) .EQ. '') E_HR_11(i__) = rNULL__
        if(ar__6(104:108) .EQ. '') C21_5(i__) = rNULL__
        if(ar__6(110:113) .EQ. '') e_C21_10(i__) = rNULL__
        if(ar__6(115:118) .EQ. '') E_C21_11(i__) = rNULL__
        if(ar__6(120:124) .EQ. '') C32_5(i__) = rNULL__
        if(ar__6(126:129) .EQ. '') e_C32_10(i__) = rNULL__
        if(ar__6(131:134) .EQ. '') E_C32_11(i__) = rNULL__
        if(ar__6(140:143) .EQ. '') Sig_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  v_BFK2008__6(i__),l_srcB_5(i__),srcB_5(i__),e_srcB_5(i__),
     +  srcS1_5(i__),e_srcS1_5(i__),srcS2_5(i__),e_srcS2_5(i__),
     +  srcS_5(i__),e_srcS_5(i__),srcH_5(i__),e_srcH_5(i__),
     +  srcSc_5(i__),e_srcSc_5(i__),srcHc_5(i__),e_srcHc_5(i__),
     +  HR_5(i__),e_HR_10(i__),E_HR_11(i__),C21_5(i__),e_C21_10(i__),
     +  E_C21_11(i__),C32_5(i__),e_C32_10(i__),E_C32_11(i__),
     +  BB_5(i__),K_S_5(i__),Sig_4(i__),l_logLX_6(i__),logLX_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table15.dat'	! Raw source and background counts from the
*                             co-added observation

C  Format for file interpretation

    8 format(
     +  I3,1X,F6.1,1X,I4,1X,F4.1,1X,I4,1X,F4.1,1X,I4,1X,F4.1,1X,I4,1X,
     +  F4.1,1X,I4,1X,F4.1,1X,I4,1X,F4.1,1X,I4,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table15.dat')
      write(6,*) '....Loading file: table15.dat'
      do i__=1,132
        read(1,'(A80)')ar__7
        read(ar__7,8)
     +  v_BFK2008__7(i__),netB(i__),rawB(i__),bkgB(i__),rawS1(i__),
     +  bkgS1(i__),rawS2(i__),bkgS2(i__),rawH(i__),bkgH(i__),
     +  rawBc(i__),bkgBc(i__),rawSc(i__),bkgSc(i__),rawHc(i__),
     +  bkgHc(i__)
        if(ar__7(5:10) .EQ. '') netB(i__) = rNULL__
        if(ar__7(12:15) .EQ. '') rawB(i__) = iNULL__
        if(ar__7(17:20) .EQ. '') bkgB(i__) = rNULL__
        if(ar__7(22:25) .EQ. '') rawS1(i__) = iNULL__
        if(ar__7(27:30) .EQ. '') bkgS1(i__) = rNULL__
        if(ar__7(32:35) .EQ. '') rawS2(i__) = iNULL__
        if(ar__7(37:40) .EQ. '') bkgS2(i__) = rNULL__
        if(ar__7(42:45) .EQ. '') rawH(i__) = iNULL__
        if(ar__7(47:50) .EQ. '') bkgH(i__) = rNULL__
        if(ar__7(52:55) .EQ. '') rawBc(i__) = iNULL__
        if(ar__7(57:60) .EQ. '') bkgBc(i__) = rNULL__
        if(ar__7(62:65) .EQ. '') rawSc(i__) = iNULL__
        if(ar__7(67:70) .EQ. '') bkgSc(i__) = rNULL__
        if(ar__7(72:75) .EQ. '') rawHc(i__) = iNULL__
        if(ar__7(77:80) .EQ. '') bkgHc(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  v_BFK2008__7(i__),netB(i__),rawB(i__),bkgB(i__),rawS1(i__),
     +  bkgS1(i__),rawS2(i__),bkgS2(i__),rawH(i__),bkgH(i__),
     +  rawBc(i__),bkgBc(i__),rawSc(i__),bkgSc(i__),rawHc(i__),
     +  bkgHc(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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