FORTRAN Generation
(/./ftp/cats/J/ApJ/708/560)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/708/560 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/ApJ/708/560           Spectroscopy of UMa II and Coma Ber     (Frebel+, 2010)
*================================================================================
*High-resolution spectroscopy of extremely metal-poor stars in the least evolved
*galaxies: Ursa Major II and Coma Berenices.
*    Frebel A., Simon J.D., Geha M., Willman B.
*   <Astrophys. J., 708, 560-583 (2010)>
*   =2010ApJ...708..560F
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'stars.dat'	! Photometry and other stellar parameters (table 1
                             and 5 of the paper)

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

      character*14  Gal        (nr__) ! Dwarf galaxy name (Coma Berenices or
*                                   Ursa Major II)
      character*14  SDSS       (nr__) ! Star SDSS designation (JHHMMSS+DDMMSS)
      character*9   Name       (nr__) ! Star name
      character*1   f_Name     (nr__) ! [f] Non-member star (1)
      real*4        Vmag       (nr__) ! (mag) ? V-band magnitude (2)
      real*4        rmag       (nr__) ! (mag) ? SDSS-DR7 r-band magnitude (2)
      real*4        g_r        (nr__) ! (mag) ? SDSS-DR7 g-r color index (2)
      real*4        rMag_1     (nr__) ! (mag) ? Absolute r-band magnitude (2)
      integer*4     Teff       (nr__) ! (K) Effective temperature (3)
      real*4        logg       (nr__) ! ([cm/s2]) Surface gravity
      real*4        v_Fe_H_    (nr__) ! ([-]) Fe/H abundance
      real*4        Vt         (nr__) ! (km/s) Microturbulent velocity {xi}
*Note (1):
*  f = UMa II-NM was observed but was later determined not to be a member of
*      UMa II.
*Note (2): The r magnitudes and g-r colors were obtained from the updated SDSS
*          DR7 (Abazajian et al. 2009, Cat. II/294), while the V magnitudes were
*          determined from g and r using the conversions given by Smith et al.
*          (2002AJ....123.2121S).
*Note (3): Temperature are rounded to the nearest 10K.

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

C  Declarations for 'table2.dat'	! Observing details

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__1) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__1) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*9   Name_1     (nr__1) ! Star name
      integer*4     RAh        (nr__1) ! (h) Hour of right ascension (J2000)
      integer*4     RAm        (nr__1) ! (min) Minute of right ascension (J2000)
      real*4        RAs        (nr__1) ! (s) Second of right ascension (J2000)
      character*1   DE_        (nr__1) ! Sign of declination (J2000)
      integer*4     DEd        (nr__1) ! (deg) Degree of declination (J2000)
      integer*4     DEm        (nr__1) ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs        (nr__1) ! (arcsec) Arcsecond of declination (J2000)
      real*8        JD         (nr__1) ! (d) Julian Date of observation
      real*4        Exp        (nr__1) ! (h) Exposure time
      integer*4     S_N5       (nr__1) ! 5000{AA} signal to noise ratio (4)
      integer*4     S_N6       (nr__1) ! 6000{AA} signal to noise ratio (4)
      integer*4     S_N65      (nr__1) ! 6500{AA} signal to noise ratio (4)
      real*4        RVH        (nr__1) ! (km/s) HIRES radial velocity
      real*4        e_RVH      (nr__1) ! (km/s) RVH uncertainty
      real*4        RVD        (nr__1) ! (km/s) ? DEIMOS radial velocity
      real*4        e_RVD      (nr__1) ! (km/s) ? RVD uncertainty
      character*8   Comm       (nr__1) ! Comment
*Note (4): The S/N measurements are for ~22m{AA} ("blue" CCD), ~26m{AA}
*     ("green" CCD), and ~28m{AA} ("green" CCD) pixel sizes, respectively.

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

C  Declarations for 'ew.dat'	! Equivalent width measurements (table 4)

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

      character*9   Name_2     (nr__2) ! Star name
      character*4   Ion        (nr__2) ! Species identification
      real*8        lambda     (nr__2) ! (0.1nm) Wavelength in Angstroms
      real*4        ExPot      (nr__2) ! (eV) ? Excitation potential
      real*8        log_gf     (nr__2) ! ([-]) ? Log of the oscillator strength
      character*1   l_EW       (nr__2) ! Limit flag on EW
      real*4        EW         (nr__2) ! (0.1pm) ? Equivalent width in milli-Angstroms
      character*1   f_EW       (nr__2) ! [s] s: indicates a spectrum synthesis measurement
      character*1   l_log_e    (nr__2) ! Limit flag on log(e)
      real*4        log_e      (nr__2) ! ([-]) ? log{epsilon} abundance (H=12)

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

