FORTRAN Generation
(/./ftp/cats/J/A_A/609/A26)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/609/A26 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/A+A/609/A26       Galactic Center old stars distribution (Gallego-Cano+, 2018)
*================================================================================
*The distribution of old stars around the Milky Way's central black hole.
*I. Star counts.
*    Gallego-Cano E., Schoedel R., Nogueras-Lara F., Gallego-Calvente A.T.,
*    Amaro-Seoane P., Baumgardt H.
*    <Astron. Astrophys. 609, A26 (2018)>
*    =2018A&A...609A..26G        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 1

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

      real*8        offX       (nr__) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY       (nr__) ! (arcsec) Y position offset (1)
      real*4        e_offX     (nr__) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY     (nr__) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag      (nr__) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag    (nr__) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag       (nr__) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag     (nr__) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK         (nr__) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK       (nr__) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p          (nr__) ! (%) Percentage reduction in observed stars(5)
      real*4        comp       (nr__) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag    (nr__) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table2.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 2

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

      real*8        offX_1     (nr__1) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_1     (nr__1) ! (arcsec) Y position offset (1)
      real*4        e_offX_1   (nr__1) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_1   (nr__1) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_1    (nr__1) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_1  (nr__1) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_1     (nr__1) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_1   (nr__1) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_1       (nr__1) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_1     (nr__1) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_1        (nr__1) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_1     (nr__1) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_1  (nr__1) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table3.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 3

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

      real*8        offX_2     (nr__2) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_2     (nr__2) ! (arcsec) Y position offset (1)
      real*4        e_offX_2   (nr__2) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_2   (nr__2) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_2    (nr__2) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_2  (nr__2) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_2     (nr__2) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_2   (nr__2) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_2       (nr__2) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_2     (nr__2) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_2        (nr__2) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_2     (nr__2) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_2  (nr__2) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table4.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 4

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

      real*8        offX_3     (nr__3) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_3     (nr__3) ! (arcsec) Y position offset (1)
      real*4        e_offX_3   (nr__3) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_3   (nr__3) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_3    (nr__3) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_3  (nr__3) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_3     (nr__3) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_3   (nr__3) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_3       (nr__3) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_3     (nr__3) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_3        (nr__3) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_3     (nr__3) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_3  (nr__3) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table5.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 1

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

      real*8        offX_4     (nr__4) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_4     (nr__4) ! (arcsec) Y position offset (1)
      real*4        e_offX_4   (nr__4) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_4   (nr__4) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_4    (nr__4) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_4  (nr__4) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_4     (nr__4) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_4   (nr__4) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_4       (nr__4) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_4     (nr__4) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_4        (nr__4) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_4     (nr__4) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_4  (nr__4) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table6.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 2

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

      real*8        offX_5     (nr__5) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_5     (nr__5) ! (arcsec) Y position offset (1)
      real*4        e_offX_5   (nr__5) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_5   (nr__5) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_5    (nr__5) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_5  (nr__5) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_5     (nr__5) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_5   (nr__5) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_5       (nr__5) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_5     (nr__5) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_5        (nr__5) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_5     (nr__5) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_5  (nr__5) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table7.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 3

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

      real*8        offX_6     (nr__6) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_6     (nr__6) ! (arcsec) Y position offset (1)
      real*4        e_offX_6   (nr__6) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_6   (nr__6) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_6    (nr__6) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_6  (nr__6) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_6     (nr__6) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_6   (nr__6) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_6       (nr__6) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_6     (nr__6) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_6        (nr__6) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_6     (nr__6) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_6  (nr__6) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table8.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 4

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

      real*8        offX_7     (nr__7) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_7     (nr__7) ! (arcsec) Y position offset (1)
      real*4        e_offX_7   (nr__7) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_7   (nr__7) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_7    (nr__7) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_7  (nr__7) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_7     (nr__7) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_7   (nr__7) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_7       (nr__7) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_7     (nr__7) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_7        (nr__7) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_7     (nr__7) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_7  (nr__7) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table9.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 1

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

      real*8        offX_8     (nr__8) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_8     (nr__8) ! (arcsec) Y position offset (1)
      real*4        e_offX_8   (nr__8) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_8   (nr__8) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_8    (nr__8) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_8  (nr__8) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_8     (nr__8) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_8   (nr__8) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_8       (nr__8) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_8     (nr__8) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_8        (nr__8) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_8     (nr__8) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_8  (nr__8) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table10.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 2

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

      real*8        offX_9     (nr__9) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_9     (nr__9) ! (arcsec) Y position offset (1)
      real*4        e_offX_9   (nr__9) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_9   (nr__9) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_9    (nr__9) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_9  (nr__9) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_9     (nr__9) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_9   (nr__9) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_9       (nr__9) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_9     (nr__9) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_9        (nr__9) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_9     (nr__9) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_9  (nr__9) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table11.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 3

      integer*4 nr__10
      parameter (nr__10=11217)	! Number of records
      character*79 ar__10 	! Full-size record

      real*8        offX_10    (nr__10) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_10    (nr__10) ! (arcsec) Y position offset (1)
      real*4        e_offX_10  (nr__10) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_10  (nr__10) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_10   (nr__10) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_10 (nr__10) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_10    (nr__10) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_10  (nr__10) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_10      (nr__10) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_10    (nr__10) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_10       (nr__10) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_10    (nr__10) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_10 (nr__10) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table12.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 4

      integer*4 nr__11
      parameter (nr__11=9189)	! Number of records
      character*79 ar__11 	! Full-size record

      real*8        offX_11    (nr__11) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_11    (nr__11) ! (arcsec) Y position offset (1)
      real*4        e_offX_11  (nr__11) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_11  (nr__11) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_11   (nr__11) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_11 (nr__11) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_11    (nr__11) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_11  (nr__11) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_11      (nr__11) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_11    (nr__11) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_11       (nr__11) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_11    (nr__11) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_11 (nr__11) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table13.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 1

      integer*4 nr__12
      parameter (nr__12=9457)	! Number of records
      character*79 ar__12 	! Full-size record

      real*8        offX_12    (nr__12) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_12    (nr__12) ! (arcsec) Y position offset (1)
      real*4        e_offX_12  (nr__12) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_12  (nr__12) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_12   (nr__12) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_12 (nr__12) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_12    (nr__12) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_12  (nr__12) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_12      (nr__12) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_12    (nr__12) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_12       (nr__12) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_12    (nr__12) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_12 (nr__12) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table14.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 2

      integer*4 nr__13
      parameter (nr__13=7868)	! Number of records
      character*79 ar__13 	! Full-size record

      real*8        offX_13    (nr__13) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_13    (nr__13) ! (arcsec) Y position offset (1)
      real*4        e_offX_13  (nr__13) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_13  (nr__13) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_13   (nr__13) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_13 (nr__13) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_13    (nr__13) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_13  (nr__13) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_13      (nr__13) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_13    (nr__13) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_13       (nr__13) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_13    (nr__13) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_13 (nr__13) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table15.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 3

      integer*4 nr__14
      parameter (nr__14=10698)	! Number of records
      character*79 ar__14 	! Full-size record

      real*8        offX_14    (nr__14) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_14    (nr__14) ! (arcsec) Y position offset (1)
      real*4        e_offX_14  (nr__14) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_14  (nr__14) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_14   (nr__14) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_14 (nr__14) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_14    (nr__14) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_14  (nr__14) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_14      (nr__14) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_14    (nr__14) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_14       (nr__14) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_14    (nr__14) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_14 (nr__14) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'table16.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 4

      integer*4 nr__15
      parameter (nr__15=8628)	! Number of records
      character*79 ar__15 	! Full-size record

      real*8        offX_15    (nr__15) ! (arcsec) X position offset(increasing to the East) (1)
      real*8        offY_15    (nr__15) ! (arcsec) Y position offset (1)
      real*4        e_offX_15  (nr__15) ! (arcsec) Uncertainty on X position (2)
      real*4        e_offY_15  (nr__15) ! (arcsec) Uncertainty on Y position (2)
      real*4        Ksmag_15   (nr__15) ! (mag) Ks-band magnitude (3)
      real*4        e_Ksmag_15 (nr__15) ! (mag) Uncertainty on Ks-band magnitude (2)
      real*4        Hmag_15    (nr__15) ! (mag) ?=0 H-band magnitude (3)
      real*4        e_Hmag_15  (nr__15) ! (mag) ?=0 Uncertainty on H-band magnitude (2)
      real*4        AK_15      (nr__15) ! (mag) Extinction in the Ks band (4)
      real*4        e_AK_15    (nr__15) ! (mag) Uncertainty on extinction in the Ks band
      real*4        p_15       (nr__15) ! (%) Percentage reduction in observed stars(5)
      real*4        comp_15    (nr__15) ! (mag) Source detection completeness in the Ks-images
      integer*4     ageFlag_15 (nr__15) ! [0,1] Age Flag (0: old star, 1: young star) (6)
