FORTRAN Generation
(/./ftp/cats/J/AJ/152/180)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/152/180 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/152/180 Bolometric fluxes of eclipsing binaries in Tycho-2 (Stassun+, 2016)
*================================================================================
*Eclipsing binaries as benchmarks for trigonometric parallaxes in the Gaia era.
*    Stassun K.G., Torres G.
*   <Astron. J., 152, 180-180 (2016)>
*   =2016AJ....152..180S    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Stellar properties of our sample of eclipsing binaries

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

      character*26  Name       (nr__) ! System name (G1)
      integer*4     TYC1       (nr__) ! TYC1 of Tycho name
      integer*4     TYC2       (nr__) ! TYC2 of Tycho name
      integer*4     TYC3       (nr__) ! TYC3 of Tycho name
      integer*4     HIP        (nr__) ! [817/114484]? Hipparcos identifier
      real*8        Rad1       (nr__) ! (Rsun) [0.4323/52.12] Radius of primary
      real*4        e_Rad1     (nr__) ! (Rsun) [0.002/1.38] Uncertainty in Rad1
      real*8        Rad2       (nr__) ! (Rsun) [0.3916/67.63] Radius of secondary
      real*4        e_Rad2     (nr__) ! (Rsun) [0.001/1.4] Uncertainty in Rad2
      integer*4     Teff1      (nr__) ! (K) [3160/38000] Effective temperature of primary
      integer*4     e_Teff1    (nr__) ! (K) [45/2500] Uncertainty in Teff1
      integer*4     Teff2      (nr__) ! (K) [3125/34100] Effective temperature of secondary
      integer*4     e_Teff2    (nr__) ! (K) [50/2500] Uncertainty in Teff2
      real*4        E_B_V      (nr__) ! (mag) [0/0.85]? Reddening
      real*4        e_E_B_V    (nr__) ! (mag) [0.001/0.08]? Uncertainty in E(B-V)
      character*7   v_Fe_H_    (nr__) ! ([Sun]) ? Metallicity of the binary (1)
      real*4        e__Fe_H_   (nr__) ! ([Sun]) [0.01/0.21]? Uncertainty in [Fe/H]
      character*5   f__Fe_H_   (nr__) ! Metallic-line star flag (Am/Am, Am/-, -/Am) (2)
      character*5   Ref        (nr__) ! Literature reference; in refs.dat file
*Note (1):
*     Am systems having components with individually measured [Fe/H]:
*     YZ Cas    (+0.54+/-0.11, Am)+(+0.01+/-0.11, normal);
*     V501 Mon  (+0.33+/-0.08, Am)+(+0.01+/-0.06, normal);
*     SW CMa    (+0.49+/-0.15, Am)+(+0.61+/-0.15, Am);
*     HW CMa    (+0.28+/-0.10, Am)+(+0.33+/-0.10, Am).
*     All uncertainties correspond to 1{sigma} errors.
*Note (2): Metallicity flags are defined as follows:
*     Am/Am = Both components are metallic-line stars;
*      Am/- = Only the primary is known to have that anomaly.

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

C  Declarations for 'table2.dat'	! Results

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

      character*26  Name_1     (nr__1) ! System name (G1)
      integer*4     TYC1_1     (nr__1) ! TYC1 of Tycho name
      integer*4     TYC2_1     (nr__1) ! TYC2 of Tycho name
      integer*4     TYC3_1     (nr__1) ! TYC3 of Tycho name
      real*4        Vmag       (nr__1) ! (mag) [1.9/12.43] V band magnitude
      real*4        chi2       (nr__1) ! [0.25/85.13] SED fit reduced {chi}_{nu}_^2^ value
      character*1   fit        (nr__1) ! [*X] Fit flag(s) (* or X) (3)
      real*4        Fbol       (nr__1) ! (mW/m2) Integrated bolometric flux F_bol_ (in erg/s/cm^2^)
      real*4        E_Fbol     (nr__1) ! (mW/m2) Upper uncertainty on Fbol
      real*4        e_Fbol_1   (nr__1) ! (mW/m2) Lower uncertainty on Fbol
      real*4        AV         (nr__1) ! (mag) [0/2.29] Extinction A_V_
      real*4        E_AV       (nr__1) ! (mag) [0/0.19] Upper uncertainty on AV
      real*4        e_AV_1     (nr__1) ! (mag) [0/0.31] Lower uncertainty on AV
      real*4        Dist       (nr__1) ! (pc) [14.2/5231.2] Predicted distance (4)
      real*4        E_Dist     (nr__1) ! (pc) [0.9/769.2] Upper uncertainty on Dist
      real*4        e_Dist_1   (nr__1) ! (pc) [0.9/712.1] Lower uncertainty on Dist
      real*4        plx        (nr__1) ! (mas) [0.19/70.58] Predicted parallax
      real*4        E_plx      (nr__1) ! (mas) [0.02/8.49] Upper uncertainty on plx
      real*4        e_plx_1    (nr__1) ! (mas) [0.02/7.38] Lower uncertainty on plx