C  Declarations for 'table6.dat'	! HIRES abundances of the UMa II stars

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

      character*9   Name_3     (nr__3) ! Star name
      character*5   Ion_1      (nr__3) ! Species identification
      character*1   l_log_e_1  (nr__3) ! Limit flag on loge
      real*4        log_e_1    (nr__3) ! ([-]) log{epsilon} abundance (X)
      character*1   f_log_e    (nr__3) ! [a] distorted lines (1)
      character*1   l__X_H_    (nr__3) ! Limit flag on [X/H]
      real*4        v_X_H_     (nr__3) ! ([Sun]) X/H abundance
      character*1   l__X_Fe_   (nr__3) ! Limit flag on [X/Fe]
      real*4        v_X_Fe_    (nr__3) ! ([Sun]) ? X/Fe abundance (G1)
      integer*4     N          (nr__3) ! Number of spectra used
      real*4        Sig        (nr__3) ! ([-]) ? Abundance uncertainty
*Note (1):
*   a = The Zn lines are somewhat distorted; this may lead to an overestimated
*       Zn abundance.

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

C  Declarations for 'table7.dat'	! HIRES abundances of the Coma Ber stars

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

      character*9   Name_4     (nr__4) ! Star name
      character*5   Ion_2      (nr__4) ! Species identification
      character*1   l_log_e_1  (nr__4) ! Limit flag on loge
      real*4        log_e_1    (nr__4) ! ([-]) log{epsilon} abundance (X)
      character*1   f_log_e_1  (nr__4) ! [a] distorted lines (1)
      character*1   l__X_H__1  (nr__4) ! Limit flag on [X/H]
      real*4        v_X_H__1   (nr__4) ! ([Sun]) X/H abundance
      character*1   l__X_Fe__1 (nr__4) ! Limit flag on [X/Fe]
      real*4        v_X_Fe__1  (nr__4) ! ([Sun]) ? X/Fe abundance (G1)
      integer*4     N_1        (nr__4) ! Number of spectra used
      real*4        Sig_1      (nr__4) ! ([-]) ? Abundance uncertainty
*Note (1):
*   a = The Zn lines are somewhat distorted; this may lead to an overestimated
*       Zn abundance.

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

C  Declarations for 'table9.dat'	! HIRES abundances of standard star HD 122563

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

      character*9   Name_5     (nr__5) ! Star name
      character*5   Ion_3      (nr__5) ! Species identification
      real*4        log_e_X0   (nr__5) ! ([-]) ? Log {epsilon} abundance (X)_0_
      character*1   l_log_e_1  (nr__5) ! Limit flag on loge
      real*4        log_e_1    (nr__5) ! ([-]) ? Log {epsilon} abundance (X)
      character*1   f_log_e_1  (nr__5) ! [b] neutron capture element (1)
      character*1   l__X_H__2  (nr__5) ! Limit flag on [X/H]
      real*4        v_X_H__2   (nr__5) ! ([Sun]) ? X/H abundance
      character*1   l__X_Fe__2 (nr__5) ! Limit flag on [X/Fe]
      real*4        v_X_Fe__2  (nr__5) ! ([Sun]) ? X/Fe abundance (G1)
      integer*4     N_2        (nr__5) ! ? Number of spectra used
      real*4        Sig_2      (nr__5) ! ([-]) ? Abundance uncertainty
      character*1   ref        (nr__5) ! [cd] source of the data: c=this paper,
*                                  d=Aoki et al. (2007ApJ...660..747A)
*Note (1): b = Neutron-capture abundances are taken from Honda et al.
*       (2006ApJ...643.1180H).

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

C  Loading file 'stars.dat'	! Photometry and other stellar parameters (table 1
*                             and 5 of the paper)