*Note (1): Offset is from Sagittarius A*:
*    tables  1-4  : Pointing 1: x0= 672 and y0=1448
*    tables  5-8  : Pointing 2: x0= 671 and y0= 431
*    tables  9-12 : Pointing 3: x0=1692 and y0= 425
*    tables 13-16 : Pointing 4: x0=1694 and y0=1446
*Note (2): Uncertainties on  positions and magnitudes are the formal
*    uncertainties from StarFinder.
*Note (3): The photometry was calibrated with the stars IRS 16C, IRS 16NW, and
*    IRS 33N. See Schoedel et al. (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (4): The extinction maps for each of the four pointings were created with
*    the same method as applied in Schoedel et al.
*    (2010A&A...511A..18S, Cat. J/A+A/511/A18)
*Note (5): Correction factor that we apply to each detected star to correct the
*    individual stellar magnitudes for differential extinction by applying the
*    methodology of Chatzopoulos et al. (2015MNRAS.453..939C)
*Note (6): We use the list of spectroscopically identified young stars
*    from Do et al. (2009ApJ...703.1323D, Cat. J/ApJ/703/1323)

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

C  Declarations for 'tablea1.dat'	! *Surface density for brighter stars

      integer*4 nr__16
      parameter (nr__16=43)	! Number of records
      character*29 ar__16 	! Full-size record

      real*8        Rad        (nr__16) ! (arcsec) Projected radius from Saggitarius A*
      real*8        Den        (nr__16) ! (arcsec-2) Number stellar density for late-type stars
      real*8        e_Den      (nr__16) ! (arcsec-2) Uncertainty on number stellar density

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

C  Declarations for 'tablea2.dat'	! *Surface density for faint stars

      integer*4 nr__17
      parameter (nr__17=43)	! Number of records
      character*29 ar__17 	! Full-size record

      real*8        Rad_1      (nr__17) ! (arcsec) Projected radius from Saggitarius A*
      real*8        Den_1      (nr__17) ! (arcsec-2) Number stellar density for late-type stars
      real*8        e_Den_1    (nr__17) ! (arcsec-2) Uncertainty on number stellar density

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

C  Declarations for 'tablea3.dat'	! *Surface density for fainter stars

      integer*4 nr__18
      parameter (nr__18=43)	! Number of records
      character*29 ar__18 	! Full-size record

      real*8        Rad_2      (nr__18) ! (arcsec) Projected radius from Saggitarius A*
      real*8        Den_2      (nr__18) ! (arcsec-2) Number stellar density for late-type stars
      real*8        e_Den_2    (nr__18) ! (arcsec-2) Uncertainty on number stellar density

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

C  Loading file 'table1.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 1

C  Format for file interpretation

    1 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,9237
        read(1,'(A79)')ar__
        read(ar__,1)
     +  offX(i__),offY(i__),e_offX(i__),e_offY(i__),Ksmag(i__),
     +  e_Ksmag(i__),Hmag(i__),e_Hmag(i__),AK(i__),e_AK(i__),p(i__),
     +  comp(i__),ageFlag(i__)
c    ..............Just test output...........
        write(6,1)
     +  offX(i__),offY(i__),e_offX(i__),e_offY(i__),Ksmag(i__),
     +  e_Ksmag(i__),Hmag(i__),e_Hmag(i__),AK(i__),e_AK(i__),p(i__),
     +  comp(i__),ageFlag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 2

C  Format for file interpretation

    2 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,7497
        read(1,'(A79)')ar__1
        read(ar__1,2)
     +  offX_1(i__),offY_1(i__),e_offX_1(i__),e_offY_1(i__),
     +  Ksmag_1(i__),e_Ksmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  AK_1(i__),e_AK_1(i__),p_1(i__),comp_1(i__),ageFlag_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  offX_1(i__),offY_1(i__),e_offX_1(i__),e_offY_1(i__),
     +  Ksmag_1(i__),e_Ksmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  AK_1(i__),e_AK_1(i__),p_1(i__),comp_1(i__),ageFlag_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 3

C  Format for file interpretation

    3 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,10551
        read(1,'(A79)')ar__2
        read(ar__2,3)
     +  offX_2(i__),offY_2(i__),e_offX_2(i__),e_offY_2(i__),
     +  Ksmag_2(i__),e_Ksmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),
     +  AK_2(i__),e_AK_2(i__),p_2(i__),comp_2(i__),ageFlag_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  offX_2(i__),offY_2(i__),e_offX_2(i__),e_offY_2(i__),
     +  Ksmag_2(i__),e_Ksmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),
     +  AK_2(i__),e_AK_2(i__),p_2(i__),comp_2(i__),ageFlag_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! *Galactic Centre stars. Pointing 1 and SF-set 4