*Note (3): SED Fit flags are defined as follows:
*     * = Systems with Teff>15000~K and a large fraction of their flux
*         extrapolated from a blackbody distribution on the blue side;
*     X = Considered to have poor SED fits.
*Note (4): With F_bol_ and L_bol_ in hand for each eclipsing binary, we can
*     calculate the predicted distance to each eclipsing binary according to
*     Equation (1): Dist=(L_bol_/4{pi}F_bol_)^1/2^.

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref_1      (nr__2) ! Reference identifier
      character*19  BibCode    (nr__2) ! Bibliographic Code
      character*21  Aut        (nr__2) ! Author's name
      character*24  Cat        (nr__2) ! VizieR catalog

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

C  Loading file 'table1.dat'	! Stellar properties of our sample of eclipsing binaries

C  Format for file interpretation

    1 format(
     +  A26,1X,I4,1X,I4,1X,I1,1X,I6,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,
     +  1X,I5,1X,I4,1X,I5,1X,I4,1X,F5.3,1X,F5.3,1X,A7,1X,F4.2,1X,A5,
     +  1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,158
        read(1,'(A134)')ar__
        read(ar__,1)
     +  Name(i__),TYC1(i__),TYC2(i__),TYC3(i__),HIP(i__),Rad1(i__),
     +  e_Rad1(i__),Rad2(i__),e_Rad2(i__),Teff1(i__),e_Teff1(i__),
     +  Teff2(i__),e_Teff2(i__),E_B_V(i__),e_E_B_V(i__),v_Fe_H_(i__),
     +  e__Fe_H_(i__),f__Fe_H_(i__),Ref(i__)
        if(ar__(40:45) .EQ. '') HIP(i__) = iNULL__
        if(ar__(99:103) .EQ. '') E_B_V(i__) = rNULL__
        if(ar__(105:109) .EQ. '') e_E_B_V(i__) = rNULL__
        if(ar__(119:122) .EQ. '') e__Fe_H_(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),TYC1(i__),TYC2(i__),TYC3(i__),HIP(i__),Rad1(i__),
     +  e_Rad1(i__),Rad2(i__),e_Rad2(i__),Teff1(i__),e_Teff1(i__),
     +  Teff2(i__),e_Teff2(i__),E_B_V(i__),e_E_B_V(i__),v_Fe_H_(i__),
     +  e__Fe_H_(i__),f__Fe_H_(i__),Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Results

C  Format for file interpretation

    2 format(
     +  A26,1X,I4,1X,I4,1X,I1,1X,F5.2,1X,F5.2,1X,A1,1X,E9.3,1X,E7.1,
     +  1X,E7.1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.1,1X,F5.1,1X,F5.1,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,158
        read(1,'(A128)')ar__1
        read(ar__1,2)
     +  Name_1(i__),TYC1_1(i__),TYC2_1(i__),TYC3_1(i__),Vmag(i__),
     +  chi2(i__),fit(i__),Fbol(i__),E_Fbol(i__),e_Fbol_1(i__),
     +  AV(i__),E_AV(i__),e_AV_1(i__),Dist(i__),E_Dist(i__),
     +  e_Dist_1(i__),plx(i__),E_plx(i__),e_plx_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),TYC1_1(i__),TYC2_1(i__),TYC3_1(i__),Vmag(i__),
     +  chi2(i__),fit(i__),Fbol(i__),E_Fbol(i__),e_Fbol_1(i__),
     +  AV(i__),E_AV(i__),e_AV_1(i__),Dist(i__),E_Dist(i__),
     +  e_Dist_1(i__),plx(i__),E_plx(i__),e_plx_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(I2,1X,A19,1X,A21,1X,A24)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,73
        read(1,'(A69)')ar__2
        read(ar__2,3)Ref_1(i__),BibCode(i__),Aut(i__),Cat(i__)
c    ..............Just test output...........
        write(6,3)Ref_1(i__),BibCode(i__),Aut(i__),Cat(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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