FORTRAN Generation
(/./ftp/cats/J/ApJ/801/L10)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/801/L10 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 2026-Oct-15
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/801/L10          XO-2N and XO-2S equivalent widths          (Teske+, 2015)
*================================================================================
*Abundance differences between exoplanet binary host stars XO-2N and XO-2S -
*Dependence on stellar parameters.
*    Teske J.K., Ghezzi L., Cunha K., Smith V.V., Schuler S.C., Bergemann M.
*   <Astrophys. J., 801, L10 (2015)>
*   =2015ApJ...801L..10T
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Derived stellar parameters and {Delta}(N-S)
                              abundances

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

      character*4   Param       ! "Orig" or "Alt" set
      integer*4     Set         ! ? Set number for "Alt" params (1)
      character*5   Name        ! Star name (XO-2N or XO-2S)
      integer*4     Teff        ! (K) [5343/5547] Effective temperature
      integer*4     e_Teff      ! (K) [32/78] Teff uncertainty
      real*4        logg        ! ([cm/s2]) [4.2/4.5] Log of surface gravity
      real*4        e_logg      ! ([cm/s2]) [0.08/0.3] logg uncertainty
      real*4        xi          ! (km/s) [1/1.4] Microturbulent velocity {xi}
      real*4        e_xi        ! (km/s) [0.05/0.2] xi uncertainty
      real*4        v_Fe_H_     ! ([-]) [0.2/0.5] Metallicity
      real*4        e__Fe_H_    ! ([-]) [0.04/0.2] [Fe/H] uncertainty
      integer*4     N1          ! [23/72] Number of FeI lines measured
      integer*4     N2          ! [8/10] Number of FeII lines measured
      character*3   Sol_        ! [yes no] Relative to solar?
      character*4   IRAF_       ! [IRAF ARES] IRAF/ARES reduction
*Note (1): Alternative Parameter as follows:
*  1 = We repeated our analysis relative to solar, using by-hand IRAF
*      measurements of the lines in Ghezzi et al. (2010, J/ApJ/720/1290). 
*  2 = We adopted the full analysis approach of Ghezzi et al. (2010,
*      J/ApJ/720/1290) and derived parameters based on absolute abundances, not
*      relative to solar. The Ghezzi et al. (2010, J/ApJ/720/1290) line list,
*      with lab-measured gf-values, was compiled specifically for an absolute
*      analysis. 
*  3 = We used ODFNEW models with a different and longer line list of 72/67 FeI
*      and 9/8 FeII lines for N/S, compiled from Sousa et al. (2008,
*      J/A+A/487/373). The equivalent widths were measured with ARES
*      (Sousa et al. 2007, J/A+A/469/783); lines >0.30dex from the average in
*      the Sun were removed and an interative 2{sigma} clipping was applied to
*      the stars during the parameter determination. 
*      See section 2.2 for full details.

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

C  Declarations for 'table2.dat'	! Lines and equivalent widths measured

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

      character*5   Ion         ! Ion designation
      real*8        lambda      ! (0.1nm) [4465.8/7114.6] Wavelength, in Angstroms
      real*4        ExPot       ! (eV) [0/8] Excitation potential
      real*4        log_gf      ! ([-]) [-9.8/0.6] Log of oscillator strength
      real*4        EW0         ! (10-13m) [4/123]? Solar equivalent width (m{AA})
      real*4        EW2S        ! (10-13m) [11/170]? XO-2S equivalent width (m{AA})
      character*1   f_EW2S      ! [a-d] Flag on XO-2S (1)
*                                             (for TiI only)
      real*4        dA2S        ! ([-]) [0.03/0.1]? NLTE abundance correction for XO-2S
      real*4        EW2N        ! (10-13m) [8.2/176.3]? XO-2N equivalent width (m{AA})
      character*1   f_EW2N      ! [b-d] Flag on XO2N (1)
*                                             (for TiI only)
      real*4        dA2N        ! ([-]) [0.02/0.1]? NLTE abundance correction for XO-2N
*Note (1): Flag as follows:
*  a = abundance derived through synthesis;
*  b = outlier, removed from final abundance;
*  c = measured for alt param sets 1, 2, 3 (see section 2.2 and table1);
*  d = measured for alt param set 4

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

C  Loading file 'table1.dat'	! Derived stellar parameters and {Delta}(N-S)
*                              abundances

C  Format for file interpretation

    1 format(
     +  A4,8X,I1,1X,A5,1X,I4,1X,I2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  F4.2,1X,F4.2,1X,I2,6X,I2,7X,A3,1X,A4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,8
        read(1,'(A83)')ar__
        read(ar__,1)
     +  Param,Set,Name,Teff,e_Teff,logg,e_logg,xi,e_xi,v_Fe_H_,
     +  e__Fe_H_,N1,N2,Sol_,IRAF_
        if(ar__(13:13) .EQ. '') Set = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Param,Set,Name,Teff,e_Teff,logg,e_logg,xi,e_xi,v_Fe_H_,
     +  e__Fe_H_,N1,N2,Sol_,IRAF_
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Lines and equivalent widths measured

C  Format for file interpretation

    2 format(
     +  A5,1X,F7.2,1X,F4.2,1X,F5.2,1X,F5.1,1X,F5.1,1X,A1,1X,F4.2,1X,
     +  F5.1,1X,A1,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,342
        read(1,'(A56)')ar__1
        read(ar__1,2)
     +  Ion,lambda,ExPot,log_gf,EW0,EW2S,f_EW2S,dA2S,EW2N,f_EW2N,dA2N
        if(ar__1(26:30) .EQ. '') EW0 = rNULL__
        if(ar__1(32:36) .EQ. '') EW2S = rNULL__
        if(ar__1(40:43) .EQ. '') dA2S = rNULL__
        if(ar__1(45:49) .EQ. '') EW2N = rNULL__
        if(ar__1(53:56) .EQ. '') dA2N = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Ion,lambda,ExPot,log_gf,EW0,EW2S,f_EW2S,dA2S,EW2N,f_EW2N,dA2N
c    .......End.of.Just test output...........
      end do
      close(1)

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