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

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/517/A81 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.5, on 2013-May-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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/A+A/517/A81  2nd and 3d parameters of HB of globular clusters (Gratton+, 2010)
*================================================================================
*The second and third parameters of the horizontal branch of globular clusters.
*    Gratton R.G., Carretta E., Bragaglia A., Lucatello S., S'Orazii V.
*   <Astron. Astrophys. 517, A81 (2010)>
*   =2010A&A...517A..81G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Basic parameters for globular clusters considered
                           in this paper

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

      character*7   NGC_IC     (nr__) ! Cluster name
      character*7   Name       (nr__) ! Other name
      real*4        v_Fe_H_    (nr__) ! ([Sun]) Metallicity
      real*8        v_m_M_V    (nr__) ! (mag) Distance modulus calculated from V magnitude
      real*8        E_B_V      (nr__) ! (mag) Redenning
      real*8        VMAG       (nr__) ! (mag) Absolute V magnitude
      real*8        HBR        (nr__) ! ? Horizontal-branch ratio = (B-R)/(B+V+R)
      real*8        Age        (nr__) ! (Gyr) ? Estimated cluster age
      real*8        M_TO       (nr__) ! (solMass) ? Mass at turn-off
      real*8        M_RGB      (nr__) ! (solMass) ? Mass at the tip of the red-giant branch
      real*8        IQR1       (nr__) ! ? Interquartile of the [Na/O]
*                                     (anti-)correlations
      real*8        IQR2       (nr__) ! ? Interquartile of the [Mg/Al]
*                                     (anti-)correlations
      real*8        T_HB       (nr__) ! ([K]) ? Horizontal-branch effective temperature
*                                     (log scale)

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

C  Declarations for 'table2.dat'	! Photometric data from the HST-snapshot
                           (Piotto et al., 2002, Cat. J/A+A/391/945)

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

      character*7   NGC_IC_1   (nr__1) ! Cluster name
      integer*4     n_HB       (nr__1) ! Number of Horizon Branch stars
      integer*4     n_fld      (nr__1) ! [0/125] Number of field stars (G1)
      real*4        F555       (nr__1) ! (mag) Median F555W magnitude
      real*4        e_F555     (nr__1) ! (mag) rms uncertainty on F555
      real*4        v439_555   (nr__1) ! (mag) Median F439W-F555W colour index
      real*4        e_439_555  (nr__1) ! (mag) rms uncertainty on F439-F555
      real*4        F555m      (nr__1) ! (mag) Minimum F555W magnitude
      real*4        e_F555m    (nr__1) ! (mag) rms uncertainty on F555m
      real*4        v439_555m  (nr__1) ! (mag) Minimum F439W-F555W colour index
      real*4        e_439_555m (nr__1) ! (mag) rms uncertainty on (F439-F555)m
      real*4        v439_555M_1(nr__1) ! (mag) Maximum F439W-F555W colour index
      real*4        e_439_555M_1(nr__1) ! (mag) rms uncertainty on (F439-F555)M
      integer*4     n_RGB      (nr__1) ! Number of RGB stars
      real*4        v439_555G  (nr__1) ! (mag) F439W-F555W colour index for RGB
      real*4        e_439_555G (nr__1) ! (mag) rms uncertainty on (F439-F555)RGB

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

C  Declarations for 'table3.dat'	! Photometric data from ground-based data (Rosenberg
                           et al. 2000, Cat. J/A+AS/144/5 and J/A+AS/145/451)

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

      character*7   NGC_IC_2   (nr__2) ! Cluster name
      integer*4     n_HB_1     (nr__2) ! Number of HB stars in the cluster
      integer*4     n_fld_1    (nr__2) ! [0/47] Number of field stars (G1)
      real*4        Vmag_1     (nr__2) ! (mag) Median V magnitude
      real*4        e_Vmag     (nr__2) ! (mag) rms uncertainty on
      real*4        V_I        (nr__2) ! (mag) Median V-I colour index
      real*4        e_V_I      (nr__2) ! (mag) rms uncertainty on V-I
      real*4        Vmagm      (nr__2) ! (mag) Minimum V magnitude
      real*4        e_Vmagm    (nr__2) ! (mag) rms uncertainty on
      real*4        v_V_I_m    (nr__2) ! (mag) Minimum V magnitude
      real*4        e__V_I_m   (nr__2) ! (mag) rms uncertainty on (V-I)m
      real*4        v_V_I_M_1  (nr__2) ! (mag) Maximum V-I colour index
      real*4        e__V_I_M_1 (nr__2) ! (mag) rms uncertainty on (V-I)M
      integer*4     n_RGB_1    (nr__2) ! Number of RGB stars
      real*4        v_V_I_G    (nr__2) ! (mag) V-I colour index for RGB
      real*4        e__V_I_G   (nr__2) ! (mag) rms uncertainty on (V-I)RGB

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

