FORTRAN Generation
(/./ftp/cats/J/A_A/480/599)

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

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-28
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/480/599       XMM-Newton survey of M31 centre field   (Stiele+, 2008)
*================================================================================
*Time variability of X-ray sources in the M31 centre field.
*    Stiele H., Pietsch W., Haberl F., Freyberg M.
*   <Astron. Astrophys. 480, 599 (2008)>
*   =2008A&A...480..599S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! X-ray source catalogue extension of the
                               XMM-Newton EPIC M31 central observations

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

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

      integer*4     Seq         ! [857/895] Source number (1)
      character*6   M31ID       ! M31 field identification
      integer*4     RAh         ! (h) Right Ascension J2000 (hours)
      integer*4     RAm         ! (min) Right Ascension J2000 (minutes)
      real*4        RAs         ! (s) Right Ascension J2000 (seconds)
      character*1   DE_         ! Declination J2000 (sign)
      integer*4     DEd         ! (deg) Declination J2000 (degrees)
      integer*4     DEm         ! (arcmin) Declination J2000 (minutes)
      real*4        DEs         ! (arcsec) Declination J2000 (seconds)
      real*4        e_Pos       ! (arcsec) ? Positional error (2)
      real*4        LH          ! Source detection likelihood
*                                          (combined EPIC)
      real*4        CRate       ! (ct/s) Combined EPIC count rate (G2)
      real*4        e_CRate     ! (ct/s) Error in CRate
      real*4        CFlux       ! (mW/m2) Combined EPIC flux in erg/cm2/s (G2)
      real*4        e_CFlux     ! (mW/m2) Error in CFlux
      real*4        CHR1        ! ? Combined EPIC hardness ratio 1 (G3)
      real*4        e_CHR1      ! ? Error in CHR1
      real*4        CHR2        ! ? Combined EPIC hardness ratio 2 (G3)
      real*4        e_CHR2      ! ? Error in CHR2
      real*4        CHR3        ! ? Combined EPIC hardness ratio 3 (G3)
      real*4        e_CHR3      ! ? Error in CHR3
      real*4        CHR4        ! ? Combined EPIC hardness ratio 4 (G3)
      real*4        e_CHR4      ! ? Error in CHR4
      character*3   VAL         ! EPIC instruments contributing to source
*                                          detection (G4)
      real*4        PNExpo      ! (ks) ? Exposure of EPIC pn instrument
      real*4        PNLH        ! ? EPIC pn source detection likelihood
      real*4        PNRate      ! (ct/s) ? EPIC pn count rate (G2)
      real*4        e_PNRate    ! (ct/s) ? Error in PNRate
      real*4        PNFlux      ! (mW/m2) ? EPIC pn flux in erg/cm^2^/s (G2)
      real*4        e_PNFlux    ! (mW/m2) ? Error in PNFlux
      real*4        PNHR1       ! ? EPIC pn hardness ratio 1 (G3)
      real*4        e_PNHR1     ! ? Error in PNHR1
      real*4        PNHR2       ! ? EPIC pn hardness ratio 2 (G3)
      real*4        e_PNHR2     ! ? Error in PNHR2
      real*4        PNHR3       ! ? EPIC pn hardness ratio 3 (G3)
      real*4        e_PNHR3     ! ? Error in PNHR3
      real*4        PNHR4       ! ? EPIC pn hardness ratio 4 (G3)
      real*4        e_PNHR4     ! ? Error in PNHR4
      real*4        M1Expo      ! (ks) ? Exposure of EPIC MOS1 instrument
      real*4        M1LH        ! ? EPIC MOS1 source detection likelihood
      real*4        M1Rate      ! (ct/s) ? EPIC MOS1 count rate (G2)
      real*4        e_M1Rate    ! (ct/s) ? Error in M1Rate
      real*4        M1Flux      ! (mW/m2) ? EPIC MOS1 flux in erg/cm^2^/s (G2)
      real*4        e_M1Flux    ! (mW/m2) ? Error in M1Flux
      real*4        M1HR1       ! ? EPIC MOS1 hardness ratio 1 (G3)
      real*4        e_M1HR1     ! ? Error in M1HR1
      real*4        M1HR2       ! ? EPIC MOS1 hardness ratio 2 (G3)
      real*4        e_M1HR2     ! ? Error in M1HR2
      real*4        M1HR3       ! ? EPIC MOS1 hardness ratio 3 (G3)
      real*4        e_M1HR3     ! ? Error in M1HR3
      real*4        M1HR4       ! ? EPIC MOS1 hardness ratio 4 (G3)
      real*4        e_M1HR4     ! ? Error in M1HR4
      real*4        M2Expo      ! (ks) ? Exposure of EPIC MOS2 instrument
      real*4        M2LH        ! ? EPIC MOS2 source detection likelihood
      real*4        M2Rate      ! (ct/s) ? EPIC MOS2 count rate (G2)
      real*4        e_M2Rate    ! (ct/s) ? Error in M2Rate
      real*4        M2Flux      ! (mW/m2) ? EPIC MOS2 flux in erg/cm^2^/s (G2)
      real*4        e_M2Flux    ! (mW/m2) ? Error in M2Flux
      real*4        M2HR1       ! ? EPIC MOS2 hardness ratio 1 (G3)
      real*4        e_M2HR1     ! ? Error in M2HR1
      real*4        M2HR2       ! ? EPIC MOS2 hardness ratio 2 (G3)
      real*4        e_M2HR2     ! ? Error in M2HR2
      real*4        M2HR3       ! ? EPIC MOS2 hardness ratio 3 (G3)
      real*4        e_M2HR3     ! ? Error in M2HR3
      real*4        M2HR4       ! ? EPIC MOS2 hardness ratio 4 (G3)
      real*4        e_M2HR4     ! ? Error in M2HR4
      character*100 XID         ! ? X-ray identification (G5)
      character*12  USNOB1_0    ! ? Name of brightest USNO-B1.0
