FORTRAN Generation
(/./ftp/cats/J/A_A/661/A5)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/661/A5 into FORTRAN code for reading data files line by line.

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-May-20
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/A+A/661/A5               The eFEDS AGN catalog                    (Liu+, 2022)
*================================================================================
*The eROSITA Final Equatorial-Depth Survey (eFEDS).
*The AGN catalog and its X-ray spectral properties.
*    Liu T., Buchner J., Nandra K., Merloni A., Dwelly T., Sanders J.S.,
*    Salvato M., Arcodia R., Brusa M., Wolf J., Georgakakis A., Boller T.,
*    Krumpe M., Lamer G., Waddell S., Urrutia T., Schwope A., Robrade J.,
*    Wilms J., Dauser T., Comparat J., Toba Y., Ichikawa K., Iwasawa K.,
*    Shen Y., Medel H.I.
*   <Astron. Astrophys. 661, A5 (2022)>
*   =2022A&A...661A...5L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! The eFEDS AGN catalog

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

C  J2000 position composed of: RAcdeg DEcdeg
      integer*4     ID_SRC      ! [1/32683] ID of the sources in the eFEDS
*                                         main X-ray catalog (Paper I) (ID_SRC)
      real*8        RAcdeg      ! (deg) [126.24/145.59] Right ascension (J2000),
*                                         astrometric corrected (Paper I)
*                                         (RA_CORR)
      real*4        DEcdeg      ! (deg) [-2.91/5.97] Declination (J2000),
*                                         astrometric corrected (Paper I)
*                                         (DEC_CORR)
      real*8        DetLike     ! 0.2-2.3 keV source detection likelihood
*                                         (Paper I) (DET_LIKE)
      real*4        inArea90    ! Whether located inside the inner
*                                         90%-area region of eFEDS (Paper I)
*                                         (inArea90)
      character*11  CTPLS8      ! ID of the best LS8 counterpart
*                                         (Paper II) (CTP_LS8_UNIQUE_OBJID)
      real*8        RALdeg      ! (deg) [126.23/145.59] Right ascension (J2000)
*                                         of the best LS8 counterpart (Paper II)
*                                         (CTP_LS8_RA)
      real*4        DELdeg      ! (deg) [-2.91/5.97] Declination (J2000) of the
*                                         best LS8 counterpart (Paper II)
*                                         (CTP_LS8_DEC)
      integer*4     CTPqual     ! [2/4] Counterpart quality (Paper II)
*                                         A value >=2 is recommended
*                                         (CTP_quality)
      integer*4     CTPClass    ! [2/3] Classification of the optical
*                                         counterpart (Paper II) (CTP_CLASS) (1)
      real*8        CTPz        ! [0.0/8.0] Redshift of the optical
*                                         counterpart (Paper II) (CTP_REDSHIFT)
      integer*4     CTPzgrade   ! [2/5] Redshift Grade (Paper II)
*                                         A value >=3 is recommended
*                                         (CTP_REDSHIFT_GRADE)
      integer*4     inKiDSFlag  ! [0/1]  Whether located inside the region
*                                         of the KiDS survey (Paper II)
*                                         (in_KiDS_flag)
      integer*4     LxModel     ! [1/5] Index of selected model for
*                                         luminosity measurement (LxModel) (2)
      integer*4     NHclass     ! [1/4] Class of AGN NH measurement with
*                                         model 1 (single-powerlaw) (NHclass) (3)
      integer*4     FSModel     ! [3/6] Index of selected model for
*                                         0.5-2keV flux (FSModel) (4)
      integer*4     FHModel     ! [5/7] Index of selected model for
*                                         2.3-5keV flux (FHModel) (5)
      real*4        galNH       ! (cm-2) Total column density of Galactic
*                                         absorption (galNH)
      real*4        galNHI      ! (cm-2) HI column density from HI4PI (galNHI)
      real*8        SrcCts      ! (ct) Source net counts in the 0.2-5 keV band
*                                         (SrcCts)
      real*4        FluxcMedS   ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, posterior median
*                                         (FluxCorr_Med_s)
      real*4        b_FluxcS    ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, 1sigma lower limit
*                                         (FluxCorr_Lo1_s)
      real*4        B_FluxcS_1  ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, 1sigma upper limit
*                                         (FluxCorr_Up1_s)
      real*4        FluxcMedT   ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, posterior median
*                                         (FluxCorr_Med_t)
      real*4        b_FluxcT    ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, 1sigma lower limit
*                                         (FluxCorr_Lo1_t)
      real*4        B_FluxcT_1  ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, 1sigma upper limit
*                                         (FluxCorr_Up1_t)
      real*4        FluxIMedS   ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2 keV, posterior median
*                                         (FluxIntr_Med_s)
      real*4        b_FluxIS    ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2keV, 1sigma lower limit
*                                         (FluxIntr_Lo1_s)
      real*4        B_FluxIS_1  ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2keV, 1sigma upper limit
