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 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 2026-Oct-16
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         (nr__) ! Identification (G1)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension in decimal degrees (J2000) (1)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000) (1)
      character*3   Mask       (nr__) ! [1,4] Mask(s) (2)
      real*4        Vmag       (nr__) ! (mag) Subaru/Suprime-Cam V band magnitude (1)
      real*4        B_V        (nr__) ! (mag) Subaru/Suprime-Cam (B-V) color index (1)
      real*4        V_I        (nr__) ! (mag) Subaru/Suprime-Cam (V-I) color index (1)
      character*4   Class      (nr__) ! Classification (G2)
      integer*4     HRV        (nr__) ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV      (nr__) ! (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       (nr__1) ! Identification (G1)
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_1     (nr__1) ! [1,4]? Mask (for table 5 only)
      real*4        imag       (nr__1) ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r        (nr__1) ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i        (nr__1) ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_1    (nr__1) ! Classification (G2)
      integer*4     HRV_1      (nr__1) ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_1    (nr__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       (nr__2) ! Identification (G1)
      real*8        RAdeg_2    (nr__2) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_2     (nr__2) ! [1,4]? Mask (for table 5 only)
      real*4        imag_1     (nr__2) ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r_1      (nr__2) ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i_1      (nr__2) ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_2    (nr__2) ! Classification (G2)
      integer*4     HRV_2      (nr__2) ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_2    (nr__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       (nr__3) ! Identification (G1)
      real*8        RAdeg_3    (nr__3) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination in decimal degrees (J2000)
      integer*4     Mask_3     (nr__3) ! [1,4]? Mask (for table 5 only)
      real*4        imag_2     (nr__3) ! (mag) ? The CFHT/Megacam i band magnitude (1)
      real*4        g_r_2      (nr__3) ! (mag) ? The CFHT/Megacam (g-r) color index (1)
      real*4        g_i_2      (nr__3) ! (mag) ? The CFHT/Megacam (g-i) color index (1)
      character*4   Class_3    (nr__3) ! Classification (G2)
      integer*4     HRV_3      (nr__3) ! (km/s) ? Heliocentric radial velocity
      integer*4     e_HRV_3    (nr__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       (nr__4) ! ID Identification (G1)
      character*1   f_ID       (nr__4) ! [b] Flag on ID (1)
      integer*4     H01        (nr__4) ! (km/s) ? Hanes et al. 2001, Cat .J/ApJ/559/812
*                                heliocentric radial velocity
      integer*4     e_H01      (nr__4) ! (km/s) ? Uncertainty in H01
      integer*4     K97        (nr__4) ! (km/s) ? Kohen & Ryzhov, 1997ApJ...486..230C,
*                                heliocentric radial velocity
      integer*4     e_K97      (nr__4) ! (km/s) ? Uncertainty in K97
      integer*4     K00        (nr__4) ! (km/s) ? Cohen, 2000AJ....119..162C, heliocentric
*                                radial velocity
      integer*4     e_K00      (nr__4) ! (km/s) ? Uncertainty in K00
      integer*4     M90        (nr__4) ! (km/s) ? Mould et al., 1990AJ.....99.1823M,
*                                heliocentric radial velocity
      integer*4     e_M90      (nr__4) ! (km/s) ? Uncertainty in M90
      integer*4     M87        (nr__4) ! (km/s) ? Mould et al., 1987AJ.....93...53M,
*                                heliocentric radial velocity
      integer*4     e_M87      (nr__4) ! (km/s) ? Uncertainty in M87
      integer*4     HB87       (nr__4) ! (km/s) ? Huchra & Brodie, 1987AJ.....93..779H,
*                                heliocentric radial velocity
      integer*4     e_HB87     (nr__4) ! (km/s) ? Uncertainty in HB87
      integer*4     MMT        (nr__4) ! (km/s) ? This paper, MMT/Hectospec heliocentric radial
*                                velocity
      integer*4     e_MMT      (nr__4) ! (km/s) ? Uncertainty in MMT
      integer*4     LRIS       (nr__4) ! (km/s) ? This paper, Keck/LRIS heliocentric radial
*                                velocity
      integer*4     e_LRIS     (nr__4) ! (km/s) ? Uncertainty in LRIS
      integer*4     D1         (nr__4) ! (km/s) ? This paper, Keck/DEIMOS run 1 heliocentric
*                                radial velocity
      integer*4     e_D1       (nr__4) ! (km/s) ? Uncertainty in D1
      integer*4     D2         (nr__4) ! (km/s) ? This paper, Keck/DEIMOS run 2 heliocentric
*                                radial velocity
      integer*4     e_D2       (nr__4) ! (km/s) ? Uncertainty in D2
      integer*4     H05_H07    (nr__4) ! (km/s) ? Hasegan et al., 2005ApJ...627..203H, or
*                                Hasegan 2007PhDT.........4H, heliocentric radial
*                                velocity
      integer*4     e_H05_H07  (nr__4) ! (km/s) ? Uncertainty in H05/H07
      integer*4     J06        (nr__4) ! (km/s) ? Jones et al., 2006AJ....131..312J,
*                                heliocentric radial velocity
      integer*4     e_J06      (nr__4) ! (km/s) ? Uncertainty in J06
      integer*4     E07        (nr__4) ! (km/s) ? Evstigneeva et al., 2007AJ....133.1722E,
*                                heliocentric radial velocity
      integer*4     e_E07      (nr__4) ! (km/s) ? Uncertainty in E07
      integer*4     F08        (nr__4) ! (km/s) ? Firth et al., 2008MNRAS.389.1539F,
*                                heliocentric radial velocity
      integer*4     e_F08      (nr__4) ! (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       (nr__5) ! Star name (G1)
      real*4        rh         (nr__5) ! (pc) ? Half-light radius assuming a distance of 16.5Mpc
      character*1   f_rh       (nr__5) ! [d] d: unresolved; the upper limit is r_h_~1pc
      character*5   Inst       (nr__5) ! Instrument
      character*6   Filt       (nr__5) ! 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       (nr__6) ! Star name (G1)
      character*1   f_ID_1     (nr__6) ! [f] Flag on ID (1)
      character*13  AID        (nr__6) ! Alternative identification(s) (G1)
      real*8        RAdeg_4    (nr__6) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_4    (nr__6) ! (deg) Declination in decimal degrees (J2000)
      character*1   f_imag     (nr__6) ! [i] Flag on imag (2)
      real*4        imag_3     (nr__6) ! (mag) ? The CFHT i-band magnitude
      real*4        g_r_3      (nr__6) ! (mag) ? The CFHT (g-r) color index
      real*4        g_i_3      (nr__6) ! (mag) ? The CFHT (g-i) color index
      integer*4     HRV_4      (nr__6) ! (km/s) Heliocentric radial velocity of the object
      integer*4     e_HRV_4    (nr__6) ! (km/s) ? HRV uncertainty
      character*23  r_HRV      (nr__6) ! Source(s) (3)
      character*1   HST        (nr__6) ! [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       (nr__7) ! Identification (G1)
      character*1   f_ID_2     (nr__7) ! [i] Flag on ID (1)
      character*13  AID_1      (nr__7) ! Alternate identification (G1)
      real*8        RAdeg_5    (nr__7) ! (deg) Right Ascension in decimal degrees (J2000) (2)
      real*8        DEdeg_5    (nr__7) ! (deg) Declination in decimal degrees (J2000) (2)
      real*4        imag_4     (nr__7) ! (mag) The CFHT i band magnitude
      character*1   f_imag_1   (nr__7) ! [i] Flag on imag (3)
      real*4        g_r_4      (nr__7) ! (mag) ? The CFHT (g-r) color index
      real*4        g_i_4      (nr__7) ! (mag) ? The CFHT (g-i) color index
      integer*4     HRV_5      (nr__7) ! (km/s) Heliocenric radial velocity
      integer*4     e_HRV_5    (nr__7) ! (km/s) Error in HRV
      character*19  r_HRV_1    (nr__7) ! Source of HRV (4)
      real*4        rh_1       (nr__7) ! (pc) ? Half-light radius; assuming D=16.5Mpc
      character*1   f_rh_1     (nr__7) ! [u] Indicates unresolved source for rh
      character*3   r_rh       (nr__7) ! Source of rh (5)
      character*4   spop       (nr__7) ! Cluster subpopulation (blue or red) (6)
      character*5   Class_4    (nr__7) ! 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       (nr__8) ! Identification (G1)
      real*4        T1mag      (nr__8) ! (mag) ? The Washington T_1_ band magnitude (1)
      real*4        C_T1       (nr__8) ! (mag) ? The Washington (C-T_1_) color index (1)
      real*4        gmag       (nr__8) ! (mag) ? HST/ACS g band magnitude (2)
      real*4        g_z        (nr__8) ! (mag) ? HST/ACS (g-z) color index (2)
      real*4        Vmag_1     (nr__8) ! (mag) ? The Subaru V band magnitude
      real*4        B_V_1      (nr__8) ! (mag) ? The Subaru (B-V) color index
      real*4        V_I_1      (nr__8) ! (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(i__),RAdeg(i__),DEdeg(i__),Mask(i__),Vmag(i__),B_V(i__),
     +  V_I(i__),Class(i__),HRV(i__),e_HRV(i__)
        if(ar__(52:55) .EQ. '') HRV(i__) = iNULL__
        if(ar__(57:58) .EQ. '') e_HRV(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),RAdeg(i__),DEdeg(i__),Mask(i__),Vmag(i__),B_V(i__),
     +  V_I(i__),Class(i__),HRV(i__),e_HRV(i__)
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(i__),RAdeg_1(i__),DEdeg_1(i__),Mask_1(i__),imag(i__),
     +  g_r(i__),g_i(i__),Class_1(i__),HRV_1(i__),e_HRV_1(i__)
        if(ar__1(27:27) .EQ. '') Mask_1(i__) = iNULL__
        if(ar__1(29:33) .EQ. '') imag(i__) = rNULL__
        if(ar__1(35:38) .EQ. '') g_r(i__) = rNULL__
        if(ar__1(40:43) .EQ. '') g_i(i__) = rNULL__
        if(ar__1(50:53) .EQ. '') HRV_1(i__) = iNULL__
        if(ar__1(55:56) .EQ. '') e_HRV_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),Mask_1(i__),imag(i__),
     +  g_r(i__),g_i(i__),Class_1(i__),HRV_1(i__),e_HRV_1(i__)
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(i__),RAdeg_2(i__),DEdeg_2(i__),Mask_2(i__),imag_1(i__),
     +  g_r_1(i__),g_i_1(i__),Class_2(i__),HRV_2(i__),e_HRV_2(i__)
        if(ar__2(27:27) .EQ. '') Mask_2(i__) = iNULL__
        if(ar__2(29:33) .EQ. '') imag_1(i__) = rNULL__
        if(ar__2(35:38) .EQ. '') g_r_1(i__) = rNULL__
        if(ar__2(40:43) .EQ. '') g_i_1(i__) = rNULL__
        if(ar__2(50:53) .EQ. '') HRV_2(i__) = iNULL__
        if(ar__2(55:56) .EQ. '') e_HRV_2(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),RAdeg_2(i__),DEdeg_2(i__),Mask_2(i__),imag_1(i__),
     +  g_r_1(i__),g_i_1(i__),Class_2(i__),HRV_2(i__),e_HRV_2(i__)
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(i__),RAdeg_3(i__),DEdeg_3(i__),Mask_3(i__),imag_2(i__),
     +  g_r_2(i__),g_i_2(i__),Class_3(i__),HRV_3(i__),e_HRV_3(i__)
        if(ar__3(27:27) .EQ. '') Mask_3(i__) = iNULL__
        if(ar__3(29:33) .EQ. '') imag_2(i__) = rNULL__
        if(ar__3(35:38) .EQ. '') g_r_2(i__) = rNULL__
        if(ar__3(40:43) .EQ. '') g_i_2(i__) = rNULL__
        if(ar__3(50:53) .EQ. '') HRV_3(i__) = iNULL__
        if(ar__3(55:56) .EQ. '') e_HRV_3(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  ID_3(i__),RAdeg_3(i__),DEdeg_3(i__),Mask_3(i__),imag_2(i__),
     +  g_r_2(i__),g_i_2(i__),Class_3(i__),HRV_3(i__),e_HRV_3(i__)
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(i__),f_ID(i__),H01(i__),e_H01(i__),K97(i__),e_K97(i__),
     +  K00(i__),e_K00(i__),M90(i__),e_M90(i__),M87(i__),e_M87(i__),
     +  HB87(i__),e_HB87(i__),MMT(i__),e_MMT(i__),LRIS(i__),
     +  e_LRIS(i__),D1(i__),e_D1(i__),D2(i__),e_D2(i__),H05_H07(i__),
     +  e_H05_H07(i__),J06(i__),e_J06(i__),E07(i__),e_E07(i__),
     +  F08(i__),e_F08(i__)
        if(ar__4(10:13) .EQ. '') H01(i__) = iNULL__
        if(ar__4(15:17) .EQ. '') e_H01(i__) = iNULL__
        if(ar__4(19:22) .EQ. '') K97(i__) = iNULL__
        if(ar__4(24:26) .EQ. '') e_K97(i__) = iNULL__
        if(ar__4(28:31) .EQ. '') K00(i__) = iNULL__
        if(ar__4(33:34) .EQ. '') e_K00(i__) = iNULL__
        if(ar__4(36:39) .EQ. '') M90(i__) = iNULL__
        if(ar__4(41:43) .EQ. '') e_M90(i__) = iNULL__
        if(ar__4(45:48) .EQ. '') M87(i__) = iNULL__
        if(ar__4(50:52) .EQ. '') e_M87(i__) = iNULL__
        if(ar__4(54:57) .EQ. '') HB87(i__) = iNULL__
        if(ar__4(59:61) .EQ. '') e_HB87(i__) = iNULL__
        if(ar__4(63:66) .EQ. '') MMT(i__) = iNULL__
        if(ar__4(68:69) .EQ. '') e_MMT(i__) = iNULL__
        if(ar__4(71:74) .EQ. '') LRIS(i__) = iNULL__
        if(ar__4(76:77) .EQ. '') e_LRIS(i__) = iNULL__
        if(ar__4(79:82) .EQ. '') D1(i__) = iNULL__
        if(ar__4(84:85) .EQ. '') e_D1(i__) = iNULL__
        if(ar__4(87:90) .EQ. '') D2(i__) = iNULL__
        if(ar__4(92:93) .EQ. '') e_D2(i__) = iNULL__
        if(ar__4(95:98) .EQ. '') H05_H07(i__) = iNULL__
        if(ar__4(100:101) .EQ. '') e_H05_H07(i__) = iNULL__
        if(ar__4(103:106) .EQ. '') J06(i__) = iNULL__
        if(ar__4(108:110) .EQ. '') e_J06(i__) = iNULL__
        if(ar__4(112:115) .EQ. '') E07(i__) = iNULL__
        if(ar__4(117:117) .EQ. '') e_E07(i__) = iNULL__
        if(ar__4(119:122) .EQ. '') F08(i__) = iNULL__
        if(ar__4(124:125) .EQ. '') e_F08(i__) = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  ID_4(i__),f_ID(i__),H01(i__),e_H01(i__),K97(i__),e_K97(i__),
     +  K00(i__),e_K00(i__),M90(i__),e_M90(i__),M87(i__),e_M87(i__),
     +  HB87(i__),e_HB87(i__),MMT(i__),e_MMT(i__),LRIS(i__),
     +  e_LRIS(i__),D1(i__),e_D1(i__),D2(i__),e_D2(i__),H05_H07(i__),
     +  e_H05_H07(i__),J06(i__),e_J06(i__),E07(i__),e_E07(i__),
     +  F08(i__),e_F08(i__)
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(i__),rh(i__),f_rh(i__),Inst(i__),Filt(i__)
        if(ar__5(8:11) .EQ. '') rh(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)ID_5(i__),rh(i__),f_rh(i__),Inst(i__),Filt(i__)
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(i__),f_ID_1(i__),AID(i__),RAdeg_4(i__),DEdeg_4(i__),
     +  f_imag(i__),imag_3(i__),g_r_3(i__),g_i_3(i__),HRV_4(i__),
     +  e_HRV_4(i__),r_HRV(i__),HST(i__)
        if(ar__6(44:48) .EQ. '') imag_3(i__) = rNULL__
        if(ar__6(50:53) .EQ. '') g_r_3(i__) = rNULL__
        if(ar__6(55:58) .EQ. '') g_i_3(i__) = rNULL__
        if(ar__6(65:67) .EQ. '') e_HRV_4(i__) = iNULL__
c    ..............Just test output...........
        write(6,7)
     +  ID_6(i__),f_ID_1(i__),AID(i__),RAdeg_4(i__),DEdeg_4(i__),
     +  f_imag(i__),imag_3(i__),g_r_3(i__),g_i_3(i__),HRV_4(i__),
     +  e_HRV_4(i__),r_HRV(i__),HST(i__)
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(i__),f_ID_2(i__),AID_1(i__),RAdeg_5(i__),DEdeg_5(i__),
     +  imag_4(i__),f_imag_1(i__),g_r_4(i__),g_i_4(i__),HRV_5(i__),
     +  e_HRV_5(i__),r_HRV_1(i__),rh_1(i__),f_rh_1(i__),r_rh(i__),
     +  spop(i__),Class_4(i__)
        if(ar__7(51:54) .EQ. '') g_r_4(i__) = rNULL__
        if(ar__7(56:59) .EQ. '') g_i_4(i__) = rNULL__
        if(ar__7(90:94) .EQ. '') rh_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  ID_7(i__),f_ID_2(i__),AID_1(i__),RAdeg_5(i__),DEdeg_5(i__),
     +  imag_4(i__),f_imag_1(i__),g_r_4(i__),g_i_4(i__),HRV_5(i__),
     +  e_HRV_5(i__),r_HRV_1(i__),rh_1(i__),f_rh_1(i__),r_rh(i__),
     +  spop(i__),Class_4(i__)
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(i__),T1mag(i__),C_T1(i__),gmag(i__),g_z(i__),Vmag_1(i__),
     +  B_V_1(i__),V_I_1(i__)
        if(ar__8(8:12) .EQ. '') T1mag(i__) = rNULL__
        if(ar__8(14:17) .EQ. '') C_T1(i__) = rNULL__
        if(ar__8(19:23) .EQ. '') gmag(i__) = rNULL__
        if(ar__8(25:28) .EQ. '') g_z(i__) = rNULL__
        if(ar__8(30:34) .EQ. '') Vmag_1(i__) = rNULL__
        if(ar__8(36:39) .EQ. '') B_V_1(i__) = rNULL__
        if(ar__8(41:44) .EQ. '') V_I_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  ID_8(i__),T1mag(i__),C_T1(i__),gmag(i__),g_z(i__),Vmag_1(i__),
     +  B_V_1(i__),V_I_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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