C  Declarations for 'table4.dat'	! Intrinsic HB photometric parameters for GCs
                           with HST data

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

      character*7   NGC_IC_3   (nr__3) ! Cluster name
      real*4        v_B_V_0G   (nr__3) ! (mag) Intrinsic B-V colour index value, for RGB
      real*4        v_B_V_0m   (nr__3) ! (mag) Minimim B-V colour index value
      real*4        v_B_V_0    (nr__3) ! (mag) Median B-V colour index value
      real*4        v_B_V_0M_1 (nr__3) ! (mag) Maximum B-V colour index value
      real*4        VMAGM_1    (nr__3) ! (mag) Maximum absolute V magnitude
      real*4        VMAG_2     (nr__3) ! (mag) Median absolute V magnitude

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

C  Declarations for 'table5.dat'	! Intrinsic HB photometric parameters for GCs
                           with ground-based data

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

      character*7   NGC_IC_4   (nr__4) ! Cluster name
      real*4        v_V_I_0G   (nr__4) ! (mag) Intrinsic V-I colour index value, for RGB
      real*4        v_V_I_0m   (nr__4) ! (mag) Mimimum V-I colour index value
      real*4        v_V_I_0    (nr__4) ! (mag) Median V-I colour index value
      real*4        v_V_I_0M_1 (nr__4) ! (mag) Maximum V-I colour index value
      real*4        VMAGM_2    (nr__4) ! (mag) Maximum absolute V magnitude
      real*4        VMAG_3     (nr__4) ! (mag) Median absolute V magnitude

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

C  Declarations for 'table6.dat'	! Masses on HB for GCs with HST data

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

      character*7   NGC_IC_5   (nr__5) ! Cluster name
      real*4        Mmin       (nr__5) ! (solMass) Mimimum mass
      real*4        e_Mmin     (nr__5) ! (solMass) ? rms uncertainty on Mmin
      real*4        Mmed       (nr__5) ! (solMass) Median mass
      real*4        e_Mmed     (nr__5) ! (solMass) ? rms uncertainty on Mmed
      real*4        Mmax       (nr__5) ! (solMass) ? Maximum mass
      real*4        e_Mmax     (nr__5) ! (solMass) ? rms uncertainty on Mmax
      real*4        dMmed      (nr__5) ! (solMass) ? Difference between median and minimum masses

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

C  Declarations for 'table7.dat'	! Masses on HB for GCs with ground-based data

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

      character*7   NGC_IC_6   (nr__6) ! Cluster name
      real*4        Mmin_1     (nr__6) ! (solMass) Mimimum mass
      real*4        e_Mmin_1   (nr__6) ! (solMass) ? rms uncertainty on Mmin
      real*4        Mmed_1     (nr__6) ! (solMass) Median mass
      real*4        e_Mmed_1   (nr__6) ! (solMass) ? rms uncertainty on Mmed
      real*4        Mmax_1     (nr__6) ! (solMass) ? Maximum mass
      real*4        e_Mmax_1   (nr__6) ! (solMass) ? rms uncertainty on Mmax
      real*4        dMmed_1    (nr__6) ! (solMass) ? Difference between median and minimum masses

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

