FORTRAN Generation
(/./ftp/cats/J/ApJS/197/33)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/197/33 into FORTRAN code for reading data files line by line.

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-05
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/ApJS/197/33          The M87 globular cluster system          (Strader+, 2011)
*================================================================================
*Wide-field precision kinematics of the M87 globular cluster system.
*    Strader J., Romanowsky A.J., Brodie J.P., Spitler L.R., Beasley M.A.,
*    Arnold J.A., Tamura N., Sharples R.M., Arimoto N.
*   <Astrophys. J. Suppl. Ser., 197, 33 (2011)>
*   =2011ApJS..197...33S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Data for Keck/DEIMOS spectroscopic candidates

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

C  J2000 position composed of: RAdeg DEdeg
      character*6   ID          ! Identification (G1)
      real*8        RAdeg       ! (deg) Right Ascension in decimal degrees (J2000) (1)
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000) (1)
      character*3   Mask        ! [1,4] Mask(s) (2)
      real*4        Vmag        ! (mag) Subaru/Suprime-Cam V band magnitude (1)
      real*4        B_V         ! (mag) Subaru/Suprime-Cam (B-V) color index (1)
      real*4        V_I         ! (mag) Subaru/Suprime-Cam (V-I) color index (1)
      character*4   Class       ! Classification (G2)
      integer*4     HRV         ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV       ! (km/s) ? Error in HRV
*Note (1): The astrometry and photometry in this table are taken from Tamura
*     et al. (2006MNRAS.373..588T); the photometry is corrected for
*     foreground reddening.
*Note (2): We designed four Keck/DEIMOS slit masks, each covering an area of
*     ~16'x5'. These were arranged with minimal overlap in a 2x2 grid
*     covering a total area of ~32'x10', starting ~8' from the center of M87.

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

C  Declarations for 'table3.dat'	! Data for supplementary Keck/DEIMOS spectroscopy

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

C  J2000 position composed of: RAdeg DEdeg
      character*6   ID_1        ! Identification (G1)
      real*8        RAdeg_1     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_1     ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_1      ! [1,4]? Mask (for table 5 only)
      real*4        imag        ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r         ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i         ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_1     ! Classification (G2)
      integer*4     HRV_1       ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_1     ! (km/s) ? Error in HRV
*Note (1): The photometry is corrected for foreground reddening.

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

C  Declarations for 'table4.dat'	! Data for MMT/Hectospec spectroscopic candidates

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

C  J2000 position composed of: RAdeg DEdeg
      character*6   ID_2        ! Identification (G1)
      real*8        RAdeg_2     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_2     ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_2      ! [1,4]? Mask (for table 5 only)
      real*4        imag_1      ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r_1       ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i_1       ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_2     ! Classification (G2)
      integer*4     HRV_2       ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_2     ! (km/s) ? Error in HRV
*Note (1): The photometry is corrected for foreground reddening.

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

C  Declarations for 'table5.dat'	! Data for Keck/LRIS spectroscopic candidates

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

C  J2000 position composed of: RAdeg DEdeg
      character*6   ID_3        ! Identification (G1)
      real*8        RAdeg_3     ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_3     ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_3      ! [1,4]? Mask (for table 5 only)
      real*4        imag_2      ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r_2       ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i_2       ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_3     ! Classification (G2)
      integer*4     HRV_3       ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_3     ! (km/s) ? Error in HRV
*Note (1): The photometry is corrected for foreground reddening.

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

C  Declarations for 'table6.dat'	! Globular cluster/Ultracompact dwarfs (UCDs)
                             radial velocities with multiple measurements

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

      character*6   ID_4        ! ID Identification (G1)
      character*1   f_ID        ! [b] Flag on ID (1)
      integer*4     H01         ! (km/s) ? Hanes et al. 2001, Cat .J/ApJ/559/812
*                                heliocentric radial velocity
      integer*4     e_H01       ! (km/s) ? Uncertainty in H01
      integer*4     K97         ! (km/s) ? Kohen & Ryzhov, 1997ApJ...486..230C,
*                                heliocentric radial velocity
      integer*4     e_K97       ! (km/s) ? Uncertainty in K97
      integer*4     K00         ! (km/s) ? Cohen, 2000AJ....119..162C, heliocentric