*                                       source within search area
      integer*4     USNOMul     ! ? Number of USNO-B1.0 correlations
      real*4        USNOdist    ! (arcsec) ? Distance of source USNO-B1.0
      real*4        USNOB2      ! (mag) ? B2 magnitude of source USNO-B1.0
      real*4        USNOR2      ! (mag) ? R2 magnitude of source USNO-B1.0
      real*4        USNOI       ! (mag) ? I magnitude of source USNO-B1.0
      real*4        log_fX_fopt ! ? Logarithm of fX/fopt
      character*15  Class       ! ? Classification of X-ray source (G6)
      character*100 Remarks     ! ? Additional remarks, in refs.dat files
      character*16  XMMM31      ! Source Name XMMM31 (Jhhmmss.s+ddmmss)
*Note (1): The numbering is a continuation of Pietsch et al., 2005,
*          Cat. J/A+A/434/483, [PFH2005] NNN in Simbad)
*Note (2): 1{sigma} statistical error, 0.5" to be added as 1{sigma}
*          systematic error

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

C  Declarations for 'table5.dat'	! X-ray variability of the central sources in M31
                               detected in XMM-Newton EPIC observations

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

      integer*4     SRCID       ! [136/891] Source Number (1)
      real*4        CFlux_1     ! (mW/m2) Combined EPIC flux in erg/cm2/s (G2)
      real*4        e_CFlux_1   ! (mW/m2) Error in CFlux
      real*4        CHR1_1      ! ? Combined EPIC hardness ratio 1 (G3)
      real*4        e_CHR1_1    ! ? Error in CHR1
      real*4        CHR2_1      ! ? Combined EPIC hardness ratio 2 (G3)
      real*4        e_CHR2_1    ! ? Error in CHR2
      real*4        CHR3_1      ! ? Combined EPIC hardness ratio 3 (G3)
      real*4        e_CHR3_1    ! ? Error in CHR3
      real*4        CHR4_1      ! ? Combined EPIC hardness ratio 4 (G3)
      real*4        e_CHR4_1    ! ? Error in CHR4
      character*156 XID_1       ! ? X-ray identification (G5)
      character*15  Class_1     ! ? Classification of X-ray source (G6)
      character*100 Remarks_1   ! ? Additional remarks, in refs.dat file
      character*3   c1val       ! ? EPIC instruments contributing to source
*                                       detection in observation c1 (G4)
      real*4        c1Rate      ! (ct/s) ? EPIC count rate in observation c1 (G2)
      real*4        e_c1Rate    ! (ct/s) ? Error in c1Rate
      real*4        c1Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation c1 (G2)
      real*4        e_c1Flux    ! (mW/m2) ? Error in c1Flux
      real*4        c1HR1       ! ? EPIC hardness ratio 1 in observation c1 (G3)
      real*4        e_c1HR1     ! ? Error in c1HR1
      real*4        c1HR2       ! ? EPIC hardness ratio 2 in observation c1 (G3)
      real*4        e_c1HR2     ! ? Error in c1HR2
      real*4        c1HR3       ! ? EPIC hardness ratio 3 in observation c1 (G3)
      real*4        e_c1HR3     ! ? Error in c1HR3
      real*4        c1HR4       ! ? EPIC hardness ratio 4 in observation c1 (G3)
      real*4        e_c1HR4     ! ? Error in c1HR4
      character*3   c2val       ! ? EPIC instruments contributing to source
