FORTRAN Generation
(/./ftp/cats/J/MNRAS/416/1844)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/416/1844 into FORTRAN code for loading all data files into arrays.

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-20
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

      implicit none
*  Unspecified or NULL values, generally corresponding to blank columns,
*  are assigned one of the following special values:
*     rNULL__    for unknown or NULL floating-point values
*     iNULL__    for unknown or NULL   integer      values
      real*4     rNULL__
      integer*4  iNULL__
      parameter  (rNULL__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. J/MNRAS/416/1844    2XMM ultraluminous X-ray source candidates  (Walton+, 2011)
*================================================================================
*2XMM ultraluminous X-ray source candidates in nearby galaxies.
*    Walton D.J., Roberts T.P., Mateos S., Heard V.
*   <Mon. Not. R. Astron. Soc., 416, 1844-1861 (2011)>
*   =2011MNRAS.416.1844W    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! catalogue of ultraluminous X-ray source candidates

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

C  J2000 position composed of: RAdeg DEdeg
      character*12  name       (nr__) ! RC3 Name
      character*14  altname    (nr__) ! RC3 Alternate name
      character*14  desig      (nr__) ! RC3 Other designation
      character*9   PGC        (nr__) ! RC3 PGC number
      character*7   type       (nr__) ! RC3 Type
      character*3   typesr     (nr__) ! RC3 Source of revised type
      real*4        T          (nr__) ! ? RC3 T (Hubble stage)
      real*4        e_T        (nr__) ! ? RC3 Mean error of T
      real*4        lumcl      (nr__) ! ? RC3 Luminosity class
      real*4        e_lumcl    (nr__) ! ? RC3 Mean error on luminosity class
      integer*4     o_lumcl    (nr__) ! ? RC3 Num. of luminosity class
*                                          estimates
      real*4        D25        (nr__) ! ([0.1arcmin]) ? RC3 Log D25
      character*1   u_D25      (nr__) ! [?] RC3 indicates Log D25 very uncer.
      real*4        e_D25      (nr__) ! ([0.1arcmin]) ? RC3 Mean error on log D25
      real*4        R25        (nr__) ! ? RC3 Log R25 (isophotal diameter
*                                          ratio)
      character*1   u_R25      (nr__) ! [?] RC3 indicates Log R25 very
*                                          uncertain
      real*4        e_R25      (nr__) ! ? RC3 Mean error on log R25
      real*4        Do         (nr__) ! ([0.1arcmin]) ? RC3 Log Do
      real*4        Ae         (nr__) ! ([0.1arcmin]) ? RC3 Log Ae (effective aperture)
      real*4        e_Ae       (nr__) ! ([0.1arcmin]) ? RC3 Mean error on log Ae
      integer*4     PA         (nr__) ! (deg) ? RC3 Position angle of the major axis
      real*4        BT         (nr__) ! (mag) ? RC3 BT (total B magnitude)
      character*1   BT_code    (nr__) ! ? RC3 BT code
      real*4        e_BT       (nr__) ! (mag) ? RC3 Mean error on BT
      real*4        Bmag       (nr__) ! (mag) ? RC3 mB (photographic magnitude)
      real*4        e_Bmag     (nr__) ! (mag) ? RC3 Mean error on Bmag
      real*4        BoT        (nr__) ! (mag) ? RC3 BoT
      real*4        m_25       (nr__) ! (mag/arcmin2) ? RC3 m'25
      real*4        e_m_25     (nr__) ! (mag/arcmin2) ? RC3 Mean error m'25
      real*4        m_e        (nr__) ! (mag/arcmin2) ? RC3 m'e
      real*4        e_m_e      (nr__) ! (mag/arcmin2) ? RC3 Mean error m'e
      real*4        mFIR       (nr__) ! (mag) ? RC3 mFIR
      real*4        m21        (nr__) ! (mag) ? RC3 m21
      real*4        e_m21      (nr__) ! (mag) ? RC3 Mean error on m21
      real*4        B_VT       (nr__) ! (mag) ? RC3 (B-V)T (total (B-V))
      real*4        e_B_VT     (nr__) ! (mag) ? RC3 Mean error on (B-V)T
      real*4        B_Ve       (nr__) ! (mag) ? RC3 (B-V)e (mean B-V within Ae)
      real*4        e_B_Ve     (nr__) ! (mag) ? RC3 Mean error on (B-V)e
      real*4        B_VoT      (nr__) ! (mag) ? RC3 (B-V)oT
      real*4        U_BT       (nr__) ! (mag) ? RC3 (U-B)T (total (U-B))
      real*4        e_U_BT     (nr__) ! (mag) ? RC3 Mean error on (U-B)T
      real*4        U_Be       (nr__) ! (mag) ? RC3 (U-B)e
      real*4        e_U_Be     (nr__) ! (mag) ? RC3 Mean error on (U-B)e
      real*4        U_BoT      (nr__) ! (mag) ? RC3 (U-B)To
      real*4        HI         (nr__) ! (mag) ? RC3 HI (neutral hydrogen index)
      real*4        Ai         (nr__) ! (mag) ? RC3 Ai (internal B extinction)
      real*4        A21        (nr__) ! (mag) ? RC3 A21 (HI self absorption)
      real*4        Ag         (nr__) ! (mag) ? Ag (galactic extinction in B)
      integer*4     W20        (nr__) ! (km/s) ? RC3 W20
      integer*4     e_W20      (nr__) ! (km/s) ? RC3 Mean error on W20
      integer*4     W50        (nr__) ! (km/s) ? RC3 W50
      integer*4     e_W50      (nr__) ! (km/s) ? RC3 Mean error of W50
      integer*4     V21        (nr__) ! (km/s) ? RC3 V21 (heliocentric radial
*                                          velocity)
      integer*4     e_V21      (nr__) ! (km/s) ? RC3 Mean error on V21
      integer*4     cz         (nr__) ! (km/s) ? RC3 Vopt= cz
      integer*4     e_cz       (nr__) ! (km/s) ? RC3 Mean error on Vopt
      integer*4     VGSR       (nr__) ! (km/s) ? RC3 VGSR
      integer*4     V3K        (nr__) ! (km/s) ? RC3 V3K
      real*8        RAdeg      (nr__) ! (deg) NED right ascension (J2000) (Gal_RA)
      real*8        DEdeg      (nr__) ! (deg) NED declination (J2000) (Gal_DEC)
      real*8        ePos       (nr__) ! (arcsec) NWD position uncertainty (Gal_PosErr)
      real*4        Dist       (nr__) ! (Mpc) Galaxy distance, calculated assuming
*                                          the Hubble flow for galaxies with
*                                          cz>1000km/s, or as in Tully (1998,
*                                          Cat. VII/145/) otherwise (Distance)
      integer*4     DETID      (nr__) ! ?=- 2XMM Detection ID of the
*                                          identified target (DETID)
      integer*4     SRCID      (nr__) ! 2XMM Source ID (SRCID)
      character*21  IAUName    (nr__) ! 2XMM IAU name (IAUName)
      integer*4     SRC        (nr__) ! ?=- 2XMM Source number per observation
*                                          of the identified target (SRC)
      character*16  v1XMM      (nr__) ! 1XMM match, JHHMMSS.s+DDMMSS 
*                                          (MATCH_1XMM)
      real*8        Sep1XMM    (nr__) ! (arcsec) ?=- Separation to 1XMM match (SEP_1XMM)
      integer*4     v2XMMpID   (nr__) ! ?=- 2XMMp source ID (SRCID_2XMMP)
      character*16  v2XMMp     (nr__) ! 2XMMp match (MATCH_2XMMP)
      real*8        Sep2XMMP   (nr__) ! (arcsec) ?=- Separation to 2XMMP source
*                                          (SEP_2XMMP)
*----------------------------- 2XMM data ----------------------------------------
      integer*8     ObsID      (nr__) ! 2XMM Observation number (ObsID)
      integer*4     Rev        (nr__) ! 2XMM Satellite revolution number
*                                          (consecutive in time) (REVOLUT)
      real*8        MJD0       (nr__) ! (d) 2XMM Date of first observation
*                                          (MJD_START)
      real*8        MJD1       (nr__) ! (d) 2XMM Date of last observation
*                                          (MJD_STOP)
      integer*4     Q          (nr__) ! 2XMM  Observations quality (OBS_CLASS)
      character*6   f_PN       (nr__) ! 2XMM (PN_FILT) Filter of the PN
*                                          exposure(s) (PN_FILTER)
      character*6   f_M1       (nr__) ! 2XMM (M1_FILT) Filter of the M1
*                                          exposure(s) (M1_FILTER)
      character*6   f_M2       (nr__) ! 2XMM (M2_FILT) Filter of the M2
*                                          exposure(s) (M2_FILTER)
      character*23  m_PN       (nr__) ! 2XMM (PN_MODE) Mode of the PN
*                                          exposure(s) (PN_SUBMODE)
      character*15  m_M1       (nr__) ! 2XMM (M1_MODE) Mode of the M1
*                                          exposures (M1_SUBMODE)
      character*15  m_M2       (nr__) ! 2XMM (M2_MODE) Mode of the M2
*                                          exposures (M2_SUBMODE)
      real*8        RA1deg     (nr__) ! (deg) 2XMM Detection right ascension (J2000)
*                                          (Det_RA)
      real*8        DE1deg     (nr__) ! (deg) 2XMM Detection declination (J2000)
*                                          (Det_DEC)
      real*8        ePos1      (nr__) ! (arcsec) 2XMM Total radial detection position
*                                          uncertainty (POSERR)
      real*8        GLON       (nr__) ! (deg) 2XMM XMM1 Galactic longitude (LII)
      real*8        GLAT       (nr__) ! (deg) 2XMM XMM1 Galactic latitude (BII)
      real*8        ePos2      (nr__) ! (arcsec) 2XMM Mean error on position (RADEC_ERR)
      real*4        SYSERR     (nr__) ! (arcsec) 2XMM estimated error arising from the
*                                          field rectification process on
*                                          position (SYSERR)
      real*8        RAudeg     (nr__) ! (deg) 2XMM Right Ascension as determined by
*                                          the SAS task (J2000) (RA_UNC)
      real*8        DEudeg     (nr__) ! (deg) 2XMM DEclination as determined by the
*                                          SAS task (J2000)  (DECA_UNC)
      real*4        FluxEP1    (nr__) ! (mW/m2) 2XMM Mean flux in 0.2-0.5keV band
*                                          (EP_1_FLUX)
      real*4        e_FluxEP1  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP1
*                                          (EP_1_FLUX_ERR)
      real*4        FluxEP2    (nr__) ! (mW/m2) 2XMM Mean flux in 0.5-1.0keV band
*                                          (EP_2_FLUX)
      real*4        e_FluxEP2  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP2
*                                          (EP_2_FLUX_ERR)
      real*4        FluxEP3    (nr__) ! (mW/m2) 2XMM Mean flux in 1.0-2.0keV band
*                                          (EP_3_FLUX)
      real*4        e_FluxEP3  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP3
*                                          (EP_3_FLUX_ERR)
      real*4        FluxEP4    (nr__) ! (mW/m2) 2XMM Mean flux in 2.0-4.5keV band
*                                          (EP_4_FLUX)
      real*4        e_FluxEP4  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP4
*                                          (EP_4_FLUX_ERR)
      real*4        FluxEP5    (nr__) ! (mW/m2) 2XMM Mean flux in 4.5-12keV band
*                                          (EP_5_FLUX)
      real*4        e_FluxEP5  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP5
*                                          (EP_5_FLUX_ERR)
      real*4        FluxEP8    (nr__) ! (mW/m2) 2XMM Mean flux in 0.2-12keV band
*                                          (EP_8_FLUX)
      real*4        e_FluxEP8  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP8
*                                          (EP_8_FLUX_ERR)
      real*4        FluxEP9    (nr__) ! (mW/m2) 2XMM Mean flux in 0.5-4.5keV band
*                                          (EP_9_FLUX)
      real*4        e_FluxEP9  (nr__) ! (mW/m2) 2XMM Mean error on FluxEP9
*                                          (EP_9_FLUX_ERR)
      real*4        FluxPN1    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-0.5keV band
*                                          (SC_PN_1_FLUX)
      real*4        e_FluxPN1  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN1
*                                          (SC_PN_1_FLUX_ERR)
      real*4        FluxPN2    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-1.0keV band
*                                          (SC_PN_2_FLUX)
      real*4        e_FluxPN2  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN2
*                                          (SC_PN_2_FLUX_ERR)
      real*4        FluxPN3    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 1.0-2.0keV band
*                                          (SC_PN_3_FLUX)
      real*4        e_FluxPN3  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN3
*                                          (SC_PN_3_FLUX_ERR)
      real*4        FluxPN4    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 2.0-4.5keV band
*                                          (SC_PN_4_FLUX)
      real*4        e_FluxPN4  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN4
*                                          (SC_PN_4_FLUX_ERR)
      real*4        FluxPN5    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 4.5-12keV band
*                                          (SC_PN_5_FLUX)
      real*4        e_FluxPN5  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN5
*                                          (SC_PN_5_FLUX_ERR)
      real*4        FluxPN8    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-12keV band
*                                          (SC_PN_8_FLUX)
      real*4        e_FluxPN8  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN8
*                                          (SC_PN_8_FLUX_ERR)
      real*4        FluxPN9    (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-4.5keV band
*                                          (SC_PN_9_FLUX)
      real*4        e_FluxPN9  (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxPN9
*                                          (SC_PN_9_FLUX_ERR)
      real*4        FluxM1_1   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-0.5keV band
*                                          (SC_M1_1_FLUX)
      real*4        e_FluxM1_1 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-1
*                                          (SC_M1_1_FLUX_ERR)
      real*4        FluxM1_2   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-1.0keV band
*                                          (SC_M1_2_FLUX)
      real*4        e_FluxM1_2 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-2
*                                          (SC_M1_2_FLUX_ERR)
      real*4        FluxM1_3   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 1.0-2.0keV band
*                                          (SC_M1_3_FLUX)
      real*4        e_FluxM1_3 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-3
*                                          (SC_M1_3_FLUX_ERR)
      real*4        FluxM1_4   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 2.0-4.5keV band
*                                          (SC_M1_4_FLUX)
      real*4        e_FluxM1_4 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-4
*                                          (SC_M1_4_FLUX_ERR)
      real*4        FluxM1_5   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 4.5-12keV band
*                                          (SC_M1_5_FLUX)
      real*4        e_FluxM1_5 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-5
*                                          (SC_M1_5_FLUX_ERR)
      real*4        FluxM1_8   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-12keV band
*                                          (SC_M1_8_FLUX)
      real*4        e_FluxM1_8 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-8
*                                          (SC_M1_8_FLUX_ERR)
      real*4        FluxM1_9   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-4.5keV band
*                                          (SC_M1_9_FLUX)
      real*4        e_FluxM1_9 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM1-9
*                                          (SC_M1_9_FLUX_ERR)
      real*4        FluxM2_1   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-0.5keV band
*                                          (SC_M2_1_FLUX)
      real*4        e_FluxM2_1 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-1
*                                          (SC_M2_1_FLUX_ERR)
      real*4        FluxM2_2   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-1.0keV band
*                                          (SC_M2_2_FLUX)
      real*4        e_FluxM2_2 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-2
*                                          (SC_M2_2_FLUX_ERR)
      real*4        FluxM2_3   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 1.0-2.0keV band
*                                          (SC_M2_3_FLUX)
      real*4        e_FluxM2_3 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-3
*                                          (SC_M2_3_FLUX_ERR)
      real*4        FluxM2_4   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 2.0-4.5keV band
*                                          (SC_M2_4_FLUX)
      real*4        e_FluxM2_4 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-4
*                                          (SC_M2_4_FLUX_ERR)
      real*4        FluxM2_5   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 4.5-12keV band
*                                          (SC_M2_5_FLUX)
      real*4        e_FluxM2_5 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-5
*                                          (SC_M1_5_FLUX_ERR)
      real*4        FluxM2_8   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.2-12keV band
*                                          (SC_M2_8_FLUX)
      real*4        e_FluxM2_8 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-8
*                                          (SC_M2_8_FLUX_ERR)
      real*4        FluxM2_9   (nr__) ! (mW/m2) ?=- 2XMM Mean flux in 0.5-4.5keV band
*                                          (SC_M2_9_FLUX)
      real*4        e_FluxM2_9 (nr__) ! (mW/m2) ?=- 2XMM Mean error on FluxM2-9
*                                          (SC_M2_9_FLUX_ERR)
      real*8        EP8rate    (nr__) ! (ct/s) 2XMM EPIC 8 count rate (EP_8_RATE)
      real*8        e_EP8rate  (nr__) ! (ct/s) 2XMM rms uncertainty on EP8rate
*                                          (EP_8_RATE_ERR)
      real*8        EP9rate    (nr__) ! (ct/s) 2XMM EPIC 9 count rate (EP_9_RATE)
      real*8        e_EP9rate  (nr__) ! (ct/s) 2XMM rms uncertainty on EP9rate
*                                          (EP_9_RATE_ERR)
      real*8        PN1rate    (nr__) ! (ct/s) 2XMM PN1 count rate (PN_1_RATE)
      real*8        e_PN1rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN1rate
*                                          (PN_1_RATE_ERR)
      real*8        PN2rate    (nr__) ! (ct/s) 2XMM PN2 count rate (PN_2_RATE)
      real*8        e_PN2rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN2rate
*                                          (PN_2_RATE_ERR)
      real*8        PN3rate    (nr__) ! (ct/s) 2XMM PN3 count rate (PN_3_RATE)
      real*8        e_PN3rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN3rate
*                                          (PN_3_RATE_ERR)
      real*8        PN4rate    (nr__) ! (ct/s) 2XMM PN4 count rate (PN_4_RATE)
      real*8        e_PN4rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN4rate
*                                          (PN_4_RATE_ERR)
      real*8        PN5rate    (nr__) ! (ct/s) 2XMM PN5 count rate (PN_5_RATE)
      real*8        e_PN5rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN5rate
*                                          (PN_5_RATE_ERR)
      real*8        PN8rate    (nr__) ! (ct/s) 2XMM PN8 count rate (PN_8_RATE)
      real*8        e_PN8rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN8rate
*                                          (PN_8_RATE_ERR)
      real*8        PN9rate    (nr__) ! (ct/s) 2XMM PN9 count rate (PN_9_RATE)
      real*8        e_PN9rate  (nr__) ! (ct/s) 2XMM rms uncertainty on PN9rate
*                                          (PN_9_RATE_ERR)
      real*8        M1_1rate   (nr__) ! (ct/s) 2XMM M1-1 count rate (M1_1_RATE)
      real*8        e_M1_1rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-1rate
*                                          (M1_1_RATE_ERR)
      real*8        M1_2rate   (nr__) ! (ct/s) 2XMM M1-2 count rate (M1_2_RATE)
      real*8        e_M1_2rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-2rate
*                                          (M1_2_RATE_ERR)
      real*8        M1_3rate   (nr__) ! (ct/s) 2XMM M1-3 count rate (M1_3_RATE)
      real*8        e_M1_3rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-3rate
*                                          (M1_3_RATE_ERR)
      real*8        M1_4rate   (nr__) ! (ct/s) 2XMM M1-4 count rate (M1_4_RATE)
      real*8        e_M1_4rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-4rate
*                                          (M1_4_RATE_ERR)
      real*8        M1_5rate   (nr__) ! (ct/s) 2XMM M1-5 count rate (M1_5_RATE)
      real*8        e_M1_5rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-5rate
*                                          (M1_5_RATE_ERR)
      real*8        M1_8rate   (nr__) ! (ct/s) 2XMM M1-8 count rate (M1_8_RATE)
      real*8        e_M1_8rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-8rate
*                                          (M1_8_RATE_ERR)
      real*8        M1_9rate   (nr__) ! (ct/s) 2XMM M1-9 count rate (M1_9_RATE)
      real*8        e_M1_9rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-9rate
*                                          (M1_9_RATE_ERR)
      real*8        M2_1rate   (nr__) ! (ct/s) 2XMM M2-1 count rate (M2_1_RATE)
      real*8        e_M2_1rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-1rate
*                                          (M2_1_RATE_ERR)
      real*8        M2_2rate   (nr__) ! (ct/s) 2XMM M2-2 count rate (M2_2_RATE)
      real*8        e_M2_2rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-2rate
*                                          (M2_2_RATE_ERR)
      real*8        M2_3rate   (nr__) ! (ct/s) 2XMM M2-3 count rate (M2_3_RATE)
      real*8        e_M2_3rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-3rate
*                                          (M2_3_RATE_ERR)
      real*8        M2_4rate   (nr__) ! (ct/s) 2XMM M2-4 count rate (M2_4_RATE)
      real*8        e_M2_4rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-4rate
*                                          (M2_4_RATE_ERR)
      real*8        M2_5rate   (nr__) ! (ct/s) 2XMM M2-5 count rate (M2_5_RATE)
      real*8        e_M2_5rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-5rate
*                                          (M2_5_RATE_ERR)
      real*8        M2_8rate   (nr__) ! (ct/s) 2XMM M2-8 count rate (M2_8_RATE)
      real*8        e_M2_8rate (nr__) ! (ct/s) 2XMM rms uncertainty on M2-8rate
*                                          (M2_8_RATE_ERR)
      real*8        M2_9rate   (nr__) ! (ct/s) 2XMM M2-9 count rate (M2_9_RATE)
      real*8        e_M2_9rate (nr__) ! (ct/s) 2XMM rms uncertainty on M1-9rate
*                                          (M2_9_RATE_ERR)
      real*8        EP8cts     (nr__) ! (ct) 2XMM EP 8 counts (EP_8_CTS)
      real*8        e_EP8cts   (nr__) ! (ct) 2XMM rms uncertainty on EP8cts
*                                          (EP_8_CTS_ERR)
      real*8        PN8cts     (nr__) ! (ct) ?=- 2XMM PN 8 counts (PN_8_CTS)
      real*8        e_PN8cts   (nr__) ! (ct) ?=- 2XMM rms uncertainty on PN8cts
*                                          (PN_8_CTS_ERR)
      real*8        M1_8cts    (nr__) ! (ct) ?=- 2XMM M1-8 counts (M1_8_CTS)
      real*8        e_M1_8cts  (nr__) ! (ct) ?=- 2XMM rms uncertainty on M1-8cts
*                                          (M1_8_CTS_ERR)
      real*8        M2_8cts    (nr__) ! (ct) ?=- 2XMM M2-8 counts (M2_8_CTS)
      real*8        e_M2_8cts  (nr__) ! (ct) ?=- 2XMM rms uncertainty on M2-8cts
*                                          (M2_8_CTS_ERR)
      real*4        EP8detML   (nr__) ! 2XMM EP8 detection likelihood
*                                          (EP_8_DET_ML)
      real*4        EP9detML   (nr__) ! 2XMM EP9 detection likelihood
*                                          (EP_9_DET_ML)
      real*4        PN1detML   (nr__) ! ?=- 2XMM PN1 detection likelihood
*                                          (PN_1_DET_ML)
      real*4        PN2detML   (nr__) ! ?=- 2XMM PN2 detection likelihood
*                                          (PN_2_DET_ML)
      real*4        PN3detML   (nr__) ! ?=- 2XMM PN2 detection likelihood
*                                          (PN_3_DET_ML)
      real*4        PN4detML   (nr__) ! ?=- 2XMM PN4 detection likelihood
*                                          (PN_4_DET_ML)
      real*4        PN5detML   (nr__) ! ?=- 2XMM PN5 detection likelihood
*                                          (PN_5_DET_ML)
      real*4        PN8detML   (nr__) ! ?=- 2XMM PN8 detection likelihood
*                                          (PN_8_DET_ML)
      real*4        PN9detML   (nr__) ! ?=- 2XMM PN9 detection likelihood
*                                          (PN_9_DET_ML)
      real*4        M1_1detML  (nr__) ! ?=- 2XMM M1-1 detection likelihood
*                                          (M1_1_DET_ML)
      real*4        M1_2detML  (nr__) ! ?=- 2XMM M1-2 detection likelihood
*                                          (M1_2_DET_ML)
      real*4        M1_3detML  (nr__) ! ?=- 2XMM M1-2 detection likelihood
*                                          (M1_3_DET_ML)
      real*4        M1_4detML  (nr__) ! ?=- 2XMM M1-4 detection likelihood
*                                          (M1_4_DET_ML)
      real*4        M1_5detML  (nr__) ! ?=- 2XMM M1-5 detection likelihood
*                                          (M1_5_DET_ML)
      real*4        M1_8detML  (nr__) ! ?=- 2XMM M1-8 detection likelihood
*                                          (M1_8_DET_ML)
      real*4        M1_9detML  (nr__) ! ?=- 2XMM M1-9 detection likelihood
*                                          (M1_9_DET_ML)
      real*4        M2_1detML  (nr__) ! ?=- 2XMM M2-1 detection likelihood
*                                          (M2_1_DET_ML)
      real*4        M2_2detML  (nr__) ! ?=- 2XMM M2-2 detection likelihood
*                                          (M2_2_DET_ML)
      real*4        M2_3detML  (nr__) ! ?=- 2XMM M2-2 detection likelihood
*                                          (M2_3_DET_ML)
      real*4        M2_4detML  (nr__) ! ?=- 2XMM M2-4 detection likelihood
*                                          (M2_4_DET_ML)
      real*4        M2_5detML  (nr__) ! ?=- 2XMM M2-5 detection likelihood
*                                          (M2_5_DET_ML)
      real*4        M2_8detML  (nr__) ! ?=- 2XMM M2-8 detection likelihood
*                                          (M2_8_DET_ML)
      real*4        M2_9detML  (nr__) ! ?=- 2XMM M2-9 detection likelihood
*                                          (M2_9_DET_ML)
      real*8        HREP1      (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 1,2 (EP_HR1)
      real*8        e_HREP1    (nr__) ! ?=- 2XMM Mean error on HREP1
*                                          (EP_HR1_ERR)
      real*8        HREP2      (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 2,3 (EP_HR2)
      real*8        e_HREP2    (nr__) ! ?=- 2XMM Mean error on HREP2
*                                          (EP_HR2_ERR)
      real*8        HREP3      (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 3,4 (EP_HR3)
      real*8        e_HREP3    (nr__) ! ?=- 2XMM Mean error on HREP3
*                                          (EP_HR3_ERR)
      real*8        HREP4      (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 4,5 (EP_HR4)
      real*8        e_HREP4    (nr__) ! ?=- 2XMM Mean error on HREP4
*                                          (EP_HR4_ERR)
      real*8        HRPN1      (nr__) ! [-1/1]?=- 2XMM PN Hardness ratio,
*                                          bands 1,2 (PN_HR1)
      real*8        e_HRPN1    (nr__) ! ?=- 2XMM Mean error on HRPN1
*                                          (PN_HR1_ERR)
      real*8        HRPN2      (nr__) ! [-1/1]?=- 2XMM PN Hardness ratio,
*                                          bands 2,3 (PN_HR2)
      real*8        e_HRPN2    (nr__) ! ?=- 2XMM Mean error on HRPN2
*                                          (PN_HR2_ERR)
      real*8        HRPN3      (nr__) ! [-1/1]?=- 2XMM PN Hardness ratioPN,
*                                          bands 3,4 (PN_HR3)
      real*8        e_HRPN3    (nr__) ! ?=- 2XMM Mean error on HRPN3
*                                          (PN_HR3_ERR)
      real*8        HRPN4      (nr__) ! [-1/1]?=- 2XMM PN Hardness ratio,
*                                          bands 4,5 (PN_HR4)
      real*8        e_HRPN4    (nr__) ! ?=- 2XMM Mean error on HRPN4
*                                          (PN_HR4_ERR)
      real*8        HRM1_1     (nr__) ! [-1/1]?=- 2XMM M1 Hardness ratio,
*                                          bands 1,2 (M1_HR1)
      real*8        e_HRM1_1   (nr__) ! ?=- 2XMM Mean error on HRM1-1
*                                          (M1_HR1_ERR)
      real*8        HRM1_2     (nr__) ! [-1/1]?=- 2XMM M1  Hardness ratio,
*                                          bands 2,3 (M1-_HR2)
      real*8        e_HRM1_2   (nr__) ! ?=- 2XMM Mean error on HRM1-2
*                                          (M1_HR2_ERR)
      real*8        HRM1_3     (nr__) ! [-1/1]?=- 2XMM M1 Hardness ratio,
*                                          bands 3,4 (M1-_HR3)
      real*8        e_HRM1_3   (nr__) ! ?=- 2XMM Mean error on HRM1-3
*                                          (M1_HR3_ERR)
      real*8        HRM1_4     (nr__) ! [-1/1]?=- 2XMM M1  Hardness ratio,
*                                          bands 4,5 (M1-_HR4)
      real*8        e_HRM1_4   (nr__) ! ?=- 2XMM Mean error on HRM1-4
*                                          (M1_HR4_ERR)
      real*8        HRM2_1     (nr__) ! [-1/1]?=- 2XMM M1 Hardness ratio,
*                                          bands 1,2 (M2_HR1)
      real*8        e_HRM2_1   (nr__) ! ?=- 2XMM Mean error on HRM2-1
*                                          (M2_HR1_ERR)
      real*8        HRM2_2     (nr__) ! [-1/1]?=- 2XMM M2 Hardness ratio,
*                                          bands 2,3 (M2-_HR2)
      real*8        e_HRM2_2   (nr__) ! ?=- 2XMM Mean error on HRM2-2
*                                          (M2_HR2_ERR)
      real*8        HRM2_3     (nr__) ! [-1/1]?=- 2XMM M2 Hardness ratio,
*                                          bands 3,4 (M2-_HR3)
      real*8        e_HRM2_3   (nr__) ! ?=- 2XMM Mean error on HRM2-3
*                                          (M2_HR3_ERR)
      real*8        HRM2_4     (nr__) ! [-1/1]?=- 2XMM M2 Hardness ratio,
*                                          bands 4,5 (M2-_HR4)
      real*8        e_HRM2_4   (nr__) ! ?=- 2XMM Mean error on HRM2-4
*                                          (M2_HR4_ERR)
      real*4        PN1ExpT    (nr__) ! (s) ?=- 2XMM PN1 exposure time (PN_1_EXP)
      real*4        PN2ExpT    (nr__) ! (s) ?=- 2XMM PN2 exposure time (PN_2_EXP)
      real*4        PN3ExpT    (nr__) ! (s) ?=- 2XMM PN3 exposure time (PN_3_EXP)
      real*4        PN4ExpT    (nr__) ! (s) ?=- 2XMM PN4 exposure time (PN_4_EXP)
      real*4        PN5ExpT    (nr__) ! (s) ?=- 2XMM PN5 exposure time (PN_5_EXP)
      real*4        M1_1ExpT   (nr__) ! (s) ?=- 2XMM M1-1 exposure time (M1_1_EXP)
      real*4        M1_2ExpT   (nr__) ! (s) ?=- 2XMM M1-2 exposure time (M1_2_EXP)
      real*4        M1_3ExpT   (nr__) ! (s) ?=- 2XMM M1-3 exposure time (M1_3_EXP)
      real*4        M1_4ExpT   (nr__) ! (s) ?=- 2XMM M1-4 exposure time (M1_4_EXP)
      real*4        M1_5ExpT   (nr__) ! (s) ?=- 2XMM M1-5 exposure time (M1_5_EXP)
      real*4        M2_1ExpT   (nr__) ! (s) ?=- 2XMM M2-1 exposure time (M2_1_EXP)
      real*4        M2_2ExpT   (nr__) ! (s) ?=- 2XMM M2-2 exposure time (M2_2_EXP)
      real*4        M2_3ExpT   (nr__) ! (s) ?=- 2XMM M2-3 exposure time (M2_3_EXP)
      real*4        M2_4ExpT   (nr__) ! (s) ?=- 2XMM M2-4 exposure time (M2_4_EXP)
      real*4        M2_5ExpT   (nr__) ! (s) ?=- 2XMM M2-5 exposure time (M2_5_EXP)
      real*8        PN1BG      (nr__) ! (ct/pix) ?=- 2XMM PN1 background map value
*                                          (PN_1_BG)
      real*8        PN2BG      (nr__) ! (ct/pix) ?=- 2XMM PN2 background map value
*                                          (PN_2_BG)
      real*8        PN3BG      (nr__) ! (ct/pix) ?=- 2XMM PN3 background map value
*                                          (PN_3_BG)
      real*8        PN4BG      (nr__) ! (ct/pix) ?=- 2XMM PN4 background map value
*                                          (PN_4_BG)
      real*8        PN5BG      (nr__) ! (ct/pix) ?=- 2XMM PN5 background map value
*                                          (PN_5_BG)
      real*8        M1_1BG     (nr__) ! (ct/pix) ?=- 2XMM M1-1 background map value
*                                          (M1_1_BG)
      real*8        M1_2BG     (nr__) ! (ct/pix) ?=- 2XMM M1-2 background map value
*                                          (M1_2_BG)
      real*8        M1_3BG     (nr__) ! (ct/pix) ?=- 2XMM M1-3 background map value
*                                          (M1_3_BG)
      real*8        M1_4BG     (nr__) ! (ct/pix) ?=- 2XMM M1-4 background map value
*                                          (M1_4_BG)
      real*8        M1_5BG     (nr__) ! (ct/pix) ?=- 2XMM M1-5 background map value
*                                          (M1_5_BG)
      real*8        M2_1BG     (nr__) ! (ct/pix) ?=- 2XMM M2-1 background map value
*                                          (M2_1_BG)
      real*8        M2_2BG     (nr__) ! (ct/pix) ?=- 2XMM M2-2 background map value
*                                          (M2_2_BG)
      real*8        M2_3BG     (nr__) ! (ct/pix) ?=- 2XMM M2-3 background map value
*                                          (M2_3_BG)
      real*8        M2_4BG     (nr__) ! (ct/pix) ?=- 2XMM M2-4 background map value
*                                          (M2_4_BG)
      real*8        M2_5BG     (nr__) ! (ct/pix) ?=- 2XMM M2-5 background map value
*                                          (M2_5_BG)
      real*8        PN1Vig     (nr__) ! ?=- 2XMM PN1 vignetting value
*                                          (PN_1_VIG)
      real*8        PN2Vig     (nr__) ! ?=- 2XMM PN2 vignetting value
*                                          (PN_2_VIG)
      real*8        PN3Vig     (nr__) ! ?=- 2XMM PN3 vignetting value
*                                          (PN_3_VIG)
      real*8        PN4Vig     (nr__) ! ?=- 2XMM PN4 vignetting value
*                                          (PN_4_VIG)
      real*8        PN5Vig     (nr__) ! ?=- 2XMM PN5 vignetting value
*                                          (PN_5_VIG)
      real*8        M1_1Vig    (nr__) ! ?=- 2XMM M1-1 vignetting value
*                                          (M1_1_VIG)
      real*8        M1_2Vig    (nr__) ! ?=- 2XMM M1-2 vignetting value
*                                          (M1_2_VIG)
      real*8        M1_3Vig    (nr__) ! ?=- 2XMM M1-3 vignetting value
*                                          (M1_3_VIG)
      real*8        M1_4Vig    (nr__) ! ?=- 2XMM M1-4 vignetting value
*                                          (M1_3_VIG)
      real*8        M1_5Vig    (nr__) ! ?=- 2XMM M1-5 vignetting value
*                                          (M1_5_VIG)
      real*8        M2_1Vig    (nr__) ! ?=- 2XMM M2-1 vignetting value
*                                          (M2_1_VIG)
      real*8        M2_2Vig    (nr__) ! ?=- 2XMM M2-2 vignetting value
*                                          (M2_3_VIG)
      real*8        M2_3Vig    (nr__) ! ?=- 2XMM M2-3 vignetting value
*                                          (M2_3_VIG)
      real*8        M2_4Vig    (nr__) ! ?=- 2XMM M2-4 vignetting value
*                                          (M2_4_VIG)
      real*8        M2_5Vig    (nr__) ! ?=- 2XMM M2-5 vignetting value
*                                          (M2_5_VIG)
      integer*4     PNonTime   (nr__) ! (s) ?=- 2XMM PN total good exposure time
*                                          (PN_ONTIME)
      integer*4     M1onTime   (nr__) ! (s) ?=- 2XMM M& total good exposure time
*                                          (M1_ONTIME)
      integer*4     M2onTime   (nr__) ! (s) ?=- 2XMM M2 total good exposure time
*                                          (M2_ONTIME)
      real*8        PNoffax    (nr__) ! (arcmin) ?=- 2XMM PN offaxis angle (PN_OFFAX)
      real*8        M1offax    (nr__) ! (arcmin) ?=- 2XMM M1 offaxis angle (M1_OFFAX)
      real*8        M2offax    (nr__) ! (arcmin) ?=- 2XMM M2 offaxis angle (M2_OFFAX)
      real*8        PNMaskFr   (nr__) ! ?=- 2XMM PN PSF weighted mean of the
*                                          detector coverage of a detection
*                                          (PN_MASKFRAC)
      real*8        M1MaskFr   (nr__) ! ?=- 2XMM M1 PSF weighted mean of the
*                                          detector coverage of a detection
*                                          (M1_MASKFRAC)
      real*8        M2MaskFr   (nr__) ! ?=- 2XMM M2 PSF weighted mean of the
*                                          detector coverage of a detection
*                                          (M2_MASKFRAC)
      real*8        DistNN     (nr__) ! (arcsec) 2XMM distance to the nearest
*                                          neighbouring detection (DIST_NN)
      integer*4     SumFlag    (nr__) ! 2XMM summary flag of the detection
*                                          (SUM_FLAG)
      character*12  EPFlag     (nr__) ! 2XMM EPIC flag string (EP_FLAG)
      character*12  PNFlag     (nr__) ! 2XMM PN flag string (PN_FLAG)
      character*12  M1Flag     (nr__) ! 2XMM M1 flag string (M1_FLAG)
      character*12  M2Flag     (nr__) ! 2XMM M2 flag string (M2_FLAG)
      character*1   Tseries    (nr__) ! 2XMM flag set to True if this detection
*                                          has a time series made in at least one
*                                           exposure (TSERIES)
      character*1   spectra    (nr__) ! 2XMM flag set to True if this detection
*                                          has a spectrum made in at least one
*                                          exposure (SPECTRA)
      real*8        EPchi2     (nr__) ! ?=- 2XMM EPIC chi^2^ probability that
*                                          the source is constant (EP_CHI2PROB)
      real*8        PNchi2     (nr__) ! ?=- 2XMM PN chi^2^ probability that the
*                                           source is constant (PN_CHI2PROB)
      real*8        M1chi2     (nr__) ! ?=- 2XMM M1 chi^2^ probability that the
*                                           source is constant (M1_CHI2PROB)
      real*8        M2chi2     (nr__) ! ?=- 2XMM M2 chi^2^ probability that the
*                                           source is constant (M2_CHI2PROB)
      integer*4     VarFlag    (nr__) ! 2XMM flag set to True if this source
*                                          was detected as variable (VAR_FLAG)
      character*4   VarExpID   (nr__) ! 2XMM exposure ID if the source is
*                                          detected as variable (VAR_EXP_ID)
      character*2   VarInstID  (nr__) ! 2XMM instrument ID if the source is
*                                          detected as variable  (VAR_INST_ID)
      real*8        RAsdeg     (nr__) ! (deg) 2XMM source right ascension (J2000)
*                                          (SC_RA)
      real*8        DEsdeg     (nr__) ! (deg) 2XMM source declination (J2000)
*                                          (SC_DEC)
      real*8        ePoss      (nr__) ! (arcsec) 2XMM source total radial detection
*                                          position uncertainty (SC_POSERR)
      real*4        Flux1      (nr__) ! (mW/m2) 2XMM Mean flux in 0.2-0.5keV band
*                                          (SC_EP_1_FLUX)
      real*4        e_Flux1    (nr__) ! (mW/m2) 2XMM Mean error on Flux1
*                                          (SC_EP_1_FLUX_ERR)
      real*4        Flux2      (nr__) ! (mW/m2) 2XMM Mean flux in 0.5-1.0keV band
*                                           (SC_EP_2_FLUX)
      real*4        e_Flux2    (nr__) ! (mW/m2) 2XMM Mean error on Flux2
*                                           (SC_EP_2_FLUX_ERR)
      real*4        Flux3      (nr__) ! (mW/m2) 2XMM Mean flux in 1.0-2.0keV band
*                                           (SC_EP_3_FLUX)
      real*4        e_Flux3    (nr__) ! (mW/m2) 2XMM Mean error on Flux3
*                                           (SC_EP_3_FLUX_ERR)
      real*4        Flux4      (nr__) ! (mW/m2) 2XMM Mean flux in 2.0-4.5keV band
*                                           (SC_EP_4_FLUX)
      real*4        e_Flux4    (nr__) ! (mW/m2) 2XMM Mean error on Flux4
*                                           (SC_EP_4_FLUX_ERR)
      real*4        Flux5      (nr__) ! (mW/m2) 2XMM Mean flux in 4.5-12keV band
*                                           (SC_EP_5_FLUX)
      real*4        e_Flux5    (nr__) ! (mW/m2) 2XMM Mean error on Flux5
*                                           (SC_EP_5_FLUX_ERR)
      real*4        Flux8      (nr__) ! (mW/m2) 2XMM Mean flux in 0.2-12keV band
*                                           (SC_EP_8_FLUX)
      real*4        e_Flux8    (nr__) ! (mW/m2) 2XMM Mean error on Flux8
*                                           (SC_EP_8_FLUX_ERR)
      real*4        Flux9      (nr__) ! (mW/m2) 2XMM Mean flux in 0.5-4.5keV band
*                                           (SC_EP_9_FLUX)
      real*4        e_Flux9    (nr__) ! (mW/m2) 2XMM Mean error on Flux9
*                                           (SC_EP_9_FLUX_ERR)
      real*8        HR1        (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 1,2 (SC_HR1)
      real*8        e_HR1      (nr__) ! ?=- 2XMM Mean error on HR1 (SC_HR1_ERR)
      real*8        HR2        (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 2,3 (SC_HR2)
      real*8        e_HR2      (nr__) ! ?=- 2XMM Mean error on HR2 (SC_HR2_ERR)
      real*8        HR3        (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 3,4 (SC_HR3)
      real*8        e_HR3      (nr__) ! ?=- 2XMM Mean error on HR3 (SC_HR3_ERR)
      real*8        HR4        (nr__) ! [-1/1]?=- 2XMM Hardness ratio,
*                                          bands 4,5 (SC_HR4)
      real*8        e_HR4      (nr__) ! ?=- 2XMM  Mean error on HR4
*                                          (SC_HR4_ERR)
      real*4        srcML      (nr__) ! 2XMM Source detection likelihood
*                                           (SC_DET_ML)
      real*8        extML      (nr__) ! 2XMM Detection likelihood of the
*                                          extended source (SC_EXT_ML)
      real*8        Cst        (nr__) ! [0/1]? 2XMM Constant probability
*                                           (SC_CHI2PROB)
      integer*4     V          (nr__) ! [0/1] 2XMM Variability flag (V=1)
*                                           (SC_VAR_FLAG)
      integer*4     S          (nr__) ! [0/4] 2XMM summary flag, 0=best
*                                           (SC_SUM_FLAG)
      integer*4     Ndet       (nr__) ! 2XMM Number of detection (N_DETECTIONS)
*--------------------------------------------------------------------------------
      real*8        Sep        (nr__) ! (arcsec) ?=- Separation between the 2XMM
*                                          detection and the galaxy position
*                                          (Galactic Separation)
      real*4        LumDet     (nr__) ! (10-7W) Detection 0.2-12.0 keV luminosity and
*                                          its uncertainty (Det_Luminosity)
      real*4        e_LumDet   (nr__) ! (10-7W) rms uncertainty on LumDet
*                                          (Det_Luminosity Error)
      character*6   CP02       (nr__) ! Source names of any ULXs in the
*                                          Colbert & Ptak (2002ApJS..143...25C,
*                                          Cat. J/ApJS/143/25) catalogue likely
*                                          to be matched with the 2XMM detection
*                                          (CP02)
      character*14  LB05       (nr__) ! Source names of any ULXs in the
*                                          Liu & Bregman (2005ApJS..157...59L,
*                                          Cat. J/ApJS/157/59) catalogue likely
*                                          to be matched with the 2XMM detection
*                                          (LB05)
      character*8   SW04       (nr__) ! Record number of any ULXs in the
*                                          Swartz et al. (2009ApJ...703..159S)
*                                          catalogue likely to be matched with
*                                          the 2XMM detection (SW04)
      character*21  LM05       (nr__) ! Source names of any ULXs in the
*                                          Liu & Mirabel (2005A&A...429.1125L,
*                                          Cat. J/A+A/429/1125) catalogue likely
*                                          to be matched with the 2XMM detection
*                                          (LM05)
      character*13  LM05O      (nr__) ! Alternative source names of any ULXs in
*                                          the Liu & Mirabel (2005A&A...429.1125L
*                                          Cat. J/A+A/429/1125) catalogue likely
*                                          to be matched with the 2XMM detection
*                                          (LM05O)
      character*3   New_       (nr__) ! Flag stating whether the source is a
*                                          'new' ULX candidate (i.e. not present
*                                          in the above catalogues) (New?)
      character*3   Comp       (nr__) ! Flag stating whether the detection of a
*                                          source is included in the sub-sample
*                                          compiled from complete observations
*                                          (Complete_subsample?)

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

C  Loading file 'catalog.dat'	! catalogue of ultraluminous X-ray source candidates

C  Format for file interpretation

    1 format(
     +  A12,1X,A14,1X,A14,1X,A9,A7,1X,A3,1X,F4.1,1X,F3.1,1X,F4.1,1X,
     +  F3.1,1X,I1,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,A1,1X,F4.2,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,I3,1X,F5.2,1X,A1,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,
     +  1X,F3.1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,
     +  1X,I3,1X,I2,1X,I3,1X,I2,1X,I5,1X,I2,1X,I5,1X,I3,1X,I5,1X,I5,
     +  1X,F10.6,1X,F10.6,1X,F8.5,1X,F6.2,1X,I6,1X,I6,1X,A21,1X,I3,6X,
     +  A16,1X,F8.6,1X,I6,7X,A16,1X,F8.6,1X,I10,1X,I4,1X,F7.1,1X,F7.1,
     +  1X,I1,1X,A6,1X,A6,1X,A6,1X,A23,1X,A15,1X,A15,1X,F10.6,1X,
     +  F10.6,1X,F8.6,1X,F10.6,1X,F10.6,1X,F8.6,1X,F4.2,1X,F10.6,1X,
     +  F10.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F8.6,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
     +  1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
     +  1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
     +  1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F12.5,1X,
     +  F9.5,1X,F12.5,1X,F9.5,1X,F12.5,1X,F9.5,1X,F12.5,1X,F9.5,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,7X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,
     +  1X,F9.6,1X,F9.6,1X,F9.6,1X,F11.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,
     +  1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F10.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F9.6,
     +  1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F8.6,
     +  1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,I6,1X,
     +  I6,1X,I6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
     +  F9.5,1X,I1,1X,A12,1X,A12,1X,A12,1X,A12,1X,A1,1X,A1,1X,F8.6,1X,
     +  F8.6,1X,F8.6,1X,F8.6,1X,I1,1X,A4,1X,A2,1X,F10.6,1X,F10.6,1X,
     +  F8.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F9.6,1X,E12.6,1X,F8.6,1X,F8.6,1X,I1,1X,I1,1X,I2,1X,
     +  F10.6,1X,E22.15,1X,E22.15,1X,A6,1X,A14,1X,A8,1X,A21,1X,A13,1X,
     +  A3,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,650
        read(1,'(A3658)')ar__
        read(ar__,1)
     +  name(i__),altname(i__),desig(i__),PGC(i__),type(i__),
     +  typesr(i__),T(i__),e_T(i__),lumcl(i__),e_lumcl(i__),
     +  o_lumcl(i__),D25(i__),u_D25(i__),e_D25(i__),R25(i__),
     +  u_R25(i__),e_R25(i__),Do(i__),Ae(i__),e_Ae(i__),PA(i__),
     +  BT(i__),BT_code(i__),e_BT(i__),Bmag(i__),e_Bmag(i__),BoT(i__),
     +  m_25(i__),e_m_25(i__),m_e(i__),e_m_e(i__),mFIR(i__),m21(i__),
     +  e_m21(i__),B_VT(i__),e_B_VT(i__),B_Ve(i__),e_B_Ve(i__),
     +  B_VoT(i__),U_BT(i__),e_U_BT(i__),U_Be(i__),e_U_Be(i__),
     +  U_BoT(i__),HI(i__),Ai(i__),A21(i__),Ag(i__),W20(i__),
     +  e_W20(i__),W50(i__),e_W50(i__),V21(i__),e_V21(i__),cz(i__),
     +  e_cz(i__),VGSR(i__),V3K(i__),RAdeg(i__),DEdeg(i__),ePos(i__),
     +  Dist(i__),DETID(i__),SRCID(i__),IAUName(i__),SRC(i__),
     +  v1XMM(i__),Sep1XMM(i__),v2XMMpID(i__),v2XMMp(i__),
     +  Sep2XMMP(i__),ObsID(i__),Rev(i__),MJD0(i__),MJD1(i__),Q(i__),
     +  f_PN(i__),f_M1(i__),f_M2(i__),m_PN(i__),m_M1(i__),m_M2(i__),
     +  RA1deg(i__),DE1deg(i__),ePos1(i__),GLON(i__),GLAT(i__),
     +  ePos2(i__),SYSERR(i__),RAudeg(i__),DEudeg(i__),FluxEP1(i__),
     +  e_FluxEP1(i__),FluxEP2(i__),e_FluxEP2(i__),FluxEP3(i__),
     +  e_FluxEP3(i__),FluxEP4(i__),e_FluxEP4(i__),FluxEP5(i__),
     +  e_FluxEP5(i__),FluxEP8(i__),e_FluxEP8(i__),FluxEP9(i__),
     +  e_FluxEP9(i__),FluxPN1(i__),e_FluxPN1(i__),FluxPN2(i__),
     +  e_FluxPN2(i__),FluxPN3(i__),e_FluxPN3(i__),FluxPN4(i__),
     +  e_FluxPN4(i__),FluxPN5(i__),e_FluxPN5(i__),FluxPN8(i__),
     +  e_FluxPN8(i__),FluxPN9(i__),e_FluxPN9(i__),FluxM1_1(i__),
     +  e_FluxM1_1(i__),FluxM1_2(i__),e_FluxM1_2(i__),FluxM1_3(i__),
     +  e_FluxM1_3(i__),FluxM1_4(i__),e_FluxM1_4(i__),FluxM1_5(i__),
     +  e_FluxM1_5(i__),FluxM1_8(i__),e_FluxM1_8(i__),FluxM1_9(i__),
     +  e_FluxM1_9(i__),FluxM2_1(i__),e_FluxM2_1(i__),FluxM2_2(i__),
     +  e_FluxM2_2(i__),FluxM2_3(i__),e_FluxM2_3(i__),FluxM2_4(i__),
     +  e_FluxM2_4(i__),FluxM2_5(i__),e_FluxM2_5(i__),FluxM2_8(i__),
     +  e_FluxM2_8(i__),FluxM2_9(i__),e_FluxM2_9(i__),EP8rate(i__),
     +  e_EP8rate(i__),EP9rate(i__),e_EP9rate(i__),PN1rate(i__),
     +  e_PN1rate(i__),PN2rate(i__),e_PN2rate(i__),PN3rate(i__),
     +  e_PN3rate(i__),PN4rate(i__),e_PN4rate(i__),PN5rate(i__),
     +  e_PN5rate(i__),PN8rate(i__),e_PN8rate(i__),PN9rate(i__),
     +  e_PN9rate(i__),M1_1rate(i__),e_M1_1rate(i__),M1_2rate(i__),
     +  e_M1_2rate(i__),M1_3rate(i__),e_M1_3rate(i__),M1_4rate(i__),
     +  e_M1_4rate(i__),M1_5rate(i__),e_M1_5rate(i__),M1_8rate(i__),
     +  e_M1_8rate(i__),M1_9rate(i__),e_M1_9rate(i__),M2_1rate(i__),
     +  e_M2_1rate(i__),M2_2rate(i__),e_M2_2rate(i__),M2_3rate(i__),
     +  e_M2_3rate(i__),M2_4rate(i__),e_M2_4rate(i__),M2_5rate(i__),
     +  e_M2_5rate(i__),M2_8rate(i__),e_M2_8rate(i__),M2_9rate(i__),
     +  e_M2_9rate(i__),EP8cts(i__),e_EP8cts(i__),PN8cts(i__),
     +  e_PN8cts(i__),M1_8cts(i__),e_M1_8cts(i__),M2_8cts(i__),
     +  e_M2_8cts(i__),EP8detML(i__),EP9detML(i__),PN1detML(i__),
     +  PN2detML(i__),PN3detML(i__),PN4detML(i__),PN5detML(i__),
     +  PN8detML(i__),PN9detML(i__),M1_1detML(i__),M1_2detML(i__),
     +  M1_3detML(i__),M1_4detML(i__),M1_5detML(i__),M1_8detML(i__),
     +  M1_9detML(i__),M2_1detML(i__),M2_2detML(i__),M2_3detML(i__),
     +  M2_4detML(i__),M2_5detML(i__),M2_8detML(i__),M2_9detML(i__),
     +  HREP1(i__),e_HREP1(i__),HREP2(i__),e_HREP2(i__),HREP3(i__),
     +  e_HREP3(i__),HREP4(i__),e_HREP4(i__),HRPN1(i__),e_HRPN1(i__),
     +  HRPN2(i__),e_HRPN2(i__),HRPN3(i__),e_HRPN3(i__),HRPN4(i__),
     +  e_HRPN4(i__),HRM1_1(i__),e_HRM1_1(i__),HRM1_2(i__),
     +  e_HRM1_2(i__),HRM1_3(i__),e_HRM1_3(i__),HRM1_4(i__),
     +  e_HRM1_4(i__),HRM2_1(i__),e_HRM2_1(i__),HRM2_2(i__),
     +  e_HRM2_2(i__),HRM2_3(i__),e_HRM2_3(i__),HRM2_4(i__),
     +  e_HRM2_4(i__),PN1ExpT(i__),PN2ExpT(i__),PN3ExpT(i__),
     +  PN4ExpT(i__),PN5ExpT(i__),M1_1ExpT(i__),M1_2ExpT(i__),
     +  M1_3ExpT(i__),M1_4ExpT(i__),M1_5ExpT(i__),M2_1ExpT(i__),
     +  M2_2ExpT(i__),M2_3ExpT(i__),M2_4ExpT(i__),M2_5ExpT(i__),
     +  PN1BG(i__),PN2BG(i__),PN3BG(i__),PN4BG(i__),PN5BG(i__),
     +  M1_1BG(i__),M1_2BG(i__),M1_3BG(i__),M1_4BG(i__),M1_5BG(i__),
     +  M2_1BG(i__),M2_2BG(i__),M2_3BG(i__),M2_4BG(i__),M2_5BG(i__),
     +  PN1Vig(i__),PN2Vig(i__),PN3Vig(i__),PN4Vig(i__),PN5Vig(i__),
     +  M1_1Vig(i__),M1_2Vig(i__),M1_3Vig(i__),M1_4Vig(i__),
     +  M1_5Vig(i__),M2_1Vig(i__),M2_2Vig(i__),M2_3Vig(i__),
     +  M2_4Vig(i__),M2_5Vig(i__),PNonTime(i__),M1onTime(i__),
     +  M2onTime(i__),PNoffax(i__),M1offax(i__),M2offax(i__),
     +  PNMaskFr(i__),M1MaskFr(i__),M2MaskFr(i__),DistNN(i__),
     +  SumFlag(i__),EPFlag(i__),PNFlag(i__),M1Flag(i__),M2Flag(i__),
     +  Tseries(i__),spectra(i__),EPchi2(i__),PNchi2(i__),M1chi2(i__),
     +  M2chi2(i__),VarFlag(i__),VarExpID(i__),VarInstID(i__),
     +  RAsdeg(i__),DEsdeg(i__),ePoss(i__),Flux1(i__),e_Flux1(i__),
     +  Flux2(i__),e_Flux2(i__),Flux3(i__),e_Flux3(i__),Flux4(i__),
     +  e_Flux4(i__),Flux5(i__),e_Flux5(i__),Flux8(i__),e_Flux8(i__),
     +  Flux9(i__),e_Flux9(i__),HR1(i__),e_HR1(i__),HR2(i__),
     +  e_HR2(i__),HR3(i__),e_HR3(i__),HR4(i__),e_HR4(i__),srcML(i__),
     +  extML(i__),Cst(i__),V(i__),S(i__),Ndet(i__),Sep(i__),
     +  LumDet(i__),e_LumDet(i__),CP02(i__),LB05(i__),SW04(i__),
     +  LM05(i__),LM05O(i__),New_(i__),Comp(i__)
        if(ar__(65:68) .EQ. '') T(i__) = rNULL__
        if(ar__(70:72) .EQ. '') e_T(i__) = rNULL__
        if(ar__(74:77) .EQ. '') lumcl(i__) = rNULL__
        if(ar__(79:81) .EQ. '') e_lumcl(i__) = rNULL__
        if(ar__(83:83) .EQ. '') o_lumcl(i__) = iNULL__
        if(ar__(85:88) .EQ. '') D25(i__) = rNULL__
        if(ar__(92:95) .EQ. '') e_D25(i__) = rNULL__
        if(ar__(97:100) .EQ. '') R25(i__) = rNULL__
        if(ar__(104:107) .EQ. '') e_R25(i__) = rNULL__
        if(ar__(109:112) .EQ. '') Do(i__) = rNULL__
        if(ar__(114:117) .EQ. '') Ae(i__) = rNULL__
        if(ar__(119:122) .EQ. '') e_Ae(i__) = rNULL__
        if(ar__(124:126) .EQ. '') PA(i__) = iNULL__
        if(ar__(128:132) .EQ. '') BT(i__) = rNULL__
        if(ar__(136:139) .EQ. '') e_BT(i__) = rNULL__
        if(ar__(141:145) .EQ. '') Bmag(i__) = rNULL__
        if(ar__(147:150) .EQ. '') e_Bmag(i__) = rNULL__
        if(ar__(152:156) .EQ. '') BoT(i__) = rNULL__
        if(ar__(158:162) .EQ. '') m_25(i__) = rNULL__
        if(ar__(164:167) .EQ. '') e_m_25(i__) = rNULL__
        if(ar__(169:173) .EQ. '') m_e(i__) = rNULL__
        if(ar__(175:178) .EQ. '') e_m_e(i__) = rNULL__
        if(ar__(180:184) .EQ. '') mFIR(i__) = rNULL__
        if(ar__(186:190) .EQ. '') m21(i__) = rNULL__
        if(ar__(192:194) .EQ. '') e_m21(i__) = rNULL__
        if(ar__(196:199) .EQ. '') B_VT(i__) = rNULL__
        if(ar__(201:204) .EQ. '') e_B_VT(i__) = rNULL__
        if(ar__(206:209) .EQ. '') B_Ve(i__) = rNULL__
        if(ar__(211:214) .EQ. '') e_B_Ve(i__) = rNULL__
        if(ar__(216:219) .EQ. '') B_VoT(i__) = rNULL__
        if(ar__(221:225) .EQ. '') U_BT(i__) = rNULL__
        if(ar__(227:230) .EQ. '') e_U_BT(i__) = rNULL__
        if(ar__(232:236) .EQ. '') U_Be(i__) = rNULL__
        if(ar__(238:241) .EQ. '') e_U_Be(i__) = rNULL__
        if(ar__(243:247) .EQ. '') U_BoT(i__) = rNULL__
        if(ar__(249:253) .EQ. '') HI(i__) = rNULL__
        if(ar__(255:258) .EQ. '') Ai(i__) = rNULL__
        if(ar__(260:263) .EQ. '') A21(i__) = rNULL__
        if(ar__(265:268) .EQ. '') Ag(i__) = rNULL__
        if(ar__(270:272) .EQ. '') W20(i__) = iNULL__
        if(ar__(274:275) .EQ. '') e_W20(i__) = iNULL__
        if(ar__(277:279) .EQ. '') W50(i__) = iNULL__
        if(ar__(281:282) .EQ. '') e_W50(i__) = iNULL__
        if(ar__(284:288) .EQ. '') V21(i__) = iNULL__
        if(ar__(290:291) .EQ. '') e_V21(i__) = iNULL__
        if(ar__(293:297) .EQ. '') cz(i__) = iNULL__
        if(ar__(299:301) .EQ. '') e_cz(i__) = iNULL__
        if(ar__(303:307) .EQ. '') VGSR(i__) = iNULL__
        if(ar__(309:313) .EQ. '') V3K(i__) = iNULL__
        if (DETID(i__) .EQ. 45) DETID(i__) =  iNULL__
        if (SRC(i__) .EQ. 45) SRC(i__) =  iNULL__
        if (idig(ar__(415:422)).EQ.0) Sep1XMM(i__) =  rNULL__
        if (v2XMMpID(i__) .EQ. 45) v2XMMpID(i__) =  iNULL__
        if (idig(ar__(454:461)).EQ.0) Sep2XMMP(i__) =  rNULL__
        if (idig(ar__(845:856)).EQ.0) FluxPN1(i__) =  rNULL__
        if (idig(ar__(858:869)).EQ.0) e_FluxPN1(i__) =  rNULL__
        if (idig(ar__(871:882)).EQ.0) FluxPN2(i__) =  rNULL__
        if (idig(ar__(884:895)).EQ.0) e_FluxPN2(i__) =  rNULL__
        if (idig(ar__(897:908)).EQ.0) FluxPN3(i__) =  rNULL__
        if (idig(ar__(910:921)).EQ.0) e_FluxPN3(i__) =  rNULL__
        if (idig(ar__(923:934)).EQ.0) FluxPN4(i__) =  rNULL__
        if (idig(ar__(936:947)).EQ.0) e_FluxPN4(i__) =  rNULL__
        if (idig(ar__(949:960)).EQ.0) FluxPN5(i__) =  rNULL__
        if (idig(ar__(962:973)).EQ.0) e_FluxPN5(i__) =  rNULL__
        if (idig(ar__(975:986)).EQ.0) FluxPN8(i__) =  rNULL__
        if (idig(ar__(988:999)).EQ.0) e_FluxPN8(i__) =  rNULL__
        if (idig(ar__(1001:1012)).EQ.0) FluxPN9(i__) =  rNULL__
        if (idig(ar__(1014:1025)).EQ.0) e_FluxPN9(i__) =  rNULL__
        if (idig(ar__(1027:1038)).EQ.0) FluxM1_1(i__) =  rNULL__
        if (idig(ar__(1040:1051)).EQ.0) e_FluxM1_1(i__) =  rNULL__
        if (idig(ar__(1053:1064)).EQ.0) FluxM1_2(i__) =  rNULL__
        if (idig(ar__(1066:1077)).EQ.0) e_FluxM1_2(i__) =  rNULL__
        if (idig(ar__(1079:1090)).EQ.0) FluxM1_3(i__) =  rNULL__
        if (idig(ar__(1092:1103)).EQ.0) e_FluxM1_3(i__) =  rNULL__
        if (idig(ar__(1105:1116)).EQ.0) FluxM1_4(i__) =  rNULL__
        if (idig(ar__(1118:1129)).EQ.0) e_FluxM1_4(i__) =  rNULL__
        if (idig(ar__(1131:1142)).EQ.0) FluxM1_5(i__) =  rNULL__
        if (idig(ar__(1144:1155)).EQ.0) e_FluxM1_5(i__) =  rNULL__
        if (idig(ar__(1157:1168)).EQ.0) FluxM1_8(i__) =  rNULL__
        if (idig(ar__(1170:1181)).EQ.0) e_FluxM1_8(i__) =  rNULL__
        if (idig(ar__(1183:1194)).EQ.0) FluxM1_9(i__) =  rNULL__
        if (idig(ar__(1196:1207)).EQ.0) e_FluxM1_9(i__) =  rNULL__
        if (idig(ar__(1209:1220)).EQ.0) FluxM2_1(i__) =  rNULL__
        if (idig(ar__(1222:1233)).EQ.0) e_FluxM2_1(i__) =  rNULL__
        if (idig(ar__(1235:1246)).EQ.0) FluxM2_2(i__) =  rNULL__
        if (idig(ar__(1248:1259)).EQ.0) e_FluxM2_2(i__) =  rNULL__
        if (idig(ar__(1261:1272)).EQ.0) FluxM2_3(i__) =  rNULL__
        if (idig(ar__(1274:1285)).EQ.0) e_FluxM2_3(i__) =  rNULL__
        if (idig(ar__(1287:1298)).EQ.0) FluxM2_4(i__) =  rNULL__
        if (idig(ar__(1300:1311)).EQ.0) e_FluxM2_4(i__) =  rNULL__
        if (idig(ar__(1313:1324)).EQ.0) FluxM2_5(i__) =  rNULL__
        if (idig(ar__(1326:1337)).EQ.0) e_FluxM2_5(i__) =  rNULL__
        if (idig(ar__(1339:1350)).EQ.0) FluxM2_8(i__) =  rNULL__
        if (idig(ar__(1352:1363)).EQ.0) e_FluxM2_8(i__) =  rNULL__
        if (idig(ar__(1365:1376)).EQ.0) FluxM2_9(i__) =  rNULL__
        if (idig(ar__(1378:1389)).EQ.0) e_FluxM2_9(i__) =  rNULL__
        if (idig(ar__(1828:1839)).EQ.0) PN8cts(i__) =  rNULL__
        if (idig(ar__(1841:1849)).EQ.0) e_PN8cts(i__) =  rNULL__
        if (idig(ar__(1851:1862)).EQ.0) M1_8cts(i__) =  rNULL__
        if (idig(ar__(1864:1872)).EQ.0) e_M1_8cts(i__) =  rNULL__
        if (idig(ar__(1874:1885)).EQ.0) M2_8cts(i__) =  rNULL__
        if (idig(ar__(1887:1895)).EQ.0) e_M2_8cts(i__) =  rNULL__
        if (idig(ar__(1923:1934)).EQ.0) PN1detML(i__) =  rNULL__
        if (idig(ar__(1936:1947)).EQ.0) PN2detML(i__) =  rNULL__
        if (idig(ar__(1949:1960)).EQ.0) PN3detML(i__) =  rNULL__
        if (idig(ar__(1962:1973)).EQ.0) PN4detML(i__) =  rNULL__
        if (idig(ar__(1975:1986)).EQ.0) PN5detML(i__) =  rNULL__
        if (idig(ar__(1988:1999)).EQ.0) PN8detML(i__) =  rNULL__
        if (idig(ar__(2001:2012)).EQ.0) PN9detML(i__) =  rNULL__
        if (idig(ar__(2014:2025)).EQ.0) M1_1detML(i__) =  rNULL__
        if (idig(ar__(2027:2038)).EQ.0) M1_2detML(i__) =  rNULL__
        if (idig(ar__(2040:2051)).EQ.0) M1_3detML(i__) =  rNULL__
        if (idig(ar__(2053:2064)).EQ.0) M1_4detML(i__) =  rNULL__
        if (idig(ar__(2066:2077)).EQ.0) M1_5detML(i__) =  rNULL__
        if (idig(ar__(2079:2090)).EQ.0) M1_8detML(i__) =  rNULL__
        if (idig(ar__(2092:2103)).EQ.0) M1_9detML(i__) =  rNULL__
        if (idig(ar__(2105:2116)).EQ.0) M2_1detML(i__) =  rNULL__
        if (idig(ar__(2118:2129)).EQ.0) M2_2detML(i__) =  rNULL__
        if (idig(ar__(2131:2142)).EQ.0) M2_3detML(i__) =  rNULL__
        if (idig(ar__(2144:2155)).EQ.0) M2_4detML(i__) =  rNULL__
        if (idig(ar__(2157:2168)).EQ.0) M2_5detML(i__) =  rNULL__
        if (idig(ar__(2170:2181)).EQ.0) M2_8detML(i__) =  rNULL__
        if (idig(ar__(2183:2194)).EQ.0) M2_9detML(i__) =  rNULL__
        if (idig(ar__(2202:2210)).EQ.0) HREP1(i__) =  rNULL__
        if (idig(ar__(2212:2220)).EQ.0) e_HREP1(i__) =  rNULL__
        if (idig(ar__(2222:2230)).EQ.0) HREP2(i__) =  rNULL__
        if (idig(ar__(2232:2240)).EQ.0) e_HREP2(i__) =  rNULL__
        if (idig(ar__(2242:2250)).EQ.0) HREP3(i__) =  rNULL__
        if (idig(ar__(2252:2260)).EQ.0) e_HREP3(i__) =  rNULL__
        if (idig(ar__(2262:2270)).EQ.0) HREP4(i__) =  rNULL__
        if (idig(ar__(2272:2280)).EQ.0) e_HREP4(i__) =  rNULL__
        if (idig(ar__(2282:2290)).EQ.0) HRPN1(i__) =  rNULL__
        if (idig(ar__(2292:2300)).EQ.0) e_HRPN1(i__) =  rNULL__
        if (idig(ar__(2302:2310)).EQ.0) HRPN2(i__) =  rNULL__
        if (idig(ar__(2312:2320)).EQ.0) e_HRPN2(i__) =  rNULL__
        if (idig(ar__(2322:2330)).EQ.0) HRPN3(i__) =  rNULL__
        if (idig(ar__(2332:2340)).EQ.0) e_HRPN3(i__) =  rNULL__
        if (idig(ar__(2342:2350)).EQ.0) HRPN4(i__) =  rNULL__
        if (idig(ar__(2352:2362)).EQ.0) e_HRPN4(i__) =  rNULL__
        if (idig(ar__(2364:2372)).EQ.0) HRM1_1(i__) =  rNULL__
        if (idig(ar__(2374:2382)).EQ.0) e_HRM1_1(i__) =  rNULL__
        if (idig(ar__(2384:2392)).EQ.0) HRM1_2(i__) =  rNULL__
        if (idig(ar__(2394:2402)).EQ.0) e_HRM1_2(i__) =  rNULL__
        if (idig(ar__(2404:2412)).EQ.0) HRM1_3(i__) =  rNULL__
        if (idig(ar__(2414:2422)).EQ.0) e_HRM1_3(i__) =  rNULL__
        if (idig(ar__(2424:2432)).EQ.0) HRM1_4(i__) =  rNULL__
        if (idig(ar__(2434:2442)).EQ.0) e_HRM1_4(i__) =  rNULL__
        if (idig(ar__(2444:2452)).EQ.0) HRM2_1(i__) =  rNULL__
        if (idig(ar__(2454:2462)).EQ.0) e_HRM2_1(i__) =  rNULL__
        if (idig(ar__(2464:2472)).EQ.0) HRM2_2(i__) =  rNULL__
        if (idig(ar__(2474:2482)).EQ.0) e_HRM2_2(i__) =  rNULL__
        if (idig(ar__(2484:2492)).EQ.0) HRM2_3(i__) =  rNULL__
        if (idig(ar__(2494:2502)).EQ.0) e_HRM2_3(i__) =  rNULL__
        if (idig(ar__(2504:2512)).EQ.0) HRM2_4(i__) =  rNULL__
        if (idig(ar__(2514:2523)).EQ.0) e_HRM2_4(i__) =  rNULL__
        if (idig(ar__(2525:2536)).EQ.0) PN1ExpT(i__) =  rNULL__
        if (idig(ar__(2538:2549)).EQ.0) PN2ExpT(i__) =  rNULL__
        if (idig(ar__(2551:2562)).EQ.0) PN3ExpT(i__) =  rNULL__
        if (idig(ar__(2564:2575)).EQ.0) PN4ExpT(i__) =  rNULL__
        if (idig(ar__(2577:2588)).EQ.0) PN5ExpT(i__) =  rNULL__
        if (idig(ar__(2590:2601)).EQ.0) M1_1ExpT(i__) =  rNULL__
        if (idig(ar__(2603:2614)).EQ.0) M1_2ExpT(i__) =  rNULL__
        if (idig(ar__(2616:2627)).EQ.0) M1_3ExpT(i__) =  rNULL__
        if (idig(ar__(2629:2640)).EQ.0) M1_4ExpT(i__) =  rNULL__
        if (idig(ar__(2642:2653)).EQ.0) M1_5ExpT(i__) =  rNULL__
        if (idig(ar__(2655:2666)).EQ.0) M2_1ExpT(i__) =  rNULL__
        if (idig(ar__(2668:2679)).EQ.0) M2_2ExpT(i__) =  rNULL__
        if (idig(ar__(2681:2692)).EQ.0) M2_3ExpT(i__) =  rNULL__
        if (idig(ar__(2694:2705)).EQ.0) M2_4ExpT(i__) =  rNULL__
        if (idig(ar__(2707:2718)).EQ.0) M2_5ExpT(i__) =  rNULL__
        if (idig(ar__(2720:2728)).EQ.0) PN1BG(i__) =  rNULL__
        if (idig(ar__(2730:2738)).EQ.0) PN2BG(i__) =  rNULL__
        if (idig(ar__(2740:2748)).EQ.0) PN3BG(i__) =  rNULL__
        if (idig(ar__(2750:2758)).EQ.0) PN4BG(i__) =  rNULL__
        if (idig(ar__(2760:2768)).EQ.0) PN5BG(i__) =  rNULL__
        if (idig(ar__(2770:2778)).EQ.0) M1_1BG(i__) =  rNULL__
        if (idig(ar__(2780:2788)).EQ.0) M1_2BG(i__) =  rNULL__
        if (idig(ar__(2790:2798)).EQ.0) M1_3BG(i__) =  rNULL__
        if (idig(ar__(2800:2808)).EQ.0) M1_4BG(i__) =  rNULL__
        if (idig(ar__(2810:2818)).EQ.0) M1_5BG(i__) =  rNULL__
        if (idig(ar__(2820:2828)).EQ.0) M2_1BG(i__) =  rNULL__
        if (idig(ar__(2830:2838)).EQ.0) M2_2BG(i__) =  rNULL__
        if (idig(ar__(2840:2848)).EQ.0) M2_3BG(i__) =  rNULL__
        if (idig(ar__(2850:2858)).EQ.0) M2_4BG(i__) =  rNULL__
        if (idig(ar__(2860:2868)).EQ.0) M2_5BG(i__) =  rNULL__
        if (idig(ar__(2870:2877)).EQ.0) PN1Vig(i__) =  rNULL__
        if (idig(ar__(2879:2886)).EQ.0) PN2Vig(i__) =  rNULL__
        if (idig(ar__(2888:2895)).EQ.0) PN3Vig(i__) =  rNULL__
        if (idig(ar__(2897:2904)).EQ.0) PN4Vig(i__) =  rNULL__
        if (idig(ar__(2906:2913)).EQ.0) PN5Vig(i__) =  rNULL__
        if (idig(ar__(2915:2922)).EQ.0) M1_1Vig(i__) =  rNULL__
        if (idig(ar__(2924:2931)).EQ.0) M1_2Vig(i__) =  rNULL__
        if (idig(ar__(2933:2940)).EQ.0) M1_3Vig(i__) =  rNULL__
        if (idig(ar__(2942:2949)).EQ.0) M1_4Vig(i__) =  rNULL__
        if (idig(ar__(2951:2958)).EQ.0) M1_5Vig(i__) =  rNULL__
        if (idig(ar__(2960:2967)).EQ.0) M2_1Vig(i__) =  rNULL__
        if (idig(ar__(2969:2976)).EQ.0) M2_2Vig(i__) =  rNULL__
        if (idig(ar__(2978:2985)).EQ.0) M2_3Vig(i__) =  rNULL__
        if (idig(ar__(2987:2994)).EQ.0) M2_4Vig(i__) =  rNULL__
        if (idig(ar__(2996:3003)).EQ.0) M2_5Vig(i__) =  rNULL__
        if (PNonTime(i__) .EQ. 45) PNonTime(i__) =  iNULL__
        if (M1onTime(i__) .EQ. 45) M1onTime(i__) =  iNULL__
        if (M2onTime(i__) .EQ. 45) M2onTime(i__) =  iNULL__
        if (idig(ar__(3026:3034)).EQ.0) PNoffax(i__) =  rNULL__
        if (idig(ar__(3036:3044)).EQ.0) M1offax(i__) =  rNULL__
        if (idig(ar__(3046:3054)).EQ.0) M2offax(i__) =  rNULL__
        if (idig(ar__(3056:3063)).EQ.0) PNMaskFr(i__) =  rNULL__
        if (idig(ar__(3065:3072)).EQ.0) M1MaskFr(i__) =  rNULL__
        if (idig(ar__(3074:3081)).EQ.0) M2MaskFr(i__) =  rNULL__
        if (idig(ar__(3151:3158)).EQ.0) EPchi2(i__) =  rNULL__
        if (idig(ar__(3160:3167)).EQ.0) PNchi2(i__) =  rNULL__
        if (idig(ar__(3169:3176)).EQ.0) M1chi2(i__) =  rNULL__
        if (idig(ar__(3178:3185)).EQ.0) M2chi2(i__) =  rNULL__
        if (idig(ar__(3410:3418)).EQ.0) HR1(i__) =  rNULL__
        if (idig(ar__(3420:3428)).EQ.0) e_HR1(i__) =  rNULL__
        if (idig(ar__(3430:3438)).EQ.0) HR2(i__) =  rNULL__
        if (idig(ar__(3440:3448)).EQ.0) e_HR2(i__) =  rNULL__
        if (idig(ar__(3450:3458)).EQ.0) HR3(i__) =  rNULL__
        if (idig(ar__(3460:3468)).EQ.0) e_HR3(i__) =  rNULL__
        if (idig(ar__(3470:3478)).EQ.0) HR4(i__) =  rNULL__
        if (idig(ar__(3480:3488)).EQ.0) e_HR4(i__) =  rNULL__
        if(ar__(3512:3519) .EQ. '') Cst(i__) = rNULL__
        if (idig(ar__(3528:3537)).EQ.0) Sep(i__) =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  name(i__),altname(i__),desig(i__),PGC(i__),type(i__),
     +  typesr(i__),T(i__),e_T(i__),lumcl(i__),e_lumcl(i__),
     +  o_lumcl(i__),D25(i__),u_D25(i__),e_D25(i__),R25(i__),
     +  u_R25(i__),e_R25(i__),Do(i__),Ae(i__),e_Ae(i__),PA(i__),
     +  BT(i__),BT_code(i__),e_BT(i__),Bmag(i__),e_Bmag(i__),BoT(i__),
     +  m_25(i__),e_m_25(i__),m_e(i__),e_m_e(i__),mFIR(i__),m21(i__),
     +  e_m21(i__),B_VT(i__),e_B_VT(i__),B_Ve(i__),e_B_Ve(i__),
     +  B_VoT(i__),U_BT(i__),e_U_BT(i__),U_Be(i__),e_U_Be(i__),
     +  U_BoT(i__),HI(i__),Ai(i__),A21(i__),Ag(i__),W20(i__),
     +  e_W20(i__),W50(i__),e_W50(i__),V21(i__),e_V21(i__),cz(i__),
     +  e_cz(i__),VGSR(i__),V3K(i__),RAdeg(i__),DEdeg(i__),ePos(i__),
     +  Dist(i__),DETID(i__),SRCID(i__),IAUName(i__),SRC(i__),
     +  v1XMM(i__),Sep1XMM(i__),v2XMMpID(i__),v2XMMp(i__),
     +  Sep2XMMP(i__),ObsID(i__),Rev(i__),MJD0(i__),MJD1(i__),Q(i__),
     +  f_PN(i__),f_M1(i__),f_M2(i__),m_PN(i__),m_M1(i__),m_M2(i__),
     +  RA1deg(i__),DE1deg(i__),ePos1(i__),GLON(i__),GLAT(i__),
     +  ePos2(i__),SYSERR(i__),RAudeg(i__),DEudeg(i__),FluxEP1(i__),
     +  e_FluxEP1(i__),FluxEP2(i__),e_FluxEP2(i__),FluxEP3(i__),
     +  e_FluxEP3(i__),FluxEP4(i__),e_FluxEP4(i__),FluxEP5(i__),
     +  e_FluxEP5(i__),FluxEP8(i__),e_FluxEP8(i__),FluxEP9(i__),
     +  e_FluxEP9(i__),FluxPN1(i__),e_FluxPN1(i__),FluxPN2(i__),
     +  e_FluxPN2(i__),FluxPN3(i__),e_FluxPN3(i__),FluxPN4(i__),
     +  e_FluxPN4(i__),FluxPN5(i__),e_FluxPN5(i__),FluxPN8(i__),
     +  e_FluxPN8(i__),FluxPN9(i__),e_FluxPN9(i__),FluxM1_1(i__),
     +  e_FluxM1_1(i__),FluxM1_2(i__),e_FluxM1_2(i__),FluxM1_3(i__),
     +  e_FluxM1_3(i__),FluxM1_4(i__),e_FluxM1_4(i__),FluxM1_5(i__),
     +  e_FluxM1_5(i__),FluxM1_8(i__),e_FluxM1_8(i__),FluxM1_9(i__),
     +  e_FluxM1_9(i__),FluxM2_1(i__),e_FluxM2_1(i__),FluxM2_2(i__),
     +  e_FluxM2_2(i__),FluxM2_3(i__),e_FluxM2_3(i__),FluxM2_4(i__),
     +  e_FluxM2_4(i__),FluxM2_5(i__),e_FluxM2_5(i__),FluxM2_8(i__),
     +  e_FluxM2_8(i__),FluxM2_9(i__),e_FluxM2_9(i__),EP8rate(i__),
     +  e_EP8rate(i__),EP9rate(i__),e_EP9rate(i__),PN1rate(i__),
     +  e_PN1rate(i__),PN2rate(i__),e_PN2rate(i__),PN3rate(i__),
     +  e_PN3rate(i__),PN4rate(i__),e_PN4rate(i__),PN5rate(i__),
     +  e_PN5rate(i__),PN8rate(i__),e_PN8rate(i__),PN9rate(i__),
     +  e_PN9rate(i__),M1_1rate(i__),e_M1_1rate(i__),M1_2rate(i__),
     +  e_M1_2rate(i__),M1_3rate(i__),e_M1_3rate(i__),M1_4rate(i__),
     +  e_M1_4rate(i__),M1_5rate(i__),e_M1_5rate(i__),M1_8rate(i__),
     +  e_M1_8rate(i__),M1_9rate(i__),e_M1_9rate(i__),M2_1rate(i__),
     +  e_M2_1rate(i__),M2_2rate(i__),e_M2_2rate(i__),M2_3rate(i__),
     +  e_M2_3rate(i__),M2_4rate(i__),e_M2_4rate(i__),M2_5rate(i__),
     +  e_M2_5rate(i__),M2_8rate(i__),e_M2_8rate(i__),M2_9rate(i__),
     +  e_M2_9rate(i__),EP8cts(i__),e_EP8cts(i__),PN8cts(i__),
     +  e_PN8cts(i__),M1_8cts(i__),e_M1_8cts(i__),M2_8cts(i__),
     +  e_M2_8cts(i__),EP8detML(i__),EP9detML(i__),PN1detML(i__),
     +  PN2detML(i__),PN3detML(i__),PN4detML(i__),PN5detML(i__),
     +  PN8detML(i__),PN9detML(i__),M1_1detML(i__),M1_2detML(i__),
     +  M1_3detML(i__),M1_4detML(i__),M1_5detML(i__),M1_8detML(i__),
     +  M1_9detML(i__),M2_1detML(i__),M2_2detML(i__),M2_3detML(i__),
     +  M2_4detML(i__),M2_5detML(i__),M2_8detML(i__),M2_9detML(i__),
     +  HREP1(i__),e_HREP1(i__),HREP2(i__),e_HREP2(i__),HREP3(i__),
     +  e_HREP3(i__),HREP4(i__),e_HREP4(i__),HRPN1(i__),e_HRPN1(i__),
     +  HRPN2(i__),e_HRPN2(i__),HRPN3(i__),e_HRPN3(i__),HRPN4(i__),
     +  e_HRPN4(i__),HRM1_1(i__),e_HRM1_1(i__),HRM1_2(i__),
     +  e_HRM1_2(i__),HRM1_3(i__),e_HRM1_3(i__),HRM1_4(i__),
     +  e_HRM1_4(i__),HRM2_1(i__),e_HRM2_1(i__),HRM2_2(i__),
     +  e_HRM2_2(i__),HRM2_3(i__),e_HRM2_3(i__),HRM2_4(i__),
     +  e_HRM2_4(i__),PN1ExpT(i__),PN2ExpT(i__),PN3ExpT(i__),
     +  PN4ExpT(i__),PN5ExpT(i__),M1_1ExpT(i__),M1_2ExpT(i__),
     +  M1_3ExpT(i__),M1_4ExpT(i__),M1_5ExpT(i__),M2_1ExpT(i__),
     +  M2_2ExpT(i__),M2_3ExpT(i__),M2_4ExpT(i__),M2_5ExpT(i__),
     +  PN1BG(i__),PN2BG(i__),PN3BG(i__),PN4BG(i__),PN5BG(i__),
     +  M1_1BG(i__),M1_2BG(i__),M1_3BG(i__),M1_4BG(i__),M1_5BG(i__),
     +  M2_1BG(i__),M2_2BG(i__),M2_3BG(i__),M2_4BG(i__),M2_5BG(i__),
     +  PN1Vig(i__),PN2Vig(i__),PN3Vig(i__),PN4Vig(i__),PN5Vig(i__),
     +  M1_1Vig(i__),M1_2Vig(i__),M1_3Vig(i__),M1_4Vig(i__),
     +  M1_5Vig(i__),M2_1Vig(i__),M2_2Vig(i__),M2_3Vig(i__),
     +  M2_4Vig(i__),M2_5Vig(i__),PNonTime(i__),M1onTime(i__),
     +  M2onTime(i__),PNoffax(i__),M1offax(i__),M2offax(i__),
     +  PNMaskFr(i__),M1MaskFr(i__),M2MaskFr(i__),DistNN(i__),
     +  SumFlag(i__),EPFlag(i__),PNFlag(i__),M1Flag(i__),M2Flag(i__),
     +  Tseries(i__),spectra(i__),EPchi2(i__),PNchi2(i__),M1chi2(i__),
     +  M2chi2(i__),VarFlag(i__),VarExpID(i__),VarInstID(i__),
     +  RAsdeg(i__),DEsdeg(i__),ePoss(i__),Flux1(i__),e_Flux1(i__),
     +  Flux2(i__),e_Flux2(i__),Flux3(i__),e_Flux3(i__),Flux4(i__),
     +  e_Flux4(i__),Flux5(i__),e_Flux5(i__),Flux8(i__),e_Flux8(i__),
     +  Flux9(i__),e_Flux9(i__),HR1(i__),e_HR1(i__),HR2(i__),
     +  e_HR2(i__),HR3(i__),e_HR3(i__),HR4(i__),e_HR4(i__),srcML(i__),
     +  extML(i__),Cst(i__),V(i__),S(i__),Ndet(i__),Sep(i__),
     +  LumDet(i__),e_LumDet(i__),CP02(i__),LB05(i__),SW04(i__),
     +  LM05(i__),LM05O(i__),New_(i__),Comp(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end