*                                radial velocity
      integer*4     e_K00       ! (km/s) ? Uncertainty in K00
      integer*4     M90         ! (km/s) ? Mould et al., 1990AJ.....99.1823M,
*                                heliocentric radial velocity
      integer*4     e_M90       ! (km/s) ? Uncertainty in M90
      integer*4     M87         ! (km/s) ? Mould et al., 1987AJ.....93...53M,
*                                heliocentric radial velocity
      integer*4     e_M87       ! (km/s) ? Uncertainty in M87
      integer*4     HB87        ! (km/s) ? Huchra & Brodie, 1987AJ.....93..779H,
*                                heliocentric radial velocity
      integer*4     e_HB87      ! (km/s) ? Uncertainty in HB87
      integer*4     MMT         ! (km/s) ? This paper, MMT/Hectospec heliocentric radial
*                                velocity
      integer*4     e_MMT       ! (km/s) ? Uncertainty in MMT
      integer*4     LRIS        ! (km/s) ? This paper, Keck/LRIS heliocentric radial
*                                velocity
      integer*4     e_LRIS      ! (km/s) ? Uncertainty in LRIS
      integer*4     D1          ! (km/s) ? This paper, Keck/DEIMOS run 1 heliocentric
*                                radial velocity
      integer*4     e_D1        ! (km/s) ? Uncertainty in D1
      integer*4     D2          ! (km/s) ? This paper, Keck/DEIMOS run 2 heliocentric
*                                radial velocity
      integer*4     e_D2        ! (km/s) ? Uncertainty in D2
      integer*4     H05_H07     ! (km/s) ? Hasegan et al., 2005ApJ...627..203H, or
*                                Hasegan 2007PhDT.........4H, heliocentric radial
*                                velocity
      integer*4     e_H05_H07   ! (km/s) ? Uncertainty in H05/H07
      integer*4     J06         ! (km/s) ? Jones et al., 2006AJ....131..312J,
*                                heliocentric radial velocity
      integer*4     e_J06       ! (km/s) ? Uncertainty in J06
      integer*4     E07         ! (km/s) ? Evstigneeva et al., 2007AJ....133.1722E,
*                                heliocentric radial velocity
      integer*4     e_E07       ! (km/s) ? Uncertainty in E07
      integer*4     F08         ! (km/s) ? Firth et al., 2008MNRAS.389.1539F,
*                                heliocentric radial velocity
      integer*4     e_F08       ! (km/s) ? Uncertainty in F08
*Note (1):
*    b = "catastrophic" outlier, described in the text.

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

C  Declarations for 'table7.dat'	! New half-light radii

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

      character*6   ID_5        ! Star name (G1)
      real*4        rh          ! (pc) ? Half-light radius assuming a distance of 16.5Mpc
      character*1   f_rh        ! [d] d: unresolved; the upper limit is r_h_~1pc
      character*5   Inst        ! Instrument
      character*6   Filt        ! Filter (1)
*Note (1): The STIS "50CCD" denotes a clear, unfiltered aperture.

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

C  Declarations for 'table8.dat'	! Spectroscopic data for all foreground stars

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

C  J2000 position composed of: RAdeg DEdeg
      character*7   ID_6        ! Star name (G1)
      character*1   f_ID_1      ! [f] Flag on ID (1)
      character*13  AID         ! Alternative identification(s) (G1)
      real*8        RAdeg_4     ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_4     ! (deg) Declination in decimal degrees (J2000)
      character*1   f_imag      ! [i] Flag on imag (2)
      real*4        imag_3      ! (mag) ? The CFHT i-band magnitude
      real*4        g_r_3       ! (mag) ? The CFHT (g-r) color index
      real*4        g_i_3       ! (mag) ? The CFHT (g-i) color index
      integer*4     HRV_4       ! (km/s) Heliocentric radial velocity of the object
      integer*4     e_HRV_4     ! (km/s) ? HRV uncertainty
      character*23  r_HRV       ! Source(s) (3)
      character*1   HST         ! [Y/N] Confirmed by high-resolution HST imaging?