*                                       detection in observation c2 (G4)
      real*4        c2Rate      ! (ct/s) ? EPIC count rate in observation c2 (G2)
      real*4        e_c2Rate    ! (ct/s) ? Error in c2Rate
      real*4        c2Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation c2 (G2)
      real*4        e_c2Flux    ! (mW/m2) ? Error in c2Flux
      real*4        c2HR1       ! ? EPIC hardness ratio 1 in observation c2 (G3)
      real*4        e_c2HR1     ! ? Error in c2HR1
      real*4        c2HR2       ! ? EPIC hardness ratio 2 in observation c2 (G3)
      real*4        e_c2HR2     ! ? Error in c2HR2
      real*4        c2HR3       ! ? EPIC hardness ratio 3 in observation c2 (G3)
      real*4        e_c2HR3     ! ? Error in c2HR3
      real*4        c2HR4       ! ? EPIC hardness ratio 4 in observation c2 (G3)
      real*4        e_c2HR4     ! ? Error in c2HR4
      character*3   c3val       ! ? EPIC instruments contributing to source
*                                        detection in observation c3 (G4)
      real*4        c3Rate      ! (ct/s) ? EPIC count rate in observation c3 (G2)
      real*4        e_c3Rate    ! (ct/s) ? Error in c3Rate
      real*4        c3Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation c3 (G2)
      real*4        e_c3Flux    ! (mW/m2) ? Error in c3Flux
      real*4        c3HR1       ! ? EPIC hardness ratio 1 in observation c3 (G3)
      real*4        e_c3HR1     ! ? Error in c3HR1
      real*4        c3HR2       ! ? EPIC hardness ratio 2 in observation c3 (G3)
      real*4        e_c3HR2     ! ? Error in c3HR2
      real*4        c3HR3       ! ? EPIC hardness ratio 3 in observation c3 (G3)
      real*4        e_c3HR3     ! ? Error in c3HR3
      real*4        c3HR4       ! ? EPIC hardness ratio 4 in observation c3 (G3)
      real*4        e_c3HR4     ! ? Error in c3HR4
      character*3   n1val       ! ? EPIC instruments contributing to source
*                                        detection in observation n1 (G4)
      real*4        n1Rate      ! (ct/s) ? EPIC count rate in observation n1 (G2)
      real*4        e_n1Rate    ! (ct/s) ? Error in n1Rate
      real*4        n1Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation n1 (G2)
      real*4        e_n1Flux    ! (mW/m2) ? Error in n1Flux
      real*4        n1HR1       ! ? EPIC hardness ratio 1 in observation n1 (G3)
      real*4        e_n1HR1     ! ? Error in n1HR1
      real*4        n1HR2       ! ? EPIC hardness ratio 2 in observation n1 (G3)
      real*4        e_n1HR2     ! ? Error in n1HR2
      real*4        n1HR3       ! ? EPIC hardness ratio 3 in observation n1 (G3)
      real*4        e_n1HR3     ! ? Error in n1HR3
      real*4        n1HR4       ! ? EPIC hardness ratio 4 in observation n1 (G3)
      real*4        e_n1HR4     ! ? Error in n1HR4
      character*3   c4val       ! ? EPIC instruments contributing to source
*                                        detection in observation c4 (G4)
      real*4        c4Rate      ! (ct/s) ? EPIC count rate in observation c4 (G2)
      real*4        e_c4Rate    ! (ct/s) ? Error in c4Rate
      real*4        c4Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation c4 (G2)
      real*4        e_c4Flux    ! (mW/m2) ? Error in c4Flux
      real*4        c4HR1       ! ? EPIC hardness ratio 1 in observation c4 (G3)
      real*4        e_c4HR1     ! ? Error in c4HR1
      real*4        c4HR2       ! ? EPIC hardness ratio 2 in observation c4 (G3)
      real*4        e_c4HR2     ! ? Error in c4HR2
      real*4        c4HR3       ! ? EPIC hardness ratio 3 in observation c4 (G3)
      real*4        e_c4HR3     ! ? Error in c4HR3
      real*4        c4HR4       ! ? EPIC hardness ratio 4 in observation c4 (G3)
      real*4        e_c4HR4     ! ? Error in c4HR4
      character*3   s1val       ! ? EPIC instruments contributing to source
