FORTRAN Generation
(/./ftp/cats/J/AJ/154/155)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/154/155 into FORTRAN code for loading all data files into arrays.

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-29
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/AJ/154/155 Abundance variations in the outer halo GC NGC 6229 (Johnson+, 2017)
*================================================================================
*Light and heavy element abundance variations in the outer halo globular cluster
*NGC 6229.
*    Johnson C.I., Caldwell N., Rich R.M., Walker M.G.
*   <Astron. J., 154, 155-155 (2017)>
*   =2017AJ....154..155J    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Star Identifiers, Coordinates, Photometry,
                                and Velocities

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

C  J2000 position composed of: RAdeg DEdeg
      character*1   Memb       (nr__) ! [MN] Cluster membership flag (M=Cluster member,
*                                 N=Non-member)
      character*16  v2MASS     (nr__) ! 2MASS star name (HHMMSSss+DDMMSSs)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000)
      real*4        Jmag       (nr__) ! (mag) [12.315/14.893] 2MASS J band magnitude
      real*4        Hmag       (nr__) ! (mag) [11.463/14.527] 2MASS H band magnitude
      real*4        Ksmag      (nr__) ! (mag) [11.283/14.158] 2MASS K_s_ band magnitude
      real*4        HRV        (nr__) ! (km/s) [-162.7/7.9] Heliocentric radial velocity
      real*4        e_HRV      (nr__) ! (km/s) [0.2/33.3] Uncertainty in HRV

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

C  Declarations for 'table2.dat'	! Model Atmosphere Parameters

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

      character*16  v2MASS_1   (nr__1) ! 2MASS star name (HHMMSSss+DDMMSSs)
      integer*4     Teff       (nr__1) ! (K) [3925/4450]? Effective temperature
      real*4        logg       (nr__1) ! ([cm/s2]) [0.7/1.6]? Surface gravity (in cgs units)
      real*4        v_Fe_H_    (nr__1) ! ([-]) [-1.23/-1.04]? Metallicity (in dex)
      real*4        e__Fe_H_   (nr__1) ! ([-]) [0.04/0.05]? Uncertainty in [Fe/H]
      real*4        ksi        (nr__1) ! (km/s) [1.4/2.2]? Microturbulence velocity {ksi}_mic_

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

C  Declarations for 'table3.dat'	! Line List and Adopted Reference Abundances

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

      real*8        lambda     (nr__2) ! (0.1nm) [5731.76/6862.48] Wavelength ({AA})
      character*5   Species    (nr__2) ! Species identifier
      real*4        EP         (nr__2) ! (eV) [0/5.86] Lower excitation potential
      real*4        log_gf     (nr__2) ! ([-]) [-9.75/-0.183]? Log of the oscillator strength
      character*3   n_log_gf   (nr__2) ! [hfs ] Note on log(gf) (1)
      real*4        logEsun    (nr__2) ! ([-]) [1.13/8.69] Adopted solar abundance
*                                  log{epsilon}(X)_{sun}_ (in dex)
      real*4        logeEarc   (nr__2) ! ([-]) [0.57/8.63] Adopted Arcturus abundance
*                                  log{epsilon}(X)_Arcturus_ (in dex)
*Note (1): Note as follows:
*  hfs = We adopted the hyperfine structure and isotope line lists of
*        Cunha et al. (2002AJ....124..379C) and Lawler et al.
*        (2001ApJ...556..452L) for Cu I and La II, respectively.

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

C  Declarations for 'table4.dat'	! Chemical Abundances and Uncertainties: Oxygen
                                to Calcium

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

      character*16  v2MASS_2   (nr__3) ! 2MASS star name (HHMMSSss+DDMMSSs)
      real*4        v_O_Fe_    (nr__3) ! ([-]) [-0.11/0.71]? [O/Fe] abundance (in dex)
      real*4        e__O_Fe_   (nr__3) ! ([-]) [0.07]? Uncertainty in [O/Fe] (in dex)
      real*4        v_Na_Fe_   (nr__3) ! ([-]) [-0.58/0.25]? [Na/Fe] abundance (in dex)
      real*4        e__Na_Fe_  (nr__3) ! ([-]) [0.11]? Uncertainty in [Na/Fe] (in dex)
      real*4        v_Mg_Fe_   (nr__3) ! ([-]) [0.05/0.34]? [Mg/Fe] abundance (in dex)
      real*4        e__Mg_Fe_  (nr__3) ! ([-]) [0.06]? Uncertainty in [Mg/Fe] (in dex)
      real*4        v_Al_Fe_   (nr__3) ! ([-]) [-0.45/0.39]? [Al/Fe] abundance (in dex)
      real*4        e__Al_Fe_  (nr__3) ! ([-]) [0.1]? Uncertainty in [Al/Fe] (in dex)
      real*4        v_Si_Fe_   (nr__3) ! ([-]) [0.18/0.38]? [Si/Fe] abundance (in dex)
      real*4        e__Si_Fe_  (nr__3) ! ([-]) [0.08/0.14]? Uncertainty in [Si/Fe] (in dex)
      real*4        v_Ca_Fe_   (nr__3) ! ([-]) [0.06/0.34]? [Ca/Fe] abundance (in dex)
      real*4        e__Ca_Fe_  (nr__3) ! ([-]) [0.12]? Uncertainty in [Ca/Fe] (in dex)

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