C  Declarations for 'table8.dat'	! Properties of the RR Lyrae variables for each cluster

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

      character*7   NGC_IC_7   (nr__7) ! Cluster name
      real*4        S_RR       (nr__7) ! ? RR Lyrae richness parameter (Harris, 1996,
*                                    Cat. VII/195) (1)
      integer*4     n_RR_H     (nr__7) ! ? Number of RR Lyrae stars for HST sample
      real*4        f_RR_H     (nr__7) ! ? Fraction of RR Lyrae stars for HST sample
      integer*4     n_RR_g     (nr__7) ! ? Number of RR Lyrae for ground-based sample
      real*4        f_RR_g     (nr__7) ! ? Fraction of RR Lyrae for ground-based sample
      real*4        f_RR       (nr__7) ! ? Fraction of RR Lyrae stars
      real*4        v_Pab_     (nr__7) ! (d) ? Mean period
      integer*4     n_Pab      (nr__7) ! ? Number of stars used to determine the period
*Note (1): Number of RR Lyrae per unit cluster luminosity

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

C  Declarations for 'table10.dat'	! Clusters ages from main sequence fitting and
                           from HB both for HST and ground based data-set

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

      character*7   NGC_IC_8   (nr__8) ! Cluster name
      real*8        Age_1      (nr__8) ! (Gyr) ? Age from main-sequence fitting
      real*8        Age_HB_H   (nr__8) ! (Gyr) ? Age from horizontal branch fitting,
*                                      for HST data-set
      real*8        Age_HB_g   (nr__8) ! (Gyr) ? Age from horizontal branch fitting,
*                                      for ground-based data-set
      real*8        Age_Y      (nr__8) ! (Gyr) ? Age value, determined with Ymed

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

C  Declarations for 'table11.dat'	! R' parameter as derived from both HST and
                           ground based data

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

      character*7   NGC_IC_9   (nr__9) ! Cluster name
      real*8        R__HST     (nr__9) ! ? R' parameter as derived from HST data (1)
      real*4        e_R__HST   (nr__9) ! ? rms uncertainty on R'(HST)
      real*8        R__GR      (nr__9) ! ? R' parameter as derived from ground-based