*                                        detection in observation s1 (G4)
      real*4        s1Rate      ! (ct/s) ? EPIC count rate in observation s1 (G2)
      real*4        e_s1Rate    ! (ct/s) ? Error in s1Rate
      real*4        s1Flux      ! (mW/m2) ? EPIC flux in erg/cm2/s in observation s1 (G2)
      real*4        e_s1Flux    ! (mW/m2) ? Error in s1Flux
      real*4        s1HR1       ! ? EPIC hardness ratio 1 in observation s1 (G3)
      real*4        e_s1HR1     ! ? Error in s1HR1
      real*4        s1HR2       ! ? EPIC hardness ratio 2 in observation s1 (G3)
      real*4        e_s1HR2     ! ? Error in s1HR2
      real*4        s1HR3       ! ? EPIC hardness ratio 3 in observation s1 (G3)
      real*4        e_s1HR3     ! ? Error in s1HR3
      real*4        s1HR4       ! ? EPIC hardness ratio 4 in observation s1 (G3)
      real*4        e_s1HR4     ! ? Error in s1HR4
      character*3   bval        ! ? EPIC instruments contributing to source
*                                        detection in observation b (G4)
      real*4        bRate       ! (ct/s) ? EPIC count rate in observation b (G2)
      real*4        e_bRate     ! (ct/s) ? Error in bRate
      real*4        bFlux       ! (mW/m2) ? EPIC flux in erg/cm2/s in observation b (G2)
      real*4        e_bFlux     ! (mW/m2) ? Error in bFlux
      real*4        bHR1        ! ? EPIC hardness ratio 1 in observation b (G3)
      real*4        e_bHR1      ! ? Error in bHR1
      real*4        bHR2        ! ? EPIC hardness ratio 2 in observation b (G3)
      real*4        e_bHR2      ! ? Error in bHR2
      real*4        bHR3        ! ? EPIC hardness ratio 3 in observation b (G3)
      real*4        e_bHR3      ! ? Error in bHR3
      real*4        bHR4        ! ? EPIC hardness ratio 4 in observation b (G3)
      real*4        e_bHR4      ! ? Error in bHR4
      integer*4     ndet        ! Number of observations covering the source
      real*4        svarmax     ! Maximum of significance of variation
      real*4        fvarmax     ! Maximum of flux ratio
      integer*4     nupper      ! Number of upper limits
      real*4        fmax        ! (mW/m2) ? Maximum flux
      real*4        e_fmax      ! (mW/m2) ? Error in fmax
      character*16  XMMM31_1    ! Source Name XMMM31 (Jhhmmss.s+ddmmss)
*Note (1): For numbers 136 to 856, source number from Table 2 in Pietsch et
*     al., 2005, Cat. J/A+A/434/483, [PFH2005] NNN in Simbad.
*     For numbers 857 to 891, source number of Table 3 in this paper

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

C  Declarations for 'refs.dat'	! References

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

      character*8   Ref         ! Reference code
      character*19  BibCode     ! BibCode
      character*27  Name        ! Author's name
      character*81  Com         ! Comments

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

C  Loading file 'table3.dat'	! X-ray source catalogue extension of the
*                               XMM-Newton EPIC M31 central observations