C  Declarations for 'table5.dat'	! Chemical Abundances and Uncertainties: Chromium
                                to Neodymium

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

      character*16  v2MASS_3   (nr__4) ! 2MASS star name (HHMMSSss+DDMMSSs)
      real*4        v_Cr_Fe_   (nr__4) ! ([-]) [-0.13/0.06]? [Cr/Fe] abundance (in dex)
      real*4        e__Cr_Fe_  (nr__4) ! ([-]) [0.12/0.14]? Uncertainty in [Cr/Fe] (in dex)
      real*4        v_Ni_Fe_   (nr__4) ! ([-]) [-0.19/-0.05]? [Ni/Fe] abundance (in dex)
      real*4        e__Ni_Fe_  (nr__4) ! ([-]) [0.03/0.06]? Uncertainty in [Ni/Fe] (in dex)
      real*4        v_Cu_Fe_   (nr__4) ! ([-]) [-0.53/-0.39]? [Cu/Fe] abundance (in dex)
      real*4        e__Cu_Fe_  (nr__4) ! ([-]) [0.09]? Uncertainty in [Cu/Fe] (in dex)
      real*4        v_Zr_Fe_   (nr__4) ! ([-]) [0.21/0.87]? [Zr/Fe] abundance (in dex)
      real*4        e__Zr_Fe_  (nr__4) ! ([-]) [0.12/0.17]? Uncertainty in [Zr/Fe] (in dex)
      real*4        v_La_Fe_   (nr__4) ! ([-]) [0.22/0.99]? [La/Fe] abundance (in dex)
      real*4        e__La_Fe_  (nr__4) ! ([-]) [0.05/0.08]? Uncertainty in [La/Fe] (in dex)
      real*4        v_Nd_Fe_   (nr__4) ! ([-]) [0.2/0.86]? [Nd/Fe] abundance (in dex)
      real*4        e__Nd_Fe_  (nr__4) ! ([-]) [0.05/0.14]? Uncertainty in [Nd/Fe] (in dex)

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

C  Loading file 'table1.dat'	! Star Identifiers, Coordinates, Photometry,
*                                and Velocities