C  Format for file interpretation

    1 format(
     +  A14,6X,A14,1X,A9,1X,A1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,I4,
     +  1X,F3.1,1X,F5.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'stars.dat')
      write(6,*) '....Loading file: stars.dat'
      do i__=1,8
        read(1,'(A88)')ar__
        read(ar__,1)
     +  Gal(i__),SDSS(i__),Name(i__),f_Name(i__),Vmag(i__),rmag(i__),
     +  g_r(i__),rMag_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),Vt(i__)
        if(ar__(48:52) .EQ. '') Vmag(i__) = rNULL__
        if(ar__(54:58) .EQ. '') rmag(i__) = rNULL__
        if(ar__(60:63) .EQ. '') g_r(i__) = rNULL__
        if(ar__(65:69) .EQ. '') rMag_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Gal(i__),SDSS(i__),Name(i__),f_Name(i__),Vmag(i__),rmag(i__),
     +  g_r(i__),rMag_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),Vt(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Observing details

C  Format for file interpretation

    2 format(
     +  A9,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F9.1,1X,F4.2,
     +  1X,I3,1X,I3,1X,I3,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,A8)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,8
        read(1,'(A92)')ar__1
        read(ar__1,2)
     +  Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),JD(i__),Exp(i__),S_N5(i__),S_N6(i__),
     +  S_N65(i__),RVH(i__),e_RVH(i__),RVD(i__),e_RVD(i__),Comm(i__)
        if(ar__1(74:79) .EQ. '') RVD(i__) = rNULL__
        if(ar__1(81:83) .EQ. '') e_RVD(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
     +  DEm(i__),DEs(i__),JD(i__),Exp(i__),S_N5(i__),S_N6(i__),
     +  S_N65(i__),RVH(i__),e_RVH(i__),RVD(i__),e_RVD(i__),Comm(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ew.dat'	! Equivalent width measurements (table 4)

C  Format for file interpretation

    3 format(
     +  A9,1X,A4,1X,F7.2,1X,F4.2,1X,F7.3,1X,A1,F5.1,1X,A1,1X,A1,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ew.dat')
      write(6,*) '....Loading file: ew.dat'
      do i__=1,2191
        read(1,'(A52)')ar__2
        read(ar__2,3)
     +  Name_2(i__),Ion(i__),lambda(i__),ExPot(i__),log_gf(i__),
     +  l_EW(i__),EW(i__),f_EW(i__),l_log_e(i__),log_e(i__)
        if(ar__2(24:27) .EQ. '') ExPot(i__) = rNULL__
        if(ar__2(29:35) .EQ. '') log_gf(i__) = rNULL__
        if(ar__2(38:42) .EQ. '') EW(i__) = rNULL__
        if(ar__2(47:52) .EQ. '') log_e(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Name_2(i__),Ion(i__),lambda(i__),ExPot(i__),log_gf(i__),
     +  l_EW(i__),EW(i__),f_EW(i__),l_log_e(i__),log_e(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! HIRES abundances of the UMa II stars

C  Format for file interpretation

    4 format(
     +  A9,1X,A5,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,84
        read(1,'(A47)')ar__3
        read(ar__3,4)
     +  Name_3(i__),Ion_1(i__),l_log_e_1(i__),log_e_1(i__),
     +  f_log_e(i__),l__X_H_(i__),v_X_H_(i__),l__X_Fe_(i__),
     +  v_X_Fe_(i__),N(i__),Sig(i__)
        if(ar__3(34:38) .EQ. '') v_X_Fe_(i__) = rNULL__
        if(ar__3(44:47) .EQ. '') Sig(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Name_3(i__),Ion_1(i__),l_log_e_1(i__),log_e_1(i__),
     +  f_log_e(i__),l__X_H_(i__),v_X_H_(i__),l__X_Fe_(i__),
     +  v_X_Fe_(i__),N(i__),Sig(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! HIRES abundances of the Coma Ber stars

C  Format for file interpretation

    5 format(
     +  A9,1X,A5,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,84
        read(1,'(A47)')ar__4
        read(ar__4,5)
     +  Name_4(i__),Ion_2(i__),l_log_e_1(i__),log_e_1(i__),
     +  f_log_e_1(i__),l__X_H__1(i__),v_X_H__1(i__),l__X_Fe__1(i__),
     +  v_X_Fe__1(i__),N_1(i__),Sig_1(i__)
        if(ar__4(34:38) .EQ. '') v_X_Fe__1(i__) = rNULL__
        if(ar__4(44:47) .EQ. '') Sig_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Name_4(i__),Ion_2(i__),l_log_e_1(i__),log_e_1(i__),
     +  f_log_e_1(i__),l__X_H__1(i__),v_X_H__1(i__),l__X_Fe__1(i__),
     +  v_X_Fe__1(i__),N_1(i__),Sig_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! HIRES abundances of standard star HD 122563

C  Format for file interpretation

    6 format(
     +  A9,1X,A5,1X,F4.2,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3,
     +  1X,F4.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,54
        read(1,'(A54)')ar__5
        read(ar__5,6)
     +  Name_5(i__),Ion_3(i__),log_e_X0(i__),l_log_e_1(i__),
     +  log_e_1(i__),f_log_e_1(i__),l__X_H__2(i__),v_X_H__2(i__),
     +  l__X_Fe__2(i__),v_X_Fe__2(i__),N_2(i__),Sig_2(i__),ref(i__)
        if(ar__5(17:20) .EQ. '') log_e_X0(i__) = rNULL__
        if(ar__5(23:27) .EQ. '') log_e_1(i__) = rNULL__
        if(ar__5(32:36) .EQ. '') v_X_H__2(i__) = rNULL__
        if(ar__5(39:43) .EQ. '') v_X_Fe__2(i__) = rNULL__
        if(ar__5(45:47) .EQ. '') N_2(i__) = iNULL__
        if(ar__5(49:52) .EQ. '') Sig_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Name_5(i__),Ion_3(i__),log_e_X0(i__),l_log_e_1(i__),
     +  log_e_1(i__),f_log_e_1(i__),l__X_H__2(i__),v_X_H__2(i__),
     +  l__X_Fe__2(i__),v_X_Fe__2(i__),N_2(i__),Sig_2(i__),ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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