C  Format for file interpretation

    1 format(
     +  I3,1X,A6,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,
     +  E8.2,1X,E8.2,1X,E7.1,1X,E8.2,1X,E7.1,1X,F5.2,1X,F6.2,1X,F5.2,
     +  1X,F6.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3,1X,F5.1,1X,E8.2,
     +  1X,E8.2,1X,E7.1,1X,E8.2,1X,E7.1,1X,F5.2,1X,F6.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,1X,E8.2,1X,E8.2,
     +  1X,E7.1,1X,E8.2,1X,E7.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,1X,E8.2,1X,E8.2,1X,E7.1,
     +  1X,E8.2,1X,E7.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,A100,A12,1X,I1,1X,F3.1,1X,F4.1,1X,
     +  F4.1,1X,F4.1,1X,F4.1,1X,A15,1X,A100,8X,A16)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,39
        read(1,'(A696)')ar__
        read(ar__,1)
     +  Seq,M31ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,e_Pos,LH,CRate,e_CRate,
     +  CFlux,e_CFlux,CHR1,e_CHR1,CHR2,e_CHR2,CHR3,e_CHR3,CHR4,e_CHR4,
     +  VAL,PNExpo,PNLH,PNRate,e_PNRate,PNFlux,e_PNFlux,PNHR1,e_PNHR1,
     +  PNHR2,e_PNHR2,PNHR3,e_PNHR3,PNHR4,e_PNHR4,M1Expo,M1LH,M1Rate,
     +  e_M1Rate,M1Flux,e_M1Flux,M1HR1,e_M1HR1,M1HR2,e_M1HR2,M1HR3,
     +  e_M1HR3,M1HR4,e_M1HR4,M2Expo,M2LH,M2Rate,e_M2Rate,M2Flux,
     +  e_M2Flux,M2HR1,e_M2HR1,M2HR2,e_M2HR2,M2HR3,e_M2HR3,M2HR4,
     +  e_M2HR4,XID,USNOB1_0,USNOMul,USNOdist,USNOB2,USNOR2,USNOI,
     +  log_fX_fopt,Class,Remarks,XMMM31
        if(ar__(36:40) .EQ. '') e_Pos = rNULL__
        if(ar__(85:89) .EQ. '') CHR1 = rNULL__
        if(ar__(91:96) .EQ. '') e_CHR1 = rNULL__
        if(ar__(98:102) .EQ. '') CHR2 = rNULL__
        if(ar__(104:109) .EQ. '') e_CHR2 = rNULL__
        if(ar__(111:115) .EQ. '') CHR3 = rNULL__
        if(ar__(117:120) .EQ. '') e_CHR3 = rNULL__
        if(ar__(122:126) .EQ. '') CHR4 = rNULL__
        if(ar__(128:131) .EQ. '') e_CHR4 = rNULL__
        if(ar__(137:141) .EQ. '') PNExpo = rNULL__
        if(ar__(143:150) .EQ. '') PNLH = rNULL__
        if(ar__(152:159) .EQ. '') PNRate = rNULL__
        if(ar__(161:167) .EQ. '') e_PNRate = rNULL__
        if(ar__(169:176) .EQ. '') PNFlux = rNULL__
        if(ar__(178:184) .EQ. '') e_PNFlux = rNULL__
        if(ar__(186:190) .EQ. '') PNHR1 = rNULL__
        if(ar__(192:197) .EQ. '') e_PNHR1 = rNULL__
        if(ar__(199:203) .EQ. '') PNHR2 = rNULL__
        if(ar__(205:208) .EQ. '') e_PNHR2 = rNULL__
        if(ar__(210:214) .EQ. '') PNHR3 = rNULL__
        if(ar__(216:219) .EQ. '') e_PNHR3 = rNULL__
        if(ar__(221:225) .EQ. '') PNHR4 = rNULL__
        if(ar__(227:230) .EQ. '') e_PNHR4 = rNULL__
        if(ar__(232:236) .EQ. '') M1Expo = rNULL__
        if(ar__(238:245) .EQ. '') M1LH = rNULL__
        if(ar__(247:254) .EQ. '') M1Rate = rNULL__
        if(ar__(256:262) .EQ. '') e_M1Rate = rNULL__
        if(ar__(264:271) .EQ. '') M1Flux = rNULL__
        if(ar__(273:279) .EQ. '') e_M1Flux = rNULL__
        if(ar__(281:285) .EQ. '') M1HR1 = rNULL__
        if(ar__(287:290) .EQ. '') e_M1HR1 = rNULL__
        if(ar__(292:296) .EQ. '') M1HR2 = rNULL__
        if(ar__(298:301) .EQ. '') e_M1HR2 = rNULL__
        if(ar__(303:307) .EQ. '') M1HR3 = rNULL__
        if(ar__(309:312) .EQ. '') e_M1HR3 = rNULL__
        if(ar__(314:318) .EQ. '') M1HR4 = rNULL__
        if(ar__(320:323) .EQ. '') e_M1HR4 = rNULL__
        if(ar__(325:329) .EQ. '') M2Expo = rNULL__
        if(ar__(331:338) .EQ. '') M2LH = rNULL__
        if(ar__(340:347) .EQ. '') M2Rate = rNULL__
        if(ar__(349:355) .EQ. '') e_M2Rate = rNULL__
        if(ar__(357:364) .EQ. '') M2Flux = rNULL__
        if(ar__(366:372) .EQ. '') e_M2Flux = rNULL__
        if(ar__(374:378) .EQ. '') M2HR1 = rNULL__
        if(ar__(380:383) .EQ. '') e_M2HR1 = rNULL__
        if(ar__(385:389) .EQ. '') M2HR2 = rNULL__
        if(ar__(391:394) .EQ. '') e_M2HR2 = rNULL__
        if(ar__(396:400) .EQ. '') M2HR3 = rNULL__
        if(ar__(402:405) .EQ. '') e_M2HR3 = rNULL__
        if(ar__(407:411) .EQ. '') M2HR4 = rNULL__
        if(ar__(413:416) .EQ. '') e_M2HR4 = rNULL__
        if(ar__(531:531) .EQ. '') USNOMul = iNULL__
        if(ar__(533:535) .EQ. '') USNOdist = rNULL__
        if(ar__(537:540) .EQ. '') USNOB2 = rNULL__
        if(ar__(542:545) .EQ. '') USNOR2 = rNULL__
        if(ar__(547:550) .EQ. '') USNOI = rNULL__
        if(ar__(552:555) .EQ. '') log_fX_fopt = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Seq,M31ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,e_Pos,LH,CRate,e_CRate,
     +  CFlux,e_CFlux,CHR1,e_CHR1,CHR2,e_CHR2,CHR3,e_CHR3,CHR4,e_CHR4,
     +  VAL,PNExpo,PNLH,PNRate,e_PNRate,PNFlux,e_PNFlux,PNHR1,e_PNHR1,
     +  PNHR2,e_PNHR2,PNHR3,e_PNHR3,PNHR4,e_PNHR4,M1Expo,M1LH,M1Rate,
     +  e_M1Rate,M1Flux,e_M1Flux,M1HR1,e_M1HR1,M1HR2,e_M1HR2,M1HR3,
     +  e_M1HR3,M1HR4,e_M1HR4,M2Expo,M2LH,M2Rate,e_M2Rate,M2Flux,
     +  e_M2Flux,M2HR1,e_M2HR1,M2HR2,e_M2HR2,M2HR3,e_M2HR3,M2HR4,
     +  e_M2HR4,XID,USNOB1_0,USNOMul,USNOdist,USNOB2,USNOR2,USNOI,
     +  log_fX_fopt,Class,Remarks,XMMM31
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! X-ray variability of the central sources in M31
*                               detected in XMM-Newton EPIC observations

