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 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 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/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      ! Cluster name
      character*7   Name        ! Other name
      real*4        v_Fe_H_     ! ([Sun]) Metallicity
      real*8        v_m_M_V     ! (mag) Distance modulus calculated from V magnitude
      real*8        E_B_V       ! (mag) Redenning
      real*8        VMAG        ! (mag) Absolute V magnitude
      real*8        HBR         ! ? Horizontal-branch ratio = (B-R)/(B+V+R)
      real*8        Age         ! (Gyr) ? Estimated cluster age
      real*8        M_TO        ! (solMass) ? Mass at turn-off
      real*8        M_RGB       ! (solMass) ? Mass at the tip of the red-giant branch
      real*8        IQR1        ! ? Interquartile of the [Na/O]
*                                     (anti-)correlations
      real*8        IQR2        ! ? Interquartile of the [Mg/Al]
*                                     (anti-)correlations
      real*8        T_HB        ! ([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    ! Cluster name
      integer*4     n_HB        ! Number of Horizon Branch stars
      integer*4     n_fld       ! [0/125] Number of field stars (G1)
      real*4        F555        ! (mag) Median F555W magnitude
      real*4        e_F555      ! (mag) rms uncertainty on F555
      real*4        v439_555    ! (mag) Median F439W-F555W colour index
      real*4        e_439_555   ! (mag) rms uncertainty on F439-F555
      real*4        F555m       ! (mag) Minimum F555W magnitude
      real*4        e_F555m     ! (mag) rms uncertainty on F555m
      real*4        v439_555m   ! (mag) Minimum F439W-F555W colour index
      real*4        e_439_555m  ! (mag) rms uncertainty on (F439-F555)m
      real*4        v439_555M_1 ! (mag) Maximum F439W-F555W colour index
      real*4        e_439_555M_1 ! (mag) rms uncertainty on (F439-F555)M
      integer*4     n_RGB       ! Number of RGB stars
      real*4        v439_555G   ! (mag) F439W-F555W colour index for RGB
      real*4        e_439_555G  ! (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    ! Cluster name
      integer*4     n_HB_1      ! Number of HB stars in the cluster
      integer*4     n_fld_1     ! [0/47] Number of field stars (G1)
      real*4        Vmag_1      ! (mag) Median V magnitude
      real*4        e_Vmag      ! (mag) rms uncertainty on
      real*4        V_I         ! (mag) Median V-I colour index
      real*4        e_V_I       ! (mag) rms uncertainty on V-I
      real*4        Vmagm       ! (mag) Minimum V magnitude
      real*4        e_Vmagm     ! (mag) rms uncertainty on
      real*4        v_V_I_m     ! (mag) Minimum V magnitude
      real*4        e__V_I_m    ! (mag) rms uncertainty on (V-I)m
      real*4        v_V_I_M_1   ! (mag) Maximum V-I colour index
      real*4        e__V_I_M_1  ! (mag) rms uncertainty on (V-I)M
      integer*4     n_RGB_1     ! Number of RGB stars
      real*4        v_V_I_G     ! (mag) V-I colour index for RGB
      real*4        e__V_I_G    ! (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    ! Cluster name
      real*4        v_B_V_0G    ! (mag) Intrinsic B-V colour index value, for RGB
      real*4        v_B_V_0m    ! (mag) Minimim B-V colour index value
      real*4        v_B_V_0     ! (mag) Median B-V colour index value
      real*4        v_B_V_0M_1  ! (mag) Maximum B-V colour index value
      real*4        VMAGM_1     ! (mag) Maximum absolute V magnitude
      real*4        VMAG_2      ! (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    ! Cluster name
      real*4        v_V_I_0G    ! (mag) Intrinsic V-I colour index value, for RGB
      real*4        v_V_I_0m    ! (mag) Mimimum V-I colour index value
      real*4        v_V_I_0     ! (mag) Median V-I colour index value
      real*4        v_V_I_0M_1  ! (mag) Maximum V-I colour index value
      real*4        VMAGM_2     ! (mag) Maximum absolute V magnitude
      real*4        VMAG_3      ! (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    ! Cluster name
      real*4        Mmin        ! (solMass) Mimimum mass
      real*4        e_Mmin      ! (solMass) ? rms uncertainty on Mmin
      real*4        Mmed        ! (solMass) Median mass
      real*4        e_Mmed      ! (solMass) ? rms uncertainty on Mmed
      real*4        Mmax        ! (solMass) ? Maximum mass
      real*4        e_Mmax      ! (solMass) ? rms uncertainty on Mmax
      real*4        dMmed       ! (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    ! Cluster name
      real*4        Mmin_1      ! (solMass) Mimimum mass
      real*4        e_Mmin_1    ! (solMass) ? rms uncertainty on Mmin
      real*4        Mmed_1      ! (solMass) Median mass
      real*4        e_Mmed_1    ! (solMass) ? rms uncertainty on Mmed
      real*4        Mmax_1      ! (solMass) ? Maximum mass
      real*4        e_Mmax_1    ! (solMass) ? rms uncertainty on Mmax
      real*4        dMmed_1     ! (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    ! Cluster name
      real*4        S_RR        ! ? RR Lyrae richness parameter (Harris, 1996,
*                                    Cat. VII/195) (1)
      integer*4     n_RR_H      ! ? Number of RR Lyrae stars for HST sample
      real*4        f_RR_H      ! ? Fraction of RR Lyrae stars for HST sample
      integer*4     n_RR_g      ! ? Number of RR Lyrae for ground-based sample
      real*4        f_RR_g      ! ? Fraction of RR Lyrae for ground-based sample
      real*4        f_RR        ! ? Fraction of RR Lyrae stars
      real*4        v_Pab_      ! (d) ? Mean period
      integer*4     n_Pab       ! ? 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    ! Cluster name
      real*8        Age_1       ! (Gyr) ? Age from main-sequence fitting
      real*8        Age_HB_H    ! (Gyr) ? Age from horizontal branch fitting,
*                                      for HST data-set
      real*8        Age_HB_g    ! (Gyr) ? Age from horizontal branch fitting,
*                                      for ground-based data-set
      real*8        Age_Y       ! (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    ! Cluster name
      real*8        R__HST      ! ? R' parameter as derived from HST data (1)
      real*4        e_R__HST    ! ? rms uncertainty on R'(HST)
      real*8        R__GR       ! ? R' parameter as derived from ground-based
*                                    data (1)
      real*4        e_R__GR     ! ? rms uncertainty on R'(GR)
      real*8        v_R__       ! ? Mean R' parameter (1) (2)
      real*4        e__R__      ! ? rms uncertainty on <R'>
      real*8        Y__R__      ! ? He abundance for R'=<R'>
      real*8        e_Y__R__    ! ? 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   ! Cluster name
      real*8        dM_Y        ! (solMass) Spread in mass {delta}M(Y)
      real*8        Ymed        ! Median He abundance value
      real*8        Ymax        ! Maximum He abundance value
      real*8        dY          ! Difference between maximum and median Y values
      real*8        d_B_V       ! (mag) Colour spread in {delta}(B-V)
      real*8        d_V_I       ! (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,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,94
        read(1,'(A98)')ar__
        read(ar__,1)
     +  NGC_IC,Name,v_Fe_H_,v_m_M_V,E_B_V,VMAG,HBR,Age,M_TO,M_RGB,
     +  IQR1,IQR2,T_HB
        if(ar__(44:51) .EQ. '') HBR = rNULL__
        if(ar__(52:58) .EQ. '') Age = rNULL__
        if(ar__(59:66) .EQ. '') M_TO = rNULL__
        if(ar__(67:74) .EQ. '') M_RGB = rNULL__
        if(ar__(75:81) .EQ. '') IQR1 = rNULL__
        if(ar__(82:89) .EQ. '') IQR2 = rNULL__
        if(ar__(90:98) .EQ. '') T_HB = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  NGC_IC,Name,v_Fe_H_,v_m_M_V,E_B_V,VMAG,HBR,Age,M_TO,M_RGB,
     +  IQR1,IQR2,T_HB
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,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,74
        read(1,'(A108)')ar__1
        read(ar__1,2)
     +  NGC_IC_1,n_HB,n_fld,F555,e_F555,v439_555,e_439_555,F555m,
     +  e_F555m,v439_555m,e_439_555m,v439_555M_1,e_439_555M_1,n_RGB,
     +  v439_555G,e_439_555G
c    ..............Just test output...........
        write(6,2)
     +  NGC_IC_1,n_HB,n_fld,F555,e_F555,v439_555,e_439_555,F555m,
     +  e_F555m,v439_555m,e_439_555m,v439_555M_1,e_439_555M_1,n_RGB,
     +  v439_555G,e_439_555G
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,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,49
        read(1,'(A108)')ar__2
        read(ar__2,3)
     +  NGC_IC_2,n_HB_1,n_fld_1,Vmag_1,e_Vmag,V_I,e_V_I,Vmagm,e_Vmagm,
     +  v_V_I_m,e__V_I_m,v_V_I_M_1,e__V_I_M_1,n_RGB_1,v_V_I_G,e__V_I_G
c    ..............Just test output...........
        write(6,3)
     +  NGC_IC_2,n_HB_1,n_fld_1,Vmag_1,e_Vmag,V_I,e_V_I,Vmagm,e_Vmagm,
     +  v_V_I_m,e__V_I_m,v_V_I_M_1,e__V_I_M_1,n_RGB_1,v_V_I_G,e__V_I_G
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,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,74
        read(1,'(A57)')ar__3
        read(ar__3,4)
     +  NGC_IC_3,v_B_V_0G,v_B_V_0m,v_B_V_0,v_B_V_0M_1,VMAGM_1,VMAG_2
c    ..............Just test output...........
        write(6,4)
     +  NGC_IC_3,v_B_V_0G,v_B_V_0m,v_B_V_0,v_B_V_0M_1,VMAGM_1,VMAG_2
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,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,49
        read(1,'(A57)')ar__4
        read(ar__4,5)
     +  NGC_IC_4,v_V_I_0G,v_V_I_0m,v_V_I_0,v_V_I_0M_1,VMAGM_2,VMAG_3
c    ..............Just test output...........
        write(6,5)
     +  NGC_IC_4,v_V_I_0G,v_V_I_0m,v_V_I_0,v_V_I_0M_1,VMAGM_2,VMAG_3
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,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,74
        read(1,'(A61)')ar__5
        read(ar__5,6)
     +  NGC_IC_5,Mmin,e_Mmin,Mmed,e_Mmed,Mmax,e_Mmax,dMmed
        if(ar__5(18:22) .EQ. '') e_Mmin = rNULL__
        if(ar__5(33:37) .EQ. '') e_Mmed = rNULL__
        if(ar__5(41:45) .EQ. '') Mmax = rNULL__
        if(ar__5(49:53) .EQ. '') e_Mmax = rNULL__
        if(ar__5(57:61) .EQ. '') dMmed = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  NGC_IC_5,Mmin,e_Mmin,Mmed,e_Mmed,Mmax,e_Mmax,dMmed
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,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,49
        read(1,'(A61)')ar__6
        read(ar__6,7)
     +  NGC_IC_6,Mmin_1,e_Mmin_1,Mmed_1,e_Mmed_1,Mmax_1,e_Mmax_1,
     +  dMmed_1
        if(ar__6(18:22) .EQ. '') e_Mmin_1 = rNULL__
        if(ar__6(33:37) .EQ. '') e_Mmed_1 = rNULL__
        if(ar__6(41:45) .EQ. '') Mmax_1 = rNULL__
        if(ar__6(49:53) .EQ. '') e_Mmax_1 = rNULL__
        if(ar__6(57:61) .EQ. '') dMmed_1 = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  NGC_IC_6,Mmin_1,e_Mmin_1,Mmed_1,e_Mmed_1,Mmax_1,e_Mmax_1,
     +  dMmed_1
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,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,66
        read(1,'(A73)')ar__7
        read(ar__7,8)
     +  NGC_IC_7,S_RR,n_RR_H,f_RR_H,n_RR_g,f_RR_g,f_RR,v_Pab_,n_Pab
        if(ar__7(10:14) .EQ. '') S_RR = rNULL__
        if(ar__7(18:19) .EQ. '') n_RR_H = iNULL__
        if(ar__7(27:30) .EQ. '') f_RR_H = rNULL__
        if(ar__7(37:39) .EQ. '') n_RR_g = iNULL__
        if(ar__7(48:51) .EQ. '') f_RR_g = rNULL__
        if(ar__7(56:59) .EQ. '') f_RR = rNULL__
        if(ar__7(62:66) .EQ. '') v_Pab_ = rNULL__
        if(ar__7(70:73) .EQ. '') n_Pab = iNULL__
c    ..............Just test output...........
        write(6,8)
     +  NGC_IC_7,S_RR,n_RR_H,f_RR_H,n_RR_g,f_RR_g,f_RR,v_Pab_,n_Pab
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,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,94
        read(1,'(A40)')ar__8
        read(ar__8,9)NGC_IC_8,Age_1,Age_HB_H,Age_HB_g,Age_Y
        if(ar__8(10:16) .EQ. '') Age_1 = rNULL__
        if(ar__8(17:24) .EQ. '') Age_HB_H = rNULL__
        if(ar__8(25:32) .EQ. '') Age_HB_g = rNULL__
        if(ar__8(33:40) .EQ. '') Age_Y = rNULL__
c    ..............Just test output...........
        write(6,9)NGC_IC_8,Age_1,Age_HB_H,Age_HB_g,Age_Y
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,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,94
        read(1,'(A61)')ar__9
        read(ar__9,10)
     +  NGC_IC_9,R__HST,e_R__HST,R__GR,e_R__GR,v_R__,e__R__,Y__R__,
     +  e_Y__R__
        if(ar__9(8:14) .EQ. '') R__HST = rNULL__
        if(ar__9(15:20) .EQ. '') e_R__HST = rNULL__
        if(ar__9(21:27) .EQ. '') R__GR = rNULL__
        if(ar__9(28:33) .EQ. '') e_R__GR = rNULL__
        if(ar__9(34:40) .EQ. '') v_R__ = rNULL__
        if(ar__9(41:46) .EQ. '') e__R__ = rNULL__
        if(ar__9(47:54) .EQ. '') Y__R__ = rNULL__
        if(ar__9(55:61) .EQ. '') e_Y__R__ = rNULL__
c    ..............Just test output...........
        write(6,10)
     +  NGC_IC_9,R__HST,e_R__HST,R__GR,e_R__GR,v_R__,e__R__,Y__R__,
     +  e_Y__R__
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,status='old',file=
     +'table12.dat')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,64
        read(1,'(A55)')ar__10
        read(ar__10,11)NGC_IC_10,dM_Y,Ymed,Ymax,dY,d_B_V,d_V_I
c    ..............Just test output...........
        write(6,11)NGC_IC_10,dM_Y,Ymed,Ymax,dY,d_B_V,d_V_I
c    .......End.of.Just test output...........
      end do
      close(1)

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