*                                    data (1)
      real*4        e_R__GR    (nr__9) ! ? rms uncertainty on R'(GR)
      real*8        v_R__      (nr__9) ! ? Mean R' parameter (1) (2)
      real*4        e__R__     (nr__9) ! ? rms uncertainty on <R'>
      real*8        Y__R__     (nr__9) ! ? He abundance for R'=<R'>
      real*8        e_Y__R__   (nr__9) ! ? rms uncertainty on Y(<R'>)
*Note (1): R' = N_HB_/N'_RGB_, where N'_RGB_ is the number of stars on the
*     RGB brighter than V(HB)+1.
*Note (2): Final adopted value <R'> is the weighted average between the two
*     estimates, when available.

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

C  Declarations for 'table12.dat'	! Spread in He and mass required explaining the HB

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

      character*7   NGC_IC_10  (nr__10) ! Cluster name
      real*8        dM_Y       (nr__10) ! (solMass) Spread in mass {delta}M(Y)
      real*8        Ymed       (nr__10) ! Median He abundance value
      real*8        Ymax       (nr__10) ! Maximum He abundance value
      real*8        dY         (nr__10) ! Difference between maximum and median Y values
      real*8        d_B_V      (nr__10) ! (mag) Colour spread in {delta}(B-V)
      real*8        d_V_I      (nr__10) ! (mag) Colour spread in {delta}(V-I)

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

C  Loading file 'table1.dat'	! Basic parameters for globular clusters considered
*                           in this paper

C  Format for file interpretation

    1 format(
     +  A7,2X,A7,F6.2,F7.2,F7.2,F7.2,F8.2,F7.2,F8.3,F8.3,F7.3,F8.3,
     +  F9.3)

C  Effective file loading

      open(unit=1,file='table1.dat', status='old')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,94
        read(1,'(A98)')ar__
        read(ar__,1)
     +  NGC_IC(i__),Name(i__),v_Fe_H_(i__),v_m_M_V(i__),E_B_V(i__),
     +  VMAG(i__),HBR(i__),Age(i__),M_TO(i__),M_RGB(i__),IQR1(i__),
     +  IQR2(i__),T_HB(i__)
        if(ar__(44:51) .EQ. '') HBR(i__) = rNULL__
        if(ar__(52:58) .EQ. '') Age(i__) = rNULL__
        if(ar__(59:66) .EQ. '') M_TO(i__) = rNULL__
        if(ar__(67:74) .EQ. '') M_RGB(i__) = rNULL__
        if(ar__(75:81) .EQ. '') IQR1(i__) = rNULL__
        if(ar__(82:89) .EQ. '') IQR2(i__) = rNULL__
        if(ar__(90:98) .EQ. '') T_HB(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  NGC_IC(i__),Name(i__),v_Fe_H_(i__),v_m_M_V(i__),E_B_V(i__),
     +  VMAG(i__),HBR(i__),Age(i__),M_TO(i__),M_RGB(i__),IQR1(i__),
     +  IQR2(i__),T_HB(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Photometric data from the HST-snapshot
*                           (Piotto et al., 2002, Cat. J/A+A/391/945)

C  Format for file interpretation

    2 format(
     +  A7,1X,I4,1X,I3,3X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,2X,F6.3,2X,
     +  F5.3,1X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,1X,I4,1X,F6.3,2X,F5.3)

C  Effective file loading

      open(unit=1,file='table2.dat', status='old')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,74
        read(1,'(A108)')ar__1
        read(ar__1,2)
     +  NGC_IC_1(i__),n_HB(i__),n_fld(i__),F555(i__),e_F555(i__),
     +  v439_555(i__),e_439_555(i__),F555m(i__),e_F555m(i__),
     +  v439_555m(i__),e_439_555m(i__),v439_555M_1(i__),
     +  e_439_555M_1(i__),n_RGB(i__),v439_555G(i__),e_439_555G(i__)
c    ..............Just test output...........
        write(6,2)
     +  NGC_IC_1(i__),n_HB(i__),n_fld(i__),F555(i__),e_F555(i__),
     +  v439_555(i__),e_439_555(i__),F555m(i__),e_F555m(i__),
     +  v439_555m(i__),e_439_555m(i__),v439_555M_1(i__),
     +  e_439_555M_1(i__),n_RGB(i__),v439_555G(i__),e_439_555G(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Photometric data from ground-based data (Rosenberg
*                           et al. 2000, Cat. J/A+AS/144/5 and J/A+AS/145/451)

C  Format for file interpretation

    3 format(
     +  A7,I4,1X,I3,4X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,2X,F6.3,2X,F5.3,
     +  1X,F6.3,2X,F5.3,1X,F6.3,2X,F5.3,1X,I4,1X,F6.3,2X,F5.3)

C  Effective file loading

      open(unit=1,file='table3.dat', status='old')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,49
        read(1,'(A108)')ar__2
        read(ar__2,3)
     +  NGC_IC_2(i__),n_HB_1(i__),n_fld_1(i__),Vmag_1(i__),
     +  e_Vmag(i__),V_I(i__),e_V_I(i__),Vmagm(i__),e_Vmagm(i__),
     +  v_V_I_m(i__),e__V_I_m(i__),v_V_I_M_1(i__),e__V_I_M_1(i__),
     +  n_RGB_1(i__),v_V_I_G(i__),e__V_I_G(i__)
c    ..............Just test output...........
        write(6,3)
     +  NGC_IC_2(i__),n_HB_1(i__),n_fld_1(i__),Vmag_1(i__),
     +  e_Vmag(i__),V_I(i__),e_V_I(i__),Vmagm(i__),e_Vmagm(i__),
     +  v_V_I_m(i__),e__V_I_m(i__),v_V_I_M_1(i__),e__V_I_M_1(i__),
     +  n_RGB_1(i__),v_V_I_G(i__),e__V_I_G(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Intrinsic HB photometric parameters for GCs
*                           with HST data

C  Format for file interpretation

    4 format(A7,3X,F5.3,3X,F6.3,3X,F6.3,2X,F6.3,3X,F5.3,3X,F5.3)

C  Effective file loading

      open(unit=1,file='table4.dat', status='old')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,74
        read(1,'(A57)')ar__3
        read(ar__3,4)
     +  NGC_IC_3(i__),v_B_V_0G(i__),v_B_V_0m(i__),v_B_V_0(i__),
     +  v_B_V_0M_1(i__),VMAGM_1(i__),VMAG_2(i__)
c    ..............Just test output...........
        write(6,4)
     +  NGC_IC_3(i__),v_B_V_0G(i__),v_B_V_0m(i__),v_B_V_0(i__),
     +  v_B_V_0M_1(i__),VMAGM_1(i__),VMAG_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Intrinsic HB photometric parameters for GCs
*                           with ground-based data

C  Format for file interpretation

    5 format(A7,2X,F6.3,3X,F6.3,3X,F6.3,2X,F6.3,3X,F5.3,3X,F5.3)

C  Effective file loading

      open(unit=1,file='table5.dat', status='old')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,49
        read(1,'(A57)')ar__4
        read(ar__4,5)
     +  NGC_IC_4(i__),v_V_I_0G(i__),v_V_I_0m(i__),v_V_I_0(i__),
     +  v_V_I_0M_1(i__),VMAGM_2(i__),VMAG_3(i__)
c    ..............Just test output...........
        write(6,5)
     +  NGC_IC_4(i__),v_V_I_0G(i__),v_V_I_0m(i__),v_V_I_0(i__),
     +  v_V_I_0M_1(i__),VMAGM_2(i__),VMAG_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Masses on HB for GCs with HST data

C  Format for file interpretation

    6 format(
     +  A7,3X,F5.3,2X,F5.3,3X,F5.3,2X,F5.3,3X,F5.3,3X,F5.3,3X,F5.3)

C  Effective file loading

      open(unit=1,file='table6.dat', status='old')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,74
        read(1,'(A61)')ar__5
        read(ar__5,6)
     +  NGC_IC_5(i__),Mmin(i__),e_Mmin(i__),Mmed(i__),e_Mmed(i__),
     +  Mmax(i__),e_Mmax(i__),dMmed(i__)
        if(ar__5(18:22) .EQ. '') e_Mmin(i__) = rNULL__
        if(ar__5(33:37) .EQ. '') e_Mmed(i__) = rNULL__
        if(ar__5(41:45) .EQ. '') Mmax(i__) = rNULL__
        if(ar__5(49:53) .EQ. '') e_Mmax(i__) = rNULL__
        if(ar__5(57:61) .EQ. '') dMmed(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  NGC_IC_5(i__),Mmin(i__),e_Mmin(i__),Mmed(i__),e_Mmed(i__),
     +  Mmax(i__),e_Mmax(i__),dMmed(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Masses on HB for GCs with ground-based data

C  Format for file interpretation

    7 format(
     +  A7,3X,F5.3,2X,F5.3,3X,F5.3,2X,F5.3,3X,F5.3,3X,F5.3,3X,F5.3)

C  Effective file loading

      open(unit=1,file='table7.dat', status='old')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,49
        read(1,'(A61)')ar__6
        read(ar__6,7)
     +  NGC_IC_6(i__),Mmin_1(i__),e_Mmin_1(i__),Mmed_1(i__),
     +  e_Mmed_1(i__),Mmax_1(i__),e_Mmax_1(i__),dMmed_1(i__)
        if(ar__6(18:22) .EQ. '') e_Mmin_1(i__) = rNULL__
        if(ar__6(33:37) .EQ. '') e_Mmed_1(i__) = rNULL__
        if(ar__6(41:45) .EQ. '') Mmax_1(i__) = rNULL__
        if(ar__6(49:53) .EQ. '') e_Mmax_1(i__) = rNULL__
        if(ar__6(57:61) .EQ. '') dMmed_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  NGC_IC_6(i__),Mmin_1(i__),e_Mmin_1(i__),Mmed_1(i__),
     +  e_Mmed_1(i__),Mmax_1(i__),e_Mmax_1(i__),dMmed_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Properties of the RR Lyrae variables for each cluster

C  Format for file interpretation

    8 format(
     +  A7,2X,F5.1,3X,I2,7X,F4.2,6X,I3,8X,F4.2,4X,F4.2,2X,F5.3,3X,I4)

C  Effective file loading

      open(unit=1,file='table8.dat', status='old')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,66
        read(1,'(A73)')ar__7
        read(ar__7,8)
     +  NGC_IC_7(i__),S_RR(i__),n_RR_H(i__),f_RR_H(i__),n_RR_g(i__),
     +  f_RR_g(i__),f_RR(i__),v_Pab_(i__),n_Pab(i__)
        if(ar__7(10:14) .EQ. '') S_RR(i__) = rNULL__
        if(ar__7(18:19) .EQ. '') n_RR_H(i__) = iNULL__
        if(ar__7(27:30) .EQ. '') f_RR_H(i__) = rNULL__
        if(ar__7(37:39) .EQ. '') n_RR_g(i__) = iNULL__
        if(ar__7(48:51) .EQ. '') f_RR_g(i__) = rNULL__
        if(ar__7(56:59) .EQ. '') f_RR(i__) = rNULL__
        if(ar__7(62:66) .EQ. '') v_Pab_(i__) = rNULL__
        if(ar__7(70:73) .EQ. '') n_Pab(i__) = iNULL__
c    ..............Just test output...........
        write(6,8)
     +  NGC_IC_7(i__),S_RR(i__),n_RR_H(i__),f_RR_H(i__),n_RR_g(i__),
     +  f_RR_g(i__),f_RR(i__),v_Pab_(i__),n_Pab(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Clusters ages from main sequence fitting and
*                           from HB both for HST and ground based data-set

C  Format for file interpretation

    9 format(A7,2X,F7.2,F8.3,F8.3,F8.3)

C  Effective file loading

      open(unit=1,file='table10.dat', status='old')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,94
        read(1,'(A40)')ar__8
        read(ar__8,9)
     +  NGC_IC_8(i__),Age_1(i__),Age_HB_H(i__),Age_HB_g(i__),
     +  Age_Y(i__)
        if(ar__8(10:16) .EQ. '') Age_1(i__) = rNULL__
        if(ar__8(17:24) .EQ. '') Age_HB_H(i__) = rNULL__
        if(ar__8(25:32) .EQ. '') Age_HB_g(i__) = rNULL__
        if(ar__8(33:40) .EQ. '') Age_Y(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  NGC_IC_8(i__),Age_1(i__),Age_HB_H(i__),Age_HB_g(i__),
     +  Age_Y(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table11.dat'	! R' parameter as derived from both HST and
*                           ground based data

C  Format for file interpretation

   10 format(A7,F7.2,F6.2,F7.2,F6.2,F7.2,F6.2,F8.3,F7.3)

C  Effective file loading

      open(unit=1,file='table11.dat', status='old')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,94
        read(1,'(A61)')ar__9
        read(ar__9,10)
     +  NGC_IC_9(i__),R__HST(i__),e_R__HST(i__),R__GR(i__),
     +  e_R__GR(i__),v_R__(i__),e__R__(i__),Y__R__(i__),e_Y__R__(i__)
        if(ar__9(8:14) .EQ. '') R__HST(i__) = rNULL__
        if(ar__9(15:20) .EQ. '') e_R__HST(i__) = rNULL__
        if(ar__9(21:27) .EQ. '') R__GR(i__) = rNULL__
        if(ar__9(28:33) .EQ. '') e_R__GR(i__) = rNULL__
        if(ar__9(34:40) .EQ. '') v_R__(i__) = rNULL__
        if(ar__9(41:46) .EQ. '') e__R__(i__) = rNULL__
        if(ar__9(47:54) .EQ. '') Y__R__(i__) = rNULL__
        if(ar__9(55:61) .EQ. '') e_Y__R__(i__) = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  NGC_IC_9(i__),R__HST(i__),e_R__HST(i__),R__GR(i__),
     +  e_R__GR(i__),v_R__(i__),e__R__(i__),Y__R__(i__),e_Y__R__(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table12.dat'	! Spread in He and mass required explaining the HB

C  Format for file interpretation

   11 format(A7,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3)

C  Effective file loading

      open(unit=1,file='table12.dat', status='old')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,64
        read(1,'(A55)')ar__10
        read(ar__10,11)
     +  NGC_IC_10(i__),dM_Y(i__),Ymed(i__),Ymax(i__),dY(i__),
     +  d_B_V(i__),d_V_I(i__)
c    ..............Just test output...........
        write(6,11)
     +  NGC_IC_10(i__),dM_Y(i__),Ymed(i__),Ymax(i__),dY(i__),
     +  d_B_V(i__),d_V_I(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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