*Note (1):
*   f = This radial velocity is associated with an unidentified member of a
*       close pair. The values listed are for the source closer to the reported
*       position; the second object is H51472, coordinates: 187.75358 +12.45462,
*       i=21.50, g-i=0.87.
*Note (2):
*   i = For objects without CFHT imaging, an extinction-corrected i-band
*       magnitude has been estimated from the other available photometry.
*Note (3): Source of HRV as follows:
*  H01  = Hanes et al. (2001, Cat. J/ApJ/559/812);
*  K97  = Cohen & Ryzhov (1997ApJ...486..230C);
*  M90  = Mould et al. (1990AJ.....99.1823M);
*  HB87 = Huchra & Brodie (1987AJ.....93..779H);
*  MMT  = this paper, (MMT/Hectospec);
*  D1   = this paper, (Keck/DEIMOS run 1);
*  D2   = this paper, (Keck/DEIMOS run 2);
*  LRIS = this paper, (Keck/LRIS).

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

C  Declarations for 'table9.dat'	! Spectroscopic and photometric data for all M87
                             globular clusters and ultracompact dwarfs

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

C  J2000 position composed of: RAdeg DEdeg
      character*6   ID_7        ! Identification (G1)
      character*1   f_ID_2      ! [i] Flag on ID (1)
      character*13  AID_1       ! Alternate identification (G1)
      real*8        RAdeg_5     ! (deg) Right Ascension in decimal degrees (J2000) (2)
      real*8        DEdeg_5     ! (deg) Declination in decimal degrees (J2000) (2)
      real*4        imag_4      ! (mag) The CFHT i band magnitude
      character*1   f_imag_1    ! [i] Flag on imag (3)
      real*4        g_r_4       ! (mag) ? The CFHT (g-r) color index
      real*4        g_i_4       ! (mag) ? The CFHT (g-i) color index
      integer*4     HRV_5       ! (km/s) Heliocenric radial velocity
      integer*4     e_HRV_5     ! (km/s) Error in HRV
      character*19  r_HRV_1     ! Source of HRV (4)
      real*4        rh_1        ! (pc) ? Half-light radius; assuming D=16.5Mpc
      character*1   f_rh_1      ! [u] Indicates unresolved source for rh
      character*3   r_rh        ! Source of rh (5)
      character*4   spop        ! Cluster subpopulation (blue or red) (6)
      character*5   Class_4     ! Source classification (7)
*Note (1):
*    i = As discussed in Section 4.4.4, this GC may be a member of NGC 4478.
*Note (2): Right Ascension and Declination were re-calculated using the CFHT
*     g-band image, and so there may be differences between the listed
*     coordinates of objects in common between Table 1 and this table.
*     Several objects outside of the field of view of the CFHT images
*     retain coordinates from Tamura et al. (2006MNRAS.373..588T).
*Note (3):
*    i = an object without CFHT imaging thus an extinction-corrected i-band
*        magnitude was estimated from the other available photometry.
*Note (4): Multiple sources are listed chronologically. The sources are as
*     listed in Table 6, with the addition of P10 (Paudel et al.
*     2010ApJ...724L..64P).
*Note (5): Reference as follows:
*    E08 = Evstigneeva et al. (2008AJ....136..461E);
*    H05 = Hasegan et al. (2005ApJ...627..203H);
*    H07 = Hasegan (2007PhDT.........4H);
*    J09 = Jordan et al. (2009, Cat. J/ApJS/180/54);
*    L01 = Larsen et al. (2001AJ....121.2974L);
*    M09 = Madrid et al. (2009ApJ...705..237M);
*    S11 = This paper.
*Note (6): The subpopulation (blue metal-poor or red metal-rich) to which
*     the cluster belongs, as determined in Section 4.6.
*Note (7): Whether the object is classified as a transition object ("trans")
*     or ultra-compact dwarf ("UCD") on the basis of its half-light radius,
*     according to the boundaries of 5.25pc and 9.5pc discussed in S 4.4.1.

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

