FORTRAN Generation
(/./ftp/cats/J/A_A/517/A44)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/517/A44 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-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/517/A44       Spitzer sources in the Galactic bulge   (Uttenthaler+, 2010)
*================================================================================
*Galactic bulge giants: probing stellar and galactic evolution.
*I. Catalogue of Spitzer IRAC and MIPS sources.
*    Uttenthaler S., Stute M., Sahai R., Blommaert J.A.D.L., Schultheis M.,
*    Kraemer K.E., Groenewegen M.A.T., Price S.D.
*   <Astron. Astrophys. 517, A44 (2010)>
*   =2010A&A...517A..44U
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'fits.dat'	! Summary of FITS images

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

C  J2000 position composed of: RAdeg DEdeg
      character*11  FileName   (nr__) ! Name of FITS file, in "fig" subdirectory
      character*1   Ch         (nr__) ! [1324M] Channel
      integer*4     Xpix       (nr__) ! (pix) X image size
      integer*4     Ypix       (nr__) ! (pix) Y image size
      real*8        RAdeg      (nr__) ! (deg) Right Ascension (J2000) of image centre
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000) of image centre

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

C  Declarations for 'b1.dat'	! Catalogue of point sources in the Bulge 1 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB        (nr__1) ! Source identification (from position)
      real*8        RAdeg_1    (nr__1) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1      (nr__1) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2      (nr__1) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3      (nr__1) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4      (nr__1) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS       (nr__1) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1    (nr__1) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2    (nr__1) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3    (nr__1) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4    (nr__1) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS     (nr__1) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1    (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2    (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3    (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4    (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS     (nr__1) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1    (nr__1) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2    (nr__1) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3    (nr__1) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4    (nr__1) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS     (nr__1) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN       (nr__1) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag       (nr__1) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag       (nr__1) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag       (nr__1) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M        (nr__1) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag      (nr__1) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag      (nr__1) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag      (nr__1) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag    (nr__1) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag    (nr__1) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag    (nr__1) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS      (nr__1) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12        (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12      (nr__1) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25        (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25      (nr__1) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60        (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60      (nr__1) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100       (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100     (nr__1) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX       (nr__1) ! (arcsec) ? Distance to closest MSX source
      real*4        B1         (nr__1) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1       (nr__1) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2         (nr__1) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2       (nr__1) ! [0/4]?=- MSX B2 band quality flag
      real*4        A          (nr__1) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A        (nr__1) ! [0/4]?=- MSX A band quality flag
      real*4        C          (nr__1) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C        (nr__1) ! [0/4]?=- MSX C band quality flag
      real*4        D          (nr__1) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D        (nr__1) ! [0/4]?=- MSX D band quality flag
      real*4        E          (nr__1) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E        (nr__1) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'b2.dat'	! Catalogue of point sources in the Bulge 2 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_1      (nr__2) ! Source identification (from position)
      real*8        RAdeg_2    (nr__2) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_1    (nr__2) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_1    (nr__2) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_1    (nr__2) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_1    (nr__2) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_1     (nr__2) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_1  (nr__2) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_1  (nr__2) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_1  (nr__2) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_1  (nr__2) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_1   (nr__2) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_1  (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_1  (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_1  (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_1  (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_1   (nr__2) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_1  (nr__2) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_1  (nr__2) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_1  (nr__2) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_1  (nr__2) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_1   (nr__2) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_1     (nr__2) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_1     (nr__2) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_1     (nr__2) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_1     (nr__2) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_1      (nr__2) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_1    (nr__2) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_1    (nr__2) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_1    (nr__2) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_1  (nr__2) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_1  (nr__2) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_1  (nr__2) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_1    (nr__2) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_1      (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_1    (nr__2) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_1      (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_1    (nr__2) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_1      (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_1    (nr__2) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_1     (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_1   (nr__2) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_1     (nr__2) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_1       (nr__2) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_1     (nr__2) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_1       (nr__2) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_1     (nr__2) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_1        (nr__2) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_1      (nr__2) ! [0/4]?=- MSX A band quality flag
      real*4        C_1        (nr__2) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_1      (nr__2) ! [0/4]?=- MSX C band quality flag
      real*4        D_1        (nr__2) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_1      (nr__2) ! [0/4]?=- MSX D band quality flag
      real*4        E_1        (nr__2) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_1      (nr__2) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'b3.dat'	! Catalogue of point sources in the Bulge 3 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_2      (nr__3) ! Source identification (from position)
      real*8        RAdeg_3    (nr__3) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_2    (nr__3) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_2    (nr__3) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_2    (nr__3) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_2    (nr__3) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_2     (nr__3) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_2  (nr__3) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_2  (nr__3) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_2  (nr__3) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_2  (nr__3) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_2   (nr__3) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_2  (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_2  (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_2  (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_2  (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_2   (nr__3) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_2  (nr__3) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_2  (nr__3) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_2  (nr__3) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_2  (nr__3) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_2   (nr__3) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_2     (nr__3) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_2     (nr__3) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_2     (nr__3) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_2     (nr__3) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_2      (nr__3) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_2    (nr__3) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_2    (nr__3) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_2    (nr__3) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_2  (nr__3) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_2  (nr__3) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_2  (nr__3) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_2    (nr__3) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_2      (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_2    (nr__3) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_2      (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_2    (nr__3) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_2      (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_2    (nr__3) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_2     (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_2   (nr__3) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_2     (nr__3) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_2       (nr__3) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_2     (nr__3) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_2       (nr__3) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_2     (nr__3) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_2        (nr__3) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_2      (nr__3) ! [0/4]?=- MSX A band quality flag
      real*4        C_2        (nr__3) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_2      (nr__3) ! [0/4]?=- MSX C band quality flag
      real*4        D_2        (nr__3) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_2      (nr__3) ! [0/4]?=- MSX D band quality flag
      real*4        E_2        (nr__3) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_2      (nr__3) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'b4.dat'	! Catalogue of point sources in the Bulge 4 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_3      (nr__4) ! Source identification (from position)
      real*8        RAdeg_4    (nr__4) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_4    (nr__4) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_3    (nr__4) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_3    (nr__4) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_3    (nr__4) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_3    (nr__4) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_3     (nr__4) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_3  (nr__4) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_3  (nr__4) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_3  (nr__4) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_3  (nr__4) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_3   (nr__4) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_3  (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_3  (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_3  (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_3  (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_3   (nr__4) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_3  (nr__4) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_3  (nr__4) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_3  (nr__4) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_3  (nr__4) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_3   (nr__4) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_3     (nr__4) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_3     (nr__4) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_3     (nr__4) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_3     (nr__4) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_3      (nr__4) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_3    (nr__4) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_3    (nr__4) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_3    (nr__4) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_3  (nr__4) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_3  (nr__4) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_3  (nr__4) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_3    (nr__4) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_3      (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_3    (nr__4) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_3      (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_3    (nr__4) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_3      (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_3    (nr__4) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_3     (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_3   (nr__4) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_3     (nr__4) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_3       (nr__4) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_3     (nr__4) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_3       (nr__4) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_3     (nr__4) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_3        (nr__4) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_3      (nr__4) ! [0/4]?=- MSX A band quality flag
      real*4        C_3        (nr__4) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_3      (nr__4) ! [0/4]?=- MSX C band quality flag
      real*4        D_3        (nr__4) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_3      (nr__4) ! [0/4]?=- MSX D band quality flag
      real*4        E_3        (nr__4) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_3      (nr__4) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'bn1.dat'	! Catalogue of point sources in the Bulge N1 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_4      (nr__5) ! Source identification (from position)
      real*8        RAdeg_5    (nr__5) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_5    (nr__5) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_4    (nr__5) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_4    (nr__5) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_4    (nr__5) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_4    (nr__5) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_4     (nr__5) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_4  (nr__5) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_4  (nr__5) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_4  (nr__5) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_4  (nr__5) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_4   (nr__5) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_4  (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_4  (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_4  (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_4  (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_4   (nr__5) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_4  (nr__5) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_4  (nr__5) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_4  (nr__5) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_4  (nr__5) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_4   (nr__5) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_4     (nr__5) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_4     (nr__5) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_4     (nr__5) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_4     (nr__5) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_4      (nr__5) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_4    (nr__5) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_4    (nr__5) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_4    (nr__5) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_4  (nr__5) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_4  (nr__5) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_4  (nr__5) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_4    (nr__5) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_4      (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_4    (nr__5) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_4      (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_4    (nr__5) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_4      (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_4    (nr__5) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_4     (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_4   (nr__5) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_4     (nr__5) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_4       (nr__5) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_4     (nr__5) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_4       (nr__5) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_4     (nr__5) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_4        (nr__5) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_4      (nr__5) ! [0/4]?=- MSX A band quality flag
      real*4        C_4        (nr__5) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_4      (nr__5) ! [0/4]?=- MSX C band quality flag
      real*4        D_4        (nr__5) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_4      (nr__5) ! [0/4]?=- MSX D band quality flag
      real*4        E_4        (nr__5) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_4      (nr__5) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'bn2.dat'	! Catalogue of point sources in the Bulge N2 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_5      (nr__6) ! Source identification (from position)
      real*8        RAdeg_6    (nr__6) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_6    (nr__6) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_5    (nr__6) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_5    (nr__6) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_5    (nr__6) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_5    (nr__6) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_5     (nr__6) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_5  (nr__6) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_5  (nr__6) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_5  (nr__6) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_5  (nr__6) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_5   (nr__6) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_5  (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_5  (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_5  (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_5  (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_5   (nr__6) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_5  (nr__6) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_5  (nr__6) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_5  (nr__6) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_5  (nr__6) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_5   (nr__6) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_5     (nr__6) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_5     (nr__6) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_5     (nr__6) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_5     (nr__6) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_5      (nr__6) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_5    (nr__6) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_5    (nr__6) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_5    (nr__6) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_5  (nr__6) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_5  (nr__6) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_5  (nr__6) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_5    (nr__6) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_5      (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_5    (nr__6) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_5      (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_5    (nr__6) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_5      (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_5    (nr__6) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_5     (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_5   (nr__6) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_5     (nr__6) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_5       (nr__6) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_5     (nr__6) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_5       (nr__6) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_5     (nr__6) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_5        (nr__6) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_5      (nr__6) ! [0/4]?=- MSX A band quality flag
      real*4        C_5        (nr__6) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_5      (nr__6) ! [0/4]?=- MSX C band quality flag
      real*4        D_5        (nr__6) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_5      (nr__6) ! [0/4]?=- MSX D band quality flag
      real*4        E_5        (nr__6) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_5      (nr__6) ! [0/4]?=- MSX E band quality flag

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

C  Declarations for 'bngc6522.dat'	! Catalogue of point sources in the NGC 6522 field

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  USB_6      (nr__7) ! Source identification (from position)
      real*8        RAdeg_7    (nr__7) ! (deg) Right ascension in degrees (J2000)
      real*8        DEdeg_7    (nr__7) ! (deg) Declination in degrees (J2000)
      real*4        IRAC1_6    (nr__7) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux
      real*4        IRAC2_6    (nr__7) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux
      real*4        IRAC3_6    (nr__7) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux
      real*4        IRAC4_6    (nr__7) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux
      real*4        MIPS_6     (nr__7) ! (uJy) ?=-9.99 MIPS 24 (24um) flux
      character*1   q_IRAC1_6  (nr__7) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2)
      character*1   q_IRAC2_6  (nr__7) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2)
      character*1   q_IRAC3_6  (nr__7) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2)
      character*1   q_IRAC4_6  (nr__7) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2)
      character*1   q_MIPS_6   (nr__7) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2)
      real*4        e_IRAC1_6  (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux
      real*4        e_IRAC2_6  (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux
      real*4        e_IRAC3_6  (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux
      real*4        e_IRAC4_6  (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux
      real*4        e_MIPS_6   (nr__7) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux
      integer*4     f_IRAC1_6  (nr__7) ! [0/1] Outside/Inside IRAC 1 field of view
      integer*4     f_IRAC2_6  (nr__7) ! [0/1] Outside/Inside IRAC 2 field of view
      integer*4     f_IRAC3_6  (nr__7) ! [0/1] Outside/Inside IRAC 3 field of view
      integer*4     f_IRAC4_6  (nr__7) ! [0/1] Outside/Inside IRAC 4 field of view
      integer*4     f_MIPS_6   (nr__7) ! [0/1] Outside/Inside MIPS 24 field of view
      real*4        dDEN_6     (nr__7) ! (arcsec) ? Distance to closest DENIS source
      real*4        Imag_6     (nr__7) ! (mag) ? I-band magnitude of DENIS counterpart
      real*4        Jmag_6     (nr__7) ! (mag) ? J-band magnitude of DENIS counterpart
      real*4        Kmag_6     (nr__7) ! (mag) ? K-band magnitude of DENIS counterpart
      real*4        d2M_6      (nr__7) ! (arcsec) ? Distance to closest 2MASS source
      real*4        J2mag_6    (nr__7) ! (mag) ?=- (?=-9.999) J-band magnitude of
*                                   2MASS counterpart
      real*4        H2mag_6    (nr__7) ! (mag) ?=- (?=-9.999) H-band magnitude of
*                                   2MASS counterpart
      real*4        K2mag_6    (nr__7) ! (mag) ?=- (?=-9.999) K-band magnitude of
*                                   2MASS counterpart
      character*1   q_J2mag_6  (nr__7) ! [ABCDEFU-] 2MASS J-band quality flag
      character*1   q_H2mag_6  (nr__7) ! [ABCDEFU-] 2MASS H-band quality flag
      character*1   q_K2mag_6  (nr__7) ! [ABCDEFU-] 2MASS K-band quality flag
      real*4        dIRAS_6    (nr__7) ! (arcsec) ? Distance to closest IRAS source
      real*4        F12_6      (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux
      integer*4     q_F12_6    (nr__7) ! [1/3]?=- IRAS 12 micron quality flag
      real*4        F25_6      (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux
      integer*4     q_F25_6    (nr__7) ! [1/3]?=- IRAS 25 micron quality flag
      real*4        F60_6      (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux
      integer*4     q_F60_6    (nr__7) ! [1/3]?=- IRAS 60 micron quality flag
      real*4        F100_6     (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux
      integer*4     q_F100_6   (nr__7) ! [1/3]?=- IRAS 100 micron quality flag
      real*4        dMSX_6     (nr__7) ! (arcsec) ? Distance to closest MSX source
      real*4        B1_6       (nr__7) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux
      integer*4     q_B1_6     (nr__7) ! [0/4]?=- MSX B1 band quality flag
      real*4        B2_6       (nr__7) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux
      integer*4     q_B2_6     (nr__7) ! [0/4]?=- MSX B2 band quality flag
      real*4        A_6        (nr__7) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux
      integer*4     q_A_6      (nr__7) ! [0/4]?=- MSX A band quality flag
      real*4        C_6        (nr__7) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux
      integer*4     q_C_6      (nr__7) ! [0/4]?=- MSX C band quality flag
      real*4        D_6        (nr__7) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux
      integer*4     q_D_6      (nr__7) ! [0/4]?=- MSX D band quality flag
      real*4        E_6        (nr__7) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux
      integer*4     q_E_6      (nr__7) ! [0/4]?=- MSX E band quality flag

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

C  Loading file 'fits.dat'	! Summary of FITS images

C  Format for file interpretation

    1 format(A11,3X,A1,2X,I4,1X,I4,2X,F9.5,F9.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fits.dat')
      write(6,*) '....Loading file: fits.dat'
      do i__=1,31
        read(1,'(A46)')ar__
        read(ar__,1)
     +  FileName(i__),Ch(i__),Xpix(i__),Ypix(i__),RAdeg(i__),
     +  DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  FileName(i__),Ch(i__),Xpix(i__),Ypix(i__),RAdeg(i__),
     +  DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'b1.dat'	! Catalogue of point sources in the Bulge 1 field

C  Format for file interpretation

    2 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b1.dat')
      write(6,*) '....Loading file: b1.dat'
      do i__=1,69047
        read(1,'(A466)')ar__1
        read(ar__1,2)
     +  USB(i__),RAdeg_1(i__),DEdeg_1(i__),IRAC1(i__),IRAC2(i__),
     +  IRAC3(i__),IRAC4(i__),MIPS(i__),q_IRAC1(i__),q_IRAC2(i__),
     +  q_IRAC3(i__),q_IRAC4(i__),q_MIPS(i__),e_IRAC1(i__),
     +  e_IRAC2(i__),e_IRAC3(i__),e_IRAC4(i__),e_MIPS(i__),
     +  f_IRAC1(i__),f_IRAC2(i__),f_IRAC3(i__),f_IRAC4(i__),
     +  f_MIPS(i__),dDEN(i__),Imag(i__),Jmag(i__),Kmag(i__),d2M(i__),
     +  J2mag(i__),H2mag(i__),K2mag(i__),q_J2mag(i__),q_H2mag(i__),
     +  q_K2mag(i__),dIRAS(i__),F12(i__),q_F12(i__),F25(i__),
     +  q_F25(i__),F60(i__),q_F60(i__),F100(i__),q_F100(i__),
     +  dMSX(i__),B1(i__),q_B1(i__),B2(i__),q_B2(i__),A(i__),q_A(i__),
     +  C(i__),q_C(i__),D(i__),q_D(i__),E(i__),q_E(i__)
        if(ar__1(217:221) .EQ. '') dDEN(i__) = rNULL__
        if(ar__1(225:230) .EQ. '') Imag(i__) = rNULL__
        if(ar__1(234:239) .EQ. '') Jmag(i__) = rNULL__
        if(ar__1(243:248) .EQ. '') Kmag(i__) = rNULL__
        if(ar__1(252:256) .EQ. '') d2M(i__) = rNULL__
        if (idig(ar__1(260:265)).EQ.0) J2mag(i__) =  rNULL__
        if (idig(ar__1(269:274)).EQ.0) H2mag(i__) =  rNULL__
        if (idig(ar__1(278:283)).EQ.0) K2mag(i__) =  rNULL__
        if(ar__1(293:297) .EQ. '') dIRAS(i__) = rNULL__
        if (idig(ar__1(301:308)).EQ.0) F12(i__) =  rNULL__
        if (q_F12(i__) .EQ. 45) q_F12(i__) =  iNULL__
        if (idig(ar__1(316:323)).EQ.0) F25(i__) =  rNULL__
        if (q_F25(i__) .EQ. 45) q_F25(i__) =  iNULL__
        if (idig(ar__1(331:338)).EQ.0) F60(i__) =  rNULL__
        if (q_F60(i__) .EQ. 45) q_F60(i__) =  iNULL__
        if (idig(ar__1(346:353)).EQ.0) F100(i__) =  rNULL__
        if (q_F100(i__) .EQ. 45) q_F100(i__) =  iNULL__
        if(ar__1(361:364) .EQ. '') dMSX(i__) = rNULL__
        if (idig(ar__1(368:377)).EQ.0) B1(i__) =  rNULL__
        if (q_B1(i__) .EQ. 45) q_B1(i__) =  iNULL__
        if (idig(ar__1(385:394)).EQ.0) B2(i__) =  rNULL__
        if (q_B2(i__) .EQ. 45) q_B2(i__) =  iNULL__
        if (idig(ar__1(402:411)).EQ.0) A(i__) =  rNULL__
        if (q_A(i__) .EQ. 45) q_A(i__) =  iNULL__
        if (idig(ar__1(419:428)).EQ.0) C(i__) =  rNULL__
        if (q_C(i__) .EQ. 45) q_C(i__) =  iNULL__
        if (idig(ar__1(436:445)).EQ.0) D(i__) =  rNULL__
        if (q_D(i__) .EQ. 45) q_D(i__) =  iNULL__
        if (idig(ar__1(453:462)).EQ.0) E(i__) =  rNULL__
        if (q_E(i__) .EQ. 45) q_E(i__) =  iNULL__
c    ..............Just test output...........
        write(6,2)
     +  USB(i__),RAdeg_1(i__),DEdeg_1(i__),IRAC1(i__),IRAC2(i__),
     +  IRAC3(i__),IRAC4(i__),MIPS(i__),q_IRAC1(i__),q_IRAC2(i__),
     +  q_IRAC3(i__),q_IRAC4(i__),q_MIPS(i__),e_IRAC1(i__),
     +  e_IRAC2(i__),e_IRAC3(i__),e_IRAC4(i__),e_MIPS(i__),
     +  f_IRAC1(i__),f_IRAC2(i__),f_IRAC3(i__),f_IRAC4(i__),
     +  f_MIPS(i__),dDEN(i__),Imag(i__),Jmag(i__),Kmag(i__),d2M(i__),
     +  J2mag(i__),H2mag(i__),K2mag(i__),q_J2mag(i__),q_H2mag(i__),
     +  q_K2mag(i__),dIRAS(i__),F12(i__),q_F12(i__),F25(i__),
     +  q_F25(i__),F60(i__),q_F60(i__),F100(i__),q_F100(i__),
     +  dMSX(i__),B1(i__),q_B1(i__),B2(i__),q_B2(i__),A(i__),q_A(i__),
     +  C(i__),q_C(i__),D(i__),q_D(i__),E(i__),q_E(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'b2.dat'	! Catalogue of point sources in the Bulge 2 field

C  Format for file interpretation

    3 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b2.dat')
      write(6,*) '....Loading file: b2.dat'
      do i__=1,60262
        read(1,'(A466)')ar__2
        read(ar__2,3)
     +  USB_1(i__),RAdeg_2(i__),DEdeg_2(i__),IRAC1_1(i__),
     +  IRAC2_1(i__),IRAC3_1(i__),IRAC4_1(i__),MIPS_1(i__),
     +  q_IRAC1_1(i__),q_IRAC2_1(i__),q_IRAC3_1(i__),q_IRAC4_1(i__),
     +  q_MIPS_1(i__),e_IRAC1_1(i__),e_IRAC2_1(i__),e_IRAC3_1(i__),
     +  e_IRAC4_1(i__),e_MIPS_1(i__),f_IRAC1_1(i__),f_IRAC2_1(i__),
     +  f_IRAC3_1(i__),f_IRAC4_1(i__),f_MIPS_1(i__),dDEN_1(i__),
     +  Imag_1(i__),Jmag_1(i__),Kmag_1(i__),d2M_1(i__),J2mag_1(i__),
     +  H2mag_1(i__),K2mag_1(i__),q_J2mag_1(i__),q_H2mag_1(i__),
     +  q_K2mag_1(i__),dIRAS_1(i__),F12_1(i__),q_F12_1(i__),
     +  F25_1(i__),q_F25_1(i__),F60_1(i__),q_F60_1(i__),F100_1(i__),
     +  q_F100_1(i__),dMSX_1(i__),B1_1(i__),q_B1_1(i__),B2_1(i__),
     +  q_B2_1(i__),A_1(i__),q_A_1(i__),C_1(i__),q_C_1(i__),D_1(i__),
     +  q_D_1(i__),E_1(i__),q_E_1(i__)
        if(ar__2(217:221) .EQ. '') dDEN_1(i__) = rNULL__
        if(ar__2(225:230) .EQ. '') Imag_1(i__) = rNULL__
        if(ar__2(234:239) .EQ. '') Jmag_1(i__) = rNULL__
        if(ar__2(243:248) .EQ. '') Kmag_1(i__) = rNULL__
        if(ar__2(252:256) .EQ. '') d2M_1(i__) = rNULL__
        if (idig(ar__2(260:265)).EQ.0) J2mag_1(i__) =  rNULL__
        if (idig(ar__2(269:274)).EQ.0) H2mag_1(i__) =  rNULL__
        if (idig(ar__2(278:283)).EQ.0) K2mag_1(i__) =  rNULL__
        if(ar__2(293:297) .EQ. '') dIRAS_1(i__) = rNULL__
        if (idig(ar__2(301:308)).EQ.0) F12_1(i__) =  rNULL__
        if (q_F12_1(i__) .EQ. 45) q_F12_1(i__) =  iNULL__
        if (idig(ar__2(316:323)).EQ.0) F25_1(i__) =  rNULL__
        if (q_F25_1(i__) .EQ. 45) q_F25_1(i__) =  iNULL__
        if (idig(ar__2(331:338)).EQ.0) F60_1(i__) =  rNULL__
        if (q_F60_1(i__) .EQ. 45) q_F60_1(i__) =  iNULL__
        if (idig(ar__2(346:353)).EQ.0) F100_1(i__) =  rNULL__
        if (q_F100_1(i__) .EQ. 45) q_F100_1(i__) =  iNULL__
        if(ar__2(361:364) .EQ. '') dMSX_1(i__) = rNULL__
        if (idig(ar__2(368:377)).EQ.0) B1_1(i__) =  rNULL__
        if (q_B1_1(i__) .EQ. 45) q_B1_1(i__) =  iNULL__
        if (idig(ar__2(385:394)).EQ.0) B2_1(i__) =  rNULL__
        if (q_B2_1(i__) .EQ. 45) q_B2_1(i__) =  iNULL__
        if (idig(ar__2(402:411)).EQ.0) A_1(i__) =  rNULL__
        if (q_A_1(i__) .EQ. 45) q_A_1(i__) =  iNULL__
        if (idig(ar__2(419:428)).EQ.0) C_1(i__) =  rNULL__
        if (q_C_1(i__) .EQ. 45) q_C_1(i__) =  iNULL__
        if (idig(ar__2(436:445)).EQ.0) D_1(i__) =  rNULL__
        if (q_D_1(i__) .EQ. 45) q_D_1(i__) =  iNULL__
        if (idig(ar__2(453:462)).EQ.0) E_1(i__) =  rNULL__
        if (q_E_1(i__) .EQ. 45) q_E_1(i__) =  iNULL__
c    ..............Just test output...........
        write(6,3)
     +  USB_1(i__),RAdeg_2(i__),DEdeg_2(i__),IRAC1_1(i__),
     +  IRAC2_1(i__),IRAC3_1(i__),IRAC4_1(i__),MIPS_1(i__),
     +  q_IRAC1_1(i__),q_IRAC2_1(i__),q_IRAC3_1(i__),q_IRAC4_1(i__),
     +  q_MIPS_1(i__),e_IRAC1_1(i__),e_IRAC2_1(i__),e_IRAC3_1(i__),
     +  e_IRAC4_1(i__),e_MIPS_1(i__),f_IRAC1_1(i__),f_IRAC2_1(i__),
     +  f_IRAC3_1(i__),f_IRAC4_1(i__),f_MIPS_1(i__),dDEN_1(i__),
     +  Imag_1(i__),Jmag_1(i__),Kmag_1(i__),d2M_1(i__),J2mag_1(i__),
     +  H2mag_1(i__),K2mag_1(i__),q_J2mag_1(i__),q_H2mag_1(i__),
     +  q_K2mag_1(i__),dIRAS_1(i__),F12_1(i__),q_F12_1(i__),
     +  F25_1(i__),q_F25_1(i__),F60_1(i__),q_F60_1(i__),F100_1(i__),
     +  q_F100_1(i__),dMSX_1(i__),B1_1(i__),q_B1_1(i__),B2_1(i__),
     +  q_B2_1(i__),A_1(i__),q_A_1(i__),C_1(i__),q_C_1(i__),D_1(i__),
     +  q_D_1(i__),E_1(i__),q_E_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'b3.dat'	! Catalogue of point sources in the Bulge 3 field

C  Format for file interpretation

    4 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b3.dat')
      write(6,*) '....Loading file: b3.dat'
      do i__=1,52576
        read(1,'(A466)')ar__3
        read(ar__3,4)
     +  USB_2(i__),RAdeg_3(i__),DEdeg_3(i__),IRAC1_2(i__),
     +  IRAC2_2(i__),IRAC3_2(i__),IRAC4_2(i__),MIPS_2(i__),
     +  q_IRAC1_2(i__),q_IRAC2_2(i__),q_IRAC3_2(i__),q_IRAC4_2(i__),
     +  q_MIPS_2(i__),e_IRAC1_2(i__),e_IRAC2_2(i__),e_IRAC3_2(i__),
     +  e_IRAC4_2(i__),e_MIPS_2(i__),f_IRAC1_2(i__),f_IRAC2_2(i__),
     +  f_IRAC3_2(i__),f_IRAC4_2(i__),f_MIPS_2(i__),dDEN_2(i__),
     +  Imag_2(i__),Jmag_2(i__),Kmag_2(i__),d2M_2(i__),J2mag_2(i__),
     +  H2mag_2(i__),K2mag_2(i__),q_J2mag_2(i__),q_H2mag_2(i__),
     +  q_K2mag_2(i__),dIRAS_2(i__),F12_2(i__),q_F12_2(i__),
     +  F25_2(i__),q_F25_2(i__),F60_2(i__),q_F60_2(i__),F100_2(i__),
     +  q_F100_2(i__),dMSX_2(i__),B1_2(i__),q_B1_2(i__),B2_2(i__),
     +  q_B2_2(i__),A_2(i__),q_A_2(i__),C_2(i__),q_C_2(i__),D_2(i__),
     +  q_D_2(i__),E_2(i__),q_E_2(i__)
        if(ar__3(217:221) .EQ. '') dDEN_2(i__) = rNULL__
        if(ar__3(225:230) .EQ. '') Imag_2(i__) = rNULL__
        if(ar__3(234:239) .EQ. '') Jmag_2(i__) = rNULL__
        if(ar__3(243:248) .EQ. '') Kmag_2(i__) = rNULL__
        if(ar__3(252:256) .EQ. '') d2M_2(i__) = rNULL__
        if (idig(ar__3(260:265)).EQ.0) J2mag_2(i__) =  rNULL__
        if (idig(ar__3(269:274)).EQ.0) H2mag_2(i__) =  rNULL__
        if (idig(ar__3(278:283)).EQ.0) K2mag_2(i__) =  rNULL__
        if(ar__3(293:297) .EQ. '') dIRAS_2(i__) = rNULL__
        if (idig(ar__3(301:308)).EQ.0) F12_2(i__) =  rNULL__
        if (q_F12_2(i__) .EQ. 45) q_F12_2(i__) =  iNULL__
        if (idig(ar__3(316:323)).EQ.0) F25_2(i__) =  rNULL__
        if (q_F25_2(i__) .EQ. 45) q_F25_2(i__) =  iNULL__
        if (idig(ar__3(331:338)).EQ.0) F60_2(i__) =  rNULL__
        if (q_F60_2(i__) .EQ. 45) q_F60_2(i__) =  iNULL__
        if (idig(ar__3(346:353)).EQ.0) F100_2(i__) =  rNULL__
        if (q_F100_2(i__) .EQ. 45) q_F100_2(i__) =  iNULL__
        if(ar__3(361:364) .EQ. '') dMSX_2(i__) = rNULL__
        if (idig(ar__3(368:377)).EQ.0) B1_2(i__) =  rNULL__
        if (q_B1_2(i__) .EQ. 45) q_B1_2(i__) =  iNULL__
        if (idig(ar__3(385:394)).EQ.0) B2_2(i__) =  rNULL__
        if (q_B2_2(i__) .EQ. 45) q_B2_2(i__) =  iNULL__
        if (idig(ar__3(402:411)).EQ.0) A_2(i__) =  rNULL__
        if (q_A_2(i__) .EQ. 45) q_A_2(i__) =  iNULL__
        if (idig(ar__3(419:428)).EQ.0) C_2(i__) =  rNULL__
        if (q_C_2(i__) .EQ. 45) q_C_2(i__) =  iNULL__
        if (idig(ar__3(436:445)).EQ.0) D_2(i__) =  rNULL__
        if (q_D_2(i__) .EQ. 45) q_D_2(i__) =  iNULL__
        if (idig(ar__3(453:462)).EQ.0) E_2(i__) =  rNULL__
        if (q_E_2(i__) .EQ. 45) q_E_2(i__) =  iNULL__
c    ..............Just test output...........
        write(6,4)
     +  USB_2(i__),RAdeg_3(i__),DEdeg_3(i__),IRAC1_2(i__),
     +  IRAC2_2(i__),IRAC3_2(i__),IRAC4_2(i__),MIPS_2(i__),
     +  q_IRAC1_2(i__),q_IRAC2_2(i__),q_IRAC3_2(i__),q_IRAC4_2(i__),
     +  q_MIPS_2(i__),e_IRAC1_2(i__),e_IRAC2_2(i__),e_IRAC3_2(i__),
     +  e_IRAC4_2(i__),e_MIPS_2(i__),f_IRAC1_2(i__),f_IRAC2_2(i__),
     +  f_IRAC3_2(i__),f_IRAC4_2(i__),f_MIPS_2(i__),dDEN_2(i__),
     +  Imag_2(i__),Jmag_2(i__),Kmag_2(i__),d2M_2(i__),J2mag_2(i__),
     +  H2mag_2(i__),K2mag_2(i__),q_J2mag_2(i__),q_H2mag_2(i__),
     +  q_K2mag_2(i__),dIRAS_2(i__),F12_2(i__),q_F12_2(i__),
     +  F25_2(i__),q_F25_2(i__),F60_2(i__),q_F60_2(i__),F100_2(i__),
     +  q_F100_2(i__),dMSX_2(i__),B1_2(i__),q_B1_2(i__),B2_2(i__),
     +  q_B2_2(i__),A_2(i__),q_A_2(i__),C_2(i__),q_C_2(i__),D_2(i__),
     +  q_D_2(i__),E_2(i__),q_E_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'b4.dat'	! Catalogue of point sources in the Bulge 4 field

C  Format for file interpretation

    5 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'b4.dat')
      write(6,*) '....Loading file: b4.dat'
      do i__=1,58258
        read(1,'(A466)')ar__4
        read(ar__4,5)
     +  USB_3(i__),RAdeg_4(i__),DEdeg_4(i__),IRAC1_3(i__),
     +  IRAC2_3(i__),IRAC3_3(i__),IRAC4_3(i__),MIPS_3(i__),
     +  q_IRAC1_3(i__),q_IRAC2_3(i__),q_IRAC3_3(i__),q_IRAC4_3(i__),
     +  q_MIPS_3(i__),e_IRAC1_3(i__),e_IRAC2_3(i__),e_IRAC3_3(i__),
     +  e_IRAC4_3(i__),e_MIPS_3(i__),f_IRAC1_3(i__),f_IRAC2_3(i__),
     +  f_IRAC3_3(i__),f_IRAC4_3(i__),f_MIPS_3(i__),dDEN_3(i__),
     +  Imag_3(i__),Jmag_3(i__),Kmag_3(i__),d2M_3(i__),J2mag_3(i__),
     +  H2mag_3(i__),K2mag_3(i__),q_J2mag_3(i__),q_H2mag_3(i__),
     +  q_K2mag_3(i__),dIRAS_3(i__),F12_3(i__),q_F12_3(i__),
     +  F25_3(i__),q_F25_3(i__),F60_3(i__),q_F60_3(i__),F100_3(i__),
     +  q_F100_3(i__),dMSX_3(i__),B1_3(i__),q_B1_3(i__),B2_3(i__),
     +  q_B2_3(i__),A_3(i__),q_A_3(i__),C_3(i__),q_C_3(i__),D_3(i__),
     +  q_D_3(i__),E_3(i__),q_E_3(i__)
        if(ar__4(217:221) .EQ. '') dDEN_3(i__) = rNULL__
        if(ar__4(225:230) .EQ. '') Imag_3(i__) = rNULL__
        if(ar__4(234:239) .EQ. '') Jmag_3(i__) = rNULL__
        if(ar__4(243:248) .EQ. '') Kmag_3(i__) = rNULL__
        if(ar__4(252:256) .EQ. '') d2M_3(i__) = rNULL__
        if (idig(ar__4(260:265)).EQ.0) J2mag_3(i__) =  rNULL__
        if (idig(ar__4(269:274)).EQ.0) H2mag_3(i__) =  rNULL__
        if (idig(ar__4(278:283)).EQ.0) K2mag_3(i__) =  rNULL__
        if(ar__4(293:297) .EQ. '') dIRAS_3(i__) = rNULL__
        if (idig(ar__4(301:308)).EQ.0) F12_3(i__) =  rNULL__
        if (q_F12_3(i__) .EQ. 45) q_F12_3(i__) =  iNULL__
        if (idig(ar__4(316:323)).EQ.0) F25_3(i__) =  rNULL__
        if (q_F25_3(i__) .EQ. 45) q_F25_3(i__) =  iNULL__
        if (idig(ar__4(331:338)).EQ.0) F60_3(i__) =  rNULL__
        if (q_F60_3(i__) .EQ. 45) q_F60_3(i__) =  iNULL__
        if (idig(ar__4(346:353)).EQ.0) F100_3(i__) =  rNULL__
        if (q_F100_3(i__) .EQ. 45) q_F100_3(i__) =  iNULL__
        if(ar__4(361:364) .EQ. '') dMSX_3(i__) = rNULL__
        if (idig(ar__4(368:377)).EQ.0) B1_3(i__) =  rNULL__
        if (q_B1_3(i__) .EQ. 45) q_B1_3(i__) =  iNULL__
        if (idig(ar__4(385:394)).EQ.0) B2_3(i__) =  rNULL__
        if (q_B2_3(i__) .EQ. 45) q_B2_3(i__) =  iNULL__
        if (idig(ar__4(402:411)).EQ.0) A_3(i__) =  rNULL__
        if (q_A_3(i__) .EQ. 45) q_A_3(i__) =  iNULL__
        if (idig(ar__4(419:428)).EQ.0) C_3(i__) =  rNULL__
        if (q_C_3(i__) .EQ. 45) q_C_3(i__) =  iNULL__
        if (idig(ar__4(436:445)).EQ.0) D_3(i__) =  rNULL__
        if (q_D_3(i__) .EQ. 45) q_D_3(i__) =  iNULL__
        if (idig(ar__4(453:462)).EQ.0) E_3(i__) =  rNULL__
        if (q_E_3(i__) .EQ. 45) q_E_3(i__) =  iNULL__
c    ..............Just test output...........
        write(6,5)
     +  USB_3(i__),RAdeg_4(i__),DEdeg_4(i__),IRAC1_3(i__),
     +  IRAC2_3(i__),IRAC3_3(i__),IRAC4_3(i__),MIPS_3(i__),
     +  q_IRAC1_3(i__),q_IRAC2_3(i__),q_IRAC3_3(i__),q_IRAC4_3(i__),
     +  q_MIPS_3(i__),e_IRAC1_3(i__),e_IRAC2_3(i__),e_IRAC3_3(i__),
     +  e_IRAC4_3(i__),e_MIPS_3(i__),f_IRAC1_3(i__),f_IRAC2_3(i__),
     +  f_IRAC3_3(i__),f_IRAC4_3(i__),f_MIPS_3(i__),dDEN_3(i__),
     +  Imag_3(i__),Jmag_3(i__),Kmag_3(i__),d2M_3(i__),J2mag_3(i__),
     +  H2mag_3(i__),K2mag_3(i__),q_J2mag_3(i__),q_H2mag_3(i__),
     +  q_K2mag_3(i__),dIRAS_3(i__),F12_3(i__),q_F12_3(i__),
     +  F25_3(i__),q_F25_3(i__),F60_3(i__),q_F60_3(i__),F100_3(i__),
     +  q_F100_3(i__),dMSX_3(i__),B1_3(i__),q_B1_3(i__),B2_3(i__),
     +  q_B2_3(i__),A_3(i__),q_A_3(i__),C_3(i__),q_C_3(i__),D_3(i__),
     +  q_D_3(i__),E_3(i__),q_E_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'bn1.dat'	! Catalogue of point sources in the Bulge N1 field

C  Format for file interpretation

    6 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'bn1.dat')
      write(6,*) '....Loading file: bn1.dat'
      do i__=1,64798
        read(1,'(A466)')ar__5
        read(ar__5,6)
     +  USB_4(i__),RAdeg_5(i__),DEdeg_5(i__),IRAC1_4(i__),
     +  IRAC2_4(i__),IRAC3_4(i__),IRAC4_4(i__),MIPS_4(i__),
     +  q_IRAC1_4(i__),q_IRAC2_4(i__),q_IRAC3_4(i__),q_IRAC4_4(i__),
     +  q_MIPS_4(i__),e_IRAC1_4(i__),e_IRAC2_4(i__),e_IRAC3_4(i__),
     +  e_IRAC4_4(i__),e_MIPS_4(i__),f_IRAC1_4(i__),f_IRAC2_4(i__),
     +  f_IRAC3_4(i__),f_IRAC4_4(i__),f_MIPS_4(i__),dDEN_4(i__),
     +  Imag_4(i__),Jmag_4(i__),Kmag_4(i__),d2M_4(i__),J2mag_4(i__),
     +  H2mag_4(i__),K2mag_4(i__),q_J2mag_4(i__),q_H2mag_4(i__),
     +  q_K2mag_4(i__),dIRAS_4(i__),F12_4(i__),q_F12_4(i__),
     +  F25_4(i__),q_F25_4(i__),F60_4(i__),q_F60_4(i__),F100_4(i__),
     +  q_F100_4(i__),dMSX_4(i__),B1_4(i__),q_B1_4(i__),B2_4(i__),
     +  q_B2_4(i__),A_4(i__),q_A_4(i__),C_4(i__),q_C_4(i__),D_4(i__),
     +  q_D_4(i__),E_4(i__),q_E_4(i__)
        if(ar__5(217:221) .EQ. '') dDEN_4(i__) = rNULL__
        if(ar__5(225:230) .EQ. '') Imag_4(i__) = rNULL__
        if(ar__5(234:239) .EQ. '') Jmag_4(i__) = rNULL__
        if(ar__5(243:248) .EQ. '') Kmag_4(i__) = rNULL__
        if(ar__5(252:256) .EQ. '') d2M_4(i__) = rNULL__
        if (idig(ar__5(260:265)).EQ.0) J2mag_4(i__) =  rNULL__
        if (idig(ar__5(269:274)).EQ.0) H2mag_4(i__) =  rNULL__
        if (idig(ar__5(278:283)).EQ.0) K2mag_4(i__) =  rNULL__
        if(ar__5(293:297) .EQ. '') dIRAS_4(i__) = rNULL__
        if (idig(ar__5(301:308)).EQ.0) F12_4(i__) =  rNULL__
        if (q_F12_4(i__) .EQ. 45) q_F12_4(i__) =  iNULL__
        if (idig(ar__5(316:323)).EQ.0) F25_4(i__) =  rNULL__
        if (q_F25_4(i__) .EQ. 45) q_F25_4(i__) =  iNULL__
        if (idig(ar__5(331:338)).EQ.0) F60_4(i__) =  rNULL__
        if (q_F60_4(i__) .EQ. 45) q_F60_4(i__) =  iNULL__
        if (idig(ar__5(346:353)).EQ.0) F100_4(i__) =  rNULL__
        if (q_F100_4(i__) .EQ. 45) q_F100_4(i__) =  iNULL__
        if(ar__5(361:364) .EQ. '') dMSX_4(i__) = rNULL__
        if (idig(ar__5(368:377)).EQ.0) B1_4(i__) =  rNULL__
        if (q_B1_4(i__) .EQ. 45) q_B1_4(i__) =  iNULL__
        if (idig(ar__5(385:394)).EQ.0) B2_4(i__) =  rNULL__
        if (q_B2_4(i__) .EQ. 45) q_B2_4(i__) =  iNULL__
        if (idig(ar__5(402:411)).EQ.0) A_4(i__) =  rNULL__
        if (q_A_4(i__) .EQ. 45) q_A_4(i__) =  iNULL__
        if (idig(ar__5(419:428)).EQ.0) C_4(i__) =  rNULL__
        if (q_C_4(i__) .EQ. 45) q_C_4(i__) =  iNULL__
        if (idig(ar__5(436:445)).EQ.0) D_4(i__) =  rNULL__
        if (q_D_4(i__) .EQ. 45) q_D_4(i__) =  iNULL__
        if (idig(ar__5(453:462)).EQ.0) E_4(i__) =  rNULL__
        if (q_E_4(i__) .EQ. 45) q_E_4(i__) =  iNULL__
c    ..............Just test output...........
        write(6,6)
     +  USB_4(i__),RAdeg_5(i__),DEdeg_5(i__),IRAC1_4(i__),
     +  IRAC2_4(i__),IRAC3_4(i__),IRAC4_4(i__),MIPS_4(i__),
     +  q_IRAC1_4(i__),q_IRAC2_4(i__),q_IRAC3_4(i__),q_IRAC4_4(i__),
     +  q_MIPS_4(i__),e_IRAC1_4(i__),e_IRAC2_4(i__),e_IRAC3_4(i__),
     +  e_IRAC4_4(i__),e_MIPS_4(i__),f_IRAC1_4(i__),f_IRAC2_4(i__),
     +  f_IRAC3_4(i__),f_IRAC4_4(i__),f_MIPS_4(i__),dDEN_4(i__),
     +  Imag_4(i__),Jmag_4(i__),Kmag_4(i__),d2M_4(i__),J2mag_4(i__),
     +  H2mag_4(i__),K2mag_4(i__),q_J2mag_4(i__),q_H2mag_4(i__),
     +  q_K2mag_4(i__),dIRAS_4(i__),F12_4(i__),q_F12_4(i__),
     +  F25_4(i__),q_F25_4(i__),F60_4(i__),q_F60_4(i__),F100_4(i__),
     +  q_F100_4(i__),dMSX_4(i__),B1_4(i__),q_B1_4(i__),B2_4(i__),
     +  q_B2_4(i__),A_4(i__),q_A_4(i__),C_4(i__),q_C_4(i__),D_4(i__),
     +  q_D_4(i__),E_4(i__),q_E_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'bn2.dat'	! Catalogue of point sources in the Bulge N2 field

C  Format for file interpretation

    7 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'bn2.dat')
      write(6,*) '....Loading file: bn2.dat'
      do i__=1,428
        read(1,'(A466)')ar__6
        read(ar__6,7)
     +  USB_5(i__),RAdeg_6(i__),DEdeg_6(i__),IRAC1_5(i__),
     +  IRAC2_5(i__),IRAC3_5(i__),IRAC4_5(i__),MIPS_5(i__),
     +  q_IRAC1_5(i__),q_IRAC2_5(i__),q_IRAC3_5(i__),q_IRAC4_5(i__),
     +  q_MIPS_5(i__),e_IRAC1_5(i__),e_IRAC2_5(i__),e_IRAC3_5(i__),
     +  e_IRAC4_5(i__),e_MIPS_5(i__),f_IRAC1_5(i__),f_IRAC2_5(i__),
     +  f_IRAC3_5(i__),f_IRAC4_5(i__),f_MIPS_5(i__),dDEN_5(i__),
     +  Imag_5(i__),Jmag_5(i__),Kmag_5(i__),d2M_5(i__),J2mag_5(i__),
     +  H2mag_5(i__),K2mag_5(i__),q_J2mag_5(i__),q_H2mag_5(i__),
     +  q_K2mag_5(i__),dIRAS_5(i__),F12_5(i__),q_F12_5(i__),
     +  F25_5(i__),q_F25_5(i__),F60_5(i__),q_F60_5(i__),F100_5(i__),
     +  q_F100_5(i__),dMSX_5(i__),B1_5(i__),q_B1_5(i__),B2_5(i__),
     +  q_B2_5(i__),A_5(i__),q_A_5(i__),C_5(i__),q_C_5(i__),D_5(i__),
     +  q_D_5(i__),E_5(i__),q_E_5(i__)
        if(ar__6(217:221) .EQ. '') dDEN_5(i__) = rNULL__
        if(ar__6(225:230) .EQ. '') Imag_5(i__) = rNULL__
        if(ar__6(234:239) .EQ. '') Jmag_5(i__) = rNULL__
        if(ar__6(243:248) .EQ. '') Kmag_5(i__) = rNULL__
        if(ar__6(252:256) .EQ. '') d2M_5(i__) = rNULL__
        if (idig(ar__6(260:265)).EQ.0) J2mag_5(i__) =  rNULL__
        if (idig(ar__6(269:274)).EQ.0) H2mag_5(i__) =  rNULL__
        if (idig(ar__6(278:283)).EQ.0) K2mag_5(i__) =  rNULL__
        if(ar__6(293:297) .EQ. '') dIRAS_5(i__) = rNULL__
        if (idig(ar__6(301:308)).EQ.0) F12_5(i__) =  rNULL__
        if (q_F12_5(i__) .EQ. 45) q_F12_5(i__) =  iNULL__
        if (idig(ar__6(316:323)).EQ.0) F25_5(i__) =  rNULL__
        if (q_F25_5(i__) .EQ. 45) q_F25_5(i__) =  iNULL__
        if (idig(ar__6(331:338)).EQ.0) F60_5(i__) =  rNULL__
        if (q_F60_5(i__) .EQ. 45) q_F60_5(i__) =  iNULL__
        if (idig(ar__6(346:353)).EQ.0) F100_5(i__) =  rNULL__
        if (q_F100_5(i__) .EQ. 45) q_F100_5(i__) =  iNULL__
        if(ar__6(361:364) .EQ. '') dMSX_5(i__) = rNULL__
        if (idig(ar__6(368:377)).EQ.0) B1_5(i__) =  rNULL__
        if (q_B1_5(i__) .EQ. 45) q_B1_5(i__) =  iNULL__
        if (idig(ar__6(385:394)).EQ.0) B2_5(i__) =  rNULL__
        if (q_B2_5(i__) .EQ. 45) q_B2_5(i__) =  iNULL__
        if (idig(ar__6(402:411)).EQ.0) A_5(i__) =  rNULL__
        if (q_A_5(i__) .EQ. 45) q_A_5(i__) =  iNULL__
        if (idig(ar__6(419:428)).EQ.0) C_5(i__) =  rNULL__
        if (q_C_5(i__) .EQ. 45) q_C_5(i__) =  iNULL__
        if (idig(ar__6(436:445)).EQ.0) D_5(i__) =  rNULL__
        if (q_D_5(i__) .EQ. 45) q_D_5(i__) =  iNULL__
        if (idig(ar__6(453:462)).EQ.0) E_5(i__) =  rNULL__
        if (q_E_5(i__) .EQ. 45) q_E_5(i__) =  iNULL__
c    ..............Just test output...........
        write(6,7)
     +  USB_5(i__),RAdeg_6(i__),DEdeg_6(i__),IRAC1_5(i__),
     +  IRAC2_5(i__),IRAC3_5(i__),IRAC4_5(i__),MIPS_5(i__),
     +  q_IRAC1_5(i__),q_IRAC2_5(i__),q_IRAC3_5(i__),q_IRAC4_5(i__),
     +  q_MIPS_5(i__),e_IRAC1_5(i__),e_IRAC2_5(i__),e_IRAC3_5(i__),
     +  e_IRAC4_5(i__),e_MIPS_5(i__),f_IRAC1_5(i__),f_IRAC2_5(i__),
     +  f_IRAC3_5(i__),f_IRAC4_5(i__),f_MIPS_5(i__),dDEN_5(i__),
     +  Imag_5(i__),Jmag_5(i__),Kmag_5(i__),d2M_5(i__),J2mag_5(i__),
     +  H2mag_5(i__),K2mag_5(i__),q_J2mag_5(i__),q_H2mag_5(i__),
     +  q_K2mag_5(i__),dIRAS_5(i__),F12_5(i__),q_F12_5(i__),
     +  F25_5(i__),q_F25_5(i__),F60_5(i__),q_F60_5(i__),F100_5(i__),
     +  q_F100_5(i__),dMSX_5(i__),B1_5(i__),q_B1_5(i__),B2_5(i__),
     +  q_B2_5(i__),A_5(i__),q_A_5(i__),C_5(i__),q_C_5(i__),D_5(i__),
     +  q_D_5(i__),E_5(i__),q_E_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'bngc6522.dat'	! Catalogue of point sources in the NGC 6522 field

C  Format for file interpretation

    8 format(
     +  3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6,
     +  3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X,
     +  E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X,
     +  F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X,
     +  E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,
     +  3X,E10.3,3X,I1,3X,E10.3,3X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'bngc6522.dat')
      write(6,*) '....Loading file: bngc6522.dat'
      do i__=1,47232
        read(1,'(A466)')ar__7
        read(ar__7,8)
     +  USB_6(i__),RAdeg_7(i__),DEdeg_7(i__),IRAC1_6(i__),
     +  IRAC2_6(i__),IRAC3_6(i__),IRAC4_6(i__),MIPS_6(i__),
     +  q_IRAC1_6(i__),q_IRAC2_6(i__),q_IRAC3_6(i__),q_IRAC4_6(i__),
     +  q_MIPS_6(i__),e_IRAC1_6(i__),e_IRAC2_6(i__),e_IRAC3_6(i__),
     +  e_IRAC4_6(i__),e_MIPS_6(i__),f_IRAC1_6(i__),f_IRAC2_6(i__),
     +  f_IRAC3_6(i__),f_IRAC4_6(i__),f_MIPS_6(i__),dDEN_6(i__),
     +  Imag_6(i__),Jmag_6(i__),Kmag_6(i__),d2M_6(i__),J2mag_6(i__),
     +  H2mag_6(i__),K2mag_6(i__),q_J2mag_6(i__),q_H2mag_6(i__),
     +  q_K2mag_6(i__),dIRAS_6(i__),F12_6(i__),q_F12_6(i__),
     +  F25_6(i__),q_F25_6(i__),F60_6(i__),q_F60_6(i__),F100_6(i__),
     +  q_F100_6(i__),dMSX_6(i__),B1_6(i__),q_B1_6(i__),B2_6(i__),
     +  q_B2_6(i__),A_6(i__),q_A_6(i__),C_6(i__),q_C_6(i__),D_6(i__),
     +  q_D_6(i__),E_6(i__),q_E_6(i__)
        if(ar__7(217:221) .EQ. '') dDEN_6(i__) = rNULL__
        if(ar__7(225:230) .EQ. '') Imag_6(i__) = rNULL__
        if(ar__7(234:239) .EQ. '') Jmag_6(i__) = rNULL__
        if(ar__7(243:248) .EQ. '') Kmag_6(i__) = rNULL__
        if(ar__7(252:256) .EQ. '') d2M_6(i__) = rNULL__
        if (idig(ar__7(260:265)).EQ.0) J2mag_6(i__) =  rNULL__
        if (idig(ar__7(269:274)).EQ.0) H2mag_6(i__) =  rNULL__
        if (idig(ar__7(278:283)).EQ.0) K2mag_6(i__) =  rNULL__
        if(ar__7(293:297) .EQ. '') dIRAS_6(i__) = rNULL__
        if (idig(ar__7(301:308)).EQ.0) F12_6(i__) =  rNULL__
        if (q_F12_6(i__) .EQ. 45) q_F12_6(i__) =  iNULL__
        if (idig(ar__7(316:323)).EQ.0) F25_6(i__) =  rNULL__
        if (q_F25_6(i__) .EQ. 45) q_F25_6(i__) =  iNULL__
        if (idig(ar__7(331:338)).EQ.0) F60_6(i__) =  rNULL__
        if (q_F60_6(i__) .EQ. 45) q_F60_6(i__) =  iNULL__
        if (idig(ar__7(346:353)).EQ.0) F100_6(i__) =  rNULL__
        if (q_F100_6(i__) .EQ. 45) q_F100_6(i__) =  iNULL__
        if(ar__7(361:364) .EQ. '') dMSX_6(i__) = rNULL__
        if (idig(ar__7(368:377)).EQ.0) B1_6(i__) =  rNULL__
        if (q_B1_6(i__) .EQ. 45) q_B1_6(i__) =  iNULL__
        if (idig(ar__7(385:394)).EQ.0) B2_6(i__) =  rNULL__
        if (q_B2_6(i__) .EQ. 45) q_B2_6(i__) =  iNULL__
        if (idig(ar__7(402:411)).EQ.0) A_6(i__) =  rNULL__
        if (q_A_6(i__) .EQ. 45) q_A_6(i__) =  iNULL__
        if (idig(ar__7(419:428)).EQ.0) C_6(i__) =  rNULL__
        if (q_C_6(i__) .EQ. 45) q_C_6(i__) =  iNULL__
        if (idig(ar__7(436:445)).EQ.0) D_6(i__) =  rNULL__
        if (q_D_6(i__) .EQ. 45) q_D_6(i__) =  iNULL__
        if (idig(ar__7(453:462)).EQ.0) E_6(i__) =  rNULL__
        if (q_E_6(i__) .EQ. 45) q_E_6(i__) =  iNULL__
c    ..............Just test output...........
        write(6,8)
     +  USB_6(i__),RAdeg_7(i__),DEdeg_7(i__),IRAC1_6(i__),
     +  IRAC2_6(i__),IRAC3_6(i__),IRAC4_6(i__),MIPS_6(i__),
     +  q_IRAC1_6(i__),q_IRAC2_6(i__),q_IRAC3_6(i__),q_IRAC4_6(i__),
     +  q_MIPS_6(i__),e_IRAC1_6(i__),e_IRAC2_6(i__),e_IRAC3_6(i__),
     +  e_IRAC4_6(i__),e_MIPS_6(i__),f_IRAC1_6(i__),f_IRAC2_6(i__),
     +  f_IRAC3_6(i__),f_IRAC4_6(i__),f_MIPS_6(i__),dDEN_6(i__),
     +  Imag_6(i__),Jmag_6(i__),Kmag_6(i__),d2M_6(i__),J2mag_6(i__),
     +  H2mag_6(i__),K2mag_6(i__),q_J2mag_6(i__),q_H2mag_6(i__),
     +  q_K2mag_6(i__),dIRAS_6(i__),F12_6(i__),q_F12_6(i__),
     +  F25_6(i__),q_F25_6(i__),F60_6(i__),q_F60_6(i__),F100_6(i__),
     +  q_F100_6(i__),dMSX_6(i__),B1_6(i__),q_B1_6(i__),B2_6(i__),
     +  q_B2_6(i__),A_6(i__),q_A_6(i__),C_6(i__),q_C_6(i__),D_6(i__),
     +  q_D_6(i__),E_6(i__),q_E_6(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