C  Format for file interpretation

    4 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,8476
        read(1,'(A79)')ar__3
        read(ar__3,4)
     +  offX_3(i__),offY_3(i__),e_offX_3(i__),e_offY_3(i__),
     +  Ksmag_3(i__),e_Ksmag_3(i__),Hmag_3(i__),e_Hmag_3(i__),
     +  AK_3(i__),e_AK_3(i__),p_3(i__),comp_3(i__),ageFlag_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  offX_3(i__),offY_3(i__),e_offX_3(i__),e_offY_3(i__),
     +  Ksmag_3(i__),e_Ksmag_3(i__),Hmag_3(i__),e_Hmag_3(i__),
     +  AK_3(i__),e_AK_3(i__),p_3(i__),comp_3(i__),ageFlag_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 1

C  Format for file interpretation

    5 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,8201
        read(1,'(A79)')ar__4
        read(ar__4,5)
     +  offX_4(i__),offY_4(i__),e_offX_4(i__),e_offY_4(i__),
     +  Ksmag_4(i__),e_Ksmag_4(i__),Hmag_4(i__),e_Hmag_4(i__),
     +  AK_4(i__),e_AK_4(i__),p_4(i__),comp_4(i__),ageFlag_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  offX_4(i__),offY_4(i__),e_offX_4(i__),e_offY_4(i__),
     +  Ksmag_4(i__),e_Ksmag_4(i__),Hmag_4(i__),e_Hmag_4(i__),
     +  AK_4(i__),e_AK_4(i__),p_4(i__),comp_4(i__),ageFlag_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 2

