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 loading all data files into arrays.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-May-18
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     (nr__) ! [1/32683] ID of the sources in the eFEDS
*                                         main X-ray catalog (Paper I) (ID_SRC)
      real*8        RAcdeg     (nr__) ! (deg) [126.24/145.59] Right ascension (J2000),
*                                         astrometric corrected (Paper I)
*                                         (RA_CORR)
      real*4        DEcdeg     (nr__) ! (deg) [-2.91/5.97] Declination (J2000),
*                                         astrometric corrected (Paper I)
*                                         (DEC_CORR)
      real*8        DetLike    (nr__) ! 0.2-2.3 keV source detection likelihood
*                                         (Paper I) (DET_LIKE)
      real*4        inArea90   (nr__) ! Whether located inside the inner
*                                         90%-area region of eFEDS (Paper I)
*                                         (inArea90)
      character*11  CTPLS8     (nr__) ! ID of the best LS8 counterpart
*                                         (Paper II) (CTP_LS8_UNIQUE_OBJID)
      real*8        RALdeg     (nr__) ! (deg) [126.23/145.59] Right ascension (J2000)
*                                         of the best LS8 counterpart (Paper II)
*                                         (CTP_LS8_RA)
      real*4        DELdeg     (nr__) ! (deg) [-2.91/5.97] Declination (J2000) of the
*                                         best LS8 counterpart (Paper II)
*                                         (CTP_LS8_DEC)
      integer*4     CTPqual    (nr__) ! [2/4] Counterpart quality (Paper II)
*                                         A value >=2 is recommended
*                                         (CTP_quality)
      integer*4     CTPClass   (nr__) ! [2/3] Classification of the optical
*                                         counterpart (Paper II) (CTP_CLASS) (1)
      real*8        CTPz       (nr__) ! [0.0/8.0] Redshift of the optical
*                                         counterpart (Paper II) (CTP_REDSHIFT)
      integer*4     CTPzgrade  (nr__) ! [2/5] Redshift Grade (Paper II)
*                                         A value >=3 is recommended
*                                         (CTP_REDSHIFT_GRADE)
      integer*4     inKiDSFlag (nr__) ! [0/1]  Whether located inside the region
*                                         of the KiDS survey (Paper II)
*                                         (in_KiDS_flag)
      integer*4     LxModel    (nr__) ! [1/5] Index of selected model for
*                                         luminosity measurement (LxModel) (2)
      integer*4     NHclass    (nr__) ! [1/4] Class of AGN NH measurement with
*                                         model 1 (single-powerlaw) (NHclass) (3)
      integer*4     FSModel    (nr__) ! [3/6] Index of selected model for
*                                         0.5-2keV flux (FSModel) (4)
      integer*4     FHModel    (nr__) ! [5/7] Index of selected model for
*                                         2.3-5keV flux (FHModel) (5)
      real*4        galNH      (nr__) ! (cm-2) Total column density of Galactic
*                                         absorption (galNH)
      real*4        galNHI     (nr__) ! (cm-2) HI column density from HI4PI (galNHI)
      real*8        SrcCts     (nr__) ! (ct) Source net counts in the 0.2-5 keV band