*                                         (FluxIntr_Up1_s)
      real*4        FluxIMed2keV ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, median (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, 1sigma lower limit
*                                         (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_1 ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, 1sigma upper limit
*                                         (FluxIntr_Up1_2keV)
      real*8        LumIMedS    ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, posterior median
*                                         (LumiIntr_Med_s)
      real*8        b_LumIS     ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, 1sigma lower limit
*                                         (LumiIntr_Lo1_s)
      real*8        B_LumIS_1   ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, 1sigma upper limit
*                                         (LumiIntr_Up1_s)
      real*8        LumIMed2keV ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         posterior median in erg/s/eV
*                                         (LumiIntr_Med_2keV)
      real*8        b_LumI2keV  ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         1sigma lower limit in erg/s/eV
*                                         (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_1 ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         1sigma upper limit in erg/s/eV
*                                         (LumiIntr_Up1_2keV)
      real*4        FluxObsMedS ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS  ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_1 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         1sigma upper limit (FluxObsv_Up1_s)
      real*4        FluxObsMedT ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT  ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_1 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         1sigma upper limit (FluxObsv_Up1_t)
      real*8        lognHKLm1   ! log AGN column density, KL divergence
*                                         (in nats units) (lognH_KL_m1)
      real*8        b_lognHm1   ! ([cm-2]) log AGN column density, HDI lower limit
*                                         (lognH_HLo_m1)
      real*8        B_lognHm1_1 ! ([cm-2]) log AGN column density, HDI upper limit
*                                         (lognH_HUp_m1)
      real*8        lognHMedm1  ! ([cm-2]) log AGN column density, posterior median
*                                         (lognH_Med_m1)
      real*8        GammaKLm3   ! primary power-law slope in model 3,
*                                         KL divergence (in nats units)
*                                         (Gamma_KL_m3)
      real*8        b_GammaHm3  ! primary power-law slope in model 3,
*                                         HDI lower limit (Gamma_HLo_m3)
      real*8        B_GammaHm3_1 ! primary power-law slope in model 3,
*                                         HDI upper limit (Gamma_HUp_m3)
      real*8        GammaMedm3  ! primary power-law slope in model 3,
*                                         posterior median (Gamma_Med_m3)
      real*8        logZm0      ! ([-]) log Bayesian evidence with model 0:
*                                         APEC (logZ_m0)
      real*8        logZm1      ! ([-]) log Bayesian evidence with model 1:
*                                         single powerlaw (logZ_m1)
      real*8        logZm2      ! ([-]) log Bayesian evidence with model 2:
*                                         double powerlaw (logZ_m2)
      real*8        logZm3      ! ([-]) log Bayesian evidence with model 3:
*                                         powerlaw + blackbody (logZ_m3)
      real*8        logZm4      ! ([-]) log Bayesian evidence with model 4:
*                                         powerlaw with Gamma fixed at 2.0
*                                         (logZ_m4)
      real*8        logZm5      ! ([-]) log Bayesian evidence with model 5:
*                                         shape-fixed powerlaw (logZ_m5)
      real*8        L2500       ! (10-7W/Hz) ?=-99.0 The rest-frame 2500 luminosity
*                                          in erg/s/Hz (L2500)
      real*8        L5100       ! (10-7W/Hz) ?=-99.0 The rest-frame 5100 luminosity
*                                         in erg/s/Hz (L5100)
      real*8        W1mag       ! (mag) ? LS8-WISE W1 AB magnitude (Paper II)
*                                         (W1)
      real*4        e_W1mag     ! (mag) ? LS8-WISE W1 magnitude error (Paper II)
*                                         (W1_ERR)
      real*8        W2mag       ! (mag) ? LS8-WISE W2 AB magnitude (Paper II)
*                                         (W2)
      real*4        e_W2mag     ! (mag) ? LS8-WISE W2 magnitude error (Paper II)
*                                         (W2_ERR)
*Note (1): Classification of the optical counterpart as follows:
*           0 = likely Galactic
*           1 = secure Galactic
*           2 = likely extraGalactic
*           3 = secure extraGalactic
*Note (2): Index of selected model for luminosity measurement as follows:
*           1 = single-powerlaw
*           3 = powerlaw+bb
*           4 = powerlaw with Gamma
*Note (3): Class of AGN NH measurement with model 1 (single-powerlaw) as follows:
*           1 = uninformative
*           2 = unobscured
*           3 = mildly-measured
*           4 = well-measured
*Note (4): Index of selected model for 0.5-2keV flux as follows:
*           5 = counts-based measurement for faint sources
*           0 = robust spectra-based measurement
*           3 = robust spectra-based measurement
*           6 = robust spectra-based measurement
*Note (5): Index of selected model for 2.3-5keV flux as follows:
*           5 = counts-based measurement for faint sources
*           7 = robust spectra-based measurement

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

C  Declarations for 'table2.dat'	! Spectral properties

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID_SRC_1    ! [1/32683] ID of the sources in the eFEDS
*                                       main X-ray catalog (Paper I) (ID_SRC)
      real*8        RAcdeg_1    ! (deg) [126.24/145.66] Right ascension (J2000),
*                                       astrometric corrected (Paper I) (RA_CORR)
      real*4        DEcdeg_1    ! (deg) [-2.91/5.97] Declination (J2000),
*                                       astrometric corrected (Paper I)
*                                       (DEC_CORR)
      real*8        DetLike_1   ! [6.0/33254.59] 0.2-2.3 keV source
*                                       detection likelihood (Paper I) (DET_LIKE)
      real*8        ExtLike     ! X-ray source extent likelihood (Paper I)
*                                       (EXT_LIKE)
      real*4        inArea90_1  ! Whether located inside the inner 90%-area
*                                       region of eFEDS (Paper I) (inArea90)
      real*4        galNH_1     ! (cm-2) Total column density of Galactic
*                                       absorption (galNH)
      real*4        galNHI_1    ! (cm-2) HI column density from HI4PI (galNHI)
      real*8        ExpTime     ! (s) Spectra exposure time (Exposure)
      real*8        SrcCts_1    ! (ct) Source net counts in the 0.2-5 keV band
*                                       (SrcCts)
      real*8        RAdeg       ! (deg) [126.24/145.66] Right ascension (J2000)
*                                       before astrometric correction (Paper I)
*                                       (RA)
      real*4        DEdeg       ! (deg) [-2.91/5.97] Declination (J2000) before
*                                       astrometric correction (Paper I) (DEC)
      integer*4     Rad         ! (arcsec) [12/259] Source extraction radius (Rad)
      integer*4     Ann1        ! (arcsec) [22/515] Inner radius of background
*                                       extraction region (Ann1)
      integer*4     Ann2        ! (arcsec) [171/900] Outer radius of background
*                                       extraction region (Ann2)
      real*4        RegArea     ! (deg2) Geometry area of source extraction region
*                                       (REGAREA_s)
      real*4        Backscal    ! (deg2) Source BACKSCAL (Backscal_s)
      real*8        BackscalBkg ! (deg2) Background BACKSCAL (Backscal_b)
      integer*4     FSModel_1   ! [0/6] Index of selected model for
*                                       0.5-2keV flux (FSModel) (1)
      integer*4     FHModel_1   ! [5/7] Index of selected model for
*                                       2.3-5keV flux (FHModel) (2)
      real*4        Rate0_2_2_3 ! (ct/s) Net count rate in 0.2-2.3 keV
*                                       (Rate_d2_2d3)
      real*8        e_Rate0_2_2_3 ! (ct/s) Net count rate error in 0.2-2.3 keV
*                                       (RateErr_d2_2d3)
      real*4        Rate0_2_0_5 ! (ct/s) Net count rate in 0.2-0.5 keV (Rate_d2_d5)
      real*4        e_Rate0_2_0_5 ! (ct/s) Net count rate error in 0.2-0.5 keV
*                                       (RateErr_d2_d5)
      real*4        Rate0_5_1   ! (ct/s) Net count rate in 0.5-1 keV (Rate_d5_1)
      real*4        e_Rate0_5_1 ! (ct/s) Net count rate error in 0.5-1 keV
*                                       (RateErr_d5_1)
      real*4        Rate1_2     ! (ct/s) Net count rate in 1-2 keV (Rate_1_2)
      real*4        e_Rate1_2   ! (ct/s) Net count rate error in 1-2 keV
*                                       (RateErr_1_2)
      real*4        Rate2_4_5   ! (ct/s) Net count rate in 2-4.5 keV (Rate_2_4d5)
      real*4        e_Rate2_4_5 ! (ct/s) Net count rate error in 2-4.5 keV
*                                       (RateErr_2_4d5)
      real*4        Rate2_3_5   ! (ct/s) Net count rate in 2.3-5 keV (Rate_2d3_5)
      real*4        e_Rate2_3_5 ! (ct/s) Net count rate error in 2.3-5 keV
*                                       (RateErr_2d3_5)
      real*4        Rate5_8     ! (ct/s) Net count rate in 5-8 keV (Rate_5_8)
      real*4        e_Rate5_8   ! (ct/s) Net count rate error in 5-8 keV
*                                       (RateErr_5_8)
      real*4        BkgCts0_2_0_6 ! (ct) Background counts in 0.2-0.6 keV
*                                       (BkgCts_d2_d6)
      real*8        e_BkgCts0_2_0_6 ! (ct) Background counts error in 0.2-0.6 keV
*                                       (BkgCtsErr_d2_d6)
      real*4        BkgCts0_6_2_3 ! (ct) Background counts in 0.6-2.3 keV
*                                       (BkgCts_d6_2d3)
      real*8        e_BkgCts0_6_2_3 ! (ct) Background counts error in 0.6-2.3 keV
*                                       (BkgCtsErr_d6_2d3)
      real*4        BkgCts2_3_5 ! (ct) Background counts in 2.3-5 keV
*                                       (BkgCts_2d3_5)
      real*8        e_BkgCts2_3_5 ! (ct) Background counts error in 2.3-5 keV
*                                       (BkgCtsErr_2d3_5)
      real*4        BkgCts5_8   ! (ct) Background counts in 5-8 keV (BkgCts_5_8)
      real*8        e_BkgCts5_8 ! (ct) Background counts error in 5-8 keV
*                                       (BkgCtsErr_5_8)
      real*8        PSFCor0_2_0_5 ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 0.2-0.5 keV (PSFCor_d2_d5)
      real*8        PSFCor0_5_1 ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 0.5-1 keV (PSFCor_d5_1)
      real*8        PSFCor1_2   ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 1-2 keV (PSFCor_1_2)
      real*8        PSFCor2_3_5 ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 2.3-5 keV (PSFCor_2d3_5)
      integer*4     Nempty      ! [35/771] Number of empty channels between
*                                       channel 20 and 900 (Nempty)
      real*4        FluxObsMedS_1 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_2 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_3 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       1sigma upper limit (FluxObsv_Up1_s)
      real*4        FluxObsMedT_1 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                       posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_2 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                       1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_3 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                       1sigma upper limit (FluxObsv_Up1_t)
*Note (1): Index of selected model for 0.5-2keV flux as follows:
*           5 = counts-based measurement for faint sources
*           0 = robust spectra-based measurement
*           3 = robust spectra-based measurement
*           6 = robust spectra-based measurement
*Note (2): Index of selected model for 2.3-5keV flux as follows:
*           5 = counts-based measurement for faint sources
*           7 = robust spectra-based measurement

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

C  Declarations for 'table3.dat'	! APEC fitting results

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

      integer*4     ID_SRC_2    ! [1/32683] ID of the sources in the eFEDS
*                                        main X-ray catalog (Paper I) (ID_SRC)
      real*4        z           ! The redshift adopted for the spectral
*                                        model (Redshift)
      real*8        lognHKL     ! log AGN column density, KL divergence
*                                        (in nats units) (lognH_KL)
      real*8        b_lognHH    ! ([cm-2]) log AGN column density, HDI lower limit
*                                        (lognH_HLo)
      real*8        B_lognHH_1  ! ([cm-2]) log AGN column density, HDI upper limit
*                                        (lognH_HUp)
      real*8        logkTKL     ! log temperature, KL divergence
*                                        (in nats units) (logkT_KL)
      real*4        b_logkTH    ! ([keV]) log temperature, HDI lower limit
*                                        (logkT_HLo)
      real*4        B_logkTH_1  ! ([keV]) log temperature, HDI upper limit
*                                        (logkT_HUp)
      real*8        logAbundKL  ! log abundance, KL divergence
*                                        (in nats units) (logAbundanc_KL)
      real*8        b_logAbunH  ! log abundance, HDI lower limit
*                                        (logAbundanc_HLo)
      real*4        B_logAbunH_1 ! log abundance, HDI upper limit
*                                        (logAbundanc_HUp)
      real*8        logApecNormKL ! log APEC normalization, KL divergence
*                                        (in nats units) (logApecNorm_KL)
      real*8        b_logApecNormH ! ([cm-5]) log APEC normalization, HDI lower limit
*                                        (logApecNorm_HLo)
      real*8        B_logApecNormH_1 ! ([cm-5]) log APEC normalization, HDI upper limit
*                                        (logApecNorm_HUp)
      real*8        lognHBF     ! ([cm-2]) log AGN column density, best-fit model
*                                        (lognH_BF)
      real*4        logkTBF     ! ([keV]) log temperature, best-fit model
*                                        (logkT_BF)
      real*4        logAbundBF  ! log abundance, best-fit model
*                                        (logAbundanc_BF)
      real*8        logApecNBF  ! ([cm-5]) log APEC normalization, best-fit model
*                                        (logApecNorm_BF)
      real*4        FluxObsBFS  ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        best-fit model (FluxObsv_BF_s)
      real*4        FluxObsMedS_2 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_4 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_5 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        1sigma upper limit (FluxObsv_Up1_s)
      real*4        b2_FluxObsS ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        2sigma lower limit (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_1 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        2sigma upper limit (FluxObsv_Up2_s)
      real*4        FluxObsBFT  ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_2 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_4 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_5 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        1sigma upper limit (FluxObsv_Up1_t)
      real*4        b2_FluxObsT ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        2sigma lower limit (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_1 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        2sigma upper limit (FluxObsv_Up2_t)
      real*8        lognHMed    ! ([cm-2]) log AGN column density, posterior median
*                                        (lognH_Med)
      real*8        b_lognH     ! ([cm-2]) log AGN column density, percentile
*                                        lower limit (lognH_Lo)
      real*8        B_lognH_1   ! ([cm-2]) log AGN column density, percentile
*                                        upper limit (lognH_Up)
      real*8        lognHMean   ! ([cm-2]) log AGN column density, mean (lognH_Mean)
      real*8        s_lognH     ! ([cm-2]) log AGN column density,
*                                        standard deviation (lognH_Std)
      real*4        logkTMed    ! ([keV]) log temperature, posterior median
*                                        (logkT_Med)
      real*4        b_logkT     ! ([keV]) log temperature, percentile lower limit
*                                        (logkT_Lo)
      real*4        B_logkT_1   ! ([keV]) log temperature, percentile upper limit
*                                        (logkT_Up)
      real*4        logkTMean   ! ([keV]) log temperature, mean (logkT_Mean)
      real*8        s_logkT     ! ([keV]) log temperature, standard deviation
*                                        (logkT_Std)
      real*8        logAbundMed ! log abundance, posterior median
*                                        (logAbundanc_Med)
      real*8        b_logAbund  ! log abundance, percentile lower limit
*                                        (logAbundanc_Lo)
      real*8        B_logAbund_1 ! log abundance, percentile upper limit
*                                        (logAbundanc_Up)
      real*8        logAbundMean ! log abundance, mean (logAbundanc_Mean)
      real*8        s_logAbund  ! log abundance, standard deviation
*                                        (logAbundanc_Std)
      real*8        logApecNMed ! ([cm-5]) log APEC normalization, posterior median
*                                        (logApecNorm_Med)
      real*8        b_logApecN  ! ([cm-5]) log APEC normalization, percentile
*                                        lower limit (logApecNorm_Lo)
      real*8        B_logApecN_1 ! ([cm-5]) log APEC normalization, percentile
*                                        upper limit (logApecNorm_Up)
      real*8        logApecNMean ! ([cm-5]) log APEC normalization, mean
*                                        (logApecNorm_Mean)
      real*8        s_logApecN  ! ([cm-5]) log APEC normalization,
*                                        standard deviation (logApecNorm_Std)
      real*4        logBkgNMed  ! log background normalization,
*                                        posterior median (logBkgNorm_Med)
      real*4        b_logBkgN   ! log background normalization,
*                                        percentile lower limit (logBkgNorm_Lo)
      real*4        B_logBkgN_1 ! log background normalization,
*                                        percentile upper limit (logBkgNorm_Up)
      real*4        logBkgNMean ! log background normalization,
*                                        mean (logBkgNorm_Mean)
      real*8        s_logBkgN   ! log background normalization,
*                                        standard deviation (logBkgNorm_Std)
      real*8        chi25       ! ?=-99 Chi^2 of the best-fit model against
*                                        the rebinned data (>=25 cts/bin) (chi25)
      integer*4     dof25       ! ?=-99 DOF of the rebinned data, which
*                                        has at least 25 counts per bin (dof25)
      real*8        logZ        ! log10 Bayesian evidence (logZ)

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

C  Declarations for 'table4.dat'	! Single-powerlaw results

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

      integer*4     ID_SRC_3    ! [1/32683] ID of the sources in the
*                                             eFEDS main X-ray catalog (Paper I)
*                                             (ID_SRC)
      real*4        z_1         ! [-0.01/8.0] The redshift adopted for
*                                             the spectral model (Redshift)
      real*8        GammaKL     ! power-law slope, KL divergence
*                                             (in nats units) (Gamma_KL)
      real*8        b_GammaH    ! power-law slope, HDI lower limit
*                                             (Gamma_HLo)
      real*8        B_GammaH_1  ! power-law slope, HDI upper limit
*                                             (Gamma_HUp)
      real*8        lognHKL_1   ! log AGN column density,
*                                             KL divergence (in nats units)
*                                             (lognH_KL)
      real*8        b_lognHH_2  ! ([cm-2]) log AGN column density,
*                                             HDI lower limit (lognH_HLo)
      real*8        B_lognHH_3  ! ([cm-2]) log AGN column density,
*                                             HDI upper limit (lognH_HUp)
      real*8        logPowNKL   ! log power-law normalization,
*                                             KL divergence (in nats units)
*                                             (logPowNorm_KL)
      real*8        b_logPowNH  ! (ph/keV/cm2/s) log power-law normalization,
*                                             HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_1 ! (ph/keV/cm2/s) log power-law normalization,
*                                             HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF     ! power-law slope, best-fit model
*                                             (Gamma_BF)
      real*8        lognHBF_1   ! ([cm-2]) log AGN column density,
*                                             best-fit model (lognH_BF)
      real*8        logPowNBF   ! (ph/keV/cm2/s) log power-law normalization,
*                                             best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_1 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             best-fit model (FluxObsv_BF_s)
      real*4        FluxObsMedS_3 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_6 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_7 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             1sigma upper limit (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_2 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             2sigma lower limit (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_3 ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             2sigma upper limit (FluxObsv_Up2_s)
      real*4        FluxObsBFT_1 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_3 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_6 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_7 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             1sigma upper limit (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_2 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             2sigma lower limit (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_3 ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             2sigma upper limit (FluxObsv_Up2_t)
      real*4        FluxcBFS    ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             best-fit model (FluxCorr_BF_s)
      real*4        FluxcMedS_1 ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             posterior median (FluxCorr_Med_s)
      real*4        b_FluxcS_2  ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             1sigma lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxcS_3  ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             1sigma upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxcS   ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             2sigma lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxcS_1 ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             2sigma upper limit (FluxCorr_Up2_s)
      real*4        FluxcBFT    ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV, best-fit model
*                                             (FluxCorr_BF_t)
      real*4        FluxcMedT_1 ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             posterior median (FluxCorr_Med_t)
      real*4        b_FluxcT_2  ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             1sigma lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxcT_3  ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             1sigma upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxcT   ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             2sigma lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxcT_1 ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             2sigma upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS    ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_1 ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             posterior median (FluxIntr_Med_s)
      real*4        b_FluxIS_2  ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2keV,
*                                             1sigma lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_3  ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2keV,
*                                             1sigma upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS   ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_1 ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, best- fit
*                                             (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_1 ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, median
*                                             (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_2 ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 1sigma
*                                             lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_3 ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 1sigma
*                                             upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 2sigma
*                                             lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_1 ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 2sigma
*                                             upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS     ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             best-fit model (LumiIntr_BF_s)
      real*8        LumIMedS_1  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             posterior median (LumiIntr_Med_s)
      real*8        b_LumIS_2   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             1sigma lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_3   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             1sigma upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS    ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_1  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV  ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, best-fit model
*                                             (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_1 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, posterior median
*                                             (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_2 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 1sigma lower limit
*                                             (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_3 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 1sigma upper limit
*                                             (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 2sigma lower limit
*                                             (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_1 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 2sigma upper limit
*                                             (LumiIntr_Up2_2keV)
      real*8        GammaMed    ! power-law slope, posterior median
*                                             (Gamma_Med)
      real*8        b_Gamma     ! power-law slope, percentile
*                                             lower limit (Gamma_Lo)
      real*8        B_Gamma_1   ! power-law slope, percentile
*                                             upper limit (Gamma_Up)
      real*8        GammaMean   ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma     ! power-law slope, standard deviation
*                                             (Gamma_Std)
      real*8        lognHMed_1  ! ([cm-2]) log AGN column density,
*                                             posterior median (lognH_Med)
      real*8        b_lognH_2   ! ([cm-2]) log AGN column density,
*                                             percentile lower limit (lognH_Lo)
      real*8        B_lognH_3   ! ([cm-2]) log AGN column density,
*                                             percentile upper limit (lognH_Up)
      real*8        lognHMean_1 ! ([cm-2]) log AGN column density,
*                                             mean (lognH_Mean)
      real*8        s_lognH_1   ! ([cm-2]) log AGN column density,
*                                             standard deviation (lognH_Std)
      real*8        logPowNMed  ! (ph/keV/cm2/s) log power-law normalization,
*                                             posterior median (logPowNorm_Med)
      real*8        b_logPowN   ! (ph/keV/cm2/s) log power-law normalization,
*                                             percentile lower limit
*                                             (logPowNorm_Lo)
      real*8        B_logPowN_1 ! (ph/keV/cm2/s) log power-law normalization,
*                                             percentile upper limit
*                                             (logPowNorm_Up)
      real*8        logPowNMean ! (ph/keV/cm2/s) log power-law normalization,
*                                             mean (logPowNorm_Mean)
      real*8        s_logPowN   ! (ph/keV/cm2/s) log power-law normalization,
*                                             standard deviation (logPowNorm_Std)
      real*4        logBkgNMed_1 ! log background normalization,
*                                             posterior median (logBkgNorm_Med)
      real*8        b_logBkgN_2 ! log background normalization,
*                                             percentile lower limit
*                                             (logBkgNorm_Lo)
      real*4        B_logBkgN_3 ! log background normalization,
*                                             percentile upper limit
*                                             (logBkgNorm_Up)
      real*4        logBkgNMean_1 ! log background normalization,
*                                             mean (logBkgNorm_Mean)
      real*8        s_logBkgN_1 ! log background normalization,
*                                             standard deviation (logBkgNorm_Std)
      real*8        chi25_1     ! ?=-99 Chi^2 of the best-fit model
*                                             against the rebinned data
*                                             (>=25 cts/bin) (chi25)
      integer*4     dof25_1     ! ?=-99 DOF of the rebinned data,
*                                             which has at least 25 counts per
*                                             bin (dof25)
      real*8        logZ_1      ! log10 Bayesian evidence (logZ)
      integer*4     NHclass_1   ! [1/4] Class of AGN NH measurement
*                                             with model 1 (single-powerlaw)
*                                             (NHclass) (1)
*Note (1): Class of AGN NH measurement with model 1 (single-powerlaw)
*           as follows:
*             1 = uninformative
*             2 = unobscured
*             3 = mildly-measured
*             4 = well-measured

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

C  Declarations for 'table5.dat'	! Double-powerlaw results

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

      integer*4     ID_SRC_4    ! [1/32683] ID of the sources in
*                                                 the eFEDS main X-ray catalog
*                                                 (Paper I) (ID_SRC)
      real*4        z_2         ! [-0.01/8.0] The redshift adopted
*                                                 for the spectral model
*                                                 (Redshift)
      real*8        dGmKL       ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 KL divergence (in nats units)
*                                                 (dGm_KL)
      real*8        b_dGmH      ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 HDI lower limit (dGm_HLo)
      real*8        B_dGmH_1    ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 HDI upper limit (dGm_HUp)
      real*8        GammaKL_1   ! ? power-law slope,
*                                                 KL divergence (in nats units)
*                                                 (Gamma_KL)
      real*8        b_GammaH_2  ! power-law slope, HDI lower limit
*                                                 (Gamma_HLo)
      real*8        B_GammaH_3  ! power-law slope, HDI upper limit
*                                                 (Gamma_HUp)
      real*8        lognHKL_2   ! log AGN column density,
*                                                 KL divergence (in nats units)
*                                                 (lognH_KL)
      real*8        b_lognHH_4  ! ([cm-2]) log AGN column density,
*                                                 HDI lower limit (lognH_HLo)
      real*8        B_lognHH_5  ! ([cm-2]) log AGN column density,
*                                                 HDI upper limit (lognH_HUp)
      real*8        logPowNKL_1 ! log power-law normalization,
*                                                 KL divergence (in nats units)
*                                                 (logPowNorm_KL)
      real*8        b_logPowNH_2 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI lower limit
*                                                 (logPowNorm_HLo)
      real*8        B_logPowNH_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI upper limit
*                                                 (logPowNorm_HUp)
      real*8        logFracKL   ! ? Constant factor multiplied to
*                                                 the soft power-law,
*                                                 KL divergence (in nats units)
*                                                 (logFrac_KL)
      real*8        b_logFracH  ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 HDI lower limit (logFrac_HLo)
      real*4        B_logFracH_1 ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 HDI upper limit (logFrac_HUp)
      real*8        dGmBF       ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 best-fit model (dGm_BF)
      real*8        GammaBF_1   ! power-law slope, best-fit model
*                                                 (Gamma_BF)
      real*8        lognHBF_2   ! ([cm-2]) log AGN column density,
*                                                 best-fit model (lognH_BF)
      real*8        logPowNBF_1 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 best-fit model (logPowNorm_BF)
      real*4        logFracBF   ! Constant factor multiplied to
*                                                 the soft power-law, best-fit
*                                                 model (logFrac_BF)
      real*4        FluxObsBFS_2 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, best-fit model
*                                                 (FluxObsv_BF_s)
      real*4        FluxObsMedS_4 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, posterior median
*                                                 (FluxObsv_Med_s)
      real*4        b_FluxObsS_8 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_9 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_4 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_5 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_s)
      real*4        FluxObsBFT_2 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, best-fit model
*                                                 (FluxObsv_BF_t)
      real*4        FluxObsMedT_4 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, posterior median
*                                                 (FluxObsv_Med_t)
      real*4        b_FluxObsT_8 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_9 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_4 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_5 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_t)
      real*4        FluxCBFS_1  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV,
*                                                 best-fit model (FluxCorr_BF_s)
      real*4        FluxCMedS_2 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, posterior
*                                                 median (FluxCorr_Med_s)
      real*4        b_FluxCS_4  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_5  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_2 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_3 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_1  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV,
*                                                 best-fit model (FluxCorr_BF_t)
      real*4        FluxCMedT_2 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, posterior
*                                                 median (FluxCorr_Med_t)
      real*4        b_FluxCT_4  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_5  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCT_2 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCT_3 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS_1  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_2 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 posterior median
*                                                 (FluxIntr_Med_s)
      real*4        b_FluxIS_4  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_5  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_2 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_3 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_1 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, best-fit
*                                                 (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_2 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, median
*                                                 (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_4 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_5 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_2 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_3 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_1   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (LumiIntr_BF_s)
      real*8        LumIMedS_2  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (LumiIntr_Med_s)
      real*8        b_LumIS_4   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 1sigma lower limit
*                                                 (LumiIntr_Lo1_s)
      real*8        B_LumIS_5   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 1sigma upper limit
*                                                 (LumiIntr_Up1_s)
      real*8        b2_LumIS_2  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 2sigma lower limit
*                                                 (LumiIntr_Lo2_s)
      real*8        B2_LumIS_3  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 2sigma upper limit
*                                                 (LumiIntr_Up2_s)
      real*8        LumIBF2keV_1 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, best-fit
*                                                 model (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_2 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, posterior
*                                                 median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_4 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_5 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_2 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_3 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_2keV)
      real*8        dGmMed      ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 posterior median (dGm_Med)
      real*8        b_dGm       ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 1sigma lower limit (dGm_Lo)
      real*8        B_dGm_1     ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 1sigma upper limit (dGm_Up)
      real*8        dGmMean     ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 mean (dGm_Mean)
      real*8        s_dGm       ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 standard deviation (dGm_Std)
      real*8        GammaMed_1  ! power-law slope, posterior
*                                                 median (Gamma_Med)
      real*8        b_Gamma_2   ! power-law slope, percentile
*                                                 lower limit (Gamma_Lo)
      real*8        B_Gamma_3   ! power-law slope, percentile
*                                                 upper limit (Gamma_Up)
      real*8        GammaMean_1 ! power-law slope, mean
*                                                 (Gamma_Mean)
      real*4        s_Gamma_1   ! power-law slope,
*                                                 standard deviation (Gamma_Std)
      real*8        lognHMed_2  ! ([cm-2]) log AGN column density,
*                                                 posterior median (lognH_Med)
      real*8        b_lognH_4   ! ([cm-2]) log AGN column density,
*                                                 percentile lower limit
*                                                 (lognH_Lo)
      real*8        B_lognH_5   ! ([cm-2]) log AGN column density,
*                                                 percentile upper limit
*                                                 (lognH_Up)
      real*8        lognHMean_2 ! ([cm-2]) log AGN column density,
*                                                 mean (lognH_Mean)
      real*8        s_lognH_2   ! ([cm-2]) log AGN column density,
*                                                 standard deviation (lognH_Std)
      real*8        logPowNMed_1 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 posterior median
*                                                 (logPowNorm_Med)
      real*8        b_logPowN_2 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile lower limit
*                                                 (logPowNorm_Lo)
      real*8        B_logPowN_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile upper limit
*                                                 (logPowNorm_Up)
      real*8        logPowNMean_1 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 mean (logPowNorm_Mean)
      real*8        s_logPowN_1 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 standard deviation
*                                                 (logPowNorm_Std)
      real*8        logFracMed  ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 posterior median (logFrac_Med)
      real*8        b_logFrac   ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 1sigma lower limit (logFrac_Lo)
      real*8        B_logFrac_1 ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 1sigma upper limit (logFrac_Up)
      real*8        logFracMean ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 mean (logFrac_Mean)
      real*8        s_logFrac   ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 standard deviation
*                                                 (logFrac_Std)
      real*4        logBkgNMed_2 ! log background normalization,
*                                                 posterior median
*                                                 (logBkgNorm_Med)
      real*8        b_logBkgN_4 ! log background normalization,
*                                                 percentile lower limit
*                                                 (logBkgNorm_Lo)
      real*4        B_logBkgN_5 ! log background normalization,
*                                                 percentile upper limit
*                                                 (logBkgNorm_Up)
      real*4        logBkgNMean_2 ! log background normalization,
*                                                 mean (logBkgNorm_Mean)
      real*8        s_logBkgN_2 ! log background normalization,
*                                                 standard deviation
*                                                 (logBkgNorm_Std)
      real*8        chi25_2     ! ?=-99 Chi^2 of the best-fit
*                                                 model against the rebinned
*                                                 data (>=25 cts/bin) (chi25)
      integer*4     dof25_2     ! ?=-99 DOF of the rebinned data,
*                                                 which has at least 25 counts
*                                                 per bin (dof25)
      real*8        logZ_2      ! log10 Bayesian evidence (logZ)

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

C  Declarations for 'table6.dat'	! Powerlaw+blackbody results

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

      integer*4     ID_SRC_5    ! [1/32683] ID of the sources in
*                                                 the eFEDS main X-ray catalog
*                                                 (Paper I) (ID_SRC)
      real*4        z_3         ! [-0.01/8.0] The redshift adopted
*                                                 for the spectral model
*                                                 (Redshift)
      real*8        GammaKL_2   ! power-law slope, KL divergence
*                                                 (in nats units) (Gamma_KL)
      real*8        b_GammaH_4  ! power-law slope, HDI lower limit
*                                                 (Gamma_HLo)
      real*8        B_GammaH_5  ! power-law slope, HDI upper limit
*                                                 (Gamma_HUp)
      real*8        lognHKL_3   ! log AGN column density,
*                                                 KL divergence (in nats units)
*                                                 (lognH_KL)
      real*8        b_lognHH_6  ! ([cm-2]) log AGN column density,
*                                                 HDI lower limit (lognH_HLo)
      real*8        B_lognHH_7  ! ([cm-2]) log AGN column density,
*                                                 HDI upper limit (lognH_HUp)
      real*8        logPowNKL_2 ! log power-law normalization,
*                                                 KL divergence (in nats units)
*                                                 (logPowNorm_KL)
      real*8        b_logPowNH_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI lower limit
*                                                 (logPowNorm_HLo)
      real*8        B_logPowNH_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI upper limit (logPowNorm_HUp)
      real*8        logkTKL_1   ! log temperature,
*                                                 KL divergence (in nats units)
*                                                 (logkT_KL)
      real*8        b_logkTH_2  ! ([keV]) log temperature, HDI lower limit
*                                                 (logkT_HLo)
      real*8        B_logkTH_3  ! ([keV]) log temperature, HDI upper limit
*                                                 (logkT_HUp)
      real*8        logBBNKL    ! log blackbody normalization,
*                                                 KL divergence (in nats units)
*                                                 (logBBNorm_KL)
      real*8        b_logBBNH   ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, HDI
*                                                 lower limit (logBBNorm_HLo)
      real*8        B_logBBNH_1 ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, HDI
*                                                 upper limit (logBBNorm_HUp)
      real*8        GammaBF_2   ! power-law slope, best-fit model
*                                                 (Gamma_BF)
      real*8        lognHBF_3   ! ([cm-2]) log AGN column density,
*                                                 best-fit model (lognH_BF)
      real*8        logPowNBF_2 ! ([ph/keV/cm2/s]) log power- law normalization,
*                                                 best-fit model (logPowNorm_BF)
      real*8        logkTBF_1   ! ([keV]) log temperature, best- fit model
*                                                 (logkT_BF)
      real*8        logBBNBF    ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, best-fit
*                                                 model (logBBNorm_BF)
      real*4        FluxObsBFS_3 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, best-fit model
*                                                 (FluxObsv_BF_s)
      real*4        FluxObsMedS_5 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, posterior median
*                                                 (FluxObsv_Med_s)
      real*4        b_FluxObsS_10 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_11 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_6 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_7 ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_s)
      real*4        FluxObsBFT_3 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, best-fit model
*                                                 (FluxObsv_BF_t)
      real*4        FluxObsMedT_5 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, posterior median
*                                                 (FluxObsv_Med_t)
      real*4        b_FluxObsT_10 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_11 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_6 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_7 ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_t)
      real*4        FluxCBFS_2  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, best-fit
*                                                 model (FluxCorr_BF_s)
      real*4        FluxCMedS_3 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, posterior
*                                                 median (FluxCorr_Med_s)
      real*4        b_FluxCS_6  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_7  ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_4 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_5 ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_2  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, best-fit
*                                                 model (FluxCorr_BF_t)
      real*4        FluxCMedT_3 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, posterior
*                                                 median (FluxCorr_Med_t)
      real*4        b_FluxCT_6  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_7  ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCT_4 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCT_5 ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS_2  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_3 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (FluxIntr_Med_s)
      real*4        b_FluxIS_6  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_7  ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_4 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_5 ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_2 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, best-fit
*                                                 (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_3 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, median
*                                                 (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_6 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_7 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_4 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_5 ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_2   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, best-fit
*                                                 model (LumiIntr_BF_s)
      real*8        LumIMedS_3  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (LumiIntr_Med_s)
      real*8        b_LumIS_6   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_7   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS_4  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_5  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV_2 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, best-fit
*                                                 model (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_3 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, posterior
*                                                 median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_6 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_7 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_4 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_5 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_2keV)
      real*8        GammaMed_2  ! power-law slope, posterior
*                                                 median (Gamma_Med)
      real*8        b_Gamma_4   ! power-law slope, percentile
*                                                 lower limit (Gamma_Lo)
      real*8        B_Gamma_5   ! power-law slope, percentile
*                                                 upper limit (Gamma_Up)
      real*8        GammaMean_2 ! power-law slope, mean
*                                                 (Gamma_Mean)
      real*8        s_Gamma_2   ! power-law slope,
*                                                 standard deviation (Gamma_Std)
      real*8        lognHMed_3  ! ([cm-2]) log AGN column density,
*                                                 posterior median (lognH_Med)
      real*8        b_lognH_6   ! ([cm-2]) log AGN column density,
*                                                 percentile lower limit
*                                                 (lognH_Lo)
      real*8        B_lognH_7   ! ([cm-2]) log AGN column density,
*                                                 percentile upper limit
*                                                 (lognH_Up)
      real*8        lognHMean_3 ! ([cm-2]) log AGN column density,
*                                                 mean (lognH_Mean)
      real*8        s_lognH_3   ! ([cm-2]) log AGN column density,
*                                                 standard deviation (lognH_Std)
      real*8        logPowNMed_2 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 posterior median
*                                                 (logPowNorm_Med)
      real*8        b_logPowN_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile lower limit
*                                                 (logPowNorm_Lo)
      real*8        B_logPowN_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile upper limit
*                                                 (logPowNorm_Up)
      real*8        logPowNMean_2 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 mean (logPowNorm_Mean)
      real*8        s_logPowN_2 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 standard deviation
*                                                 (logPowNorm_Std)
      real*8        logkTMed_1  ! ([keV]) log temperature, posterior
*                                                 median (logkT_Med)
      real*8        b_logkT_2   ! ([keV]) log temperature, percentile
*                                                 lower limit (logkT_Lo)
      real*8        B_logkT_3   ! ([keV]) log temperature, percentile
*                                                 upper limit (logkT_Up)
      real*8        logkTMean_1 ! ([keV]) log temperature, mean
*                                                 (logkT_Mean)
      real*8        s_logkT_1   ! ([keV]) log temperature, standard
*                                                 deviation (logkT_Std)
      real*8        logBBNMed   ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, posterior
*                                                 median (logBBNorm_Med)
      real*8        b_logBBN    ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, percentile
*                                                 lower limit (logBBNorm_Lo)
      real*8        B_logBBN_1  ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, percentile
*                                                 upper limit (logBBNorm_Up)
      real*8        logBBNMean  ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^,
*                                                 mean (logBBNorm_Mean)
      real*8        s_logBBN    ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, standard
*                                                 deviation (logBBNorm_Std)
      real*4        logBkgNMed_3 ! log background normalization,
*                                                 posterior median
*                                                 (logBkgNorm_Med)
      real*8        b_logBkgN_6 ! log background normalization,
*                                                 percentile lower limit
*                                                 (logBkgNorm_Lo)
      real*4        B_logBkgN_7 ! log background normalization,
*                                                 percentile upper limit
*                                                 (logBkgNorm_Up)
      real*4        logBkgNMean_3 ! log background normalization,
*                                                 mean (logBkgNorm_Mean)
      real*8        s_logBkgN_3 ! log background normalization,
*                                                 standard deviation
*                                                 (logBkgNorm_Std)
      real*8        chi25_3     ! ?=-99 Chi^2 of the best-fit
*                                                 model against the rebinned
*                                                 data (>=25 cts/bin) (chi25)
      integer*4     dof25_3     ! ?=-99 DOF of the rebinned data,
*                                                 which has at least
*                                                 25 counts per bin (dof25)
      real*8        logZ_3      ! log10 Bayesian evidence (logZ)
      real*8        logRbb      ! log blackbody relative strength
*                                                 at 1keV (logRbb)

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

C  Declarations for 'table7.dat'	! Gamma-fixed-powerlaw results

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

      integer*4     ID_SRC_6    ! [1/32683] ID of the sources in the
*                                               eFEDS main X-ray catalog
*                                               (Paper I) (ID_SRC)
      real*4        z_4         ! [-0.01/8.0] The redshift adopted
*                                               for the spectral model (Redshift)
      real*8        lognHKL_4   ! log AGN column density,
*                                               KL divergence (in nats units)
*                                               (lognH_KL)
      real*8        b_lognHH_8  ! ([cm-2]) log AGN column density,
*                                               HDI lower limit (lognH_HLo)
      real*8        B_lognHH_9  ! ([cm-2]) log AGN column density,
*                                               HDI upper limit (lognH_HUp)
      real*8        logPowNKL_3 ! log power-law normalization,
*                                               KL divergence (in nats units)
*                                               (logPowNorm_KL)
      real*8        b_logPowNH_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_7 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               HDI upper limit (logPowNorm_HUp)
      real*8        lognHBF_4   ! ([cm-2]) log AGN column density,
*                                               best-fit model (lognH_BF)
      real*8        logPowNBF_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_4 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, best-fit model
*                                               (FluxObsv_BF_s)
      real*4        FluxObsMedS_6 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, posterior median
*                                               (FluxObsv_Med_s)
      real*4        b_FluxObsS_12 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 1sigma lower limit
*                                               (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_13 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 1sigma upper limit
*                                               (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_8 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 2sigma lower limit
*                                               (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_9 ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 2sigma upper limit
*                                               (FluxObsv_Up2_s)
      real*4        FluxObsBFT_4 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, best-fit model
*                                               (FluxObsv_BF_t)
      real*4        FluxObsMedT_6 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, posterior median
*                                               (FluxObsv_Med_t)
      real*4        b_FluxObsT_12 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 1sigma lower limit
*                                               (FluxObsv_Lo1)
      real*4        B_FluxObsT_13 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 1sigma upper limit
*                                               (FluxObsv_Up1)
      real*4        b2_FluxObsT_8 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 2sigma lower limit
*                                               (FluxObsv_Lo2)
      real*4        B2_FluxObsT_9 ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 2sigma upper limit
*                                               (FluxObsv_Up2)
      real*4        FluxCBFS_3  ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, best-fit
*                                               model (FluxCorr_BF_s)
      real*4        FluxCMedS_4 ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, posterior
*                                               median (FluxCorr_Med_s)
      real*4        b_FluxCS_8  ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 1sigma
*                                               lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_9  ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 1sigma
*                                               upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_6 ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 2sigma
*                                               lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_7 ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 2sigma
*                                               upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_3  ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, best-fit
*                                               model (FluxCorr_BF_t)
      real*4        FluxCMedT_4 ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, posterior
*                                               median (FluxCorr_Med_t)
      real*4        b_FluxCT_8  ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 1sigma
*                                               lower limit (FluxCorr_Lo1)
      real*4        B_FluxCT_9  ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 1sigma
*                                               upper limit (FluxCorr_Up1)
      real*4        b2_FluxCT_6 ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 2sigma
*                                               lower limit (FluxCorr_Lo2)
      real*4        B2_FluxCT_7 ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 2sigma
*                                               upper limit (FluxCorr_Up2)
      real*4        FluxIBFS_3  ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, best-fit
*                                               model (FluxIntr_BF_s)
      real*4        FluxIMedS_4 ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, posterior
*                                               median (FluxIntr_Med_s)
      real*4        b_FluxIS_8  ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2keV, 1sigma
*                                               lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_9  ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2keV, 1sigma
*                                               upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_6 ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_7 ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_3 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, best- fit
*                                               (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_4 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, median
*                                               (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_8 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 1sigma
*                                               lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_9 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 1sigma
*                                               upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_6 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 2sigma
*                                               lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_7 ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 2sigma
*                                               upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_3   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, best-fit
*                                               model (LumiIntr_BF_s)
      real*8        LumIMedS_4  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, posterior
*                                               median (LumiIntr_Med_s)
      real*8        b_LumIS_8   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 1sigma
*                                               lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_9   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 1sigma
*                                               upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS_6  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_7  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV_3 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, best-fit model
*                                               (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_4 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, posterior
*                                               median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_8 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 1sigma
*                                               lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_9 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 1sigma
*                                               upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_6 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 2sigma
*                                               lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_7 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 2sigma upper
*                                               limit (LumiIntr_Up2_2keV)
      real*8        lognHMed_4  ! ([cm-2]) log AGN column density,
*                                               posterior median (lognH_Med)
      real*8        b_lognH_8   ! ([cm-2]) log AGN column density,
*                                               percentile lower limit (lognH_Lo)
      real*8        B_lognH_9   ! ([cm-2]) log AGN column density,
*                                               percentile upper limit (lognH_Up)
      real*8        lognHMean_4 ! ([cm-2]) log AGN column density, mean
*                                               (lognH_Mean)
      real*8        s_lognH_4   ! ([cm-2]) log AGN column density,
*                                               standard deviation (lognH_Std)
      real*8        logPowNMed_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               posterior median (logPowNorm_Med)
      real*8        b_logPowN_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               percentile lower limit
*                                               (logPowNorm_Lo)
      real*8        B_logPowN_7 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               percentile upper limit
*                                               (logPowNorm_Up)
      real*8        logPowNMean_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               mean (logPowNorm_Mean)
      real*8        s_logPowN_3 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               standard deviation
*                                               (logPowNorm_Std)
      real*4        logBkgNMed_4 ! log background normalization,
*                                               posterior median (logBkgNorm_Med)
      real*4        b_logBkgN_8 ! log background normalization,
*                                               percentile lower limit
*                                               (logBkgNorm_Lo)
      real*4        B_logBkgN_9 ! log background normalization,
*                                               percentile upper limit
*                                               (logBkgNorm_Up)
      real*4        logBkgNMean_4 ! log background normalization,
*                                               mean (logBkgNorm_Mean)
      real*8        s_logBkgN_4 ! log background normalization,
*                                               standard deviation
*                                               (logBkgNorm_Std)
      real*8        logZ_4      ! log10 Bayesian evidence (logZ)

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

C  Declarations for 'table8.dat'	! Shape-fixed-powerlaw results

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

      integer*4     ID_SRC_7    ! [1/32683] ID of the sources in
*                                                the eFEDS main X-ray catalog
*                                                (Paper I) (ID_SRC)
      real*4        z_5         ! [-0.01/8.0] The redshift adopted
*                                                for the spectral model
*                                                (Redshift)
      real*8        logPowNKL_4 ! log power-law normalization,
*                                                KL divergence (in nats units)
*                                                (logPowNorm_KL)
      real*8        b_logPowNH_8 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_9 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI upper limit (logPowNorm_HUp)
      real*8        logPowNBF_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_5 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, best-fit model
*                                                (FluxObsv_BF_s)
      real*4        FluxObsMedS_7 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, posterior median
*                                                (FluxObsv_Med_s)
      real*4        b_FluxObsS_14 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 1sigma lower limit
*                                                (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_15 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 1sigma upper limit
*                                                (FluxObsvs)
      real*4        b2_FluxObsS_10 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 2sigma lower limit
*                                                (FluxObsvs)
      real*4        B2_FluxObsS_11 ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 2sigma upper limit
*                                                (FluxObsvs)
      real*4        FluxObsBFT_5 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, best-fit model
*                                                (FluxObsv_BF_t)
      real*4        FluxObsMedT_7 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, posterior median
*                                                (FluxObsv_Med_t)
      real*4        b_FluxObsT_14 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 1sigma lower limit
*                                                (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_15 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 1sigma upper limit
*                                                (FluxObsvt)
      real*4        b2_FluxObsT_10 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 2sigma lower limit
*                                                (FluxObsvt)
      real*4        B2_FluxObsT_11 ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 2sigma upper limit
*                                                (FluxObsvt)
      real*4        FluxCBFS_4  ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV,
*                                                best-fit model (FluxCorr_BF_s)
      real*4        FluxCMedS_5 ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, posterior
*                                                median (FluxCorr_Med_s)
      real*4        b_FluxCS_10 ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 1sigma
*                                                lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_11 ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 1sigma
*                                                upper limit (FluxCorrs)
      real*4        b2_FluxCS_8 ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 2sigma
*                                                lower limit (FluxCorrs)
      real*4        B2_FluxCS_9 ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 2sigma
*                                                upper limit (FluxCorrs)
      real*4        FluxCBFT_4  ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, best-fit
*                                                model (FluxCorr_BF_t)
      real*4        FluxCMedT_5 ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, posterior
*                                                median (FluxCorr_Med_t)
      real*4        b_FluxCT_10 ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 1sigma
*                                                lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_11 ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 1sigma
*                                                upper limit (FluxCorrt)
      real*4        b2_FluxCT_8 ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 2sigma
*                                                lower limit (FluxCorrt)
      real*4        B2_FluxCT_9 ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 2sigma
*                                                upper limit (FluxCorrt)
      real*4        FluxIBFS_4  ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, best-fit
*                                                model (FluxIntr_BF_s)
      real*4        FluxIMedS_5 ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, posterior
*                                                median (FluxIntr_Med_s)
      real*4        b_FluxIS_10 ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2keV, 1sigma
*                                                lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_11 ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2keV, 1sigma
*                                                upper limit (FluxIntrs)
      real*4        b2_FluxIS_8 ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                lower limit (FluxIntrs)
      real*4        B2_FluxIS_9 ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                upper limit (FluxIntrs)
      real*4        FluxIBF2keV_4 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, best- fit
*                                                (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_5 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, median
*                                                (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_10 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 1sigma
*                                                lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_11 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 1sigma
*                                                upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_8 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 2sigma
*                                                lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_9 ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 2sigma
*                                                upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_4   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, best-fit
*                                                model (LumiIntr_BF_s)
      real*8        LumIMedS_5  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, posterior
*                                                median (LumiIntr_Med_s)
      real*8        b_LumIS_10  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 1sigma
*                                                lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_11  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 1sigma
*                                                upper limit (LumiIntrs)
      real*8        b2_LumI2S   ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                lower limit (LumiIntrs)
      real*8        B2_LumIS_8  ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                upper limit (LumiIntrs)
      real*8        LumIBF2keV_4 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, best-fit model
*                                                (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_5 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, posterior
*                                                median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_10 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 1sigma lower
*                                                limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_11 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 1sigma upper
*                                                limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_8 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 2sigma lower
*                                                limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_9 ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 2sigma upper
*                                                limit (LumiIntr_Up2_2keV)
      real*8        logPowNMed_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                posterior median
*                                                (logPowNorm_Med)
      real*8        b_logPowN_8 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                percentile lower limit
*                                                (logPowNorm_Lo)
      real*8        B_logPowN_9 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                percentile upper limit
*                                                (logPowNorm_Up)
      real*8        logPowNMean_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                mean (logPowNorm_Mean)
      real*8        s_logPowN_4 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                standard deviation
*                                                (logPowNorm_Std)
      real*4        logBkgNMed_5 ! log background normalization,
*                                                posterior median
*                                                (logBkgNorm_Med)
      real*4        b_logBkgN_10 ! log background normalization,
*                                                percentile lower limit
*                                                (logBkgNorm_Lo)
      real*4        B_logBkgN_11 ! log background normalization,
*                                                percentile upper limit
*                                                (logBkgNorm_Up)
      real*4        logBkgNMean_5 ! log background normalization,
*                                                mean (logBkgNorm_Mean)
      real*8        s_logBkgN_5 ! log background normalization,
*                                                standard deviation
*                                                (logBkgNorm_Std)
      real*8        logZ_5      ! log10 Bayesian evidence (logZ)

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

C  Declarations for 'table9.dat'	! Soft-band fitting results

      integer*4 nr__8
      parameter (nr__8=27910)	! Number of records
      character*489 ar__8  	! Full-size record

      integer*4     ID_SRC_8    ! [1/32683] ID of the sources in the
*                                              eFEDS main X-ray catalog
*                                              (Paper I) (ID_SRC)
      real*4        z_6         ! [-0.01/8.0] The redshift adopted
*                                              for the spectral model (Redshift)
      real*8        GammaKL_3   ! power-law slope, KL divergence
*                                              (in nats units) (Gamma_KL)
      real*8        b_GammaH_6  ! power-law slope, HDI lower limit
*                                              (Gamma_HLo)
      real*8        B_GammaH_7  ! power-law slope, HDI upper limit
*                                              (Gamma_HUp)
      real*8        lognHKL_5   ! log AGN column density,
*                                              KL divergence (in nats units)
*                                              (lognH_KL)
      real*8        b_lognHH_10 ! ([cm-2]) log AGN column density,
*                                              HDI lower limit (lognH_HLo)
      real*8        B_lognHH_11 ! ([cm-2]) log AGN column density,
*                                              HDI upper limit (lognH_HUp)
      real*8        logPowNKL_5 ! log power-law normalization,
*                                              KL divergence (in nats units)
*                                              (logPowNorm_KL)
      real*8        b_logPowNH_10 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_11 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF_3   ! power-law slope, best-fit model
*                                              (Gamma_BF)
      real*8        lognHBF_5   ! ([cm-2]) log AGN column density,
*                                              best-fit model (lognH_BF)
      real*4        logPowNBF_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_6 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, best-fit model
*                                              (FluxObsv_BF_s)
      real*4        FluxObsMedS_8 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, posterior median
*                                              (FluxObsv_Med_s)
      real*4        b_FluxObsS_16 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 1sigma lower limit
*                                              (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_17 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 1sigma upper limit
*                                              (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_12 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 2sigma lower limit
*                                              (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_13 ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 2sigma upper limit
*                                              (FluxObsv_Up2_s)
      real*8        GammaMed_3  ! power-law slope, posterior median
*                                              (Gamma_Med)
      real*8        b_Gamma_6   ! power-law slope, percentile
*                                              lower limit (Gamma_Lo)
      real*8        B_Gamma_7   ! power-law slope, percentile
*                                              upper limit (Gamma_Up)
      real*8        GammaMean_3 ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma_3   ! power-law slope,
*                                              standard deviation (Gamma_Std)
      real*8        lognHMed_5  ! ([cm-2]) log AGN column density,
*                                              posterior median (lognH_Med)
      real*8        b_lognH_10  ! ([cm-2]) log AGN column density, percentile
*                                              lower limit (lognH_Lo)
      real*8        B_lognH_11  ! ([cm-2]) log AGN column density, percentile
*                                              upper limit (lognH_Up)
      real*8        lognHMean_5 ! ([cm-2]) log AGN column density,
*                                              mean (lognH_Mean)
      real*8        s_lognH_5   ! ([cm-2]) log AGN column density,
*                                              standard deviation (lognH_Std)
      real*8        logPowNMed_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              posterior median (logPowNorm_Med)
      real*8        b_logPowN_10 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile lower limit
*                                              (logPowNorm_Lo)
      real*8        B_logPowN_11 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile upper limit
*                                              (logPowNorm_Up)
      real*8        logPowNMean_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              mean (logPowNorm_Mean)
      real*8        s_logPowN_5 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              standard deviation
*                                              (logPowNorm_Std)
      real*4        logBkgNMed_6 ! log background normalization,
*                                              posterior median (logBkgNorm_Med)
      real*4        blogBkgN    ! log background normalization,
*                                              percentile lower limit
*                                              (logBkgNorm_Lo)
      real*4        BlogBkgN_1  ! log background normalization,
*                                              percentile upper limit
*                                              (logBkgNorm_Up)
      real*4        logBkgNMean_6 ! log background normalization,
*                                              mean (logBkgNorm_Mean)
      real*8        logBkgN     ! log background normalization,
*                                              standard deviation
*                                              (logBkgNorm_Std)
      real*8        logZ_6      ! log10 Bayesian evidence (logZ)

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

C  Declarations for 'table10.dat'	! Hard-band fitting results

      integer*4 nr__9
      parameter (nr__9=27910)	! Number of records
      character*543 ar__9  	! Full-size record

      integer*4     ID_SRC_9    ! [1/32683] ID of the sources in the
*                                              eFEDS main X-ray catalog
*                                              (Paper I) (ID_SRC)
      real*4        z_7         ! [-0.01/8.0] The redshift adopted
*                                              for the spectral model (Redshift)
      real*8        GammaKL_4   ! power-law slope, KL divergence
*                                              (in nats units) (Gamma_KL)
      real*8        e_Gamma     ! power-law slope, HDI lower limit
*                                              (Gamma_HLo)
      real*8        E_Gamma_1   ! power-law slope, HDI upper limit
*                                              (Gamma_HUp)
      real*8        lognHKL_6   ! log AGN column density,
*                                              KL divergence (in nats units)
*                                              (lognH_KL)
      real*8        b_lognH_12  ! ([cm-2]) log AGN column density,
*                                              HDI lower limit (lognH_HLo)
      real*8        B_lognH_13  ! ([cm-2]) log AGN column density,
*                                              HDI upper limit (lognH_HUp)
      real*8        logPowNKL_6 ! log power-law normalization,
*                                              KL divergence (in nats units)
*                                              (logPowNorm_KL)
      real*8        b_logPowN_12 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowN_13 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF_4   ! power-law slope, best-fit model
*                                              (Gamma_BF)
      real*8        lognHBF_6   ! ([cm-2]) log AGN column density,
*                                              best-fit model (lognH_BF)
      real*8        logPowNBF_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              best-fit model (logPowNorm_BF)
      real*4        FluxObsBFT_6 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_8 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_16 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              1sigma lower limit
*                                              (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_17 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              1sigma upper limit
*                                              (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_12 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              2sigma lower limit
*                                              (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_13 ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              2sigma upper limit
*                                              (FluxObsv_Up2_t)
      real*4        FluxCBFT_5  ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, best-fit model
*                                              (FluxCorr_BF_t)
      real*4        FluxCMedT_6 ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, posterior
*                                              median (FluxCorr_Med_t)
      real*4        b_FluxCMedT ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 1sigma
*                                              lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCMedT_1 ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 1sigma
*                                              upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCMedT ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 2sigma
*                                              lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCMedT_1 ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 2sigma
*                                              upper limit (FluxCorr_Up2_t)
      real*8        GammaMed_4  ! power-law slope, posterior median
*                                              (Gamma_Med)
      real*8        b_GammaMed  ! power-law slope, percentile
*                                              lower limit (Gamma_Lo)
      real*8        B_GammaMed_1 ! power-law slope, percentile
*                                              upper limit (Gamma_Up)
      real*8        GammaMean_4 ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma_4   ! power-law slope,
*                                              standard deviation (Gamma_Std)
      real*8        lognHMed_6  ! ([cm-2]) log AGN column density,
*                                              posterior median (lognH_Med)
      real*8        b_lognHMed  ! ([cm-2]) log AGN column density,
*                                              percentile lower limit (lognH_Lo)
      real*8        B_lognHMed_1 ! ([cm-2]) log AGN column density,
*                                              percentile upper limit (lognH_Up)
      real*8        lognHMean_6 ! ([cm-2]) log AGN column density,
*                                              mean (lognH_Mean)
      real*8        s_lognH_6   ! ([cm-2]) log AGN column density,
*                                              standard deviation (lognH_Std)
      real*8        logPowNMed_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              posterior median (logPowNorm_Med)
      real*8        b_logPowNMed ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile lower limit
*                                              (logPowNorm_Lo)
      real*8        B_logPowNMed_1 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile upper limit
*                                              (logPowNorm_Up)
      real*8        logPowNMean_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              mean (logPowNorm_Mean)
      real*8        s_logPowN_6 ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              standard deviation
*                                              (logPowNorm_Std)
      real*4        logBkgNMed_7 ! log background normalization,
*                                              posterior median (logBkgNorm_Med)
      real*4        b_logBkgNMed ! log background normalization,
*                                              percentile lower limit
*                                              (logBkgNorm_Lo)
      real*4        B_logBkgNMed_1 ! log background normalization,
*                                              percentile upper limit
*                                              (logBkgNorm_Up)
      real*4        logBkgNMean_7 ! log background normalization,
*                                              mean (logBkgNorm_Mean)
      real*8        s_logBkgN_6 ! log background normalization,
*                                              standard deviation
*                                              (logBkgNorm_Std)
      real*8        logZ_7      ! log10 Bayesian evidence (logZ)

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

C  Loading file 'table1.dat'	! The eFEDS AGN catalog

C  Format for file interpretation

    1 format(
     +  I5,1X,F18.14,1X,E24.17,1X,F13.7,1X,E5.4,1X,A11,1X,F18.14,1X,
     +  E23.17,1X,I1,1X,I1,1X,F12.10,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,
     +  I1,1X,E22.17,1X,E22.17,1X,F24.18,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E14.9,1X,E14.9,1X,E14.9,1X,F7.4,1X,F7.4,1X,F7.4,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,F11.9,1X,F9.6,1X,F9.6,1X,F9.6,1X,F11.9,1X,F11.9,1X,
     +  F9.7,1X,F10.8,1X,F11.6,1X,F11.6,1X,F11.6,1X,F11.6,1X,F11.6,1X,
     +  F11.6,1X,F10.6,1X,F10.6,1X,F11.7,1X,E13.8,1X,F11.7,1X,E13.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,22079
        read(1,'(A736)')ar__
        read(ar__,1)
     +  ID_SRC,RAcdeg,DEcdeg,DetLike,inArea90,CTPLS8,RALdeg,DELdeg,
     +  CTPqual,CTPClass,CTPz,CTPzgrade,inKiDSFlag,LxModel,NHclass,
     +  FSModel,FHModel,galNH,galNHI,SrcCts,FluxcMedS,b_FluxcS,
     +  B_FluxcS_1,FluxcMedT,b_FluxcT,B_FluxcT_1,FluxIMedS,b_FluxIS,
     +  B_FluxIS_1,FluxIMed2keV,b_FluxI2keV,B_FluxI2keV_1,LumIMedS,
     +  b_LumIS,B_LumIS_1,LumIMed2keV,b_LumI2keV,B_LumI2keV_1,
     +  FluxObsMedS,b_FluxObsS,B_FluxObsS_1,FluxObsMedT,b_FluxObsT,
     +  B_FluxObsT_1,lognHKLm1,b_lognHm1,B_lognHm1_1,lognHMedm1,
     +  GammaKLm3,b_GammaHm3,B_GammaHm3_1,GammaMedm3,logZm0,logZm1,
     +  logZm2,logZm3,logZm4,logZm5,L2500,L5100,W1mag,e_W1mag,W2mag,
     +  e_W2mag
        if(ar__(686:696) .EQ. '') W1mag = rNULL__
        if(ar__(698:710) .EQ. '') e_W1mag = rNULL__
        if(ar__(712:722) .EQ. '') W2mag = rNULL__
        if(ar__(724:736) .EQ. '') e_W2mag = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID_SRC,RAcdeg,DEcdeg,DetLike,inArea90,CTPLS8,RALdeg,DELdeg,
     +  CTPqual,CTPClass,CTPz,CTPzgrade,inKiDSFlag,LxModel,NHclass,
     +  FSModel,FHModel,galNH,galNHI,SrcCts,FluxcMedS,b_FluxcS,
     +  B_FluxcS_1,FluxcMedT,b_FluxcT,B_FluxcT_1,FluxIMedS,b_FluxIS,
     +  B_FluxIS_1,FluxIMed2keV,b_FluxI2keV,B_FluxI2keV_1,LumIMedS,
     +  b_LumIS,B_LumIS_1,LumIMed2keV,b_LumI2keV,B_LumI2keV_1,
     +  FluxObsMedS,b_FluxObsS,B_FluxObsS_1,FluxObsMedT,b_FluxObsT,
     +  B_FluxObsT_1,lognHKLm1,b_lognHm1,B_lognHm1_1,lognHMedm1,
     +  GammaKLm3,b_GammaHm3,B_GammaHm3_1,GammaMedm3,logZm0,logZm1,
     +  logZm2,logZm3,logZm4,logZm5,L2500,L5100,W1mag,e_W1mag,W2mag,
     +  e_W2mag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Spectral properties

C  Format for file interpretation

    2 format(
     +  I5,1X,F18.14,1X,E24.17,1X,F13.7,1X,F11.7,1X,E5.4,1X,E22.17,1X,
     +  E22.17,1X,F11.6,1X,F24.18,1X,F10.6,1X,E10.3,1X,I3,1X,I3,1X,I3,
     +  1X,E13.8,1X,E13.8,1X,F12.10,1X,I1,1X,I1,1X,E13.7,1X,F13.11,1X,
     +  E16.9,1X,E13.8,1X,E16.9,1X,E13.8,1X,E16.9,1X,E13.8,1X,E16.9,
     +  1X,E14.9,1X,E16.9,1X,E14.9,1X,E16.9,1X,E14.9,1X,F6.1,1X,F10.7,
     +  1X,F6.1,1X,F10.7,1X,F6.1,1X,F10.7,1X,F6.1,1X,F10.7,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F9.7,1X,I3,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,27910
        read(1,'(A645)')ar__1
        read(ar__1,2)
     +  ID_SRC_1,RAcdeg_1,DEcdeg_1,DetLike_1,ExtLike,inArea90_1,
     +  galNH_1,galNHI_1,ExpTime,SrcCts_1,RAdeg,DEdeg,Rad,Ann1,Ann2,
     +  RegArea,Backscal,BackscalBkg,FSModel_1,FHModel_1,Rate0_2_2_3,
     +  e_Rate0_2_2_3,Rate0_2_0_5,e_Rate0_2_0_5,Rate0_5_1,e_Rate0_5_1,
     +  Rate1_2,e_Rate1_2,Rate2_4_5,e_Rate2_4_5,Rate2_3_5,e_Rate2_3_5,
     +  Rate5_8,e_Rate5_8,BkgCts0_2_0_6,e_BkgCts0_2_0_6,BkgCts0_6_2_3,
     +  e_BkgCts0_6_2_3,BkgCts2_3_5,e_BkgCts2_3_5,BkgCts5_8,
     +  e_BkgCts5_8,PSFCor0_2_0_5,PSFCor0_5_1,PSFCor1_2,PSFCor2_3_5,
     +  Nempty,FluxObsMedS_1,b_FluxObsS_2,B_FluxObsS_3,FluxObsMedT_1,
     +  b_FluxObsT_2,B_FluxObsT_3
c    ..............Just test output...........
        write(6,2)
     +  ID_SRC_1,RAcdeg_1,DEcdeg_1,DetLike_1,ExtLike,inArea90_1,
     +  galNH_1,galNHI_1,ExpTime,SrcCts_1,RAdeg,DEdeg,Rad,Ann1,Ann2,
     +  RegArea,Backscal,BackscalBkg,FSModel_1,FHModel_1,Rate0_2_2_3,
     +  e_Rate0_2_2_3,Rate0_2_0_5,e_Rate0_2_0_5,Rate0_5_1,e_Rate0_5_1,
     +  Rate1_2,e_Rate1_2,Rate2_4_5,e_Rate2_4_5,Rate2_3_5,e_Rate2_3_5,
     +  Rate5_8,e_Rate5_8,BkgCts0_2_0_6,e_BkgCts0_2_0_6,BkgCts0_6_2_3,
     +  e_BkgCts0_6_2_3,BkgCts2_3_5,e_BkgCts2_3_5,BkgCts5_8,
     +  e_BkgCts5_8,PSFCor0_2_0_5,PSFCor0_5_1,PSFCor1_2,PSFCor2_3_5,
     +  Nempty,FluxObsMedS_1,b_FluxObsS_2,B_FluxObsS_3,FluxObsMedT_1,
     +  b_FluxObsT_2,B_FluxObsT_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! APEC fitting results

C  Format for file interpretation

    3 format(
     +  I5,1X,F3.1,1X,F12.10,1X,F10.7,1X,F10.7,1X,F11.9,1X,E15.8,1X,
     +  E15.8,1X,F12.10,1X,F11.8,1X,E15.8,1X,F10.8,1X,F10.7,1X,F10.7,
     +  1X,F10.7,1X,E15.8,1X,E15.8,1X,F10.7,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,F10.7,1X,F10.7,1X,F10.7,1X,
     +  F10.7,1X,F12.10,1X,E14.8,1X,E14.7,1X,E14.8,1X,E15.8,1X,F12.10,
     +  1X,F12.9,1X,F11.8,1X,F12.9,1X,F11.8,1X,F12.10,1X,F10.7,1X,
     +  F11.7,1X,F10.7,1X,F10.7,1X,F12.10,1X,E16.9,1X,E13.7,1X,E13.8,
     +  1X,E15.9,1X,F12.10,1X,F12.6,1X,I3,1X,F11.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,27910
        read(1,'(A715)')ar__2
        read(ar__2,3)
     +  ID_SRC_2,z,lognHKL,b_lognHH,B_lognHH_1,logkTKL,b_logkTH,
     +  B_logkTH_1,logAbundKL,b_logAbunH,B_logAbunH_1,logApecNormKL,
     +  b_logApecNormH,B_logApecNormH_1,lognHBF,logkTBF,logAbundBF,
     +  logApecNBF,FluxObsBFS,FluxObsMedS_2,b_FluxObsS_4,B_FluxObsS_5,
     +  b2_FluxObsS,B2_FluxObsS_1,FluxObsBFT,FluxObsMedT_2,
     +  b_FluxObsT_4,B_FluxObsT_5,b2_FluxObsT,B2_FluxObsT_1,lognHMed,
     +  b_lognH,B_lognH_1,lognHMean,s_lognH,logkTMed,b_logkT,
     +  B_logkT_1,logkTMean,s_logkT,logAbundMed,b_logAbund,
     +  B_logAbund_1,logAbundMean,s_logAbund,logApecNMed,b_logApecN,
     +  B_logApecN_1,logApecNMean,s_logApecN,logBkgNMed,b_logBkgN,
     +  B_logBkgN_1,logBkgNMean,s_logBkgN,chi25,dof25,logZ
c    ..............Just test output...........
        write(6,3)
     +  ID_SRC_2,z,lognHKL,b_lognHH,B_lognHH_1,logkTKL,b_logkTH,
     +  B_logkTH_1,logAbundKL,b_logAbunH,B_logAbunH_1,logApecNormKL,
     +  b_logApecNormH,B_logApecNormH_1,lognHBF,logkTBF,logAbundBF,
     +  logApecNBF,FluxObsBFS,FluxObsMedS_2,b_FluxObsS_4,B_FluxObsS_5,
     +  b2_FluxObsS,B2_FluxObsS_1,FluxObsBFT,FluxObsMedT_2,
     +  b_FluxObsT_4,B_FluxObsT_5,b2_FluxObsT,B2_FluxObsT_1,lognHMed,
     +  b_lognH,B_lognH_1,lognHMean,s_lognH,logkTMed,b_logkT,
     +  B_logkT_1,logkTMean,s_logkT,logAbundMed,b_logAbund,
     +  B_logAbund_1,logAbundMean,s_logAbund,logApecNMed,b_logApecN,
     +  B_logApecN_1,logApecNMean,s_logApecN,logBkgNMed,b_logBkgN,
     +  B_logBkgN_1,logBkgNMean,s_logBkgN,chi25,dof25,logZ
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Single-powerlaw results

C  Format for file interpretation

    4 format(
     +  I5,1X,E9.2,1X,F12.9,1X,F10.8,1X,F9.7,1X,F11.9,1X,F9.6,1X,F9.6,
     +  1X,F10.8,1X,F10.7,1X,F10.7,1X,F12.9,1X,F9.6,1X,F10.7,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E14.8,1X,E14.8,1X,E14.8,1X,
     +  E14.8,1X,E14.8,1X,E14.8,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,
     +  F10.6,1X,F10.8,1X,F10.8,1X,F9.7,1X,F10.8,1X,F11.9,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,F9.6,1X,F11.9,1X,F10.7,1X,F10.7,1X,F10.7,1X,
     +  F10.7,1X,F12.10,1X,E16.9,1X,F14.11,1X,E14.7,1X,E16.9,1X,
     +  F12.10,1X,F12.6,1X,I3,1X,F11.6,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,27910
        read(1,'(A990)')ar__3
        read(ar__3,4)
     +  ID_SRC_3,z_1,GammaKL,b_GammaH,B_GammaH_1,lognHKL_1,b_lognHH_2,
     +  B_lognHH_3,logPowNKL,b_logPowNH,B_logPowNH_1,GammaBF,
     +  lognHBF_1,logPowNBF,FluxObsBFS_1,FluxObsMedS_3,b_FluxObsS_6,
     +  B_FluxObsS_7,b2_FluxObsS_2,B2_FluxObsS_3,FluxObsBFT_1,
     +  FluxObsMedT_3,b_FluxObsT_6,B_FluxObsT_7,b2_FluxObsT_2,
     +  B2_FluxObsT_3,FluxcBFS,FluxcMedS_1,b_FluxcS_2,B_FluxcS_3,
     +  b2_FluxcS,B2_FluxcS_1,FluxcBFT,FluxcMedT_1,b_FluxcT_2,
     +  B_FluxcT_3,b2_FluxcT,B2_FluxcT_1,FluxIBFS,FluxIMedS_1,
     +  b_FluxIS_2,B_FluxIS_3,b2_FluxIS,B2_FluxIS_1,FluxIBF2keV,
     +  FluxIMed2keV_1,b_FluxI2keV_2,B_FluxI2keV_3,b2_FluxI2keV,
     +  B2_FluxI2keV_1,LumIBFS,LumIMedS_1,b_LumIS_2,B_LumIS_3,
     +  b2_LumIS,B2_LumIS_1,LumIBF2keV,LumIMed2keV_1,b_LumI2keV_2,
     +  B_LumI2keV_3,b2_LumI2keV,B2_LumI2keV_1,GammaMed,b_Gamma,
     +  B_Gamma_1,GammaMean,s_Gamma,lognHMed_1,b_lognH_2,B_lognH_3,
     +  lognHMean_1,s_lognH_1,logPowNMed,b_logPowN,B_logPowN_1,
     +  logPowNMean,s_logPowN,logBkgNMed_1,b_logBkgN_2,B_logBkgN_3,
     +  logBkgNMean_1,s_logBkgN_1,chi25_1,dof25_1,logZ_1,NHclass_1
c    ..............Just test output...........
        write(6,4)
     +  ID_SRC_3,z_1,GammaKL,b_GammaH,B_GammaH_1,lognHKL_1,b_lognHH_2,
     +  B_lognHH_3,logPowNKL,b_logPowNH,B_logPowNH_1,GammaBF,
     +  lognHBF_1,logPowNBF,FluxObsBFS_1,FluxObsMedS_3,b_FluxObsS_6,
     +  B_FluxObsS_7,b2_FluxObsS_2,B2_FluxObsS_3,FluxObsBFT_1,
     +  FluxObsMedT_3,b_FluxObsT_6,B_FluxObsT_7,b2_FluxObsT_2,
     +  B2_FluxObsT_3,FluxcBFS,FluxcMedS_1,b_FluxcS_2,B_FluxcS_3,
     +  b2_FluxcS,B2_FluxcS_1,FluxcBFT,FluxcMedT_1,b_FluxcT_2,
     +  B_FluxcT_3,b2_FluxcT,B2_FluxcT_1,FluxIBFS,FluxIMedS_1,
     +  b_FluxIS_2,B_FluxIS_3,b2_FluxIS,B2_FluxIS_1,FluxIBF2keV,
     +  FluxIMed2keV_1,b_FluxI2keV_2,B_FluxI2keV_3,b2_FluxI2keV,
     +  B2_FluxI2keV_1,LumIBFS,LumIMedS_1,b_LumIS_2,B_LumIS_3,
     +  b2_LumIS,B2_LumIS_1,LumIBF2keV,LumIMed2keV_1,b_LumI2keV_2,
     +  B_LumI2keV_3,b2_LumI2keV,B2_LumI2keV_1,GammaMed,b_Gamma,
     +  B_Gamma_1,GammaMean,s_Gamma,lognHMed_1,b_lognH_2,B_lognH_3,
     +  lognHMean_1,s_lognH_1,logPowNMed,b_logPowN,B_logPowN_1,
     +  logPowNMean,s_logPowN,logBkgNMed_1,b_logBkgN_2,B_logBkgN_3,
     +  logBkgNMean_1,s_logBkgN_1,chi25_1,dof25_1,logZ_1,NHclass_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Double-powerlaw results

C  Format for file interpretation

    5 format(
     +  I5,1X,E9.2,1X,F12.10,1X,F10.8,1X,F9.7,1X,F13.9,1X,F10.8,1X,
     +  F9.7,1X,F11.9,1X,F9.6,1X,F9.6,1X,F11.9,1X,F10.7,1X,F10.7,1X,
     +  F13.10,1X,F11.8,1X,E16.9,1X,F8.6,1X,F11.8,1X,F9.6,1X,F10.7,1X,
     +  E15.8,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,F8.4,1X,F8.4,1X,F8.4,
     +  1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,1X,F9.7,1X,F10.8,1X,F9.7,1X,F9.7,1X,F11.9,1X,
     +  F10.8,1X,F10.8,1X,F9.7,1X,F10.8,1X,E12.7,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,F11.9,1X,F10.7,1X,F10.7,1X,F10.7,1X,F10.7,1X,
     +  F12.10,1X,F11.8,1X,F11.8,1X,F12.9,1X,F11.8,1X,F10.8,1X,E16.9,
     +  1X,F14.11,1X,E14.7,1X,E16.9,1X,F12.10,1X,F12.6,1X,I3,1X,F11.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,27910
        read(1,'(A1175)')ar__4
        read(ar__4,5)
     +  ID_SRC_4,z_2,dGmKL,b_dGmH,B_dGmH_1,GammaKL_1,b_GammaH_2,
     +  B_GammaH_3,lognHKL_2,b_lognHH_4,B_lognHH_5,logPowNKL_1,
     +  b_logPowNH_2,B_logPowNH_3,logFracKL,b_logFracH,B_logFracH_1,
     +  dGmBF,GammaBF_1,lognHBF_2,logPowNBF_1,logFracBF,FluxObsBFS_2,
     +  FluxObsMedS_4,b_FluxObsS_8,B_FluxObsS_9,b2_FluxObsS_4,
     +  B2_FluxObsS_5,FluxObsBFT_2,FluxObsMedT_4,b_FluxObsT_8,
     +  B_FluxObsT_9,b2_FluxObsT_4,B2_FluxObsT_5,FluxCBFS_1,
     +  FluxCMedS_2,b_FluxCS_4,B_FluxCS_5,b2_FluxCS_2,B2_FluxCS_3,
     +  FluxCBFT_1,FluxCMedT_2,b_FluxCT_4,B_FluxCT_5,b2_FluxCT_2,
     +  B2_FluxCT_3,FluxIBFS_1,FluxIMedS_2,b_FluxIS_4,B_FluxIS_5,
     +  b2_FluxIS_2,B2_FluxIS_3,FluxIBF2keV_1,FluxIMed2keV_2,
     +  b_FluxI2keV_4,B_FluxI2keV_5,b2_FluxI2keV_2,B2_FluxI2keV_3,
     +  LumIBFS_1,LumIMedS_2,b_LumIS_4,B_LumIS_5,b2_LumIS_2,
     +  B2_LumIS_3,LumIBF2keV_1,LumIMed2keV_2,b_LumI2keV_4,
     +  B_LumI2keV_5,b2_LumI2keV_2,B2_LumI2keV_3,dGmMed,b_dGm,B_dGm_1,
     +  dGmMean,s_dGm,GammaMed_1,b_Gamma_2,B_Gamma_3,GammaMean_1,
     +  s_Gamma_1,lognHMed_2,b_lognH_4,B_lognH_5,lognHMean_2,
     +  s_lognH_2,logPowNMed_1,b_logPowN_2,B_logPowN_3,logPowNMean_1,
     +  s_logPowN_1,logFracMed,b_logFrac,B_logFrac_1,logFracMean,
     +  s_logFrac,logBkgNMed_2,b_logBkgN_4,B_logBkgN_5,logBkgNMean_2,
     +  s_logBkgN_2,chi25_2,dof25_2,logZ_2
        if(ar__4(51:63) .EQ. '') GammaKL_1 = rNULL__
        if(ar__4(152:164) .EQ. '') logFracKL = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  ID_SRC_4,z_2,dGmKL,b_dGmH,B_dGmH_1,GammaKL_1,b_GammaH_2,
     +  B_GammaH_3,lognHKL_2,b_lognHH_4,B_lognHH_5,logPowNKL_1,
     +  b_logPowNH_2,B_logPowNH_3,logFracKL,b_logFracH,B_logFracH_1,
     +  dGmBF,GammaBF_1,lognHBF_2,logPowNBF_1,logFracBF,FluxObsBFS_2,
     +  FluxObsMedS_4,b_FluxObsS_8,B_FluxObsS_9,b2_FluxObsS_4,
     +  B2_FluxObsS_5,FluxObsBFT_2,FluxObsMedT_4,b_FluxObsT_8,
     +  B_FluxObsT_9,b2_FluxObsT_4,B2_FluxObsT_5,FluxCBFS_1,
     +  FluxCMedS_2,b_FluxCS_4,B_FluxCS_5,b2_FluxCS_2,B2_FluxCS_3,
     +  FluxCBFT_1,FluxCMedT_2,b_FluxCT_4,B_FluxCT_5,b2_FluxCT_2,
     +  B2_FluxCT_3,FluxIBFS_1,FluxIMedS_2,b_FluxIS_4,B_FluxIS_5,
     +  b2_FluxIS_2,B2_FluxIS_3,FluxIBF2keV_1,FluxIMed2keV_2,
     +  b_FluxI2keV_4,B_FluxI2keV_5,b2_FluxI2keV_2,B2_FluxI2keV_3,
     +  LumIBFS_1,LumIMedS_2,b_LumIS_4,B_LumIS_5,b2_LumIS_2,
     +  B2_LumIS_3,LumIBF2keV_1,LumIMed2keV_2,b_LumI2keV_4,
     +  B_LumI2keV_5,b2_LumI2keV_2,B2_LumI2keV_3,dGmMed,b_dGm,B_dGm_1,
     +  dGmMean,s_dGm,GammaMed_1,b_Gamma_2,B_Gamma_3,GammaMean_1,
     +  s_Gamma_1,lognHMed_2,b_lognH_4,B_lognH_5,lognHMean_2,
     +  s_lognH_2,logPowNMed_1,b_logPowN_2,B_logPowN_3,logPowNMean_1,
     +  s_logPowN_1,logFracMed,b_logFrac,B_logFrac_1,logFracMean,
     +  s_logFrac,logBkgNMed_2,b_logBkgN_4,B_logBkgN_5,logBkgNMean_2,
     +  s_logBkgN_2,chi25_2,dof25_2,logZ_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Powerlaw+blackbody results

C  Format for file interpretation

    6 format(
     +  I5,1X,E9.2,1X,F12.9,1X,F11.9,1X,F9.7,1X,F11.9,1X,F9.6,1X,F9.6,
     +  1X,F11.9,1X,F10.7,1X,F10.7,1X,F12.10,1X,F11.8,1X,F11.8,1X,
     +  F12.10,1X,F10.7,1X,F10.7,1X,F11.8,1X,F9.6,1X,F10.7,1X,F11.8,
     +  1X,F10.7,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E14.8,1X,
     +  E14.8,1X,E14.8,1X,E14.8,1X,E14.8,1X,E14.8,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F10.6,1X,F10.6,1X,F10.6,1X,
     +  F10.6,1X,F10.6,1X,F10.6,1X,F10.8,1X,F10.8,1X,F9.7,1X,F10.8,1X,
     +  F11.9,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F11.9,1X,F10.7,1X,
     +  F10.7,1X,F10.7,1X,F10.7,1X,F12.10,1X,F11.8,1X,F11.8,1X,F11.8,
     +  1X,F11.8,1X,F12.10,1X,F10.7,1X,F10.7,1X,F10.7,1X,F10.7,1X,
     +  F11.9,1X,E16.9,1X,F14.11,1X,E15.8,1X,E16.9,1X,F12.10,1X,F11.6,
     +  1X,I3,1X,F11.6,1X,F15.11)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,27910
        read(1,'(A1217)')ar__5
        read(ar__5,6)
     +  ID_SRC_5,z_3,GammaKL_2,b_GammaH_4,B_GammaH_5,lognHKL_3,
     +  b_lognHH_6,B_lognHH_7,logPowNKL_2,b_logPowNH_4,B_logPowNH_5,
     +  logkTKL_1,b_logkTH_2,B_logkTH_3,logBBNKL,b_logBBNH,
     +  B_logBBNH_1,GammaBF_2,lognHBF_3,logPowNBF_2,logkTBF_1,
     +  logBBNBF,FluxObsBFS_3,FluxObsMedS_5,b_FluxObsS_10,
     +  B_FluxObsS_11,b2_FluxObsS_6,B2_FluxObsS_7,FluxObsBFT_3,
     +  FluxObsMedT_5,b_FluxObsT_10,B_FluxObsT_11,b2_FluxObsT_6,
     +  B2_FluxObsT_7,FluxCBFS_2,FluxCMedS_3,b_FluxCS_6,B_FluxCS_7,
     +  b2_FluxCS_4,B2_FluxCS_5,FluxCBFT_2,FluxCMedT_3,b_FluxCT_6,
     +  B_FluxCT_7,b2_FluxCT_4,B2_FluxCT_5,FluxIBFS_2,FluxIMedS_3,
     +  b_FluxIS_6,B_FluxIS_7,b2_FluxIS_4,B2_FluxIS_5,FluxIBF2keV_2,
     +  FluxIMed2keV_3,b_FluxI2keV_6,B_FluxI2keV_7,b2_FluxI2keV_4,
     +  B2_FluxI2keV_5,LumIBFS_2,LumIMedS_3,b_LumIS_6,B_LumIS_7,
     +  b2_LumIS_4,B2_LumIS_5,LumIBF2keV_2,LumIMed2keV_3,b_LumI2keV_6,
     +  B_LumI2keV_7,b2_LumI2keV_4,B2_LumI2keV_5,GammaMed_2,b_Gamma_4,
     +  B_Gamma_5,GammaMean_2,s_Gamma_2,lognHMed_3,b_lognH_6,
     +  B_lognH_7,lognHMean_3,s_lognH_3,logPowNMed_2,b_logPowN_4,
     +  B_logPowN_5,logPowNMean_2,s_logPowN_2,logkTMed_1,b_logkT_2,
     +  B_logkT_3,logkTMean_1,s_logkT_1,logBBNMed,b_logBBN,B_logBBN_1,
     +  logBBNMean,s_logBBN,logBkgNMed_3,b_logBkgN_6,B_logBkgN_7,
     +  logBkgNMean_3,s_logBkgN_3,chi25_3,dof25_3,logZ_3,logRbb
c    ..............Just test output...........
        write(6,6)
     +  ID_SRC_5,z_3,GammaKL_2,b_GammaH_4,B_GammaH_5,lognHKL_3,
     +  b_lognHH_6,B_lognHH_7,logPowNKL_2,b_logPowNH_4,B_logPowNH_5,
     +  logkTKL_1,b_logkTH_2,B_logkTH_3,logBBNKL,b_logBBNH,
     +  B_logBBNH_1,GammaBF_2,lognHBF_3,logPowNBF_2,logkTBF_1,
     +  logBBNBF,FluxObsBFS_3,FluxObsMedS_5,b_FluxObsS_10,
     +  B_FluxObsS_11,b2_FluxObsS_6,B2_FluxObsS_7,FluxObsBFT_3,
     +  FluxObsMedT_5,b_FluxObsT_10,B_FluxObsT_11,b2_FluxObsT_6,
     +  B2_FluxObsT_7,FluxCBFS_2,FluxCMedS_3,b_FluxCS_6,B_FluxCS_7,
     +  b2_FluxCS_4,B2_FluxCS_5,FluxCBFT_2,FluxCMedT_3,b_FluxCT_6,
     +  B_FluxCT_7,b2_FluxCT_4,B2_FluxCT_5,FluxIBFS_2,FluxIMedS_3,
     +  b_FluxIS_6,B_FluxIS_7,b2_FluxIS_4,B2_FluxIS_5,FluxIBF2keV_2,
     +  FluxIMed2keV_3,b_FluxI2keV_6,B_FluxI2keV_7,b2_FluxI2keV_4,
     +  B2_FluxI2keV_5,LumIBFS_2,LumIMedS_3,b_LumIS_6,B_LumIS_7,
     +  b2_LumIS_4,B2_LumIS_5,LumIBF2keV_2,LumIMed2keV_3,b_LumI2keV_6,
     +  B_LumI2keV_7,b2_LumI2keV_4,B2_LumI2keV_5,GammaMed_2,b_Gamma_4,
     +  B_Gamma_5,GammaMean_2,s_Gamma_2,lognHMed_3,b_lognH_6,
     +  B_lognH_7,lognHMean_3,s_lognH_3,logPowNMed_2,b_logPowN_4,
     +  B_logPowN_5,logPowNMean_2,s_logPowN_2,logkTMed_1,b_logkT_2,
     +  B_logkT_3,logkTMean_1,s_logkT_1,logBBNMed,b_logBBN,B_logBBN_1,
     +  logBBNMean,s_logBBN,logBkgNMed_3,b_logBkgN_6,B_logBkgN_7,
     +  logBkgNMean_3,s_logBkgN_3,chi25_3,dof25_3,logZ_3,logRbb
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Gamma-fixed-powerlaw results

C  Format for file interpretation

    7 format(
     +  I5,1X,E9.2,1X,F11.9,1X,F9.6,1X,F9.6,1X,F10.8,1X,F10.7,1X,
     +  F10.7,1X,F9.6,1X,F10.7,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,F8.4,1X,
     +  F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
     +  1X,F8.4,1X,F8.4,1X,F8.4,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F12.10,1X,F10.7,1X,F10.7,1X,F10.7,1X,F10.7,1X,F12.10,1X,E16.9,
     +  1X,E14.7,1X,E15.8,1X,E16.9,1X,F12.10,1X,F11.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,27910
        read(1,'(A841)')ar__6
        read(ar__6,7)
     +  ID_SRC_6,z_4,lognHKL_4,b_lognHH_8,B_lognHH_9,logPowNKL_3,
     +  b_logPowNH_6,B_logPowNH_7,lognHBF_4,logPowNBF_3,FluxObsBFS_4,
     +  FluxObsMedS_6,b_FluxObsS_12,B_FluxObsS_13,b2_FluxObsS_8,
     +  B2_FluxObsS_9,FluxObsBFT_4,FluxObsMedT_6,b_FluxObsT_12,
     +  B_FluxObsT_13,b2_FluxObsT_8,B2_FluxObsT_9,FluxCBFS_3,
     +  FluxCMedS_4,b_FluxCS_8,B_FluxCS_9,b2_FluxCS_6,B2_FluxCS_7,
     +  FluxCBFT_3,FluxCMedT_4,b_FluxCT_8,B_FluxCT_9,b2_FluxCT_6,
     +  B2_FluxCT_7,FluxIBFS_3,FluxIMedS_4,b_FluxIS_8,B_FluxIS_9,
     +  b2_FluxIS_6,B2_FluxIS_7,FluxIBF2keV_3,FluxIMed2keV_4,
     +  b_FluxI2keV_8,B_FluxI2keV_9,b2_FluxI2keV_6,B2_FluxI2keV_7,
     +  LumIBFS_3,LumIMedS_4,b_LumIS_8,B_LumIS_9,b2_LumIS_6,
     +  B2_LumIS_7,LumIBF2keV_3,LumIMed2keV_4,b_LumI2keV_8,
     +  B_LumI2keV_9,b2_LumI2keV_6,B2_LumI2keV_7,lognHMed_4,b_lognH_8,
     +  B_lognH_9,lognHMean_4,s_lognH_4,logPowNMed_3,b_logPowN_6,
     +  B_logPowN_7,logPowNMean_3,s_logPowN_3,logBkgNMed_4,
     +  b_logBkgN_8,B_logBkgN_9,logBkgNMean_4,s_logBkgN_4,logZ_4
c    ..............Just test output...........
        write(6,7)
     +  ID_SRC_6,z_4,lognHKL_4,b_lognHH_8,B_lognHH_9,logPowNKL_3,
     +  b_logPowNH_6,B_logPowNH_7,lognHBF_4,logPowNBF_3,FluxObsBFS_4,
     +  FluxObsMedS_6,b_FluxObsS_12,B_FluxObsS_13,b2_FluxObsS_8,
     +  B2_FluxObsS_9,FluxObsBFT_4,FluxObsMedT_6,b_FluxObsT_12,
     +  B_FluxObsT_13,b2_FluxObsT_8,B2_FluxObsT_9,FluxCBFS_3,
     +  FluxCMedS_4,b_FluxCS_8,B_FluxCS_9,b2_FluxCS_6,B2_FluxCS_7,
     +  FluxCBFT_3,FluxCMedT_4,b_FluxCT_8,B_FluxCT_9,b2_FluxCT_6,
     +  B2_FluxCT_7,FluxIBFS_3,FluxIMedS_4,b_FluxIS_8,B_FluxIS_9,
     +  b2_FluxIS_6,B2_FluxIS_7,FluxIBF2keV_3,FluxIMed2keV_4,
     +  b_FluxI2keV_8,B_FluxI2keV_9,b2_FluxI2keV_6,B2_FluxI2keV_7,
     +  LumIBFS_3,LumIMedS_4,b_LumIS_8,B_LumIS_9,b2_LumIS_6,
     +  B2_LumIS_7,LumIBF2keV_3,LumIMed2keV_4,b_LumI2keV_8,
     +  B_LumI2keV_9,b2_LumI2keV_6,B2_LumI2keV_7,lognHMed_4,b_lognH_8,
     +  B_lognH_9,lognHMean_4,s_lognH_4,logPowNMed_3,b_logPowN_6,
     +  B_logPowN_7,logPowNMean_3,s_logPowN_3,logBkgNMed_4,
     +  b_logBkgN_8,B_logBkgN_9,logBkgNMean_4,s_logBkgN_4,logZ_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Shape-fixed-powerlaw results

C  Format for file interpretation

    8 format(
     +  I5,1X,E9.2,1X,F10.8,1X,F10.7,1X,F10.7,1X,F10.7,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F10.7,
     +  1X,F10.7,1X,F10.7,1X,F10.7,1X,F12.10,1X,E15.9,1X,E15.8,1X,
     +  E15.8,1X,E16.9,1X,F12.10,1X,F11.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,27910
        read(1,'(A746)')ar__7
        read(ar__7,8)
     +  ID_SRC_7,z_5,logPowNKL_4,b_logPowNH_8,B_logPowNH_9,
     +  logPowNBF_4,FluxObsBFS_5,FluxObsMedS_7,b_FluxObsS_14,
     +  B_FluxObsS_15,b2_FluxObsS_10,B2_FluxObsS_11,FluxObsBFT_5,
     +  FluxObsMedT_7,b_FluxObsT_14,B_FluxObsT_15,b2_FluxObsT_10,
     +  B2_FluxObsT_11,FluxCBFS_4,FluxCMedS_5,b_FluxCS_10,B_FluxCS_11,
     +  b2_FluxCS_8,B2_FluxCS_9,FluxCBFT_4,FluxCMedT_5,b_FluxCT_10,
     +  B_FluxCT_11,b2_FluxCT_8,B2_FluxCT_9,FluxIBFS_4,FluxIMedS_5,
     +  b_FluxIS_10,B_FluxIS_11,b2_FluxIS_8,B2_FluxIS_9,FluxIBF2keV_4,
     +  FluxIMed2keV_5,b_FluxI2keV_10,B_FluxI2keV_11,b2_FluxI2keV_8,
     +  B2_FluxI2keV_9,LumIBFS_4,LumIMedS_5,b_LumIS_10,B_LumIS_11,
     +  b2_LumI2S,B2_LumIS_8,LumIBF2keV_4,LumIMed2keV_5,b_LumI2keV_10,
     +  B_LumI2keV_11,b2_LumI2keV_8,B2_LumI2keV_9,logPowNMed_4,
     +  b_logPowN_8,B_logPowN_9,logPowNMean_4,s_logPowN_4,
     +  logBkgNMed_5,b_logBkgN_10,B_logBkgN_11,logBkgNMean_5,
     +  s_logBkgN_5,logZ_5
c    ..............Just test output...........
        write(6,8)
     +  ID_SRC_7,z_5,logPowNKL_4,b_logPowNH_8,B_logPowNH_9,
     +  logPowNBF_4,FluxObsBFS_5,FluxObsMedS_7,b_FluxObsS_14,
     +  B_FluxObsS_15,b2_FluxObsS_10,B2_FluxObsS_11,FluxObsBFT_5,
     +  FluxObsMedT_7,b_FluxObsT_14,B_FluxObsT_15,b2_FluxObsT_10,
     +  B2_FluxObsT_11,FluxCBFS_4,FluxCMedS_5,b_FluxCS_10,B_FluxCS_11,
     +  b2_FluxCS_8,B2_FluxCS_9,FluxCBFT_4,FluxCMedT_5,b_FluxCT_10,
     +  B_FluxCT_11,b2_FluxCT_8,B2_FluxCT_9,FluxIBFS_4,FluxIMedS_5,
     +  b_FluxIS_10,B_FluxIS_11,b2_FluxIS_8,B2_FluxIS_9,FluxIBF2keV_4,
     +  FluxIMed2keV_5,b_FluxI2keV_10,B_FluxI2keV_11,b2_FluxI2keV_8,
     +  B2_FluxI2keV_9,LumIBFS_4,LumIMedS_5,b_LumIS_10,B_LumIS_11,
     +  b2_LumI2S,B2_LumIS_8,LumIBF2keV_4,LumIMed2keV_5,b_LumI2keV_10,
     +  B_LumI2keV_11,b2_LumI2keV_8,B2_LumI2keV_9,logPowNMed_4,
     +  b_logPowN_8,B_logPowN_9,logPowNMean_4,s_logPowN_4,
     +  logBkgNMed_5,b_logBkgN_10,B_logBkgN_11,logBkgNMean_5,
     +  s_logBkgN_5,logZ_5
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Soft-band fitting results

C  Format for file interpretation

    9 format(
     +  I5,1X,E9.2,1X,F12.9,1X,F10.8,1X,F9.7,1X,F11.9,1X,F9.6,1X,F9.6,
     +  1X,F11.9,1X,F11.8,1X,F13.10,1X,F11.8,1X,F9.6,1X,E15.8,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,F10.8,
     +  1X,F10.8,1X,F9.7,1X,F9.7,1X,F11.9,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F11.9,1X,F11.8,1X,F11.8,1X,F12.9,1X,F11.8,1X,F12.10,
     +  1X,E16.9,1X,E15.8,1X,E16.9,1X,E15.8,1X,F12.10,1X,F12.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,27910
        read(1,'(A489)')ar__8
        read(ar__8,9)
     +  ID_SRC_8,z_6,GammaKL_3,b_GammaH_6,B_GammaH_7,lognHKL_5,
     +  b_lognHH_10,B_lognHH_11,logPowNKL_5,b_logPowNH_10,
     +  B_logPowNH_11,GammaBF_3,lognHBF_5,logPowNBF_5,FluxObsBFS_6,
     +  FluxObsMedS_8,b_FluxObsS_16,B_FluxObsS_17,b2_FluxObsS_12,
     +  B2_FluxObsS_13,GammaMed_3,b_Gamma_6,B_Gamma_7,GammaMean_3,
     +  s_Gamma_3,lognHMed_5,b_lognH_10,B_lognH_11,lognHMean_5,
     +  s_lognH_5,logPowNMed_5,b_logPowN_10,B_logPowN_11,
     +  logPowNMean_5,s_logPowN_5,logBkgNMed_6,blogBkgN,BlogBkgN_1,
     +  logBkgNMean_6,logBkgN,logZ_6
c    ..............Just test output...........
        write(6,9)
     +  ID_SRC_8,z_6,GammaKL_3,b_GammaH_6,B_GammaH_7,lognHKL_5,
     +  b_lognHH_10,B_lognHH_11,logPowNKL_5,b_logPowNH_10,
     +  B_logPowNH_11,GammaBF_3,lognHBF_5,logPowNBF_5,FluxObsBFS_6,
     +  FluxObsMedS_8,b_FluxObsS_16,B_FluxObsS_17,b2_FluxObsS_12,
     +  B2_FluxObsS_13,GammaMed_3,b_Gamma_6,B_Gamma_7,GammaMean_3,
     +  s_Gamma_3,lognHMed_5,b_lognH_10,B_lognH_11,lognHMean_5,
     +  s_lognH_5,logPowNMed_5,b_logPowN_10,B_logPowN_11,
     +  logPowNMean_5,s_logPowN_5,logBkgNMed_6,blogBkgN,BlogBkgN_1,
     +  logBkgNMean_6,logBkgN,logZ_6
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Hard-band fitting results

C  Format for file interpretation

   10 format(
     +  I5,1X,E9.2,1X,F11.9,1X,F9.7,1X,F9.7,1X,F12.10,1X,F9.6,1X,F9.6,
     +  1X,F11.9,1X,F10.7,1X,F12.9,1X,F11.8,1X,F9.6,1X,F13.10,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,F9.7,1X,F9.7,
     +  1X,F9.7,1X,F9.7,1X,F10.8,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F10.8,1X,F10.7,1X,F10.7,1X,F11.8,1X,F10.7,1X,F10.8,1X,E15.9,
     +  1X,E15.8,1X,E15.8,1X,E15.8,1X,F12.10,1X,F11.7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,27910
        read(1,'(A543)')ar__9
        read(ar__9,10)
     +  ID_SRC_9,z_7,GammaKL_4,e_Gamma,E_Gamma_1,lognHKL_6,b_lognH_12,
     +  B_lognH_13,logPowNKL_6,b_logPowN_12,B_logPowN_13,GammaBF_4,
     +  lognHBF_6,logPowNBF_6,FluxObsBFT_6,FluxObsMedT_8,
     +  b_FluxObsT_16,B_FluxObsT_17,b2_FluxObsT_12,B2_FluxObsT_13,
     +  FluxCBFT_5,FluxCMedT_6,b_FluxCMedT,B_FluxCMedT_1,b2_FluxCMedT,
     +  B2_FluxCMedT_1,GammaMed_4,b_GammaMed,B_GammaMed_1,GammaMean_4,
     +  s_Gamma_4,lognHMed_6,b_lognHMed,B_lognHMed_1,lognHMean_6,
     +  s_lognH_6,logPowNMed_6,b_logPowNMed,B_logPowNMed_1,
     +  logPowNMean_6,s_logPowN_6,logBkgNMed_7,b_logBkgNMed,
     +  B_logBkgNMed_1,logBkgNMean_7,s_logBkgN_6,logZ_7
c    ..............Just test output...........
        write(6,10)
     +  ID_SRC_9,z_7,GammaKL_4,e_Gamma,E_Gamma_1,lognHKL_6,b_lognH_12,
     +  B_lognH_13,logPowNKL_6,b_logPowN_12,B_logPowN_13,GammaBF_4,
     +  lognHBF_6,logPowNBF_6,FluxObsBFT_6,FluxObsMedT_8,
     +  b_FluxObsT_16,B_FluxObsT_17,b2_FluxObsT_12,B2_FluxObsT_13,
     +  FluxCBFT_5,FluxCMedT_6,b_FluxCMedT,B_FluxCMedT_1,b2_FluxCMedT,
     +  B2_FluxCMedT_1,GammaMed_4,b_GammaMed,B_GammaMed_1,GammaMean_4,
     +  s_Gamma_4,lognHMed_6,b_lognHMed,B_lognHMed_1,lognHMean_6,
     +  s_lognH_6,logPowNMed_6,b_logPowNMed,B_logPowNMed_1,
     +  logPowNMean_6,s_logPowN_6,logBkgNMed_7,b_logBkgNMed,
     +  B_logBkgNMed_1,logBkgNMean_7,s_logBkgN_6,logZ_7
c    .......End.of.Just test output...........
      end do
      close(1)

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