C  Format for file interpretation

    6 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,6766
        read(1,'(A79)')ar__5
        read(ar__5,6)
     +  offX_5(i__),offY_5(i__),e_offX_5(i__),e_offY_5(i__),
     +  Ksmag_5(i__),e_Ksmag_5(i__),Hmag_5(i__),e_Hmag_5(i__),
     +  AK_5(i__),e_AK_5(i__),p_5(i__),comp_5(i__),ageFlag_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  offX_5(i__),offY_5(i__),e_offX_5(i__),e_offY_5(i__),
     +  Ksmag_5(i__),e_Ksmag_5(i__),Hmag_5(i__),e_Hmag_5(i__),
     +  AK_5(i__),e_AK_5(i__),p_5(i__),comp_5(i__),ageFlag_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 3

C  Format for file interpretation

    7 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,9321
        read(1,'(A79)')ar__6
        read(ar__6,7)
     +  offX_6(i__),offY_6(i__),e_offX_6(i__),e_offY_6(i__),
     +  Ksmag_6(i__),e_Ksmag_6(i__),Hmag_6(i__),e_Hmag_6(i__),
     +  AK_6(i__),e_AK_6(i__),p_6(i__),comp_6(i__),ageFlag_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  offX_6(i__),offY_6(i__),e_offX_6(i__),e_offY_6(i__),
     +  Ksmag_6(i__),e_Ksmag_6(i__),Hmag_6(i__),e_Hmag_6(i__),
     +  AK_6(i__),e_AK_6(i__),p_6(i__),comp_6(i__),ageFlag_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! *Galactic Centre stars. Pointing 2 and SF-set 4

C  Format for file interpretation

    8 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,7816
        read(1,'(A79)')ar__7
        read(ar__7,8)
     +  offX_7(i__),offY_7(i__),e_offX_7(i__),e_offY_7(i__),
     +  Ksmag_7(i__),e_Ksmag_7(i__),Hmag_7(i__),e_Hmag_7(i__),
     +  AK_7(i__),e_AK_7(i__),p_7(i__),comp_7(i__),ageFlag_7(i__)
c    ..............Just test output...........
        write(6,8)
     +  offX_7(i__),offY_7(i__),e_offX_7(i__),e_offY_7(i__),
     +  Ksmag_7(i__),e_Ksmag_7(i__),Hmag_7(i__),e_Hmag_7(i__),
     +  AK_7(i__),e_AK_7(i__),p_7(i__),comp_7(i__),ageFlag_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 1