*                                         (SrcCts)
      real*4        FluxcMedS  (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, posterior median
*                                         (FluxCorr_Med_s)
      real*4        b_FluxcS   (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, 1sigma lower limit
*                                         (FluxCorr_Lo1_s)
      real*4        B_FluxcS_1 (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         0.5-2 keV, 1sigma upper limit
*                                         (FluxCorr_Up1_s)
      real*4        FluxcMedT  (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, posterior median
*                                         (FluxCorr_Med_t)
      real*4        b_FluxcT   (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, 1sigma lower limit
*                                         (FluxCorr_Lo1_t)
      real*4        B_FluxcT_1 (nr__) ! (mW/m2) Absorption corrected flux in observed
*                                         2.3-5 keV, 1sigma upper limit
*                                         (FluxCorr_Up1_t)
      real*4        FluxIMedS  (nr__) ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2 keV, posterior median
*                                         (FluxIntr_Med_s)
      real*4        b_FluxIS   (nr__) ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2keV, 1sigma lower limit
*                                         (FluxIntr_Lo1_s)
      real*4        B_FluxIS_1 (nr__) ! (mW/m2) Absorption corrected flux in rest-frame
*                                         0.5-2keV, 1sigma upper limit
*                                         (FluxIntr_Up1_s)
      real*4        FluxIMed2keV(nr__) ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, median (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV(nr__) ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, 1sigma lower limit
*                                         (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_1(nr__) ! (mW/m2/eV) Absorption corrected flux at rest-frame
*                                         2keV, 1sigma upper limit
*                                         (FluxIntr_Up1_2keV)
      real*8        LumIMedS   (nr__) ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, posterior median
*                                         (LumiIntr_Med_s)
      real*8        b_LumIS    (nr__) ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, 1sigma lower limit
*                                         (LumiIntr_Lo1_s)
      real*8        B_LumIS_1  (nr__) ! (10-7W) Intrinsic luminosity in rest-frame
*                                         0.5-2 keV, 1sigma upper limit
*                                         (LumiIntr_Up1_s)
      real*8        LumIMed2keV(nr__) ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         posterior median in erg/s/eV
*                                         (LumiIntr_Med_2keV)
      real*8        b_LumI2keV (nr__) ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         1sigma lower limit in erg/s/eV
*                                         (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_1(nr__) ! (10-7W/eV) Intrinsic luminosity at rest-frame 2keV,
*                                         1sigma upper limit in erg/s/eV
*                                         (LumiIntr_Up1_2keV)
      real*4        FluxObsMedS(nr__) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS (nr__) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_1(nr__) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                         1sigma upper limit (FluxObsv_Up1_s)
      real*4        FluxObsMedT(nr__) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT (nr__) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_1(nr__) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                         1sigma upper limit (FluxObsv_Up1_t)
      real*8        lognHKLm1  (nr__) ! log AGN column density, KL divergence
*                                         (in nats units) (lognH_KL_m1)
      real*8        b_lognHm1  (nr__) ! ([cm-2]) log AGN column density, HDI lower limit
*                                         (lognH_HLo_m1)
      real*8        B_lognHm1_1(nr__) ! ([cm-2]) log AGN column density, HDI upper limit
*                                         (lognH_HUp_m1)
      real*8        lognHMedm1 (nr__) ! ([cm-2]) log AGN column density, posterior median
*                                         (lognH_Med_m1)
      real*8        GammaKLm3  (nr__) ! primary power-law slope in model 3,
*                                         KL divergence (in nats units)
*                                         (Gamma_KL_m3)
      real*8        b_GammaHm3 (nr__) ! primary power-law slope in model 3,
*                                         HDI lower limit (Gamma_HLo_m3)
      real*8        B_GammaHm3_1(nr__) ! primary power-law slope in model 3,
*                                         HDI upper limit (Gamma_HUp_m3)
      real*8        GammaMedm3 (nr__) ! primary power-law slope in model 3,
*                                         posterior median (Gamma_Med_m3)
      real*8        logZm0     (nr__) ! ([-]) log Bayesian evidence with model 0:
*                                         APEC (logZ_m0)
      real*8        logZm1     (nr__) ! ([-]) log Bayesian evidence with model 1:
*                                         single powerlaw (logZ_m1)
      real*8        logZm2     (nr__) ! ([-]) log Bayesian evidence with model 2:
*                                         double powerlaw (logZ_m2)
      real*8        logZm3     (nr__) ! ([-]) log Bayesian evidence with model 3:
*                                         powerlaw + blackbody (logZ_m3)
      real*8        logZm4     (nr__) ! ([-]) log Bayesian evidence with model 4:
*                                         powerlaw with Gamma fixed at 2.0
*                                         (logZ_m4)
      real*8        logZm5     (nr__) ! ([-]) log Bayesian evidence with model 5:
*                                         shape-fixed powerlaw (logZ_m5)
      real*8        L2500      (nr__) ! (10-7W/Hz) ?=-99.0 The rest-frame 2500 luminosity
*                                          in erg/s/Hz (L2500)
      real*8        L5100      (nr__) ! (10-7W/Hz) ?=-99.0 The rest-frame 5100 luminosity
*                                         in erg/s/Hz (L5100)
      real*8        W1mag      (nr__) ! (mag) ? LS8-WISE W1 AB magnitude (Paper II)
*                                         (W1)
      real*4        e_W1mag    (nr__) ! (mag) ? LS8-WISE W1 magnitude error (Paper II)
*                                         (W1_ERR)
      real*8        W2mag      (nr__) ! (mag) ? LS8-WISE W2 AB magnitude (Paper II)
*                                         (W2)
      real*4        e_W2mag    (nr__) ! (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   (nr__1) ! [1/32683] ID of the sources in the eFEDS
*                                       main X-ray catalog (Paper I) (ID_SRC)
      real*8        RAcdeg_1   (nr__1) ! (deg) [126.24/145.66] Right ascension (J2000),
*                                       astrometric corrected (Paper I) (RA_CORR)
      real*4        DEcdeg_1   (nr__1) ! (deg) [-2.91/5.97] Declination (J2000),
*                                       astrometric corrected (Paper I)
*                                       (DEC_CORR)
      real*8        DetLike_1  (nr__1) ! [6.0/33254.59] 0.2-2.3 keV source
*                                       detection likelihood (Paper I) (DET_LIKE)
      real*8        ExtLike    (nr__1) ! X-ray source extent likelihood (Paper I)
*                                       (EXT_LIKE)
      real*4        inArea90_1 (nr__1) ! Whether located inside the inner 90%-area
*                                       region of eFEDS (Paper I) (inArea90)
      real*4        galNH_1    (nr__1) ! (cm-2) Total column density of Galactic
*                                       absorption (galNH)
      real*4        galNHI_1   (nr__1) ! (cm-2) HI column density from HI4PI (galNHI)
      real*8        ExpTime    (nr__1) ! (s) Spectra exposure time (Exposure)
      real*8        SrcCts_1   (nr__1) ! (ct) Source net counts in the 0.2-5 keV band
*                                       (SrcCts)
      real*8        RAdeg      (nr__1) ! (deg) [126.24/145.66] Right ascension (J2000)
*                                       before astrometric correction (Paper I)
*                                       (RA)
      real*4        DEdeg      (nr__1) ! (deg) [-2.91/5.97] Declination (J2000) before
*                                       astrometric correction (Paper I) (DEC)
      integer*4     Rad        (nr__1) ! (arcsec) [12/259] Source extraction radius (Rad)
      integer*4     Ann1       (nr__1) ! (arcsec) [22/515] Inner radius of background
*                                       extraction region (Ann1)
      integer*4     Ann2       (nr__1) ! (arcsec) [171/900] Outer radius of background
*                                       extraction region (Ann2)
      real*4        RegArea    (nr__1) ! (deg2) Geometry area of source extraction region
*                                       (REGAREA_s)
      real*4        Backscal   (nr__1) ! (deg2) Source BACKSCAL (Backscal_s)
      real*8        BackscalBkg(nr__1) ! (deg2) Background BACKSCAL (Backscal_b)
      integer*4     FSModel_1  (nr__1) ! [0/6] Index of selected model for
*                                       0.5-2keV flux (FSModel) (1)
      integer*4     FHModel_1  (nr__1) ! [5/7] Index of selected model for
*                                       2.3-5keV flux (FHModel) (2)
      real*4        Rate0_2_2_3(nr__1) ! (ct/s) Net count rate in 0.2-2.3 keV
*                                       (Rate_d2_2d3)
      real*8        e_Rate0_2_2_3(nr__1) ! (ct/s) Net count rate error in 0.2-2.3 keV
*                                       (RateErr_d2_2d3)
      real*4        Rate0_2_0_5(nr__1) ! (ct/s) Net count rate in 0.2-0.5 keV (Rate_d2_d5)
      real*4        e_Rate0_2_0_5(nr__1) ! (ct/s) Net count rate error in 0.2-0.5 keV
*                                       (RateErr_d2_d5)
      real*4        Rate0_5_1  (nr__1) ! (ct/s) Net count rate in 0.5-1 keV (Rate_d5_1)
      real*4        e_Rate0_5_1(nr__1) ! (ct/s) Net count rate error in 0.5-1 keV
*                                       (RateErr_d5_1)
      real*4        Rate1_2    (nr__1) ! (ct/s) Net count rate in 1-2 keV (Rate_1_2)
      real*4        e_Rate1_2  (nr__1) ! (ct/s) Net count rate error in 1-2 keV
*                                       (RateErr_1_2)
      real*4        Rate2_4_5  (nr__1) ! (ct/s) Net count rate in 2-4.5 keV (Rate_2_4d5)
      real*4        e_Rate2_4_5(nr__1) ! (ct/s) Net count rate error in 2-4.5 keV
*                                       (RateErr_2_4d5)
      real*4        Rate2_3_5  (nr__1) ! (ct/s) Net count rate in 2.3-5 keV (Rate_2d3_5)
      real*4        e_Rate2_3_5(nr__1) ! (ct/s) Net count rate error in 2.3-5 keV
*                                       (RateErr_2d3_5)
      real*4        Rate5_8    (nr__1) ! (ct/s) Net count rate in 5-8 keV (Rate_5_8)
      real*4        e_Rate5_8  (nr__1) ! (ct/s) Net count rate error in 5-8 keV
*                                       (RateErr_5_8)
      real*4        BkgCts0_2_0_6(nr__1) ! (ct) Background counts in 0.2-0.6 keV
*                                       (BkgCts_d2_d6)
      real*8        e_BkgCts0_2_0_6(nr__1) ! (ct) Background counts error in 0.2-0.6 keV
*                                       (BkgCtsErr_d2_d6)
      real*4        BkgCts0_6_2_3(nr__1) ! (ct) Background counts in 0.6-2.3 keV
*                                       (BkgCts_d6_2d3)
      real*8        e_BkgCts0_6_2_3(nr__1) ! (ct) Background counts error in 0.6-2.3 keV
*                                       (BkgCtsErr_d6_2d3)
      real*4        BkgCts2_3_5(nr__1) ! (ct) Background counts in 2.3-5 keV
*                                       (BkgCts_2d3_5)
      real*8        e_BkgCts2_3_5(nr__1) ! (ct) Background counts error in 2.3-5 keV
*                                       (BkgCtsErr_2d3_5)
      real*4        BkgCts5_8  (nr__1) ! (ct) Background counts in 5-8 keV (BkgCts_5_8)
      real*8        e_BkgCts5_8(nr__1) ! (ct) Background counts error in 5-8 keV
*                                       (BkgCtsErr_5_8)
      real*8        PSFCor0_2_0_5(nr__1) ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 0.2-0.5 keV (PSFCor_d2_d5)
      real*8        PSFCor0_5_1(nr__1) ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 0.5-1 keV (PSFCor_d5_1)
      real*8        PSFCor1_2  (nr__1) ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 1-2 keV (PSFCor_1_2)
      real*8        PSFCor2_3_5(nr__1) ! PSF-loss correction (ARF CORRPSF) averaged
*                                       in 2.3-5 keV (PSFCor_2d3_5)
      integer*4     Nempty     (nr__1) ! [35/771] Number of empty channels between
*                                       channel 20 and 900 (Nempty)
      real*4        FluxObsMedS_1(nr__1) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_2(nr__1) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_3(nr__1) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                       1sigma upper limit (FluxObsv_Up1_s)
      real*4        FluxObsMedT_1(nr__1) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                       posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_2(nr__1) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                       1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_3(nr__1) ! (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   (nr__2) ! [1/32683] ID of the sources in the eFEDS
*                                        main X-ray catalog (Paper I) (ID_SRC)
      real*4        z          (nr__2) ! The redshift adopted for the spectral
*                                        model (Redshift)
      real*8        lognHKL    (nr__2) ! log AGN column density, KL divergence
*                                        (in nats units) (lognH_KL)
      real*8        b_lognHH   (nr__2) ! ([cm-2]) log AGN column density, HDI lower limit
*                                        (lognH_HLo)
      real*8        B_lognHH_1 (nr__2) ! ([cm-2]) log AGN column density, HDI upper limit
*                                        (lognH_HUp)
      real*8        logkTKL    (nr__2) ! log temperature, KL divergence
*                                        (in nats units) (logkT_KL)
      real*4        b_logkTH   (nr__2) ! ([keV]) log temperature, HDI lower limit
*                                        (logkT_HLo)
      real*4        B_logkTH_1 (nr__2) ! ([keV]) log temperature, HDI upper limit
*                                        (logkT_HUp)
      real*8        logAbundKL (nr__2) ! log abundance, KL divergence
*                                        (in nats units) (logAbundanc_KL)
      real*8        b_logAbunH (nr__2) ! log abundance, HDI lower limit
*                                        (logAbundanc_HLo)
      real*4        B_logAbunH_1(nr__2) ! log abundance, HDI upper limit
*                                        (logAbundanc_HUp)
      real*8        logApecNormKL(nr__2) ! log APEC normalization, KL divergence
*                                        (in nats units) (logApecNorm_KL)
      real*8        b_logApecNormH(nr__2) ! ([cm-5]) log APEC normalization, HDI lower limit
*                                        (logApecNorm_HLo)
      real*8        B_logApecNormH_1(nr__2) ! ([cm-5]) log APEC normalization, HDI upper limit
*                                        (logApecNorm_HUp)
      real*8        lognHBF    (nr__2) ! ([cm-2]) log AGN column density, best-fit model
*                                        (lognH_BF)
      real*4        logkTBF    (nr__2) ! ([keV]) log temperature, best-fit model
*                                        (logkT_BF)
      real*4        logAbundBF (nr__2) ! log abundance, best-fit model
*                                        (logAbundanc_BF)
      real*8        logApecNBF (nr__2) ! ([cm-5]) log APEC normalization, best-fit model
*                                        (logApecNorm_BF)
      real*4        FluxObsBFS (nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        best-fit model (FluxObsv_BF_s)
      real*4        FluxObsMedS_2(nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_4(nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_5(nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        1sigma upper limit (FluxObsv_Up1_s)
      real*4        b2_FluxObsS(nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        2sigma lower limit (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_1(nr__2) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                        2sigma upper limit (FluxObsv_Up2_s)
      real*4        FluxObsBFT (nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_2(nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_4(nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_5(nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        1sigma upper limit (FluxObsv_Up1_t)
      real*4        b2_FluxObsT(nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        2sigma lower limit (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_1(nr__2) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                        2sigma upper limit (FluxObsv_Up2_t)
      real*8        lognHMed   (nr__2) ! ([cm-2]) log AGN column density, posterior median
*                                        (lognH_Med)
      real*8        b_lognH    (nr__2) ! ([cm-2]) log AGN column density, percentile
*                                        lower limit (lognH_Lo)
      real*8        B_lognH_1  (nr__2) ! ([cm-2]) log AGN column density, percentile
*                                        upper limit (lognH_Up)
      real*8        lognHMean  (nr__2) ! ([cm-2]) log AGN column density, mean (lognH_Mean)
      real*8        s_lognH    (nr__2) ! ([cm-2]) log AGN column density,
*                                        standard deviation (lognH_Std)
      real*4        logkTMed   (nr__2) ! ([keV]) log temperature, posterior median
*                                        (logkT_Med)
      real*4        b_logkT    (nr__2) ! ([keV]) log temperature, percentile lower limit
*                                        (logkT_Lo)
      real*4        B_logkT_1  (nr__2) ! ([keV]) log temperature, percentile upper limit
*                                        (logkT_Up)
      real*4        logkTMean  (nr__2) ! ([keV]) log temperature, mean (logkT_Mean)
      real*8        s_logkT    (nr__2) ! ([keV]) log temperature, standard deviation
*                                        (logkT_Std)
      real*8        logAbundMed(nr__2) ! log abundance, posterior median
*                                        (logAbundanc_Med)
      real*8        b_logAbund (nr__2) ! log abundance, percentile lower limit
*                                        (logAbundanc_Lo)
      real*8        B_logAbund_1(nr__2) ! log abundance, percentile upper limit
*                                        (logAbundanc_Up)
      real*8        logAbundMean(nr__2) ! log abundance, mean (logAbundanc_Mean)
      real*8        s_logAbund (nr__2) ! log abundance, standard deviation
*                                        (logAbundanc_Std)
      real*8        logApecNMed(nr__2) ! ([cm-5]) log APEC normalization, posterior median
*                                        (logApecNorm_Med)
      real*8        b_logApecN (nr__2) ! ([cm-5]) log APEC normalization, percentile
*                                        lower limit (logApecNorm_Lo)
      real*8        B_logApecN_1(nr__2) ! ([cm-5]) log APEC normalization, percentile
*                                        upper limit (logApecNorm_Up)
      real*8        logApecNMean(nr__2) ! ([cm-5]) log APEC normalization, mean
*                                        (logApecNorm_Mean)
      real*8        s_logApecN (nr__2) ! ([cm-5]) log APEC normalization,
*                                        standard deviation (logApecNorm_Std)
      real*4        logBkgNMed (nr__2) ! log background normalization,
*                                        posterior median (logBkgNorm_Med)
      real*4        b_logBkgN  (nr__2) ! log background normalization,
*                                        percentile lower limit (logBkgNorm_Lo)
      real*4        B_logBkgN_1(nr__2) ! log background normalization,
*                                        percentile upper limit (logBkgNorm_Up)
      real*4        logBkgNMean(nr__2) ! log background normalization,
*                                        mean (logBkgNorm_Mean)
      real*8        s_logBkgN  (nr__2) ! log background normalization,
*                                        standard deviation (logBkgNorm_Std)
      real*8        chi25      (nr__2) ! ?=-99 Chi^2 of the best-fit model against
*                                        the rebinned data (>=25 cts/bin) (chi25)
      integer*4     dof25      (nr__2) ! ?=-99 DOF of the rebinned data, which
*                                        has at least 25 counts per bin (dof25)
      real*8        logZ       (nr__2) ! 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   (nr__3) ! [1/32683] ID of the sources in the
*                                             eFEDS main X-ray catalog (Paper I)
*                                             (ID_SRC)
      real*4        z_1        (nr__3) ! [-0.01/8.0] The redshift adopted for
*                                             the spectral model (Redshift)
      real*8        GammaKL    (nr__3) ! power-law slope, KL divergence
*                                             (in nats units) (Gamma_KL)
      real*8        b_GammaH   (nr__3) ! power-law slope, HDI lower limit
*                                             (Gamma_HLo)
      real*8        B_GammaH_1 (nr__3) ! power-law slope, HDI upper limit
*                                             (Gamma_HUp)
      real*8        lognHKL_1  (nr__3) ! log AGN column density,
*                                             KL divergence (in nats units)
*                                             (lognH_KL)
      real*8        b_lognHH_2 (nr__3) ! ([cm-2]) log AGN column density,
*                                             HDI lower limit (lognH_HLo)
      real*8        B_lognHH_3 (nr__3) ! ([cm-2]) log AGN column density,
*                                             HDI upper limit (lognH_HUp)
      real*8        logPowNKL  (nr__3) ! log power-law normalization,
*                                             KL divergence (in nats units)
*                                             (logPowNorm_KL)
      real*8        b_logPowNH (nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_1(nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF    (nr__3) ! power-law slope, best-fit model
*                                             (Gamma_BF)
      real*8        lognHBF_1  (nr__3) ! ([cm-2]) log AGN column density,
*                                             best-fit model (lognH_BF)
      real*8        logPowNBF  (nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_1(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             best-fit model (FluxObsv_BF_s)
      real*4        FluxObsMedS_3(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             posterior median (FluxObsv_Med_s)
      real*4        b_FluxObsS_6(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             1sigma lower limit (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_7(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             1sigma upper limit (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_2(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             2sigma lower limit (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_3(nr__3) ! (mW/m2) Observed flux in observed 0.5-2 keV,
*                                             2sigma upper limit (FluxObsv_Up2_s)
      real*4        FluxObsBFT_1(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_3(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_6(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             1sigma lower limit (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_7(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             1sigma upper limit (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_2(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             2sigma lower limit (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_3(nr__3) ! (mW/m2) Observed flux in observed 2.3-5 keV,
*                                             2sigma upper limit (FluxObsv_Up2_t)
      real*4        FluxcBFS   (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             best-fit model (FluxCorr_BF_s)
      real*4        FluxcMedS_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             posterior median (FluxCorr_Med_s)
      real*4        b_FluxcS_2 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             1sigma lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxcS_3 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             1sigma upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxcS  (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             2sigma lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxcS_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 0.5-2 keV,
*                                             2sigma upper limit (FluxCorr_Up2_s)
      real*4        FluxcBFT   (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV, best-fit model
*                                             (FluxCorr_BF_t)
      real*4        FluxcMedT_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             posterior median (FluxCorr_Med_t)
      real*4        b_FluxcT_2 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             1sigma lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxcT_3 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             1sigma upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxcT  (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             2sigma lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxcT_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             observed 2.3-5 keV,
*                                             2sigma upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS   (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             posterior median (FluxIntr_Med_s)
      real*4        b_FluxIS_2 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2keV,
*                                             1sigma lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_3 (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2keV,
*                                             1sigma upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS  (nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_1(nr__3) ! (mW/m2) Absorption corrected flux in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, best- fit
*                                             (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_1(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, median
*                                             (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_2(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 1sigma
*                                             lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_3(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 1sigma
*                                             upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 2sigma
*                                             lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_1(nr__3) ! (mW/m2/eV) Absorption corrected flux at
*                                             rest-frame 2keV, 2sigma
*                                             upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS    (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             best-fit model (LumiIntr_BF_s)
      real*8        LumIMedS_1 (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             posterior median (LumiIntr_Med_s)
      real*8        b_LumIS_2  (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             1sigma lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_3  (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             1sigma upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS   (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_1 (nr__3) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                             rest-frame 0.5-2 keV,
*                                             2sigma upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV (nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, best-fit model
*                                             (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_1(nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, posterior median
*                                             (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_2(nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 1sigma lower limit
*                                             (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_3(nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 1sigma upper limit
*                                             (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV(nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 2sigma lower limit
*                                             (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_1(nr__3) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                             rest-frame 2keV, 2sigma upper limit
*                                             (LumiIntr_Up2_2keV)
      real*8        GammaMed   (nr__3) ! power-law slope, posterior median
*                                             (Gamma_Med)
      real*8        b_Gamma    (nr__3) ! power-law slope, percentile
*                                             lower limit (Gamma_Lo)
      real*8        B_Gamma_1  (nr__3) ! power-law slope, percentile
*                                             upper limit (Gamma_Up)
      real*8        GammaMean  (nr__3) ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma    (nr__3) ! power-law slope, standard deviation
*                                             (Gamma_Std)
      real*8        lognHMed_1 (nr__3) ! ([cm-2]) log AGN column density,
*                                             posterior median (lognH_Med)
      real*8        b_lognH_2  (nr__3) ! ([cm-2]) log AGN column density,
*                                             percentile lower limit (lognH_Lo)
      real*8        B_lognH_3  (nr__3) ! ([cm-2]) log AGN column density,
*                                             percentile upper limit (lognH_Up)
      real*8        lognHMean_1(nr__3) ! ([cm-2]) log AGN column density,
*                                             mean (lognH_Mean)
      real*8        s_lognH_1  (nr__3) ! ([cm-2]) log AGN column density,
*                                             standard deviation (lognH_Std)
      real*8        logPowNMed (nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             posterior median (logPowNorm_Med)
      real*8        b_logPowN  (nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             percentile lower limit
*                                             (logPowNorm_Lo)
      real*8        B_logPowN_1(nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             percentile upper limit
*                                             (logPowNorm_Up)
      real*8        logPowNMean(nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             mean (logPowNorm_Mean)
      real*8        s_logPowN  (nr__3) ! (ph/keV/cm2/s) log power-law normalization,
*                                             standard deviation (logPowNorm_Std)
      real*4        logBkgNMed_1(nr__3) ! log background normalization,
*                                             posterior median (logBkgNorm_Med)
      real*8        b_logBkgN_2(nr__3) ! log background normalization,
*                                             percentile lower limit
*                                             (logBkgNorm_Lo)
      real*4        B_logBkgN_3(nr__3) ! log background normalization,
*                                             percentile upper limit
*                                             (logBkgNorm_Up)
      real*4        logBkgNMean_1(nr__3) ! log background normalization,
*                                             mean (logBkgNorm_Mean)
      real*8        s_logBkgN_1(nr__3) ! log background normalization,
*                                             standard deviation (logBkgNorm_Std)
      real*8        chi25_1    (nr__3) ! ?=-99 Chi^2 of the best-fit model
*                                             against the rebinned data
*                                             (>=25 cts/bin) (chi25)
      integer*4     dof25_1    (nr__3) ! ?=-99 DOF of the rebinned data,
*                                             which has at least 25 counts per
*                                             bin (dof25)
      real*8        logZ_1     (nr__3) ! log10 Bayesian evidence (logZ)
      integer*4     NHclass_1  (nr__3) ! [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   (nr__4) ! [1/32683] ID of the sources in
*                                                 the eFEDS main X-ray catalog
*                                                 (Paper I) (ID_SRC)
      real*4        z_2        (nr__4) ! [-0.01/8.0] The redshift adopted
*                                                 for the spectral model
*                                                 (Redshift)
      real*8        dGmKL      (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 KL divergence (in nats units)
*                                                 (dGm_KL)
      real*8        b_dGmH     (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 HDI lower limit (dGm_HLo)
      real*8        B_dGmH_1   (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 HDI upper limit (dGm_HUp)
      real*8        GammaKL_1  (nr__4) ! ? power-law slope,
*                                                 KL divergence (in nats units)
*                                                 (Gamma_KL)
      real*8        b_GammaH_2 (nr__4) ! power-law slope, HDI lower limit
*                                                 (Gamma_HLo)
      real*8        B_GammaH_3 (nr__4) ! power-law slope, HDI upper limit
*                                                 (Gamma_HUp)
      real*8        lognHKL_2  (nr__4) ! log AGN column density,
*                                                 KL divergence (in nats units)
*                                                 (lognH_KL)
      real*8        b_lognHH_4 (nr__4) ! ([cm-2]) log AGN column density,
*                                                 HDI lower limit (lognH_HLo)
      real*8        B_lognHH_5 (nr__4) ! ([cm-2]) log AGN column density,
*                                                 HDI upper limit (lognH_HUp)
      real*8        logPowNKL_1(nr__4) ! log power-law normalization,
*                                                 KL divergence (in nats units)
*                                                 (logPowNorm_KL)
      real*8        b_logPowNH_2(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI lower limit
*                                                 (logPowNorm_HLo)
      real*8        B_logPowNH_3(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI upper limit
*                                                 (logPowNorm_HUp)
      real*8        logFracKL  (nr__4) ! ? Constant factor multiplied to
*                                                 the soft power-law,
*                                                 KL divergence (in nats units)
*                                                 (logFrac_KL)
      real*8        b_logFracH (nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 HDI lower limit (logFrac_HLo)
      real*4        B_logFracH_1(nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 HDI upper limit (logFrac_HUp)
      real*8        dGmBF      (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 best-fit model (dGm_BF)
      real*8        GammaBF_1  (nr__4) ! power-law slope, best-fit model
*                                                 (Gamma_BF)
      real*8        lognHBF_2  (nr__4) ! ([cm-2]) log AGN column density,
*                                                 best-fit model (lognH_BF)
      real*8        logPowNBF_1(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 best-fit model (logPowNorm_BF)
      real*4        logFracBF  (nr__4) ! Constant factor multiplied to
*                                                 the soft power-law, best-fit
*                                                 model (logFrac_BF)
      real*4        FluxObsBFS_2(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, best-fit model
*                                                 (FluxObsv_BF_s)
      real*4        FluxObsMedS_4(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, posterior median
*                                                 (FluxObsv_Med_s)
      real*4        b_FluxObsS_8(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_9(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_4(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_5(nr__4) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_s)
      real*4        FluxObsBFT_2(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, best-fit model
*                                                 (FluxObsv_BF_t)
      real*4        FluxObsMedT_4(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, posterior median
*                                                 (FluxObsv_Med_t)
      real*4        b_FluxObsT_8(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_9(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_4(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_5(nr__4) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_t)
      real*4        FluxCBFS_1 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV,
*                                                 best-fit model (FluxCorr_BF_s)
      real*4        FluxCMedS_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, posterior
*                                                 median (FluxCorr_Med_s)
      real*4        b_FluxCS_4 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_5 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_3(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_1 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV,
*                                                 best-fit model (FluxCorr_BF_t)
      real*4        FluxCMedT_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, posterior
*                                                 median (FluxCorr_Med_t)
      real*4        b_FluxCT_4 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_5 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCT_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCT_3(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS_1 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 posterior median
*                                                 (FluxIntr_Med_s)
      real*4        b_FluxIS_4 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_5 (nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_2(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_3(nr__4) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_1(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, best-fit
*                                                 (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_2(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, median
*                                                 (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_4(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_5(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_2(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_3(nr__4) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_1  (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (LumiIntr_BF_s)
      real*8        LumIMedS_2 (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (LumiIntr_Med_s)
      real*8        b_LumIS_4  (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 1sigma lower limit
*                                                 (LumiIntr_Lo1_s)
      real*8        B_LumIS_5  (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 1sigma upper limit
*                                                 (LumiIntr_Up1_s)
      real*8        b2_LumIS_2 (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 2sigma lower limit
*                                                 (LumiIntr_Lo2_s)
      real*8        B2_LumIS_3 (nr__4) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV,
*                                                 2sigma upper limit
*                                                 (LumiIntr_Up2_s)
      real*8        LumIBF2keV_1(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, best-fit
*                                                 model (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_2(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, posterior
*                                                 median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_4(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_5(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_2(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_3(nr__4) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_2keV)
      real*8        dGmMed     (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 posterior median (dGm_Med)
      real*8        b_dGm      (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 1sigma lower limit (dGm_Lo)
      real*8        B_dGm_1    (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 1sigma upper limit (dGm_Up)
      real*8        dGmMean    (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 mean (dGm_Mean)
      real*8        s_dGm      (nr__4) ! Soft-powerlaw slope minus
*                                                 primary-powerlaw slope,
*                                                 standard deviation (dGm_Std)
      real*8        GammaMed_1 (nr__4) ! power-law slope, posterior
*                                                 median (Gamma_Med)
      real*8        b_Gamma_2  (nr__4) ! power-law slope, percentile
*                                                 lower limit (Gamma_Lo)
      real*8        B_Gamma_3  (nr__4) ! power-law slope, percentile
*                                                 upper limit (Gamma_Up)
      real*8        GammaMean_1(nr__4) ! power-law slope, mean
*                                                 (Gamma_Mean)
      real*4        s_Gamma_1  (nr__4) ! power-law slope,
*                                                 standard deviation (Gamma_Std)
      real*8        lognHMed_2 (nr__4) ! ([cm-2]) log AGN column density,
*                                                 posterior median (lognH_Med)
      real*8        b_lognH_4  (nr__4) ! ([cm-2]) log AGN column density,
*                                                 percentile lower limit
*                                                 (lognH_Lo)
      real*8        B_lognH_5  (nr__4) ! ([cm-2]) log AGN column density,
*                                                 percentile upper limit
*                                                 (lognH_Up)
      real*8        lognHMean_2(nr__4) ! ([cm-2]) log AGN column density,
*                                                 mean (lognH_Mean)
      real*8        s_lognH_2  (nr__4) ! ([cm-2]) log AGN column density,
*                                                 standard deviation (lognH_Std)
      real*8        logPowNMed_1(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 posterior median
*                                                 (logPowNorm_Med)
      real*8        b_logPowN_2(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile lower limit
*                                                 (logPowNorm_Lo)
      real*8        B_logPowN_3(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile upper limit
*                                                 (logPowNorm_Up)
      real*8        logPowNMean_1(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 mean (logPowNorm_Mean)
      real*8        s_logPowN_1(nr__4) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 standard deviation
*                                                 (logPowNorm_Std)
      real*8        logFracMed (nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 posterior median (logFrac_Med)
      real*8        b_logFrac  (nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 1sigma lower limit (logFrac_Lo)
      real*8        B_logFrac_1(nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 1sigma upper limit (logFrac_Up)
      real*8        logFracMean(nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 mean (logFrac_Mean)
      real*8        s_logFrac  (nr__4) ! Constant factor multiplied to
*                                                 the soft power-law,
*                                                 standard deviation
*                                                 (logFrac_Std)
      real*4        logBkgNMed_2(nr__4) ! log background normalization,
*                                                 posterior median
*                                                 (logBkgNorm_Med)
      real*8        b_logBkgN_4(nr__4) ! log background normalization,
*                                                 percentile lower limit
*                                                 (logBkgNorm_Lo)
      real*4        B_logBkgN_5(nr__4) ! log background normalization,
*                                                 percentile upper limit
*                                                 (logBkgNorm_Up)
      real*4        logBkgNMean_2(nr__4) ! log background normalization,
*                                                 mean (logBkgNorm_Mean)
      real*8        s_logBkgN_2(nr__4) ! log background normalization,
*                                                 standard deviation
*                                                 (logBkgNorm_Std)
      real*8        chi25_2    (nr__4) ! ?=-99 Chi^2 of the best-fit
*                                                 model against the rebinned
*                                                 data (>=25 cts/bin) (chi25)
      integer*4     dof25_2    (nr__4) ! ?=-99 DOF of the rebinned data,
*                                                 which has at least 25 counts
*                                                 per bin (dof25)
      real*8        logZ_2     (nr__4) ! 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   (nr__5) ! [1/32683] ID of the sources in
*                                                 the eFEDS main X-ray catalog
*                                                 (Paper I) (ID_SRC)
      real*4        z_3        (nr__5) ! [-0.01/8.0] The redshift adopted
*                                                 for the spectral model
*                                                 (Redshift)
      real*8        GammaKL_2  (nr__5) ! power-law slope, KL divergence
*                                                 (in nats units) (Gamma_KL)
      real*8        b_GammaH_4 (nr__5) ! power-law slope, HDI lower limit
*                                                 (Gamma_HLo)
      real*8        B_GammaH_5 (nr__5) ! power-law slope, HDI upper limit
*                                                 (Gamma_HUp)
      real*8        lognHKL_3  (nr__5) ! log AGN column density,
*                                                 KL divergence (in nats units)
*                                                 (lognH_KL)
      real*8        b_lognHH_6 (nr__5) ! ([cm-2]) log AGN column density,
*                                                 HDI lower limit (lognH_HLo)
      real*8        B_lognHH_7 (nr__5) ! ([cm-2]) log AGN column density,
*                                                 HDI upper limit (lognH_HUp)
      real*8        logPowNKL_2(nr__5) ! log power-law normalization,
*                                                 KL divergence (in nats units)
*                                                 (logPowNorm_KL)
      real*8        b_logPowNH_4(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 HDI lower limit
*                                                 (logPowNorm_HLo)
      real*8        B_logPowNH_5(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI upper limit (logPowNorm_HUp)
      real*8        logkTKL_1  (nr__5) ! log temperature,
*                                                 KL divergence (in nats units)
*                                                 (logkT_KL)
      real*8        b_logkTH_2 (nr__5) ! ([keV]) log temperature, HDI lower limit
*                                                 (logkT_HLo)
      real*8        B_logkTH_3 (nr__5) ! ([keV]) log temperature, HDI upper limit
*                                                 (logkT_HUp)
      real*8        logBBNKL   (nr__5) ! log blackbody normalization,
*                                                 KL divergence (in nats units)
*                                                 (logBBNorm_KL)
      real*8        b_logBBNH  (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, HDI
*                                                 lower limit (logBBNorm_HLo)
      real*8        B_logBBNH_1(nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, HDI
*                                                 upper limit (logBBNorm_HUp)
      real*8        GammaBF_2  (nr__5) ! power-law slope, best-fit model
*                                                 (Gamma_BF)
      real*8        lognHBF_3  (nr__5) ! ([cm-2]) log AGN column density,
*                                                 best-fit model (lognH_BF)
      real*8        logPowNBF_2(nr__5) ! ([ph/keV/cm2/s]) log power- law normalization,
*                                                 best-fit model (logPowNorm_BF)
      real*8        logkTBF_1  (nr__5) ! ([keV]) log temperature, best- fit model
*                                                 (logkT_BF)
      real*8        logBBNBF   (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/pc^2^, best-fit
*                                                 model (logBBNorm_BF)
      real*4        FluxObsBFS_3(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, best-fit model
*                                                 (FluxObsv_BF_s)
      real*4        FluxObsMedS_5(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, posterior median
*                                                 (FluxObsv_Med_s)
      real*4        b_FluxObsS_10(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_11(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_6(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_7(nr__5) ! (mW/m2) Observed flux in observed
*                                                 0.5-2 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_s)
      real*4        FluxObsBFT_3(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, best-fit model
*                                                 (FluxObsv_BF_t)
      real*4        FluxObsMedT_5(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, posterior median
*                                                 (FluxObsv_Med_t)
      real*4        b_FluxObsT_10(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma lower limit
*                                                 (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_11(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 1sigma upper limit
*                                                 (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_6(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma lower limit
*                                                 (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_7(nr__5) ! (mW/m2) Observed flux in observed
*                                                 2.3-5 keV, 2sigma upper limit
*                                                 (FluxObsv_Up2_t)
      real*4        FluxCBFS_2 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, best-fit
*                                                 model (FluxCorr_BF_s)
      real*4        FluxCMedS_3(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, posterior
*                                                 median (FluxCorr_Med_s)
      real*4        b_FluxCS_6 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_7 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_4(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_5(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 0.5-2 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_2 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, best-fit
*                                                 model (FluxCorr_BF_t)
      real*4        FluxCMedT_3(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, posterior
*                                                 median (FluxCorr_Med_t)
      real*4        b_FluxCT_6 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_7 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 1sigma
*                                                 upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCT_4(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCT_5(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 observed 2.3-5 keV, 2sigma
*                                                 upper limit (FluxCorr_Up2_t)
      real*4        FluxIBFS_2 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV,
*                                                 best-fit model (FluxIntr_BF_s)
      real*4        FluxIMedS_3(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (FluxIntr_Med_s)
      real*4        b_FluxIS_6 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_7 (nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_4(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_5(nr__5) ! (mW/m2) Absorption corrected flux in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_2(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, best-fit
*                                                 (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_3(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, median
*                                                 (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_6(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_7(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_4(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_5(nr__5) ! (mW/m2/eV) Absorption corrected flux at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_2  (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, best-fit
*                                                 model (LumiIntr_BF_s)
      real*8        LumIMedS_3 (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, posterior
*                                                 median (LumiIntr_Med_s)
      real*8        b_LumIS_6  (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_7  (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS_4 (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_5 (nr__5) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                 rest-frame 0.5-2 keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV_2(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, best-fit
*                                                 model (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_3(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, posterior
*                                                 median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_6(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_7(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 1sigma
*                                                 upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_4(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_5(nr__5) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                 rest-frame 2keV, 2sigma
*                                                 upper limit (LumiIntr_Up2_2keV)
      real*8        GammaMed_2 (nr__5) ! power-law slope, posterior
*                                                 median (Gamma_Med)
      real*8        b_Gamma_4  (nr__5) ! power-law slope, percentile
*                                                 lower limit (Gamma_Lo)
      real*8        B_Gamma_5  (nr__5) ! power-law slope, percentile
*                                                 upper limit (Gamma_Up)
      real*8        GammaMean_2(nr__5) ! power-law slope, mean
*                                                 (Gamma_Mean)
      real*8        s_Gamma_2  (nr__5) ! power-law slope,
*                                                 standard deviation (Gamma_Std)
      real*8        lognHMed_3 (nr__5) ! ([cm-2]) log AGN column density,
*                                                 posterior median (lognH_Med)
      real*8        b_lognH_6  (nr__5) ! ([cm-2]) log AGN column density,
*                                                 percentile lower limit
*                                                 (lognH_Lo)
      real*8        B_lognH_7  (nr__5) ! ([cm-2]) log AGN column density,
*                                                 percentile upper limit
*                                                 (lognH_Up)
      real*8        lognHMean_3(nr__5) ! ([cm-2]) log AGN column density,
*                                                 mean (lognH_Mean)
      real*8        s_lognH_3  (nr__5) ! ([cm-2]) log AGN column density,
*                                                 standard deviation (lognH_Std)
      real*8        logPowNMed_2(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 posterior median
*                                                 (logPowNorm_Med)
      real*8        b_logPowN_4(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile lower limit
*                                                 (logPowNorm_Lo)
      real*8        B_logPowN_5(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 percentile upper limit
*                                                 (logPowNorm_Up)
      real*8        logPowNMean_2(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 mean (logPowNorm_Mean)
      real*8        s_logPowN_2(nr__5) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                 standard deviation
*                                                 (logPowNorm_Std)
      real*8        logkTMed_1 (nr__5) ! ([keV]) log temperature, posterior
*                                                 median (logkT_Med)
      real*8        b_logkT_2  (nr__5) ! ([keV]) log temperature, percentile
*                                                 lower limit (logkT_Lo)
      real*8        B_logkT_3  (nr__5) ! ([keV]) log temperature, percentile
*                                                 upper limit (logkT_Up)
      real*8        logkTMean_1(nr__5) ! ([keV]) log temperature, mean
*                                                 (logkT_Mean)
      real*8        s_logkT_1  (nr__5) ! ([keV]) log temperature, standard
*                                                 deviation (logkT_Std)
      real*8        logBBNMed  (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, posterior
*                                                 median (logBBNorm_Med)
      real*8        b_logBBN   (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, percentile
*                                                 lower limit (logBBNorm_Lo)
      real*8        B_logBBN_1 (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, percentile
*                                                 upper limit (logBBNorm_Up)
      real*8        logBBNMean (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^,
*                                                 mean (logBBNorm_Mean)
      real*8        s_logBBN   (nr__5) ! ([10+24W/pc2]) log blackbody normalization in
*                                                 10^+31^erg/s/pc^2^, standard
*                                                 deviation (logBBNorm_Std)
      real*4        logBkgNMed_3(nr__5) ! log background normalization,
*                                                 posterior median
*                                                 (logBkgNorm_Med)
      real*8        b_logBkgN_6(nr__5) ! log background normalization,
*                                                 percentile lower limit
*                                                 (logBkgNorm_Lo)
      real*4        B_logBkgN_7(nr__5) ! log background normalization,
*                                                 percentile upper limit
*                                                 (logBkgNorm_Up)
      real*4        logBkgNMean_3(nr__5) ! log background normalization,
*                                                 mean (logBkgNorm_Mean)
      real*8        s_logBkgN_3(nr__5) ! log background normalization,
*                                                 standard deviation
*                                                 (logBkgNorm_Std)
      real*8        chi25_3    (nr__5) ! ?=-99 Chi^2 of the best-fit
*                                                 model against the rebinned
*                                                 data (>=25 cts/bin) (chi25)
      integer*4     dof25_3    (nr__5) ! ?=-99 DOF of the rebinned data,
*                                                 which has at least
*                                                 25 counts per bin (dof25)
      real*8        logZ_3     (nr__5) ! log10 Bayesian evidence (logZ)
      real*8        logRbb     (nr__5) ! 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   (nr__6) ! [1/32683] ID of the sources in the
*                                               eFEDS main X-ray catalog
*                                               (Paper I) (ID_SRC)
      real*4        z_4        (nr__6) ! [-0.01/8.0] The redshift adopted
*                                               for the spectral model (Redshift)
      real*8        lognHKL_4  (nr__6) ! log AGN column density,
*                                               KL divergence (in nats units)
*                                               (lognH_KL)
      real*8        b_lognHH_8 (nr__6) ! ([cm-2]) log AGN column density,
*                                               HDI lower limit (lognH_HLo)
      real*8        B_lognHH_9 (nr__6) ! ([cm-2]) log AGN column density,
*                                               HDI upper limit (lognH_HUp)
      real*8        logPowNKL_3(nr__6) ! log power-law normalization,
*                                               KL divergence (in nats units)
*                                               (logPowNorm_KL)
      real*8        b_logPowNH_6(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_7(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               HDI upper limit (logPowNorm_HUp)
      real*8        lognHBF_4  (nr__6) ! ([cm-2]) log AGN column density,
*                                               best-fit model (lognH_BF)
      real*8        logPowNBF_3(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_4(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, best-fit model
*                                               (FluxObsv_BF_s)
      real*4        FluxObsMedS_6(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, posterior median
*                                               (FluxObsv_Med_s)
      real*4        b_FluxObsS_12(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 1sigma lower limit
*                                               (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_13(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 1sigma upper limit
*                                               (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_8(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 2sigma lower limit
*                                               (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_9(nr__6) ! (mW/m2) Observed flux in observed
*                                               0.5-2 keV, 2sigma upper limit
*                                               (FluxObsv_Up2_s)
      real*4        FluxObsBFT_4(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, best-fit model
*                                               (FluxObsv_BF_t)
      real*4        FluxObsMedT_6(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, posterior median
*                                               (FluxObsv_Med_t)
      real*4        b_FluxObsT_12(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 1sigma lower limit
*                                               (FluxObsv_Lo1)
      real*4        B_FluxObsT_13(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 1sigma upper limit
*                                               (FluxObsv_Up1)
      real*4        b2_FluxObsT_8(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 2sigma lower limit
*                                               (FluxObsv_Lo2)
      real*4        B2_FluxObsT_9(nr__6) ! (mW/m2) Observed flux in observed
*                                               2.3-5 keV, 2sigma upper limit
*                                               (FluxObsv_Up2)
      real*4        FluxCBFS_3 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, best-fit
*                                               model (FluxCorr_BF_s)
      real*4        FluxCMedS_4(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, posterior
*                                               median (FluxCorr_Med_s)
      real*4        b_FluxCS_8 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 1sigma
*                                               lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_9 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 1sigma
*                                               upper limit (FluxCorr_Up1_s)
      real*4        b2_FluxCS_6(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 2sigma
*                                               lower limit (FluxCorr_Lo2_s)
      real*4        B2_FluxCS_7(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 0.5-2 keV, 2sigma
*                                               upper limit (FluxCorr_Up2_s)
      real*4        FluxCBFT_3 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, best-fit
*                                               model (FluxCorr_BF_t)
      real*4        FluxCMedT_4(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, posterior
*                                               median (FluxCorr_Med_t)
      real*4        b_FluxCT_8 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 1sigma
*                                               lower limit (FluxCorr_Lo1)
      real*4        B_FluxCT_9 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 1sigma
*                                               upper limit (FluxCorr_Up1)
      real*4        b2_FluxCT_6(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 2sigma
*                                               lower limit (FluxCorr_Lo2)
      real*4        B2_FluxCT_7(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               observed 2.3-5 keV, 2sigma
*                                               upper limit (FluxCorr_Up2)
      real*4        FluxIBFS_3 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, best-fit
*                                               model (FluxIntr_BF_s)
      real*4        FluxIMedS_4(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, posterior
*                                               median (FluxIntr_Med_s)
      real*4        b_FluxIS_8 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2keV, 1sigma
*                                               lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_9 (nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2keV, 1sigma
*                                               upper limit (FluxIntr_Up1_s)
      real*4        b2_FluxIS_6(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               lower limit (FluxIntr_Lo2_s)
      real*4        B2_FluxIS_7(nr__6) ! (mW/m2) Absorption corrected flux in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               upper limit (FluxIntr_Up2_s)
      real*4        FluxIBF2keV_3(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, best- fit
*                                               (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_4(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, median
*                                               (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_8(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 1sigma
*                                               lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_9(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 1sigma
*                                               upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_6(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 2sigma
*                                               lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_7(nr__6) ! (mW/m2/eV) Absorption corrected flux at
*                                               rest-frame 2keV, 2sigma
*                                               upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_3  (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, best-fit
*                                               model (LumiIntr_BF_s)
      real*8        LumIMedS_4 (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, posterior
*                                               median (LumiIntr_Med_s)
      real*8        b_LumIS_8  (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 1sigma
*                                               lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_9  (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 1sigma
*                                               upper limit (LumiIntr_Up1_s)
      real*8        b2_LumIS_6 (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               lower limit (LumiIntr_Lo2_s)
      real*8        B2_LumIS_7 (nr__6) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                               rest-frame 0.5-2 keV, 2sigma
*                                               upper limit (LumiIntr_Up2_s)
      real*8        LumIBF2keV_3(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, best-fit model
*                                               (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_4(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, posterior
*                                               median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_8(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 1sigma
*                                               lower limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_9(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 1sigma
*                                               upper limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_6(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 2sigma
*                                               lower limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_7(nr__6) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                               rest-frame 2keV, 2sigma upper
*                                               limit (LumiIntr_Up2_2keV)
      real*8        lognHMed_4 (nr__6) ! ([cm-2]) log AGN column density,
*                                               posterior median (lognH_Med)
      real*8        b_lognH_8  (nr__6) ! ([cm-2]) log AGN column density,
*                                               percentile lower limit (lognH_Lo)
      real*8        B_lognH_9  (nr__6) ! ([cm-2]) log AGN column density,
*                                               percentile upper limit (lognH_Up)
      real*8        lognHMean_4(nr__6) ! ([cm-2]) log AGN column density, mean
*                                               (lognH_Mean)
      real*8        s_lognH_4  (nr__6) ! ([cm-2]) log AGN column density,
*                                               standard deviation (lognH_Std)
      real*8        logPowNMed_3(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               posterior median (logPowNorm_Med)
      real*8        b_logPowN_6(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               percentile lower limit
*                                               (logPowNorm_Lo)
      real*8        B_logPowN_7(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               percentile upper limit
*                                               (logPowNorm_Up)
      real*8        logPowNMean_3(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               mean (logPowNorm_Mean)
      real*8        s_logPowN_3(nr__6) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                               standard deviation
*                                               (logPowNorm_Std)
      real*4        logBkgNMed_4(nr__6) ! log background normalization,
*                                               posterior median (logBkgNorm_Med)
      real*4        b_logBkgN_8(nr__6) ! log background normalization,
*                                               percentile lower limit
*                                               (logBkgNorm_Lo)
      real*4        B_logBkgN_9(nr__6) ! log background normalization,
*                                               percentile upper limit
*                                               (logBkgNorm_Up)
      real*4        logBkgNMean_4(nr__6) ! log background normalization,
*                                               mean (logBkgNorm_Mean)
      real*8        s_logBkgN_4(nr__6) ! log background normalization,
*                                               standard deviation
*                                               (logBkgNorm_Std)
      real*8        logZ_4     (nr__6) ! 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   (nr__7) ! [1/32683] ID of the sources in
*                                                the eFEDS main X-ray catalog
*                                                (Paper I) (ID_SRC)
      real*4        z_5        (nr__7) ! [-0.01/8.0] The redshift adopted
*                                                for the spectral model
*                                                (Redshift)
      real*8        logPowNKL_4(nr__7) ! log power-law normalization,
*                                                KL divergence (in nats units)
*                                                (logPowNorm_KL)
      real*8        b_logPowNH_8(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_9(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                HDI upper limit (logPowNorm_HUp)
      real*8        logPowNBF_4(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_5(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, best-fit model
*                                                (FluxObsv_BF_s)
      real*4        FluxObsMedS_7(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, posterior median
*                                                (FluxObsv_Med_s)
      real*4        b_FluxObsS_14(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 1sigma lower limit
*                                                (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_15(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 1sigma upper limit
*                                                (FluxObsvs)
      real*4        b2_FluxObsS_10(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 2sigma lower limit
*                                                (FluxObsvs)
      real*4        B2_FluxObsS_11(nr__7) ! (mW/m2) Observed flux in observed
*                                                0.5-2 keV, 2sigma upper limit
*                                                (FluxObsvs)
      real*4        FluxObsBFT_5(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, best-fit model
*                                                (FluxObsv_BF_t)
      real*4        FluxObsMedT_7(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, posterior median
*                                                (FluxObsv_Med_t)
      real*4        b_FluxObsT_14(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 1sigma lower limit
*                                                (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_15(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 1sigma upper limit
*                                                (FluxObsvt)
      real*4        b2_FluxObsT_10(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 2sigma lower limit
*                                                (FluxObsvt)
      real*4        B2_FluxObsT_11(nr__7) ! (mW/m2) Observed flux in observed
*                                                2.3-5 keV, 2sigma upper limit
*                                                (FluxObsvt)
      real*4        FluxCBFS_4 (nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV,
*                                                best-fit model (FluxCorr_BF_s)
      real*4        FluxCMedS_5(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, posterior
*                                                median (FluxCorr_Med_s)
      real*4        b_FluxCS_10(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 1sigma
*                                                lower limit (FluxCorr_Lo1_s)
      real*4        B_FluxCS_11(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 1sigma
*                                                upper limit (FluxCorrs)
      real*4        b2_FluxCS_8(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 2sigma
*                                                lower limit (FluxCorrs)
      real*4        B2_FluxCS_9(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 0.5-2 keV, 2sigma
*                                                upper limit (FluxCorrs)
      real*4        FluxCBFT_4 (nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, best-fit
*                                                model (FluxCorr_BF_t)
      real*4        FluxCMedT_5(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, posterior
*                                                median (FluxCorr_Med_t)
      real*4        b_FluxCT_10(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 1sigma
*                                                lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCT_11(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 1sigma
*                                                upper limit (FluxCorrt)
      real*4        b2_FluxCT_8(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 2sigma
*                                                lower limit (FluxCorrt)
      real*4        B2_FluxCT_9(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                observed 2.3-5 keV, 2sigma
*                                                upper limit (FluxCorrt)
      real*4        FluxIBFS_4 (nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, best-fit
*                                                model (FluxIntr_BF_s)
      real*4        FluxIMedS_5(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, posterior
*                                                median (FluxIntr_Med_s)
      real*4        b_FluxIS_10(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2keV, 1sigma
*                                                lower limit (FluxIntr_Lo1_s)
      real*4        B_FluxIS_11(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2keV, 1sigma
*                                                upper limit (FluxIntrs)
      real*4        b2_FluxIS_8(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                lower limit (FluxIntrs)
      real*4        B2_FluxIS_9(nr__7) ! (mW/m2) Absorption corrected flux in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                upper limit (FluxIntrs)
      real*4        FluxIBF2keV_4(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, best- fit
*                                                (FluxIntr_BF_2keV)
      real*4        FluxIMed2keV_5(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, median
*                                                (FluxIntr_Med_2keV)
      real*4        b_FluxI2keV_10(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 1sigma
*                                                lower limit (FluxIntr_Lo1_2keV)
      real*4        B_FluxI2keV_11(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 1sigma
*                                                upper limit (FluxIntr_Up1_2keV)
      real*4        b2_FluxI2keV_8(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 2sigma
*                                                lower limit (FluxIntr_Lo2_2keV)
      real*4        B2_FluxI2keV_9(nr__7) ! (mW/m2/eV) Absorption corrected flux at
*                                                rest-frame 2keV, 2sigma
*                                                upper limit (FluxIntr_Up2_2keV)
      real*8        LumIBFS_4  (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, best-fit
*                                                model (LumiIntr_BF_s)
      real*8        LumIMedS_5 (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, posterior
*                                                median (LumiIntr_Med_s)
      real*8        b_LumIS_10 (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 1sigma
*                                                lower limit (LumiIntr_Lo1_s)
      real*8        B_LumIS_11 (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 1sigma
*                                                upper limit (LumiIntrs)
      real*8        b2_LumI2S  (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                lower limit (LumiIntrs)
      real*8        B2_LumIS_8 (nr__7) ! (10-7W) ?=-99 Intrinsic luminosity in
*                                                rest-frame 0.5-2 keV, 2sigma
*                                                upper limit (LumiIntrs)
      real*8        LumIBF2keV_4(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, best-fit model
*                                                (LumiIntr_BF_2keV)
      real*8        LumIMed2keV_5(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, posterior
*                                                median (LumiIntr_Med_2keV)
      real*8        b_LumI2keV_10(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 1sigma lower
*                                                limit (LumiIntr_Lo1_2keV)
      real*8        B_LumI2keV_11(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 1sigma upper
*                                                limit (LumiIntr_Up1_2keV)
      real*8        b2_LumI2keV_8(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 2sigma lower
*                                                limit (LumiIntr_Lo2_2keV)
      real*8        B2_LumI2keV_9(nr__7) ! (10-7W/eV) ?=-99 Intrinsic luminosity at
*                                                rest-frame 2keV, 2sigma upper
*                                                limit (LumiIntr_Up2_2keV)
      real*8        logPowNMed_4(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                posterior median
*                                                (logPowNorm_Med)
      real*8        b_logPowN_8(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                percentile lower limit
*                                                (logPowNorm_Lo)
      real*8        B_logPowN_9(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                percentile upper limit
*                                                (logPowNorm_Up)
      real*8        logPowNMean_4(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                mean (logPowNorm_Mean)
      real*8        s_logPowN_4(nr__7) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                                standard deviation
*                                                (logPowNorm_Std)
      real*4        logBkgNMed_5(nr__7) ! log background normalization,
*                                                posterior median
*                                                (logBkgNorm_Med)
      real*4        b_logBkgN_10(nr__7) ! log background normalization,
*                                                percentile lower limit
*                                                (logBkgNorm_Lo)
      real*4        B_logBkgN_11(nr__7) ! log background normalization,
*                                                percentile upper limit
*                                                (logBkgNorm_Up)
      real*4        logBkgNMean_5(nr__7) ! log background normalization,
*                                                mean (logBkgNorm_Mean)
      real*8        s_logBkgN_5(nr__7) ! log background normalization,
*                                                standard deviation
*                                                (logBkgNorm_Std)
      real*8        logZ_5     (nr__7) ! 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   (nr__8) ! [1/32683] ID of the sources in the
*                                              eFEDS main X-ray catalog
*                                              (Paper I) (ID_SRC)
      real*4        z_6        (nr__8) ! [-0.01/8.0] The redshift adopted
*                                              for the spectral model (Redshift)
      real*8        GammaKL_3  (nr__8) ! power-law slope, KL divergence
*                                              (in nats units) (Gamma_KL)
      real*8        b_GammaH_6 (nr__8) ! power-law slope, HDI lower limit
*                                              (Gamma_HLo)
      real*8        B_GammaH_7 (nr__8) ! power-law slope, HDI upper limit
*                                              (Gamma_HUp)
      real*8        lognHKL_5  (nr__8) ! log AGN column density,
*                                              KL divergence (in nats units)
*                                              (lognH_KL)
      real*8        b_lognHH_10(nr__8) ! ([cm-2]) log AGN column density,
*                                              HDI lower limit (lognH_HLo)
      real*8        B_lognHH_11(nr__8) ! ([cm-2]) log AGN column density,
*                                              HDI upper limit (lognH_HUp)
      real*8        logPowNKL_5(nr__8) ! log power-law normalization,
*                                              KL divergence (in nats units)
*                                              (logPowNorm_KL)
      real*8        b_logPowNH_10(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowNH_11(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF_3  (nr__8) ! power-law slope, best-fit model
*                                              (Gamma_BF)
      real*8        lognHBF_5  (nr__8) ! ([cm-2]) log AGN column density,
*                                              best-fit model (lognH_BF)
      real*4        logPowNBF_5(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              best-fit model (logPowNorm_BF)
      real*4        FluxObsBFS_6(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, best-fit model
*                                              (FluxObsv_BF_s)
      real*4        FluxObsMedS_8(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, posterior median
*                                              (FluxObsv_Med_s)
      real*4        b_FluxObsS_16(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 1sigma lower limit
*                                              (FluxObsv_Lo1_s)
      real*4        B_FluxObsS_17(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 1sigma upper limit
*                                              (FluxObsv_Up1_s)
      real*4        b2_FluxObsS_12(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 2sigma lower limit
*                                              (FluxObsv_Lo2_s)
      real*4        B2_FluxObsS_13(nr__8) ! (mW/m2) Observed flux in observed
*                                              0.5-2 keV, 2sigma upper limit
*                                              (FluxObsv_Up2_s)
      real*8        GammaMed_3 (nr__8) ! power-law slope, posterior median
*                                              (Gamma_Med)
      real*8        b_Gamma_6  (nr__8) ! power-law slope, percentile
*                                              lower limit (Gamma_Lo)
      real*8        B_Gamma_7  (nr__8) ! power-law slope, percentile
*                                              upper limit (Gamma_Up)
      real*8        GammaMean_3(nr__8) ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma_3  (nr__8) ! power-law slope,
*                                              standard deviation (Gamma_Std)
      real*8        lognHMed_5 (nr__8) ! ([cm-2]) log AGN column density,
*                                              posterior median (lognH_Med)
      real*8        b_lognH_10 (nr__8) ! ([cm-2]) log AGN column density, percentile
*                                              lower limit (lognH_Lo)
      real*8        B_lognH_11 (nr__8) ! ([cm-2]) log AGN column density, percentile
*                                              upper limit (lognH_Up)
      real*8        lognHMean_5(nr__8) ! ([cm-2]) log AGN column density,
*                                              mean (lognH_Mean)
      real*8        s_lognH_5  (nr__8) ! ([cm-2]) log AGN column density,
*                                              standard deviation (lognH_Std)
      real*8        logPowNMed_5(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              posterior median (logPowNorm_Med)
      real*8        b_logPowN_10(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile lower limit
*                                              (logPowNorm_Lo)
      real*8        B_logPowN_11(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile upper limit
*                                              (logPowNorm_Up)
      real*8        logPowNMean_5(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              mean (logPowNorm_Mean)
      real*8        s_logPowN_5(nr__8) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              standard deviation
*                                              (logPowNorm_Std)
      real*4        logBkgNMed_6(nr__8) ! log background normalization,
*                                              posterior median (logBkgNorm_Med)
      real*4        blogBkgN   (nr__8) ! log background normalization,
*                                              percentile lower limit
*                                              (logBkgNorm_Lo)
      real*4        BlogBkgN_1 (nr__8) ! log background normalization,
*                                              percentile upper limit
*                                              (logBkgNorm_Up)
      real*4        logBkgNMean_6(nr__8) ! log background normalization,
*                                              mean (logBkgNorm_Mean)
      real*8        logBkgN    (nr__8) ! log background normalization,
*                                              standard deviation
*                                              (logBkgNorm_Std)
      real*8        logZ_6     (nr__8) ! 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   (nr__9) ! [1/32683] ID of the sources in the
*                                              eFEDS main X-ray catalog
*                                              (Paper I) (ID_SRC)
      real*4        z_7        (nr__9) ! [-0.01/8.0] The redshift adopted
*                                              for the spectral model (Redshift)
      real*8        GammaKL_4  (nr__9) ! power-law slope, KL divergence
*                                              (in nats units) (Gamma_KL)
      real*8        e_Gamma    (nr__9) ! power-law slope, HDI lower limit
*                                              (Gamma_HLo)
      real*8        E_Gamma_1  (nr__9) ! power-law slope, HDI upper limit
*                                              (Gamma_HUp)
      real*8        lognHKL_6  (nr__9) ! log AGN column density,
*                                              KL divergence (in nats units)
*                                              (lognH_KL)
      real*8        b_lognH_12 (nr__9) ! ([cm-2]) log AGN column density,
*                                              HDI lower limit (lognH_HLo)
      real*8        B_lognH_13 (nr__9) ! ([cm-2]) log AGN column density,
*                                              HDI upper limit (lognH_HUp)
      real*8        logPowNKL_6(nr__9) ! log power-law normalization,
*                                              KL divergence (in nats units)
*                                              (logPowNorm_KL)
      real*8        b_logPowN_12(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI lower limit (logPowNorm_HLo)
      real*8        B_logPowN_13(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              HDI upper limit (logPowNorm_HUp)
      real*8        GammaBF_4  (nr__9) ! power-law slope, best-fit model
*                                              (Gamma_BF)
      real*8        lognHBF_6  (nr__9) ! ([cm-2]) log AGN column density,
*                                              best-fit model (lognH_BF)
      real*8        logPowNBF_6(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              best-fit model (logPowNorm_BF)
      real*4        FluxObsBFT_6(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              best-fit model (FluxObsv_BF_t)
      real*4        FluxObsMedT_8(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              posterior median (FluxObsv_Med_t)
      real*4        b_FluxObsT_16(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              1sigma lower limit
*                                              (FluxObsv_Lo1_t)
      real*4        B_FluxObsT_17(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              1sigma upper limit
*                                              (FluxObsv_Up1_t)
      real*4        b2_FluxObsT_12(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              2sigma lower limit
*                                              (FluxObsv_Lo2_t)
      real*4        B2_FluxObsT_13(nr__9) ! (mW/m2) Observed flux in observed 2.3-5keV,
*                                              2sigma upper limit
*                                              (FluxObsv_Up2_t)
      real*4        FluxCBFT_5 (nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, best-fit model
*                                              (FluxCorr_BF_t)
      real*4        FluxCMedT_6(nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, posterior
*                                              median (FluxCorr_Med_t)
      real*4        b_FluxCMedT(nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 1sigma
*                                              lower limit (FluxCorr_Lo1_t)
      real*4        B_FluxCMedT_1(nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 1sigma
*                                              upper limit (FluxCorr_Up1_t)
      real*4        b2_FluxCMedT(nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 2sigma
*                                              lower limit (FluxCorr_Lo2_t)
      real*4        B2_FluxCMedT_1(nr__9) ! (mW/m2) Absorption corrected flux in
*                                              observed 2.3-5 keV, 2sigma
*                                              upper limit (FluxCorr_Up2_t)
      real*8        GammaMed_4 (nr__9) ! power-law slope, posterior median
*                                              (Gamma_Med)
      real*8        b_GammaMed (nr__9) ! power-law slope, percentile
*                                              lower limit (Gamma_Lo)
      real*8        B_GammaMed_1(nr__9) ! power-law slope, percentile
*                                              upper limit (Gamma_Up)
      real*8        GammaMean_4(nr__9) ! power-law slope, mean (Gamma_Mean)
      real*8        s_Gamma_4  (nr__9) ! power-law slope,
*                                              standard deviation (Gamma_Std)
      real*8        lognHMed_6 (nr__9) ! ([cm-2]) log AGN column density,
*                                              posterior median (lognH_Med)
      real*8        b_lognHMed (nr__9) ! ([cm-2]) log AGN column density,
*                                              percentile lower limit (lognH_Lo)
      real*8        B_lognHMed_1(nr__9) ! ([cm-2]) log AGN column density,
*                                              percentile upper limit (lognH_Up)
      real*8        lognHMean_6(nr__9) ! ([cm-2]) log AGN column density,
*                                              mean (lognH_Mean)
      real*8        s_lognH_6  (nr__9) ! ([cm-2]) log AGN column density,
*                                              standard deviation (lognH_Std)
      real*8        logPowNMed_6(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              posterior median (logPowNorm_Med)
      real*8        b_logPowNMed(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile lower limit
*                                              (logPowNorm_Lo)
      real*8        B_logPowNMed_1(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              percentile upper limit
*                                              (logPowNorm_Up)
      real*8        logPowNMean_6(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              mean (logPowNorm_Mean)
      real*8        s_logPowN_6(nr__9) ! ([ph/keV/cm2/s]) log power-law normalization,
*                                              standard deviation
*                                              (logPowNorm_Std)
      real*4        logBkgNMed_7(nr__9) ! log background normalization,
*                                              posterior median (logBkgNorm_Med)
      real*4        b_logBkgNMed(nr__9) ! log background normalization,
*                                              percentile lower limit
*                                              (logBkgNorm_Lo)
      real*4        B_logBkgNMed_1(nr__9) ! log background normalization,
*                                              percentile upper limit
*                                              (logBkgNorm_Up)
      real*4        logBkgNMean_7(nr__9) ! log background normalization,
*                                              mean (logBkgNorm_Mean)
      real*8        s_logBkgN_6(nr__9) ! log background normalization,
*                                              standard deviation
*                                              (logBkgNorm_Std)
      real*8        logZ_7     (nr__9) ! 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(i__),RAcdeg(i__),DEcdeg(i__),DetLike(i__),
     +  inArea90(i__),CTPLS8(i__),RALdeg(i__),DELdeg(i__),
     +  CTPqual(i__),CTPClass(i__),CTPz(i__),CTPzgrade(i__),
     +  inKiDSFlag(i__),LxModel(i__),NHclass(i__),FSModel(i__),
     +  FHModel(i__),galNH(i__),galNHI(i__),SrcCts(i__),
     +  FluxcMedS(i__),b_FluxcS(i__),B_FluxcS_1(i__),FluxcMedT(i__),
     +  b_FluxcT(i__),B_FluxcT_1(i__),FluxIMedS(i__),b_FluxIS(i__),
     +  B_FluxIS_1(i__),FluxIMed2keV(i__),b_FluxI2keV(i__),
     +  B_FluxI2keV_1(i__),LumIMedS(i__),b_LumIS(i__),B_LumIS_1(i__),
     +  LumIMed2keV(i__),b_LumI2keV(i__),B_LumI2keV_1(i__),
     +  FluxObsMedS(i__),b_FluxObsS(i__),B_FluxObsS_1(i__),
     +  FluxObsMedT(i__),b_FluxObsT(i__),B_FluxObsT_1(i__),
     +  lognHKLm1(i__),b_lognHm1(i__),B_lognHm1_1(i__),
     +  lognHMedm1(i__),GammaKLm3(i__),b_GammaHm3(i__),
     +  B_GammaHm3_1(i__),GammaMedm3(i__),logZm0(i__),logZm1(i__),
     +  logZm2(i__),logZm3(i__),logZm4(i__),logZm5(i__),L2500(i__),
     +  L5100(i__),W1mag(i__),e_W1mag(i__),W2mag(i__),e_W2mag(i__)
        if(ar__(686:696) .EQ. '') W1mag(i__) = rNULL__
        if(ar__(698:710) .EQ. '') e_W1mag(i__) = rNULL__
        if(ar__(712:722) .EQ. '') W2mag(i__) = rNULL__
        if(ar__(724:736) .EQ. '') e_W2mag(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID_SRC(i__),RAcdeg(i__),DEcdeg(i__),DetLike(i__),
     +  inArea90(i__),CTPLS8(i__),RALdeg(i__),DELdeg(i__),
     +  CTPqual(i__),CTPClass(i__),CTPz(i__),CTPzgrade(i__),
     +  inKiDSFlag(i__),LxModel(i__),NHclass(i__),FSModel(i__),
     +  FHModel(i__),galNH(i__),galNHI(i__),SrcCts(i__),
     +  FluxcMedS(i__),b_FluxcS(i__),B_FluxcS_1(i__),FluxcMedT(i__),
     +  b_FluxcT(i__),B_FluxcT_1(i__),FluxIMedS(i__),b_FluxIS(i__),
     +  B_FluxIS_1(i__),FluxIMed2keV(i__),b_FluxI2keV(i__),
     +  B_FluxI2keV_1(i__),LumIMedS(i__),b_LumIS(i__),B_LumIS_1(i__),
     +  LumIMed2keV(i__),b_LumI2keV(i__),B_LumI2keV_1(i__),
     +  FluxObsMedS(i__),b_FluxObsS(i__),B_FluxObsS_1(i__),
     +  FluxObsMedT(i__),b_FluxObsT(i__),B_FluxObsT_1(i__),
     +  lognHKLm1(i__),b_lognHm1(i__),B_lognHm1_1(i__),
     +  lognHMedm1(i__),GammaKLm3(i__),b_GammaHm3(i__),
     +  B_GammaHm3_1(i__),GammaMedm3(i__),logZm0(i__),logZm1(i__),
     +  logZm2(i__),logZm3(i__),logZm4(i__),logZm5(i__),L2500(i__),
     +  L5100(i__),W1mag(i__),e_W1mag(i__),W2mag(i__),e_W2mag(i__)
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(i__),RAcdeg_1(i__),DEcdeg_1(i__),DetLike_1(i__),
     +  ExtLike(i__),inArea90_1(i__),galNH_1(i__),galNHI_1(i__),
     +  ExpTime(i__),SrcCts_1(i__),RAdeg(i__),DEdeg(i__),Rad(i__),
     +  Ann1(i__),Ann2(i__),RegArea(i__),Backscal(i__),
     +  BackscalBkg(i__),FSModel_1(i__),FHModel_1(i__),
     +  Rate0_2_2_3(i__),e_Rate0_2_2_3(i__),Rate0_2_0_5(i__),
     +  e_Rate0_2_0_5(i__),Rate0_5_1(i__),e_Rate0_5_1(i__),
     +  Rate1_2(i__),e_Rate1_2(i__),Rate2_4_5(i__),e_Rate2_4_5(i__),
     +  Rate2_3_5(i__),e_Rate2_3_5(i__),Rate5_8(i__),e_Rate5_8(i__),
     +  BkgCts0_2_0_6(i__),e_BkgCts0_2_0_6(i__),BkgCts0_6_2_3(i__),
     +  e_BkgCts0_6_2_3(i__),BkgCts2_3_5(i__),e_BkgCts2_3_5(i__),
     +  BkgCts5_8(i__),e_BkgCts5_8(i__),PSFCor0_2_0_5(i__),
     +  PSFCor0_5_1(i__),PSFCor1_2(i__),PSFCor2_3_5(i__),Nempty(i__),
     +  FluxObsMedS_1(i__),b_FluxObsS_2(i__),B_FluxObsS_3(i__),
     +  FluxObsMedT_1(i__),b_FluxObsT_2(i__),B_FluxObsT_3(i__)
c    ..............Just test output...........
        write(6,2)
     +  ID_SRC_1(i__),RAcdeg_1(i__),DEcdeg_1(i__),DetLike_1(i__),
     +  ExtLike(i__),inArea90_1(i__),galNH_1(i__),galNHI_1(i__),
     +  ExpTime(i__),SrcCts_1(i__),RAdeg(i__),DEdeg(i__),Rad(i__),
     +  Ann1(i__),Ann2(i__),RegArea(i__),Backscal(i__),
     +  BackscalBkg(i__),FSModel_1(i__),FHModel_1(i__),
     +  Rate0_2_2_3(i__),e_Rate0_2_2_3(i__),Rate0_2_0_5(i__),
     +  e_Rate0_2_0_5(i__),Rate0_5_1(i__),e_Rate0_5_1(i__),
     +  Rate1_2(i__),e_Rate1_2(i__),Rate2_4_5(i__),e_Rate2_4_5(i__),
     +  Rate2_3_5(i__),e_Rate2_3_5(i__),Rate5_8(i__),e_Rate5_8(i__),
     +  BkgCts0_2_0_6(i__),e_BkgCts0_2_0_6(i__),BkgCts0_6_2_3(i__),
     +  e_BkgCts0_6_2_3(i__),BkgCts2_3_5(i__),e_BkgCts2_3_5(i__),
     +  BkgCts5_8(i__),e_BkgCts5_8(i__),PSFCor0_2_0_5(i__),
     +  PSFCor0_5_1(i__),PSFCor1_2(i__),PSFCor2_3_5(i__),Nempty(i__),
     +  FluxObsMedS_1(i__),b_FluxObsS_2(i__),B_FluxObsS_3(i__),
     +  FluxObsMedT_1(i__),b_FluxObsT_2(i__),B_FluxObsT_3(i__)
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(i__),z(i__),lognHKL(i__),b_lognHH(i__),
     +  B_lognHH_1(i__),logkTKL(i__),b_logkTH(i__),B_logkTH_1(i__),
     +  logAbundKL(i__),b_logAbunH(i__),B_logAbunH_1(i__),
     +  logApecNormKL(i__),b_logApecNormH(i__),B_logApecNormH_1(i__),
     +  lognHBF(i__),logkTBF(i__),logAbundBF(i__),logApecNBF(i__),
     +  FluxObsBFS(i__),FluxObsMedS_2(i__),b_FluxObsS_4(i__),
     +  B_FluxObsS_5(i__),b2_FluxObsS(i__),B2_FluxObsS_1(i__),
     +  FluxObsBFT(i__),FluxObsMedT_2(i__),b_FluxObsT_4(i__),
     +  B_FluxObsT_5(i__),b2_FluxObsT(i__),B2_FluxObsT_1(i__),
     +  lognHMed(i__),b_lognH(i__),B_lognH_1(i__),lognHMean(i__),
     +  s_lognH(i__),logkTMed(i__),b_logkT(i__),B_logkT_1(i__),
     +  logkTMean(i__),s_logkT(i__),logAbundMed(i__),b_logAbund(i__),
     +  B_logAbund_1(i__),logAbundMean(i__),s_logAbund(i__),
     +  logApecNMed(i__),b_logApecN(i__),B_logApecN_1(i__),
     +  logApecNMean(i__),s_logApecN(i__),logBkgNMed(i__),
     +  b_logBkgN(i__),B_logBkgN_1(i__),logBkgNMean(i__),
     +  s_logBkgN(i__),chi25(i__),dof25(i__),logZ(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_SRC_2(i__),z(i__),lognHKL(i__),b_lognHH(i__),
     +  B_lognHH_1(i__),logkTKL(i__),b_logkTH(i__),B_logkTH_1(i__),
     +  logAbundKL(i__),b_logAbunH(i__),B_logAbunH_1(i__),
     +  logApecNormKL(i__),b_logApecNormH(i__),B_logApecNormH_1(i__),
     +  lognHBF(i__),logkTBF(i__),logAbundBF(i__),logApecNBF(i__),
     +  FluxObsBFS(i__),FluxObsMedS_2(i__),b_FluxObsS_4(i__),
     +  B_FluxObsS_5(i__),b2_FluxObsS(i__),B2_FluxObsS_1(i__),
     +  FluxObsBFT(i__),FluxObsMedT_2(i__),b_FluxObsT_4(i__),
     +  B_FluxObsT_5(i__),b2_FluxObsT(i__),B2_FluxObsT_1(i__),
     +  lognHMed(i__),b_lognH(i__),B_lognH_1(i__),lognHMean(i__),
     +  s_lognH(i__),logkTMed(i__),b_logkT(i__),B_logkT_1(i__),
     +  logkTMean(i__),s_logkT(i__),logAbundMed(i__),b_logAbund(i__),
     +  B_logAbund_1(i__),logAbundMean(i__),s_logAbund(i__),
     +  logApecNMed(i__),b_logApecN(i__),B_logApecN_1(i__),
     +  logApecNMean(i__),s_logApecN(i__),logBkgNMed(i__),
     +  b_logBkgN(i__),B_logBkgN_1(i__),logBkgNMean(i__),
     +  s_logBkgN(i__),chi25(i__),dof25(i__),logZ(i__)
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(i__),z_1(i__),GammaKL(i__),b_GammaH(i__),
     +  B_GammaH_1(i__),lognHKL_1(i__),b_lognHH_2(i__),
     +  B_lognHH_3(i__),logPowNKL(i__),b_logPowNH(i__),
     +  B_logPowNH_1(i__),GammaBF(i__),lognHBF_1(i__),logPowNBF(i__),
     +  FluxObsBFS_1(i__),FluxObsMedS_3(i__),b_FluxObsS_6(i__),
     +  B_FluxObsS_7(i__),b2_FluxObsS_2(i__),B2_FluxObsS_3(i__),
     +  FluxObsBFT_1(i__),FluxObsMedT_3(i__),b_FluxObsT_6(i__),
     +  B_FluxObsT_7(i__),b2_FluxObsT_2(i__),B2_FluxObsT_3(i__),
     +  FluxcBFS(i__),FluxcMedS_1(i__),b_FluxcS_2(i__),
     +  B_FluxcS_3(i__),b2_FluxcS(i__),B2_FluxcS_1(i__),FluxcBFT(i__),
     +  FluxcMedT_1(i__),b_FluxcT_2(i__),B_FluxcT_3(i__),
     +  b2_FluxcT(i__),B2_FluxcT_1(i__),FluxIBFS(i__),
     +  FluxIMedS_1(i__),b_FluxIS_2(i__),B_FluxIS_3(i__),
     +  b2_FluxIS(i__),B2_FluxIS_1(i__),FluxIBF2keV(i__),
     +  FluxIMed2keV_1(i__),b_FluxI2keV_2(i__),B_FluxI2keV_3(i__),
     +  b2_FluxI2keV(i__),B2_FluxI2keV_1(i__),LumIBFS(i__),
     +  LumIMedS_1(i__),b_LumIS_2(i__),B_LumIS_3(i__),b2_LumIS(i__),
     +  B2_LumIS_1(i__),LumIBF2keV(i__),LumIMed2keV_1(i__),
     +  b_LumI2keV_2(i__),B_LumI2keV_3(i__),b2_LumI2keV(i__),
     +  B2_LumI2keV_1(i__),GammaMed(i__),b_Gamma(i__),B_Gamma_1(i__),
     +  GammaMean(i__),s_Gamma(i__),lognHMed_1(i__),b_lognH_2(i__),
     +  B_lognH_3(i__),lognHMean_1(i__),s_lognH_1(i__),
     +  logPowNMed(i__),b_logPowN(i__),B_logPowN_1(i__),
     +  logPowNMean(i__),s_logPowN(i__),logBkgNMed_1(i__),
     +  b_logBkgN_2(i__),B_logBkgN_3(i__),logBkgNMean_1(i__),
     +  s_logBkgN_1(i__),chi25_1(i__),dof25_1(i__),logZ_1(i__),
     +  NHclass_1(i__)
c    ..............Just test output...........
        write(6,4)
     +  ID_SRC_3(i__),z_1(i__),GammaKL(i__),b_GammaH(i__),
     +  B_GammaH_1(i__),lognHKL_1(i__),b_lognHH_2(i__),
     +  B_lognHH_3(i__),logPowNKL(i__),b_logPowNH(i__),
     +  B_logPowNH_1(i__),GammaBF(i__),lognHBF_1(i__),logPowNBF(i__),
     +  FluxObsBFS_1(i__),FluxObsMedS_3(i__),b_FluxObsS_6(i__),
     +  B_FluxObsS_7(i__),b2_FluxObsS_2(i__),B2_FluxObsS_3(i__),
     +  FluxObsBFT_1(i__),FluxObsMedT_3(i__),b_FluxObsT_6(i__),
     +  B_FluxObsT_7(i__),b2_FluxObsT_2(i__),B2_FluxObsT_3(i__),
     +  FluxcBFS(i__),FluxcMedS_1(i__),b_FluxcS_2(i__),
     +  B_FluxcS_3(i__),b2_FluxcS(i__),B2_FluxcS_1(i__),FluxcBFT(i__),
     +  FluxcMedT_1(i__),b_FluxcT_2(i__),B_FluxcT_3(i__),
     +  b2_FluxcT(i__),B2_FluxcT_1(i__),FluxIBFS(i__),
     +  FluxIMedS_1(i__),b_FluxIS_2(i__),B_FluxIS_3(i__),
     +  b2_FluxIS(i__),B2_FluxIS_1(i__),FluxIBF2keV(i__),
     +  FluxIMed2keV_1(i__),b_FluxI2keV_2(i__),B_FluxI2keV_3(i__),
     +  b2_FluxI2keV(i__),B2_FluxI2keV_1(i__),LumIBFS(i__),
     +  LumIMedS_1(i__),b_LumIS_2(i__),B_LumIS_3(i__),b2_LumIS(i__),
     +  B2_LumIS_1(i__),LumIBF2keV(i__),LumIMed2keV_1(i__),
     +  b_LumI2keV_2(i__),B_LumI2keV_3(i__),b2_LumI2keV(i__),
     +  B2_LumI2keV_1(i__),GammaMed(i__),b_Gamma(i__),B_Gamma_1(i__),
     +  GammaMean(i__),s_Gamma(i__),lognHMed_1(i__),b_lognH_2(i__),
     +  B_lognH_3(i__),lognHMean_1(i__),s_lognH_1(i__),
     +  logPowNMed(i__),b_logPowN(i__),B_logPowN_1(i__),
     +  logPowNMean(i__),s_logPowN(i__),logBkgNMed_1(i__),
     +  b_logBkgN_2(i__),B_logBkgN_3(i__),logBkgNMean_1(i__),
     +  s_logBkgN_1(i__),chi25_1(i__),dof25_1(i__),logZ_1(i__),
     +  NHclass_1(i__)
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(i__),z_2(i__),dGmKL(i__),b_dGmH(i__),B_dGmH_1(i__),
     +  GammaKL_1(i__),b_GammaH_2(i__),B_GammaH_3(i__),lognHKL_2(i__),
     +  b_lognHH_4(i__),B_lognHH_5(i__),logPowNKL_1(i__),
     +  b_logPowNH_2(i__),B_logPowNH_3(i__),logFracKL(i__),
     +  b_logFracH(i__),B_logFracH_1(i__),dGmBF(i__),GammaBF_1(i__),
     +  lognHBF_2(i__),logPowNBF_1(i__),logFracBF(i__),
     +  FluxObsBFS_2(i__),FluxObsMedS_4(i__),b_FluxObsS_8(i__),
     +  B_FluxObsS_9(i__),b2_FluxObsS_4(i__),B2_FluxObsS_5(i__),
     +  FluxObsBFT_2(i__),FluxObsMedT_4(i__),b_FluxObsT_8(i__),
     +  B_FluxObsT_9(i__),b2_FluxObsT_4(i__),B2_FluxObsT_5(i__),
     +  FluxCBFS_1(i__),FluxCMedS_2(i__),b_FluxCS_4(i__),
     +  B_FluxCS_5(i__),b2_FluxCS_2(i__),B2_FluxCS_3(i__),
     +  FluxCBFT_1(i__),FluxCMedT_2(i__),b_FluxCT_4(i__),
     +  B_FluxCT_5(i__),b2_FluxCT_2(i__),B2_FluxCT_3(i__),
     +  FluxIBFS_1(i__),FluxIMedS_2(i__),b_FluxIS_4(i__),
     +  B_FluxIS_5(i__),b2_FluxIS_2(i__),B2_FluxIS_3(i__),
     +  FluxIBF2keV_1(i__),FluxIMed2keV_2(i__),b_FluxI2keV_4(i__),
     +  B_FluxI2keV_5(i__),b2_FluxI2keV_2(i__),B2_FluxI2keV_3(i__),
     +  LumIBFS_1(i__),LumIMedS_2(i__),b_LumIS_4(i__),B_LumIS_5(i__),
     +  b2_LumIS_2(i__),B2_LumIS_3(i__),LumIBF2keV_1(i__),
     +  LumIMed2keV_2(i__),b_LumI2keV_4(i__),B_LumI2keV_5(i__),
     +  b2_LumI2keV_2(i__),B2_LumI2keV_3(i__),dGmMed(i__),b_dGm(i__),
     +  B_dGm_1(i__),dGmMean(i__),s_dGm(i__),GammaMed_1(i__),
     +  b_Gamma_2(i__),B_Gamma_3(i__),GammaMean_1(i__),s_Gamma_1(i__),
     +  lognHMed_2(i__),b_lognH_4(i__),B_lognH_5(i__),
     +  lognHMean_2(i__),s_lognH_2(i__),logPowNMed_1(i__),
     +  b_logPowN_2(i__),B_logPowN_3(i__),logPowNMean_1(i__),
     +  s_logPowN_1(i__),logFracMed(i__),b_logFrac(i__),
     +  B_logFrac_1(i__),logFracMean(i__),s_logFrac(i__),
     +  logBkgNMed_2(i__),b_logBkgN_4(i__),B_logBkgN_5(i__),
     +  logBkgNMean_2(i__),s_logBkgN_2(i__),chi25_2(i__),dof25_2(i__),
     +  logZ_2(i__)
        if(ar__4(51:63) .EQ. '') GammaKL_1(i__) = rNULL__
        if(ar__4(152:164) .EQ. '') logFracKL(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  ID_SRC_4(i__),z_2(i__),dGmKL(i__),b_dGmH(i__),B_dGmH_1(i__),
     +  GammaKL_1(i__),b_GammaH_2(i__),B_GammaH_3(i__),lognHKL_2(i__),
     +  b_lognHH_4(i__),B_lognHH_5(i__),logPowNKL_1(i__),
     +  b_logPowNH_2(i__),B_logPowNH_3(i__),logFracKL(i__),
     +  b_logFracH(i__),B_logFracH_1(i__),dGmBF(i__),GammaBF_1(i__),
     +  lognHBF_2(i__),logPowNBF_1(i__),logFracBF(i__),
     +  FluxObsBFS_2(i__),FluxObsMedS_4(i__),b_FluxObsS_8(i__),
     +  B_FluxObsS_9(i__),b2_FluxObsS_4(i__),B2_FluxObsS_5(i__),
     +  FluxObsBFT_2(i__),FluxObsMedT_4(i__),b_FluxObsT_8(i__),
     +  B_FluxObsT_9(i__),b2_FluxObsT_4(i__),B2_FluxObsT_5(i__),
     +  FluxCBFS_1(i__),FluxCMedS_2(i__),b_FluxCS_4(i__),
     +  B_FluxCS_5(i__),b2_FluxCS_2(i__),B2_FluxCS_3(i__),
     +  FluxCBFT_1(i__),FluxCMedT_2(i__),b_FluxCT_4(i__),
     +  B_FluxCT_5(i__),b2_FluxCT_2(i__),B2_FluxCT_3(i__),
     +  FluxIBFS_1(i__),FluxIMedS_2(i__),b_FluxIS_4(i__),
     +  B_FluxIS_5(i__),b2_FluxIS_2(i__),B2_FluxIS_3(i__),
     +  FluxIBF2keV_1(i__),FluxIMed2keV_2(i__),b_FluxI2keV_4(i__),
     +  B_FluxI2keV_5(i__),b2_FluxI2keV_2(i__),B2_FluxI2keV_3(i__),
     +  LumIBFS_1(i__),LumIMedS_2(i__),b_LumIS_4(i__),B_LumIS_5(i__),
     +  b2_LumIS_2(i__),B2_LumIS_3(i__),LumIBF2keV_1(i__),
     +  LumIMed2keV_2(i__),b_LumI2keV_4(i__),B_LumI2keV_5(i__),
     +  b2_LumI2keV_2(i__),B2_LumI2keV_3(i__),dGmMed(i__),b_dGm(i__),
     +  B_dGm_1(i__),dGmMean(i__),s_dGm(i__),GammaMed_1(i__),
     +  b_Gamma_2(i__),B_Gamma_3(i__),GammaMean_1(i__),s_Gamma_1(i__),
     +  lognHMed_2(i__),b_lognH_4(i__),B_lognH_5(i__),
     +  lognHMean_2(i__),s_lognH_2(i__),logPowNMed_1(i__),
     +  b_logPowN_2(i__),B_logPowN_3(i__),logPowNMean_1(i__),
     +  s_logPowN_1(i__),logFracMed(i__),b_logFrac(i__),
     +  B_logFrac_1(i__),logFracMean(i__),s_logFrac(i__),
     +  logBkgNMed_2(i__),b_logBkgN_4(i__),B_logBkgN_5(i__),
     +  logBkgNMean_2(i__),s_logBkgN_2(i__),chi25_2(i__),dof25_2(i__),
     +  logZ_2(i__)
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(i__),z_3(i__),GammaKL_2(i__),b_GammaH_4(i__),
     +  B_GammaH_5(i__),lognHKL_3(i__),b_lognHH_6(i__),
     +  B_lognHH_7(i__),logPowNKL_2(i__),b_logPowNH_4(i__),
     +  B_logPowNH_5(i__),logkTKL_1(i__),b_logkTH_2(i__),
     +  B_logkTH_3(i__),logBBNKL(i__),b_logBBNH(i__),B_logBBNH_1(i__),
     +  GammaBF_2(i__),lognHBF_3(i__),logPowNBF_2(i__),logkTBF_1(i__),
     +  logBBNBF(i__),FluxObsBFS_3(i__),FluxObsMedS_5(i__),
     +  b_FluxObsS_10(i__),B_FluxObsS_11(i__),b2_FluxObsS_6(i__),
     +  B2_FluxObsS_7(i__),FluxObsBFT_3(i__),FluxObsMedT_5(i__),
     +  b_FluxObsT_10(i__),B_FluxObsT_11(i__),b2_FluxObsT_6(i__),
     +  B2_FluxObsT_7(i__),FluxCBFS_2(i__),FluxCMedS_3(i__),
     +  b_FluxCS_6(i__),B_FluxCS_7(i__),b2_FluxCS_4(i__),
     +  B2_FluxCS_5(i__),FluxCBFT_2(i__),FluxCMedT_3(i__),
     +  b_FluxCT_6(i__),B_FluxCT_7(i__),b2_FluxCT_4(i__),
     +  B2_FluxCT_5(i__),FluxIBFS_2(i__),FluxIMedS_3(i__),
     +  b_FluxIS_6(i__),B_FluxIS_7(i__),b2_FluxIS_4(i__),
     +  B2_FluxIS_5(i__),FluxIBF2keV_2(i__),FluxIMed2keV_3(i__),
     +  b_FluxI2keV_6(i__),B_FluxI2keV_7(i__),b2_FluxI2keV_4(i__),
     +  B2_FluxI2keV_5(i__),LumIBFS_2(i__),LumIMedS_3(i__),
     +  b_LumIS_6(i__),B_LumIS_7(i__),b2_LumIS_4(i__),B2_LumIS_5(i__),
     +  LumIBF2keV_2(i__),LumIMed2keV_3(i__),b_LumI2keV_6(i__),
     +  B_LumI2keV_7(i__),b2_LumI2keV_4(i__),B2_LumI2keV_5(i__),
     +  GammaMed_2(i__),b_Gamma_4(i__),B_Gamma_5(i__),
     +  GammaMean_2(i__),s_Gamma_2(i__),lognHMed_3(i__),
     +  b_lognH_6(i__),B_lognH_7(i__),lognHMean_3(i__),s_lognH_3(i__),
     +  logPowNMed_2(i__),b_logPowN_4(i__),B_logPowN_5(i__),
     +  logPowNMean_2(i__),s_logPowN_2(i__),logkTMed_1(i__),
     +  b_logkT_2(i__),B_logkT_3(i__),logkTMean_1(i__),s_logkT_1(i__),
     +  logBBNMed(i__),b_logBBN(i__),B_logBBN_1(i__),logBBNMean(i__),
     +  s_logBBN(i__),logBkgNMed_3(i__),b_logBkgN_6(i__),
     +  B_logBkgN_7(i__),logBkgNMean_3(i__),s_logBkgN_3(i__),
     +  chi25_3(i__),dof25_3(i__),logZ_3(i__),logRbb(i__)
c    ..............Just test output...........
        write(6,6)
     +  ID_SRC_5(i__),z_3(i__),GammaKL_2(i__),b_GammaH_4(i__),
     +  B_GammaH_5(i__),lognHKL_3(i__),b_lognHH_6(i__),
     +  B_lognHH_7(i__),logPowNKL_2(i__),b_logPowNH_4(i__),
     +  B_logPowNH_5(i__),logkTKL_1(i__),b_logkTH_2(i__),
     +  B_logkTH_3(i__),logBBNKL(i__),b_logBBNH(i__),B_logBBNH_1(i__),
     +  GammaBF_2(i__),lognHBF_3(i__),logPowNBF_2(i__),logkTBF_1(i__),
     +  logBBNBF(i__),FluxObsBFS_3(i__),FluxObsMedS_5(i__),
     +  b_FluxObsS_10(i__),B_FluxObsS_11(i__),b2_FluxObsS_6(i__),
     +  B2_FluxObsS_7(i__),FluxObsBFT_3(i__),FluxObsMedT_5(i__),
     +  b_FluxObsT_10(i__),B_FluxObsT_11(i__),b2_FluxObsT_6(i__),
     +  B2_FluxObsT_7(i__),FluxCBFS_2(i__),FluxCMedS_3(i__),
     +  b_FluxCS_6(i__),B_FluxCS_7(i__),b2_FluxCS_4(i__),
     +  B2_FluxCS_5(i__),FluxCBFT_2(i__),FluxCMedT_3(i__),
     +  b_FluxCT_6(i__),B_FluxCT_7(i__),b2_FluxCT_4(i__),
     +  B2_FluxCT_5(i__),FluxIBFS_2(i__),FluxIMedS_3(i__),
     +  b_FluxIS_6(i__),B_FluxIS_7(i__),b2_FluxIS_4(i__),
     +  B2_FluxIS_5(i__),FluxIBF2keV_2(i__),FluxIMed2keV_3(i__),
     +  b_FluxI2keV_6(i__),B_FluxI2keV_7(i__),b2_FluxI2keV_4(i__),
     +  B2_FluxI2keV_5(i__),LumIBFS_2(i__),LumIMedS_3(i__),
     +  b_LumIS_6(i__),B_LumIS_7(i__),b2_LumIS_4(i__),B2_LumIS_5(i__),
     +  LumIBF2keV_2(i__),LumIMed2keV_3(i__),b_LumI2keV_6(i__),
     +  B_LumI2keV_7(i__),b2_LumI2keV_4(i__),B2_LumI2keV_5(i__),
     +  GammaMed_2(i__),b_Gamma_4(i__),B_Gamma_5(i__),
     +  GammaMean_2(i__),s_Gamma_2(i__),lognHMed_3(i__),
     +  b_lognH_6(i__),B_lognH_7(i__),lognHMean_3(i__),s_lognH_3(i__),
     +  logPowNMed_2(i__),b_logPowN_4(i__),B_logPowN_5(i__),
     +  logPowNMean_2(i__),s_logPowN_2(i__),logkTMed_1(i__),
     +  b_logkT_2(i__),B_logkT_3(i__),logkTMean_1(i__),s_logkT_1(i__),
     +  logBBNMed(i__),b_logBBN(i__),B_logBBN_1(i__),logBBNMean(i__),
     +  s_logBBN(i__),logBkgNMed_3(i__),b_logBkgN_6(i__),
     +  B_logBkgN_7(i__),logBkgNMean_3(i__),s_logBkgN_3(i__),
     +  chi25_3(i__),dof25_3(i__),logZ_3(i__),logRbb(i__)
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(i__),z_4(i__),lognHKL_4(i__),b_lognHH_8(i__),
     +  B_lognHH_9(i__),logPowNKL_3(i__),b_logPowNH_6(i__),
     +  B_logPowNH_7(i__),lognHBF_4(i__),logPowNBF_3(i__),
     +  FluxObsBFS_4(i__),FluxObsMedS_6(i__),b_FluxObsS_12(i__),
     +  B_FluxObsS_13(i__),b2_FluxObsS_8(i__),B2_FluxObsS_9(i__),
     +  FluxObsBFT_4(i__),FluxObsMedT_6(i__),b_FluxObsT_12(i__),
     +  B_FluxObsT_13(i__),b2_FluxObsT_8(i__),B2_FluxObsT_9(i__),
     +  FluxCBFS_3(i__),FluxCMedS_4(i__),b_FluxCS_8(i__),
     +  B_FluxCS_9(i__),b2_FluxCS_6(i__),B2_FluxCS_7(i__),
     +  FluxCBFT_3(i__),FluxCMedT_4(i__),b_FluxCT_8(i__),
     +  B_FluxCT_9(i__),b2_FluxCT_6(i__),B2_FluxCT_7(i__),
     +  FluxIBFS_3(i__),FluxIMedS_4(i__),b_FluxIS_8(i__),
     +  B_FluxIS_9(i__),b2_FluxIS_6(i__),B2_FluxIS_7(i__),
     +  FluxIBF2keV_3(i__),FluxIMed2keV_4(i__),b_FluxI2keV_8(i__),
     +  B_FluxI2keV_9(i__),b2_FluxI2keV_6(i__),B2_FluxI2keV_7(i__),
     +  LumIBFS_3(i__),LumIMedS_4(i__),b_LumIS_8(i__),B_LumIS_9(i__),
     +  b2_LumIS_6(i__),B2_LumIS_7(i__),LumIBF2keV_3(i__),
     +  LumIMed2keV_4(i__),b_LumI2keV_8(i__),B_LumI2keV_9(i__),
     +  b2_LumI2keV_6(i__),B2_LumI2keV_7(i__),lognHMed_4(i__),
     +  b_lognH_8(i__),B_lognH_9(i__),lognHMean_4(i__),s_lognH_4(i__),
     +  logPowNMed_3(i__),b_logPowN_6(i__),B_logPowN_7(i__),
     +  logPowNMean_3(i__),s_logPowN_3(i__),logBkgNMed_4(i__),
     +  b_logBkgN_8(i__),B_logBkgN_9(i__),logBkgNMean_4(i__),
     +  s_logBkgN_4(i__),logZ_4(i__)
c    ..............Just test output...........
        write(6,7)
     +  ID_SRC_6(i__),z_4(i__),lognHKL_4(i__),b_lognHH_8(i__),
     +  B_lognHH_9(i__),logPowNKL_3(i__),b_logPowNH_6(i__),
     +  B_logPowNH_7(i__),lognHBF_4(i__),logPowNBF_3(i__),
     +  FluxObsBFS_4(i__),FluxObsMedS_6(i__),b_FluxObsS_12(i__),
     +  B_FluxObsS_13(i__),b2_FluxObsS_8(i__),B2_FluxObsS_9(i__),
     +  FluxObsBFT_4(i__),FluxObsMedT_6(i__),b_FluxObsT_12(i__),
     +  B_FluxObsT_13(i__),b2_FluxObsT_8(i__),B2_FluxObsT_9(i__),
     +  FluxCBFS_3(i__),FluxCMedS_4(i__),b_FluxCS_8(i__),
     +  B_FluxCS_9(i__),b2_FluxCS_6(i__),B2_FluxCS_7(i__),
     +  FluxCBFT_3(i__),FluxCMedT_4(i__),b_FluxCT_8(i__),
     +  B_FluxCT_9(i__),b2_FluxCT_6(i__),B2_FluxCT_7(i__),
     +  FluxIBFS_3(i__),FluxIMedS_4(i__),b_FluxIS_8(i__),
     +  B_FluxIS_9(i__),b2_FluxIS_6(i__),B2_FluxIS_7(i__),
     +  FluxIBF2keV_3(i__),FluxIMed2keV_4(i__),b_FluxI2keV_8(i__),
     +  B_FluxI2keV_9(i__),b2_FluxI2keV_6(i__),B2_FluxI2keV_7(i__),
     +  LumIBFS_3(i__),LumIMedS_4(i__),b_LumIS_8(i__),B_LumIS_9(i__),
     +  b2_LumIS_6(i__),B2_LumIS_7(i__),LumIBF2keV_3(i__),
     +  LumIMed2keV_4(i__),b_LumI2keV_8(i__),B_LumI2keV_9(i__),
     +  b2_LumI2keV_6(i__),B2_LumI2keV_7(i__),lognHMed_4(i__),
     +  b_lognH_8(i__),B_lognH_9(i__),lognHMean_4(i__),s_lognH_4(i__),
     +  logPowNMed_3(i__),b_logPowN_6(i__),B_logPowN_7(i__),
     +  logPowNMean_3(i__),s_logPowN_3(i__),logBkgNMed_4(i__),
     +  b_logBkgN_8(i__),B_logBkgN_9(i__),logBkgNMean_4(i__),
     +  s_logBkgN_4(i__),logZ_4(i__)
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(i__),z_5(i__),logPowNKL_4(i__),b_logPowNH_8(i__),
     +  B_logPowNH_9(i__),logPowNBF_4(i__),FluxObsBFS_5(i__),
     +  FluxObsMedS_7(i__),b_FluxObsS_14(i__),B_FluxObsS_15(i__),
     +  b2_FluxObsS_10(i__),B2_FluxObsS_11(i__),FluxObsBFT_5(i__),
     +  FluxObsMedT_7(i__),b_FluxObsT_14(i__),B_FluxObsT_15(i__),
     +  b2_FluxObsT_10(i__),B2_FluxObsT_11(i__),FluxCBFS_4(i__),
     +  FluxCMedS_5(i__),b_FluxCS_10(i__),B_FluxCS_11(i__),
     +  b2_FluxCS_8(i__),B2_FluxCS_9(i__),FluxCBFT_4(i__),
     +  FluxCMedT_5(i__),b_FluxCT_10(i__),B_FluxCT_11(i__),
     +  b2_FluxCT_8(i__),B2_FluxCT_9(i__),FluxIBFS_4(i__),
     +  FluxIMedS_5(i__),b_FluxIS_10(i__),B_FluxIS_11(i__),
     +  b2_FluxIS_8(i__),B2_FluxIS_9(i__),FluxIBF2keV_4(i__),
     +  FluxIMed2keV_5(i__),b_FluxI2keV_10(i__),B_FluxI2keV_11(i__),
     +  b2_FluxI2keV_8(i__),B2_FluxI2keV_9(i__),LumIBFS_4(i__),
     +  LumIMedS_5(i__),b_LumIS_10(i__),B_LumIS_11(i__),
     +  b2_LumI2S(i__),B2_LumIS_8(i__),LumIBF2keV_4(i__),
     +  LumIMed2keV_5(i__),b_LumI2keV_10(i__),B_LumI2keV_11(i__),
     +  b2_LumI2keV_8(i__),B2_LumI2keV_9(i__),logPowNMed_4(i__),
     +  b_logPowN_8(i__),B_logPowN_9(i__),logPowNMean_4(i__),
     +  s_logPowN_4(i__),logBkgNMed_5(i__),b_logBkgN_10(i__),
     +  B_logBkgN_11(i__),logBkgNMean_5(i__),s_logBkgN_5(i__),
     +  logZ_5(i__)
c    ..............Just test output...........
        write(6,8)
     +  ID_SRC_7(i__),z_5(i__),logPowNKL_4(i__),b_logPowNH_8(i__),
     +  B_logPowNH_9(i__),logPowNBF_4(i__),FluxObsBFS_5(i__),
     +  FluxObsMedS_7(i__),b_FluxObsS_14(i__),B_FluxObsS_15(i__),
     +  b2_FluxObsS_10(i__),B2_FluxObsS_11(i__),FluxObsBFT_5(i__),
     +  FluxObsMedT_7(i__),b_FluxObsT_14(i__),B_FluxObsT_15(i__),
     +  b2_FluxObsT_10(i__),B2_FluxObsT_11(i__),FluxCBFS_4(i__),
     +  FluxCMedS_5(i__),b_FluxCS_10(i__),B_FluxCS_11(i__),
     +  b2_FluxCS_8(i__),B2_FluxCS_9(i__),FluxCBFT_4(i__),
     +  FluxCMedT_5(i__),b_FluxCT_10(i__),B_FluxCT_11(i__),
     +  b2_FluxCT_8(i__),B2_FluxCT_9(i__),FluxIBFS_4(i__),
     +  FluxIMedS_5(i__),b_FluxIS_10(i__),B_FluxIS_11(i__),
     +  b2_FluxIS_8(i__),B2_FluxIS_9(i__),FluxIBF2keV_4(i__),
     +  FluxIMed2keV_5(i__),b_FluxI2keV_10(i__),B_FluxI2keV_11(i__),
     +  b2_FluxI2keV_8(i__),B2_FluxI2keV_9(i__),LumIBFS_4(i__),
     +  LumIMedS_5(i__),b_LumIS_10(i__),B_LumIS_11(i__),
     +  b2_LumI2S(i__),B2_LumIS_8(i__),LumIBF2keV_4(i__),
     +  LumIMed2keV_5(i__),b_LumI2keV_10(i__),B_LumI2keV_11(i__),
     +  b2_LumI2keV_8(i__),B2_LumI2keV_9(i__),logPowNMed_4(i__),
     +  b_logPowN_8(i__),B_logPowN_9(i__),logPowNMean_4(i__),
     +  s_logPowN_4(i__),logBkgNMed_5(i__),b_logBkgN_10(i__),
     +  B_logBkgN_11(i__),logBkgNMean_5(i__),s_logBkgN_5(i__),
     +  logZ_5(i__)
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(i__),z_6(i__),GammaKL_3(i__),b_GammaH_6(i__),
     +  B_GammaH_7(i__),lognHKL_5(i__),b_lognHH_10(i__),
     +  B_lognHH_11(i__),logPowNKL_5(i__),b_logPowNH_10(i__),
     +  B_logPowNH_11(i__),GammaBF_3(i__),lognHBF_5(i__),
     +  logPowNBF_5(i__),FluxObsBFS_6(i__),FluxObsMedS_8(i__),
     +  b_FluxObsS_16(i__),B_FluxObsS_17(i__),b2_FluxObsS_12(i__),
     +  B2_FluxObsS_13(i__),GammaMed_3(i__),b_Gamma_6(i__),
     +  B_Gamma_7(i__),GammaMean_3(i__),s_Gamma_3(i__),
     +  lognHMed_5(i__),b_lognH_10(i__),B_lognH_11(i__),
     +  lognHMean_5(i__),s_lognH_5(i__),logPowNMed_5(i__),
     +  b_logPowN_10(i__),B_logPowN_11(i__),logPowNMean_5(i__),
     +  s_logPowN_5(i__),logBkgNMed_6(i__),blogBkgN(i__),
     +  BlogBkgN_1(i__),logBkgNMean_6(i__),logBkgN(i__),logZ_6(i__)
c    ..............Just test output...........
        write(6,9)
     +  ID_SRC_8(i__),z_6(i__),GammaKL_3(i__),b_GammaH_6(i__),
     +  B_GammaH_7(i__),lognHKL_5(i__),b_lognHH_10(i__),
     +  B_lognHH_11(i__),logPowNKL_5(i__),b_logPowNH_10(i__),
     +  B_logPowNH_11(i__),GammaBF_3(i__),lognHBF_5(i__),
     +  logPowNBF_5(i__),FluxObsBFS_6(i__),FluxObsMedS_8(i__),
     +  b_FluxObsS_16(i__),B_FluxObsS_17(i__),b2_FluxObsS_12(i__),
     +  B2_FluxObsS_13(i__),GammaMed_3(i__),b_Gamma_6(i__),
     +  B_Gamma_7(i__),GammaMean_3(i__),s_Gamma_3(i__),
     +  lognHMed_5(i__),b_lognH_10(i__),B_lognH_11(i__),
     +  lognHMean_5(i__),s_lognH_5(i__),logPowNMed_5(i__),
     +  b_logPowN_10(i__),B_logPowN_11(i__),logPowNMean_5(i__),
     +  s_logPowN_5(i__),logBkgNMed_6(i__),blogBkgN(i__),
     +  BlogBkgN_1(i__),logBkgNMean_6(i__),logBkgN(i__),logZ_6(i__)
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(i__),z_7(i__),GammaKL_4(i__),e_Gamma(i__),
     +  E_Gamma_1(i__),lognHKL_6(i__),b_lognH_12(i__),B_lognH_13(i__),
     +  logPowNKL_6(i__),b_logPowN_12(i__),B_logPowN_13(i__),
     +  GammaBF_4(i__),lognHBF_6(i__),logPowNBF_6(i__),
     +  FluxObsBFT_6(i__),FluxObsMedT_8(i__),b_FluxObsT_16(i__),
     +  B_FluxObsT_17(i__),b2_FluxObsT_12(i__),B2_FluxObsT_13(i__),
     +  FluxCBFT_5(i__),FluxCMedT_6(i__),b_FluxCMedT(i__),
     +  B_FluxCMedT_1(i__),b2_FluxCMedT(i__),B2_FluxCMedT_1(i__),
     +  GammaMed_4(i__),b_GammaMed(i__),B_GammaMed_1(i__),
     +  GammaMean_4(i__),s_Gamma_4(i__),lognHMed_6(i__),
     +  b_lognHMed(i__),B_lognHMed_1(i__),lognHMean_6(i__),
     +  s_lognH_6(i__),logPowNMed_6(i__),b_logPowNMed(i__),
     +  B_logPowNMed_1(i__),logPowNMean_6(i__),s_logPowN_6(i__),
     +  logBkgNMed_7(i__),b_logBkgNMed(i__),B_logBkgNMed_1(i__),
     +  logBkgNMean_7(i__),s_logBkgN_6(i__),logZ_7(i__)
c    ..............Just test output...........
        write(6,10)
     +  ID_SRC_9(i__),z_7(i__),GammaKL_4(i__),e_Gamma(i__),
     +  E_Gamma_1(i__),lognHKL_6(i__),b_lognH_12(i__),B_lognH_13(i__),
     +  logPowNKL_6(i__),b_logPowN_12(i__),B_logPowN_13(i__),
     +  GammaBF_4(i__),lognHBF_6(i__),logPowNBF_6(i__),
     +  FluxObsBFT_6(i__),FluxObsMedT_8(i__),b_FluxObsT_16(i__),
     +  B_FluxObsT_17(i__),b2_FluxObsT_12(i__),B2_FluxObsT_13(i__),
     +  FluxCBFT_5(i__),FluxCMedT_6(i__),b_FluxCMedT(i__),
     +  B_FluxCMedT_1(i__),b2_FluxCMedT(i__),B2_FluxCMedT_1(i__),
     +  GammaMed_4(i__),b_GammaMed(i__),B_GammaMed_1(i__),
     +  GammaMean_4(i__),s_Gamma_4(i__),lognHMed_6(i__),
     +  b_lognHMed(i__),B_lognHMed_1(i__),lognHMean_6(i__),
     +  s_lognH_6(i__),logPowNMed_6(i__),b_logPowNMed(i__),
     +  B_logPowNMed_1(i__),logPowNMean_6(i__),s_logPowN_6(i__),
     +  logBkgNMed_7(i__),b_logBkgNMed(i__),B_logBkgNMed_1(i__),
     +  logBkgNMean_7(i__),s_logBkgN_6(i__),logZ_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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