C  Declarations for 'table10.dat'	! Supplementary photometric data for M87 globular
                             clusters and ultracompact dwarfs

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

      character*6   ID_8        ! Identification (G1)
      real*4        T1mag       ! (mag) ? The Washington T_1_ band magnitude (1)
      real*4        C_T1        ! (mag) ? The Washington (C-T_1_) color index (1)
      real*4        gmag        ! (mag) ? HST/ACS g band magnitude (2)
      real*4        g_z         ! (mag) ? HST/ACS (g-z) color index (2)
      real*4        Vmag_1      ! (mag) ? The Subaru V band magnitude
      real*4        B_V_1       ! (mag) ? The Subaru (B-V) color index
      real*4        V_I_1       ! (mag) ? The Subaru (V-I) color index
*Note (1): From Hanes et al. 2001, Cat. J/ApJ/559/812
*Note (2): From Jordan et al., 2009, Cat. J/ApJS/180/54

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

C  Loading file 'table2.dat'	! Data for Keck/DEIMOS spectroscopic candidates

C  Format for file interpretation

    1 format(
     +  A6,1X,F9.5,1X,F8.5,1X,A3,1X,F5.2,1X,F4.2,1X,F4.2,1X,A4,1X,I4,
     +  1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,254
        read(1,'(A58)')ar__
        read(ar__,1)ID,RAdeg,DEdeg,Mask,Vmag,B_V,V_I,Class,HRV,e_HRV
        if(ar__(52:55) .EQ. '') HRV = iNULL__
        if(ar__(57:58) .EQ. '') e_HRV = iNULL__
c    ..............Just test output...........
        write(6,1)ID,RAdeg,DEdeg,Mask,Vmag,B_V,V_I,Class,HRV,e_HRV
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Data for supplementary Keck/DEIMOS spectroscopy

C  Format for file interpretation

    2 format(
     +  A6,1X,F9.5,1X,F8.5,1X,I1,1X,F5.2,1X,F4.2,1X,F4.2,1X,A4,1X,I4,
     +  1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,29
        read(1,'(A56)')ar__1
        read(ar__1,2)
     +  ID_1,RAdeg_1,DEdeg_1,Mask_1,imag,g_r,g_i,Class_1,HRV_1,e_HRV_1
        if(ar__1(27:27) .EQ. '') Mask_1 = iNULL__
        if(ar__1(29:33) .EQ. '') imag = rNULL__
        if(ar__1(35:38) .EQ. '') g_r = rNULL__
        if(ar__1(40:43) .EQ. '') g_i = rNULL__
        if(ar__1(50:53) .EQ. '') HRV_1 = iNULL__
        if(ar__1(55:56) .EQ. '') e_HRV_1 = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  ID_1,RAdeg_1,DEdeg_1,Mask_1,imag,g_r,g_i,Class_1,HRV_1,e_HRV_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Data for MMT/Hectospec spectroscopic candidates

C  Format for file interpretation

    3 format(
     +  A6,1X,F9.5,1X,F8.5,1X,I1,1X,F5.2,1X,F4.2,1X,F4.2,1X,A4,1X,I4,
     +  1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,244
        read(1,'(A56)')ar__2
        read(ar__2,3)
     +  ID_2,RAdeg_2,DEdeg_2,Mask_2,imag_1,g_r_1,g_i_1,Class_2,HRV_2,
     +  e_HRV_2
        if(ar__2(27:27) .EQ. '') Mask_2 = iNULL__
        if(ar__2(29:33) .EQ. '') imag_1 = rNULL__
        if(ar__2(35:38) .EQ. '') g_r_1 = rNULL__
        if(ar__2(40:43) .EQ. '') g_i_1 = rNULL__
        if(ar__2(50:53) .EQ. '') HRV_2 = iNULL__
        if(ar__2(55:56) .EQ. '') e_HRV_2 = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  ID_2,RAdeg_2,DEdeg_2,Mask_2,imag_1,g_r_1,g_i_1,Class_2,HRV_2,
     +  e_HRV_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Data for Keck/LRIS spectroscopic candidates

C  Format for file interpretation

    4 format(
     +  A6,1X,F9.5,1X,F8.5,1X,I1,1X,F5.2,1X,F4.2,1X,F4.2,1X,A4,1X,I4,
     +  1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,209
        read(1,'(A56)')ar__3
        read(ar__3,4)
     +  ID_3,RAdeg_3,DEdeg_3,Mask_3,imag_2,g_r_2,g_i_2,Class_3,HRV_3,
     +  e_HRV_3
        if(ar__3(27:27) .EQ. '') Mask_3 = iNULL__
        if(ar__3(29:33) .EQ. '') imag_2 = rNULL__
        if(ar__3(35:38) .EQ. '') g_r_2 = rNULL__
        if(ar__3(40:43) .EQ. '') g_i_2 = rNULL__
        if(ar__3(50:53) .EQ. '') HRV_3 = iNULL__
        if(ar__3(55:56) .EQ. '') e_HRV_3 = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  ID_3,RAdeg_3,DEdeg_3,Mask_3,imag_2,g_r_2,g_i_2,Class_3,HRV_3,
     +  e_HRV_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Globular cluster/Ultracompact dwarfs (UCDs)
*                             radial velocities with multiple measurements

C  Format for file interpretation

    5 format(
     +  A6,1X,A1,1X,I4,1X,I3,1X,I4,1X,I3,1X,I4,1X,I2,1X,I4,1X,I3,1X,
     +  I4,1X,I3,1X,I4,1X,I3,1X,I4,1X,I2,1X,I4,1X,I2,1X,I4,1X,I2,1X,
     +  I4,1X,I2,1X,I4,1X,I2,1X,I4,1X,I3,1X,I4,1X,I1,1X,I4,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,114
        read(1,'(A125)')ar__4
        read(ar__4,5)
     +  ID_4,f_ID,H01,e_H01,K97,e_K97,K00,e_K00,M90,e_M90,M87,e_M87,
     +  HB87,e_HB87,MMT,e_MMT,LRIS,e_LRIS,D1,e_D1,D2,e_D2,H05_H07,
     +  e_H05_H07,J06,e_J06,E07,e_E07,F08,e_F08
        if(ar__4(10:13) .EQ. '') H01 = iNULL__
        if(ar__4(15:17) .EQ. '') e_H01 = iNULL__
        if(ar__4(19:22) .EQ. '') K97 = iNULL__
        if(ar__4(24:26) .EQ. '') e_K97 = iNULL__
        if(ar__4(28:31) .EQ. '') K00 = iNULL__
        if(ar__4(33:34) .EQ. '') e_K00 = iNULL__
        if(ar__4(36:39) .EQ. '') M90 = iNULL__
        if(ar__4(41:43) .EQ. '') e_M90 = iNULL__
        if(ar__4(45:48) .EQ. '') M87 = iNULL__
        if(ar__4(50:52) .EQ. '') e_M87 = iNULL__
        if(ar__4(54:57) .EQ. '') HB87 = iNULL__
        if(ar__4(59:61) .EQ. '') e_HB87 = iNULL__
        if(ar__4(63:66) .EQ. '') MMT = iNULL__
        if(ar__4(68:69) .EQ. '') e_MMT = iNULL__
        if(ar__4(71:74) .EQ. '') LRIS = iNULL__
        if(ar__4(76:77) .EQ. '') e_LRIS = iNULL__
        if(ar__4(79:82) .EQ. '') D1 = iNULL__
        if(ar__4(84:85) .EQ. '') e_D1 = iNULL__
        if(ar__4(87:90) .EQ. '') D2 = iNULL__
        if(ar__4(92:93) .EQ. '') e_D2 = iNULL__
        if(ar__4(95:98) .EQ. '') H05_H07 = iNULL__
        if(ar__4(100:101) .EQ. '') e_H05_H07 = iNULL__
        if(ar__4(103:106) .EQ. '') J06 = iNULL__
        if(ar__4(108:110) .EQ. '') e_J06 = iNULL__
        if(ar__4(112:115) .EQ. '') E07 = iNULL__
        if(ar__4(117:117) .EQ. '') e_E07 = iNULL__
        if(ar__4(119:122) .EQ. '') F08 = iNULL__
        if(ar__4(124:125) .EQ. '') e_F08 = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  ID_4,f_ID,H01,e_H01,K97,e_K97,K00,e_K00,M90,e_M90,M87,e_M87,
     +  HB87,e_HB87,MMT,e_MMT,LRIS,e_LRIS,D1,e_D1,D2,e_D2,H05_H07,
     +  e_H05_H07,J06,e_J06,E07,e_E07,F08,e_F08
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! New half-light radii

C  Format for file interpretation

    6 format(A6,1X,F4.1,A1,1X,A5,1X,A6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,197
        read(1,'(A25)')ar__5
        read(ar__5,6)ID_5,rh,f_rh,Inst,Filt
        if(ar__5(8:11) .EQ. '') rh = rNULL__
c    ..............Just test output...........
        write(6,6)ID_5,rh,f_rh,Inst,Filt
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Spectroscopic data for all foreground stars

C  Format for file interpretation

    7 format(
     +  A7,A1,1X,A13,1X,F9.5,1X,F8.5,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,I4,
     +  1X,I3,1X,A23,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,189
        read(1,'(A93)')ar__6
        read(ar__6,7)
     +  ID_6,f_ID_1,AID,RAdeg_4,DEdeg_4,f_imag,imag_3,g_r_3,g_i_3,
     +  HRV_4,e_HRV_4,r_HRV,HST
        if(ar__6(44:48) .EQ. '') imag_3 = rNULL__
        if(ar__6(50:53) .EQ. '') g_r_3 = rNULL__
        if(ar__6(55:58) .EQ. '') g_i_3 = rNULL__
        if(ar__6(65:67) .EQ. '') e_HRV_4 = iNULL__
c    ..............Just test output...........
        write(6,7)
     +  ID_6,f_ID_1,AID,RAdeg_4,DEdeg_4,f_imag,imag_3,g_r_3,g_i_3,
     +  HRV_4,e_HRV_4,r_HRV,HST
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Spectroscopic and photometric data for all M87
*                             globular clusters and ultracompact dwarfs

C  Format for file interpretation

    8 format(
     +  A6,1X,A1,1X,A13,1X,F9.5,1X,F8.5,1X,F5.2,1X,A1,1X,F4.2,1X,F4.2,
     +  1X,I4,1X,I3,1X,A19,1X,F5.1,1X,A1,1X,A3,1X,A4,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,739
        read(1,'(A111)')ar__7
        read(ar__7,8)
     +  ID_7,f_ID_2,AID_1,RAdeg_5,DEdeg_5,imag_4,f_imag_1,g_r_4,g_i_4,
     +  HRV_5,e_HRV_5,r_HRV_1,rh_1,f_rh_1,r_rh,spop,Class_4
        if(ar__7(51:54) .EQ. '') g_r_4 = rNULL__
        if(ar__7(56:59) .EQ. '') g_i_4 = rNULL__
        if(ar__7(90:94) .EQ. '') rh_1 = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  ID_7,f_ID_2,AID_1,RAdeg_5,DEdeg_5,imag_4,f_imag_1,g_r_4,g_i_4,
     +  HRV_5,e_HRV_5,r_HRV_1,rh_1,f_rh_1,r_rh,spop,Class_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Supplementary photometric data for M87 globular
*                             clusters and ultracompact dwarfs

C  Format for file interpretation

    9 format(
     +  A6,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,591
        read(1,'(A44)')ar__8
        read(ar__8,9)ID_8,T1mag,C_T1,gmag,g_z,Vmag_1,B_V_1,V_I_1
        if(ar__8(8:12) .EQ. '') T1mag = rNULL__
        if(ar__8(14:17) .EQ. '') C_T1 = rNULL__
        if(ar__8(19:23) .EQ. '') gmag = rNULL__
        if(ar__8(25:28) .EQ. '') g_z = rNULL__
        if(ar__8(30:34) .EQ. '') Vmag_1 = rNULL__
        if(ar__8(36:39) .EQ. '') B_V_1 = rNULL__
        if(ar__8(41:44) .EQ. '') V_I_1 = rNULL__
c    ..............Just test output...........
        write(6,9)ID_8,T1mag,C_T1,gmag,g_z,Vmag_1,B_V_1,V_I_1
c    .......End.of.Just test output...........
      end do
      close(1)

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