C  Format for file interpretation

    9 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,10023
        read(1,'(A79)')ar__8
        read(ar__8,9)
     +  offX_8(i__),offY_8(i__),e_offX_8(i__),e_offY_8(i__),
     +  Ksmag_8(i__),e_Ksmag_8(i__),Hmag_8(i__),e_Hmag_8(i__),
     +  AK_8(i__),e_AK_8(i__),p_8(i__),comp_8(i__),ageFlag_8(i__)
c    ..............Just test output...........
        write(6,9)
     +  offX_8(i__),offY_8(i__),e_offX_8(i__),e_offY_8(i__),
     +  Ksmag_8(i__),e_Ksmag_8(i__),Hmag_8(i__),e_Hmag_8(i__),
     +  AK_8(i__),e_AK_8(i__),p_8(i__),comp_8(i__),ageFlag_8(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 2

C  Format for file interpretation

   10 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,8202
        read(1,'(A79)')ar__9
        read(ar__9,10)
     +  offX_9(i__),offY_9(i__),e_offX_9(i__),e_offY_9(i__),
     +  Ksmag_9(i__),e_Ksmag_9(i__),Hmag_9(i__),e_Hmag_9(i__),
     +  AK_9(i__),e_AK_9(i__),p_9(i__),comp_9(i__),ageFlag_9(i__)
c    ..............Just test output...........
        write(6,10)
     +  offX_9(i__),offY_9(i__),e_offX_9(i__),e_offY_9(i__),
     +  Ksmag_9(i__),e_Ksmag_9(i__),Hmag_9(i__),e_Hmag_9(i__),
     +  AK_9(i__),e_AK_9(i__),p_9(i__),comp_9(i__),ageFlag_9(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table11.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 3

C  Format for file interpretation

   11 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,11217
        read(1,'(A79)')ar__10
        read(ar__10,11)
     +  offX_10(i__),offY_10(i__),e_offX_10(i__),e_offY_10(i__),
     +  Ksmag_10(i__),e_Ksmag_10(i__),Hmag_10(i__),e_Hmag_10(i__),
     +  AK_10(i__),e_AK_10(i__),p_10(i__),comp_10(i__),ageFlag_10(i__)
c    ..............Just test output...........
        write(6,11)
     +  offX_10(i__),offY_10(i__),e_offX_10(i__),e_offY_10(i__),
     +  Ksmag_10(i__),e_Ksmag_10(i__),Hmag_10(i__),e_Hmag_10(i__),
     +  AK_10(i__),e_AK_10(i__),p_10(i__),comp_10(i__),ageFlag_10(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table12.dat'	! *Galactic Centre stars. Pointing 3 and SF-set 4

C  Format for file interpretation

   12 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table12.dat')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,9189
        read(1,'(A79)')ar__11
        read(ar__11,12)
     +  offX_11(i__),offY_11(i__),e_offX_11(i__),e_offY_11(i__),
     +  Ksmag_11(i__),e_Ksmag_11(i__),Hmag_11(i__),e_Hmag_11(i__),
     +  AK_11(i__),e_AK_11(i__),p_11(i__),comp_11(i__),ageFlag_11(i__)
c    ..............Just test output...........
        write(6,12)
     +  offX_11(i__),offY_11(i__),e_offX_11(i__),e_offY_11(i__),
     +  Ksmag_11(i__),e_Ksmag_11(i__),Hmag_11(i__),e_Hmag_11(i__),
     +  AK_11(i__),e_AK_11(i__),p_11(i__),comp_11(i__),ageFlag_11(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table13.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 1

C  Format for file interpretation

   13 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table13.dat')
      write(6,*) '....Loading file: table13.dat'
      do i__=1,9457
        read(1,'(A79)')ar__12
        read(ar__12,13)
     +  offX_12(i__),offY_12(i__),e_offX_12(i__),e_offY_12(i__),
     +  Ksmag_12(i__),e_Ksmag_12(i__),Hmag_12(i__),e_Hmag_12(i__),
     +  AK_12(i__),e_AK_12(i__),p_12(i__),comp_12(i__),ageFlag_12(i__)
c    ..............Just test output...........
        write(6,13)
     +  offX_12(i__),offY_12(i__),e_offX_12(i__),e_offY_12(i__),
     +  Ksmag_12(i__),e_Ksmag_12(i__),Hmag_12(i__),e_Hmag_12(i__),
     +  AK_12(i__),e_AK_12(i__),p_12(i__),comp_12(i__),ageFlag_12(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table14.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 2

C  Format for file interpretation

   14 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table14.dat')
      write(6,*) '....Loading file: table14.dat'
      do i__=1,7868
        read(1,'(A79)')ar__13
        read(ar__13,14)
     +  offX_13(i__),offY_13(i__),e_offX_13(i__),e_offY_13(i__),
     +  Ksmag_13(i__),e_Ksmag_13(i__),Hmag_13(i__),e_Hmag_13(i__),
     +  AK_13(i__),e_AK_13(i__),p_13(i__),comp_13(i__),ageFlag_13(i__)
c    ..............Just test output...........
        write(6,14)
     +  offX_13(i__),offY_13(i__),e_offX_13(i__),e_offY_13(i__),
     +  Ksmag_13(i__),e_Ksmag_13(i__),Hmag_13(i__),e_Hmag_13(i__),
     +  AK_13(i__),e_AK_13(i__),p_13(i__),comp_13(i__),ageFlag_13(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table15.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 3

C  Format for file interpretation

   15 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table15.dat')
      write(6,*) '....Loading file: table15.dat'
      do i__=1,10698
        read(1,'(A79)')ar__14
        read(ar__14,15)
     +  offX_14(i__),offY_14(i__),e_offX_14(i__),e_offY_14(i__),
     +  Ksmag_14(i__),e_Ksmag_14(i__),Hmag_14(i__),e_Hmag_14(i__),
     +  AK_14(i__),e_AK_14(i__),p_14(i__),comp_14(i__),ageFlag_14(i__)
c    ..............Just test output...........
        write(6,15)
     +  offX_14(i__),offY_14(i__),e_offX_14(i__),e_offY_14(i__),
     +  Ksmag_14(i__),e_Ksmag_14(i__),Hmag_14(i__),e_Hmag_14(i__),
     +  AK_14(i__),e_AK_14(i__),p_14(i__),comp_14(i__),ageFlag_14(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table16.dat'	! *Galactic Centre stars. Pointing 4 and SF-set 4

C  Format for file interpretation

   16 format(
     +  F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table16.dat')
      write(6,*) '....Loading file: table16.dat'
      do i__=1,8628
        read(1,'(A79)')ar__15
        read(ar__15,16)
     +  offX_15(i__),offY_15(i__),e_offX_15(i__),e_offY_15(i__),
     +  Ksmag_15(i__),e_Ksmag_15(i__),Hmag_15(i__),e_Hmag_15(i__),
     +  AK_15(i__),e_AK_15(i__),p_15(i__),comp_15(i__),ageFlag_15(i__)
c    ..............Just test output...........
        write(6,16)
     +  offX_15(i__),offY_15(i__),e_offX_15(i__),e_offY_15(i__),
     +  Ksmag_15(i__),e_Ksmag_15(i__),Hmag_15(i__),e_Hmag_15(i__),
     +  AK_15(i__),e_AK_15(i__),p_15(i__),comp_15(i__),ageFlag_15(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea1.dat'	! *Surface density for brighter stars

C  Format for file interpretation

   17 format(F9.6,1X,F9.6,1X,F9.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea1.dat')
      write(6,*) '....Loading file: tablea1.dat'
      do i__=1,43
        read(1,'(A29)')ar__16
        read(ar__16,17)Rad(i__),Den(i__),e_Den(i__)
c    ..............Just test output...........
        write(6,17)Rad(i__),Den(i__),e_Den(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea2.dat'	! *Surface density for faint stars

C  Format for file interpretation

   18 format(F9.6,1X,F9.6,1X,F9.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea2.dat')
      write(6,*) '....Loading file: tablea2.dat'
      do i__=1,43
        read(1,'(A29)')ar__17
        read(ar__17,18)Rad_1(i__),Den_1(i__),e_Den_1(i__)
c    ..............Just test output...........
        write(6,18)Rad_1(i__),Den_1(i__),e_Den_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablea3.dat'	! *Surface density for fainter stars

C  Format for file interpretation

   19 format(F9.6,1X,F9.6,1X,F9.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea3.dat')
      write(6,*) '....Loading file: tablea3.dat'
      do i__=1,43
        read(1,'(A29)')ar__18
        read(ar__18,19)Rad_2(i__),Den_2(i__),e_Den_2(i__)
c    ..............Just test output...........
        write(6,19)Rad_2(i__),Den_2(i__),e_Den_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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