C  Format for file interpretation

    2 format(
     +  I3,1X,E8.3,1X,E7.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,A156,1X,A15,1X,A100,1X,A3,1X,E8.3,1X,
     +  E7.2,1X,E8.3,1X,E7.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F5.2,1X,F4.2,1X,A3,1X,E8.3,1X,E7.2,1X,E8.3,1X,E7.2,
     +  1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,
     +  F4.2,1X,A3,1X,E8.3,1X,E7.2,1X,E8.3,1X,E7.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F4.2,1X,A3,1X,E8.3,1X,
     +  E7.2,1X,E8.3,1X,E7.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F5.2,1X,F4.2,1X,A3,1X,E8.3,1X,E7.2,1X,E8.3,1X,E7.2,
     +  1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,A3,1X,E8.3,1X,E7.2,1X,E8.3,1X,E7.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3,1X,E8.3,1X,
     +  E7.2,1X,E8.3,1X,E7.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F6.2,1X,I1,1X,E8.3,
     +  1X,E7.2,8X,A16)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,300
        read(1,'(A978)')ar__1
        read(ar__1,2)
     +  SRCID,CFlux_1,e_CFlux_1,CHR1_1,e_CHR1_1,CHR2_1,e_CHR2_1,
     +  CHR3_1,e_CHR3_1,CHR4_1,e_CHR4_1,XID_1,Class_1,Remarks_1,c1val,
     +  c1Rate,e_c1Rate,c1Flux,e_c1Flux,c1HR1,e_c1HR1,c1HR2,e_c1HR2,
     +  c1HR3,e_c1HR3,c1HR4,e_c1HR4,c2val,c2Rate,e_c2Rate,c2Flux,
     +  e_c2Flux,c2HR1,e_c2HR1,c2HR2,e_c2HR2,c2HR3,e_c2HR3,c2HR4,
     +  e_c2HR4,c3val,c3Rate,e_c3Rate,c3Flux,e_c3Flux,c3HR1,e_c3HR1,
     +  c3HR2,e_c3HR2,c3HR3,e_c3HR3,c3HR4,e_c3HR4,n1val,n1Rate,
     +  e_n1Rate,n1Flux,e_n1Flux,n1HR1,e_n1HR1,n1HR2,e_n1HR2,n1HR3,
     +  e_n1HR3,n1HR4,e_n1HR4,c4val,c4Rate,e_c4Rate,c4Flux,e_c4Flux,
     +  c4HR1,e_c4HR1,c4HR2,e_c4HR2,c4HR3,e_c4HR3,c4HR4,e_c4HR4,s1val,
     +  s1Rate,e_s1Rate,s1Flux,e_s1Flux,s1HR1,e_s1HR1,s1HR2,e_s1HR2,
     +  s1HR3,e_s1HR3,s1HR4,e_s1HR4,bval,bRate,e_bRate,bFlux,e_bFlux,
     +  bHR1,e_bHR1,bHR2,e_bHR2,bHR3,e_bHR3,bHR4,e_bHR4,ndet,svarmax,
     +  fvarmax,nupper,fmax,e_fmax,XMMM31_1
        if(ar__1(22:26) .EQ. '') CHR1_1 = rNULL__
        if(ar__1(28:33) .EQ. '') e_CHR1_1 = rNULL__
        if(ar__1(35:39) .EQ. '') CHR2_1 = rNULL__
        if(ar__1(41:46) .EQ. '') e_CHR2_1 = rNULL__
        if(ar__1(48:52) .EQ. '') CHR3_1 = rNULL__
        if(ar__1(54:57) .EQ. '') e_CHR3_1 = rNULL__
        if(ar__1(59:63) .EQ. '') CHR4_1 = rNULL__
        if(ar__1(65:68) .EQ. '') e_CHR4_1 = rNULL__
        if(ar__1(348:355) .EQ. '') c1Rate = rNULL__
        if(ar__1(357:363) .EQ. '') e_c1Rate = rNULL__
        if(ar__1(365:372) .EQ. '') c1Flux = rNULL__
        if(ar__1(374:380) .EQ. '') e_c1Flux = rNULL__
        if(ar__1(382:386) .EQ. '') c1HR1 = rNULL__
        if(ar__1(388:391) .EQ. '') e_c1HR1 = rNULL__
        if(ar__1(393:397) .EQ. '') c1HR2 = rNULL__
        if(ar__1(399:402) .EQ. '') e_c1HR2 = rNULL__
        if(ar__1(404:408) .EQ. '') c1HR3 = rNULL__
        if(ar__1(410:413) .EQ. '') e_c1HR3 = rNULL__
        if(ar__1(415:419) .EQ. '') c1HR4 = rNULL__
        if(ar__1(421:424) .EQ. '') e_c1HR4 = rNULL__
        if(ar__1(430:437) .EQ. '') c2Rate = rNULL__
        if(ar__1(439:445) .EQ. '') e_c2Rate = rNULL__
        if(ar__1(447:454) .EQ. '') c2Flux = rNULL__
        if(ar__1(456:462) .EQ. '') e_c2Flux = rNULL__
        if(ar__1(464:468) .EQ. '') c2HR1 = rNULL__
        if(ar__1(470:473) .EQ. '') e_c2HR1 = rNULL__
        if(ar__1(475:479) .EQ. '') c2HR2 = rNULL__
        if(ar__1(481:484) .EQ. '') e_c2HR2 = rNULL__
        if(ar__1(486:490) .EQ. '') c2HR3 = rNULL__
        if(ar__1(492:497) .EQ. '') e_c2HR3 = rNULL__
        if(ar__1(499:503) .EQ. '') c2HR4 = rNULL__
        if(ar__1(505:508) .EQ. '') e_c2HR4 = rNULL__
        if(ar__1(514:521) .EQ. '') c3Rate = rNULL__
        if(ar__1(523:529) .EQ. '') e_c3Rate = rNULL__
        if(ar__1(531:538) .EQ. '') c3Flux = rNULL__
        if(ar__1(540:546) .EQ. '') e_c3Flux = rNULL__
        if(ar__1(548:552) .EQ. '') c3HR1 = rNULL__
        if(ar__1(554:557) .EQ. '') e_c3HR1 = rNULL__
        if(ar__1(559:563) .EQ. '') c3HR2 = rNULL__
        if(ar__1(565:568) .EQ. '') e_c3HR2 = rNULL__
        if(ar__1(570:574) .EQ. '') c3HR3 = rNULL__
        if(ar__1(576:581) .EQ. '') e_c3HR3 = rNULL__
        if(ar__1(583:587) .EQ. '') c3HR4 = rNULL__
        if(ar__1(589:592) .EQ. '') e_c3HR4 = rNULL__
        if(ar__1(598:605) .EQ. '') n1Rate = rNULL__
        if(ar__1(607:613) .EQ. '') e_n1Rate = rNULL__
        if(ar__1(615:622) .EQ. '') n1Flux = rNULL__
        if(ar__1(624:630) .EQ. '') e_n1Flux = rNULL__
        if(ar__1(632:636) .EQ. '') n1HR1 = rNULL__
        if(ar__1(638:641) .EQ. '') e_n1HR1 = rNULL__
        if(ar__1(643:647) .EQ. '') n1HR2 = rNULL__
        if(ar__1(649:652) .EQ. '') e_n1HR2 = rNULL__
        if(ar__1(654:658) .EQ. '') n1HR3 = rNULL__
        if(ar__1(660:663) .EQ. '') e_n1HR3 = rNULL__
        if(ar__1(665:669) .EQ. '') n1HR4 = rNULL__
        if(ar__1(671:674) .EQ. '') e_n1HR4 = rNULL__
        if(ar__1(680:687) .EQ. '') c4Rate = rNULL__
        if(ar__1(689:695) .EQ. '') e_c4Rate = rNULL__
        if(ar__1(697:704) .EQ. '') c4Flux = rNULL__
        if(ar__1(706:712) .EQ. '') e_c4Flux = rNULL__
        if(ar__1(714:718) .EQ. '') c4HR1 = rNULL__
        if(ar__1(720:723) .EQ. '') e_c4HR1 = rNULL__
        if(ar__1(725:729) .EQ. '') c4HR2 = rNULL__
        if(ar__1(731:734) .EQ. '') e_c4HR2 = rNULL__
        if(ar__1(736:740) .EQ. '') c4HR3 = rNULL__
        if(ar__1(742:745) .EQ. '') e_c4HR3 = rNULL__
        if(ar__1(747:751) .EQ. '') c4HR4 = rNULL__
        if(ar__1(753:756) .EQ. '') e_c4HR4 = rNULL__
        if(ar__1(762:769) .EQ. '') s1Rate = rNULL__
        if(ar__1(771:777) .EQ. '') e_s1Rate = rNULL__
        if(ar__1(779:786) .EQ. '') s1Flux = rNULL__
        if(ar__1(788:794) .EQ. '') e_s1Flux = rNULL__
        if(ar__1(796:800) .EQ. '') s1HR1 = rNULL__
        if(ar__1(802:805) .EQ. '') e_s1HR1 = rNULL__
        if(ar__1(807:811) .EQ. '') s1HR2 = rNULL__
        if(ar__1(813:816) .EQ. '') e_s1HR2 = rNULL__
        if(ar__1(818:822) .EQ. '') s1HR3 = rNULL__
        if(ar__1(824:827) .EQ. '') e_s1HR3 = rNULL__
        if(ar__1(829:833) .EQ. '') s1HR4 = rNULL__
        if(ar__1(835:838) .EQ. '') e_s1HR4 = rNULL__
        if(ar__1(844:851) .EQ. '') bRate = rNULL__
        if(ar__1(853:859) .EQ. '') e_bRate = rNULL__
        if(ar__1(861:868) .EQ. '') bFlux = rNULL__
        if(ar__1(870:876) .EQ. '') e_bFlux = rNULL__
        if(ar__1(878:882) .EQ. '') bHR1 = rNULL__
        if(ar__1(884:887) .EQ. '') e_bHR1 = rNULL__
        if(ar__1(889:893) .EQ. '') bHR2 = rNULL__
        if(ar__1(895:898) .EQ. '') e_bHR2 = rNULL__
        if(ar__1(900:904) .EQ. '') bHR3 = rNULL__
        if(ar__1(906:909) .EQ. '') e_bHR3 = rNULL__
        if(ar__1(911:915) .EQ. '') bHR4 = rNULL__
        if(ar__1(917:920) .EQ. '') e_bHR4 = rNULL__
        if(ar__1(939:946) .EQ. '') fmax = rNULL__
        if(ar__1(948:954) .EQ. '') e_fmax = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  SRCID,CFlux_1,e_CFlux_1,CHR1_1,e_CHR1_1,CHR2_1,e_CHR2_1,
     +  CHR3_1,e_CHR3_1,CHR4_1,e_CHR4_1,XID_1,Class_1,Remarks_1,c1val,
     +  c1Rate,e_c1Rate,c1Flux,e_c1Flux,c1HR1,e_c1HR1,c1HR2,e_c1HR2,
     +  c1HR3,e_c1HR3,c1HR4,e_c1HR4,c2val,c2Rate,e_c2Rate,c2Flux,
     +  e_c2Flux,c2HR1,e_c2HR1,c2HR2,e_c2HR2,c2HR3,e_c2HR3,c2HR4,
     +  e_c2HR4,c3val,c3Rate,e_c3Rate,c3Flux,e_c3Flux,c3HR1,e_c3HR1,
     +  c3HR2,e_c3HR2,c3HR3,e_c3HR3,c3HR4,e_c3HR4,n1val,n1Rate,
     +  e_n1Rate,n1Flux,e_n1Flux,n1HR1,e_n1HR1,n1HR2,e_n1HR2,n1HR3,
     +  e_n1HR3,n1HR4,e_n1HR4,c4val,c4Rate,e_c4Rate,c4Flux,e_c4Flux,
     +  c4HR1,e_c4HR1,c4HR2,e_c4HR2,c4HR3,e_c4HR3,c4HR4,e_c4HR4,s1val,
     +  s1Rate,e_s1Rate,s1Flux,e_s1Flux,s1HR1,e_s1HR1,s1HR2,e_s1HR2,
     +  s1HR3,e_s1HR3,s1HR4,e_s1HR4,bval,bRate,e_bRate,bFlux,e_bFlux,
     +  bHR1,e_bHR1,bHR2,e_bHR2,bHR3,e_bHR3,bHR4,e_bHR4,ndet,svarmax,
     +  fvarmax,nupper,fmax,e_fmax,XMMM31_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(A8,1X,A19,1X,A27,1X,A81)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,52
        read(1,'(A138)')ar__2
        read(ar__2,3)Ref,BibCode,Name,Com
c    ..............Just test output...........
        write(6,3)Ref,BibCode,Name,Com
c    .......End.of.Just test output...........
      end do
      close(1)

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