C  Format for file interpretation

    1 format(
     +  A1,1X,A16,1X,F10.6,1X,F9.6,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.1,1X,
     +  F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,57
        read(1,'(A72)')ar__
        read(ar__,1)
     +  Memb(i__),v2MASS(i__),RAdeg(i__),DEdeg(i__),Jmag(i__),
     +  Hmag(i__),Ksmag(i__),HRV(i__),e_HRV(i__)
c    ..............Just test output...........
        write(6,1)
     +  Memb(i__),v2MASS(i__),RAdeg(i__),DEdeg(i__),Jmag(i__),
     +  Hmag(i__),Ksmag(i__),HRV(i__),e_HRV(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Model Atmosphere Parameters

C  Format for file interpretation

    2 format(A16,1X,I4,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,14
        read(1,'(A42)')ar__1
        read(ar__1,2)
     +  v2MASS_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),e__Fe_H_(i__),
     +  ksi(i__)
        if(ar__1(18:21) .EQ. '') Teff(i__) = iNULL__
        if(ar__1(23:26) .EQ. '') logg(i__) = rNULL__
        if(ar__1(28:32) .EQ. '') v_Fe_H_(i__) = rNULL__
        if(ar__1(34:37) .EQ. '') e__Fe_H_(i__) = rNULL__
        if(ar__1(39:42) .EQ. '') ksi(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v2MASS_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),e__Fe_H_(i__),
     +  ksi(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Line List and Adopted Reference Abundances

C  Format for file interpretation

    3 format(F7.2,1X,A5,1X,F4.2,1X,F6.3,1X,A3,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,88
        read(1,'(A39)')ar__2
        read(ar__2,3)
     +  lambda(i__),Species(i__),EP(i__),log_gf(i__),n_log_gf(i__),
     +  logEsun(i__),logeEarc(i__)
        if(ar__2(20:25) .EQ. '') log_gf(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  lambda(i__),Species(i__),EP(i__),log_gf(i__),n_log_gf(i__),
     +  logEsun(i__),logeEarc(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Chemical Abundances and Uncertainties: Oxygen
*                                to Calcium

C  Format for file interpretation

    4 format(
     +  A16,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,14
        read(1,'(A79)')ar__3
        read(ar__3,4)
     +  v2MASS_2(i__),v_O_Fe_(i__),e__O_Fe_(i__),v_Na_Fe_(i__),
     +  e__Na_Fe_(i__),v_Mg_Fe_(i__),e__Mg_Fe_(i__),v_Al_Fe_(i__),
     +  e__Al_Fe_(i__),v_Si_Fe_(i__),e__Si_Fe_(i__),v_Ca_Fe_(i__),
     +  e__Ca_Fe_(i__)
        if(ar__3(18:22) .EQ. '') v_O_Fe_(i__) = rNULL__
        if(ar__3(24:27) .EQ. '') e__O_Fe_(i__) = rNULL__
        if(ar__3(29:33) .EQ. '') v_Na_Fe_(i__) = rNULL__
        if(ar__3(35:38) .EQ. '') e__Na_Fe_(i__) = rNULL__
        if(ar__3(40:43) .EQ. '') v_Mg_Fe_(i__) = rNULL__
        if(ar__3(45:48) .EQ. '') e__Mg_Fe_(i__) = rNULL__
        if(ar__3(50:54) .EQ. '') v_Al_Fe_(i__) = rNULL__
        if(ar__3(56:59) .EQ. '') e__Al_Fe_(i__) = rNULL__
        if(ar__3(61:64) .EQ. '') v_Si_Fe_(i__) = rNULL__
        if(ar__3(66:69) .EQ. '') e__Si_Fe_(i__) = rNULL__
        if(ar__3(71:74) .EQ. '') v_Ca_Fe_(i__) = rNULL__
        if(ar__3(76:79) .EQ. '') e__Ca_Fe_(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v2MASS_2(i__),v_O_Fe_(i__),e__O_Fe_(i__),v_Na_Fe_(i__),
     +  e__Na_Fe_(i__),v_Mg_Fe_(i__),e__Mg_Fe_(i__),v_Al_Fe_(i__),
     +  e__Al_Fe_(i__),v_Si_Fe_(i__),e__Si_Fe_(i__),v_Ca_Fe_(i__),
     +  e__Ca_Fe_(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Chemical Abundances and Uncertainties: Chromium
*                                to Neodymium

C  Format for file interpretation

    5 format(
     +  A16,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,14
        read(1,'(A79)')ar__4
        read(ar__4,5)
     +  v2MASS_3(i__),v_Cr_Fe_(i__),e__Cr_Fe_(i__),v_Ni_Fe_(i__),
     +  e__Ni_Fe_(i__),v_Cu_Fe_(i__),e__Cu_Fe_(i__),v_Zr_Fe_(i__),
     +  e__Zr_Fe_(i__),v_La_Fe_(i__),e__La_Fe_(i__),v_Nd_Fe_(i__),
     +  e__Nd_Fe_(i__)
        if(ar__4(18:22) .EQ. '') v_Cr_Fe_(i__) = rNULL__
        if(ar__4(24:27) .EQ. '') e__Cr_Fe_(i__) = rNULL__
        if(ar__4(29:33) .EQ. '') v_Ni_Fe_(i__) = rNULL__
        if(ar__4(35:38) .EQ. '') e__Ni_Fe_(i__) = rNULL__
        if(ar__4(40:44) .EQ. '') v_Cu_Fe_(i__) = rNULL__
        if(ar__4(46:49) .EQ. '') e__Cu_Fe_(i__) = rNULL__
        if(ar__4(51:54) .EQ. '') v_Zr_Fe_(i__) = rNULL__
        if(ar__4(56:59) .EQ. '') e__Zr_Fe_(i__) = rNULL__
        if(ar__4(61:64) .EQ. '') v_La_Fe_(i__) = rNULL__
        if(ar__4(66:69) .EQ. '') e__La_Fe_(i__) = rNULL__
        if(ar__4(71:74) .EQ. '') v_Nd_Fe_(i__) = rNULL__
        if(ar__4(76:79) .EQ. '') e__Nd_Fe_(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v2MASS_3(i__),v_Cr_Fe_(i__),e__Cr_Fe_(i__),v_Ni_Fe_(i__),
     +  e__Ni_Fe_(i__),v_Cu_Fe_(i__),e__Cu_Fe_(i__),v_Zr_Fe_(i__),
     +  e__Zr_Fe_(i__),v_La_Fe_(i__),e__La_Fe_(i__),v_Nd_Fe_(i__),
     +  e__Nd_Fe_(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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