FORTRAN Generation
(/./ftp/cats/J/ApJS/229/30)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/229/30 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/ApJS/229/30 Revised stellar properties of Q1-17 Kepler targets (Mathur+, 2017)
*================================================================================
*Revised stellar properties of Kepler targets for the Q1-17 (DR25) transit
*detection run.
*    Mathur S., Huber D., Batalha N.M., Ciardi D.R., Bastien F.A., Bieryla A.,
*    Buchhave L.A., Cochran W.D., Endl M., Esquerdo G.A., Furlan E., Howard A.,
*    Howell S.B., Isaacson H., Latham D.W., MacQueen P.J., Silva D.R.
*   <Astrophys. J. Suppl. Ser., 229, 30-30 (2017)>
*   =2017ApJS..229...30M    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Input values of the DR25 stellar properties catalog;
                             with the updated values for 779 KOIs from erratum
                             published in 2018, ApJS, 234, 43

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

      integer*4     KIC        (nr__) ! Kepler Input Catalog identifier
      integer*4     Teffi      (nr__) ! (K) [2500/27730] Effective temperature; input
      integer*4     e_Teffi    (nr__) ! (K) [50/1283] Uncertainty in Teff
      real*4        loggi      (nr__) ! ([cm/s2]) [-0.4/6.2] Log surface gravity; input
      real*4        e_loggi    (nr__) ! ([cm/s2]) [0.007/0.7] Uncertainty in log(g)
      real*4        v_Fe_H_i   (nr__) ! ([Sun]) [-2.6/1] Metallicity; input
      real*4        e__Fe_H_i  (nr__) ! ([Sun]) [0.06/2] Uncertainty in [Fe/H]
      character*3   n_Teffi    (nr__) ! Method for Teffi estimation (G1)
      integer*4     r_Teffi    (nr__) ! ? Reference for Teffi (see refs.dat file)
      character*3   n_loggi    (nr__) ! Method for loggi estimation (G1)
      integer*4     r_loggi    (nr__) ! ? Reference for loggi (see refs.dat file)
      character*3   n__Fe_H_i  (nr__) ! Method for [Fe/H]i estimation (G1)
      integer*4     r__Fe_H_i  (nr__) ! ? Reference for [Fe/H]i (see refs.dat file)

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

C  Declarations for 'table4.dat'	! Output values of the DR25 stellar properties
                             catalog with the updated distances and extinctions;
                             with the updated values for 779 KOIs from erratum
                             published in 2018, ApJS, 234, 43

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

      integer*4     KIC_1      (nr__1) ! Kepler Input Catalog identifier
      integer*4     Teff       (nr__1) ! (K) [2500/27730] Effective temperature; output
      integer*4     E_Teff     (nr__1) ! (K) [0/1015] Upper 1{sigma} confidence interval
*                                   in Teff
      integer*4     e_Teff_1   (nr__1) ! (K) [0/1939] Lower 1{sigma} confidence interval
*                                   in Teff
      real*4        log_g      (nr__1) ! ([cm/s2]) [-0.2/5.6] Log surface gravity; output
      real*4        E_log_g    (nr__1) ! ([cm/s2]) [0/2] Upper 1{sigma} confidence interval
*                                   in log(g)
      real*4        e_log_g_1  (nr__1) ! ([cm/s2]) [0/0.9] Lower 1{sigma} confidence interval
*                                   in log(g)
      real*4        v_Fe_H_    (nr__1) ! ([Sun]) [-2/0.6] Metallicity; output
      real*4        E__Fe_H_   (nr__1) ! ([Sun]) [0/1.4] Upper 1{sigma} confidence interval
*                                   in [Fe/H]
      real*4        e__Fe_H__1 (nr__1) ! ([Sun]) [0/1.1] Lower 1{sigma} confidence interval
*                                   in [Fe/H]
      real*8        Rad        (nr__1) ! (Rsun) [0.1/301] Radius
      real*4        E_Rad      (nr__1) ! (Rsun) [0/42] Upper 1{sigma} confidence interval
*                                   in Rad
      real*8        e_Rad_1    (nr__1) ! (Rsun) [0/180] Lower 1{sigma} confidence interval
*                                   in Rad
      real*4        Mass       (nr__1) ! (Msun) [0.08/3.8] Mass
      real*4        E_Mass     (nr__1) ! (Msun) [0/1.6] Upper 1{sigma} confidence interval
*                                   in Mass
      real*4        e_Mass_1   (nr__1) ! (Msun) [0/2.6] Lower 1{sigma} confidence interval
*                                   in Mass
      real*4        rho        (nr__1) ! [/117] Density
      real*4        E_rho      (nr__1) ! [0/24] Upper 1{sigma} confidence interval
*                                   in rho
      real*4        e_rho_1    (nr__1) ! [0/100] Lower 1{sigma} confidence interval
*                                   in rho
      real*8        Dist       (nr__1) ! (kpc) [0/218828] Distance
      real*8        E_Dist     (nr__1) ! (kpc) [0/45232] Upper 1{sigma} confidence interval
*                                   in Dist
      real*8        e_Dist_1   (nr__1) ! (kpc) [0/75587] Lower 1{sigma} confidence interval
*                                   in Dist
      real*4        Av         (nr__1) ! (mag) [0/1.1] V band extinction
      real*4        E_Av       (nr__1) ! (mag) [0/0.3] Upper 1{sigma} confidence interval
*                                   in Av
      real*4        e_Av_1     (nr__1) ! (mag) [0/0.8] Lower 1{sigma} confidence interval
*                                   in Av
      character*4   Mod        (nr__1) ! Provenance of output values (1)
      integer*4     Ref        (nr__1) ! ? Provenance of output values
*                                   (see refs.dat file)
*Note (1): Provenance of the output parameters as follows:
*   DSEP = results from the Dartmouth models,
*   BTSL = results from the BT-Settl models,
*   MULT = from the reference listed in Ref column.

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

C  Declarations for 'table6.dat'	! Possible M dwarfs according to
                             Gaidos+ 2016MNRAS.457.2877G

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

      integer*4     KIC_2      (nr__2) ! Kepler Input Catalog identifier
      integer*4     Teff_      (nr__2) ! (K) [3183/4026] Effective temperature (1)
      real*4        logg_      (nr__2) ! ([cm/s2]) [4.5/5.1] Log surface gravity (1)
      real*4        Rad_       (nr__2) ! (Rsun) [0.2/0.8] Radius (1)
      integer*4     Teff_1     (nr__2) ! (K) [3271/5052] DR25 effective temperature
      real*4        logg       (nr__2) ! ([cm/s2]) [0.05/3.1] DR25 log surface gravity
      real*4        Rad_1      (nr__2) ! (Rsun) [5.3/301] DR25 radius
      character*3   n_Teff     (nr__2) ! Method for Teff estimation (G1)
      integer*4     r_Teff     (nr__2) ! ? Reference for Teff (see refs.dat file)
      character*3   n_logg     (nr__2) ! Method for logg estimation (G1)
      integer*4     r_logg     (nr__2) ! ? Reference for logg (see refs.dat file)
*Note (1): From Gaidos et al. 2016MNRAS.457.2877G

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

C  Declarations for 'refs.dat'	! References (from table 5)

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

      integer*4     Ref_1      (nr__3) ! Reference code (from table 5)
      character*19  BibCode    (nr__3) ! Bibcode of the reference
      character*23  Auth       (nr__3) ! Author's name(s)
      character*60  Comm       (nr__3) ! Comment

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

C  Loading file 'table3.dat'	! Input values of the DR25 stellar properties catalog;
*                             with the updated values for 779 KOIs from erratum
*                             published in 2018, ApJS, 234, 43

C  Format for file interpretation

    1 format(
     +  I8,1X,I5,1X,I4,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,A3,I2,1X,A3,
     +  I2,1X,A3,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,197096
        read(1,'(A63)')ar__
        read(ar__,1)
     +  KIC(i__),Teffi(i__),e_Teffi(i__),loggi(i__),e_loggi(i__),
     +  v_Fe_H_i(i__),e__Fe_H_i(i__),n_Teffi(i__),r_Teffi(i__),
     +  n_loggi(i__),r_loggi(i__),n__Fe_H_i(i__),r__Fe_H_i(i__)
        if(ar__(50:51) .EQ. '') r_Teffi(i__) = iNULL__
        if(ar__(56:57) .EQ. '') r_loggi(i__) = iNULL__
        if(ar__(62:63) .EQ. '') r__Fe_H_i(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  KIC(i__),Teffi(i__),e_Teffi(i__),loggi(i__),e_loggi(i__),
     +  v_Fe_H_i(i__),e__Fe_H_i(i__),n_Teffi(i__),r_Teffi(i__),
     +  n_loggi(i__),r_loggi(i__),n__Fe_H_i(i__),r__Fe_H_i(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Output values of the DR25 stellar properties
*                             catalog with the updated distances and extinctions;
*                             with the updated values for 779 KOIs from erratum
*                             published in 2018, ApJS, 234, 43

C  Format for file interpretation

    2 format(
     +  I8,1X,I5,1X,I4,1X,I4,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F5.3,1X,F7.3,1X,F6.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
     +  E9.3,1X,E9.3,1X,E9.3,1X,F10.3,1X,F9.3,1X,F9.3,1X,F5.3,1X,F5.3,
     +  1X,F5.3,1X,A4,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,197096
        read(1,'(A189)')ar__1
        read(ar__1,2)
     +  KIC_1(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),log_g(i__),
     +  E_log_g(i__),e_log_g_1(i__),v_Fe_H_(i__),E__Fe_H_(i__),
     +  e__Fe_H__1(i__),Rad(i__),E_Rad(i__),e_Rad_1(i__),Mass(i__),
     +  E_Mass(i__),e_Mass_1(i__),rho(i__),E_rho(i__),e_rho_1(i__),
     +  Dist(i__),E_Dist(i__),e_Dist_1(i__),Av(i__),E_Av(i__),
     +  e_Av_1(i__),Mod(i__),Ref(i__)
        if(ar__1(188:189) .EQ. '') Ref(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  KIC_1(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),log_g(i__),
     +  E_log_g(i__),e_log_g_1(i__),v_Fe_H_(i__),E__Fe_H_(i__),
     +  e__Fe_H__1(i__),Rad(i__),E_Rad(i__),e_Rad_1(i__),Mass(i__),
     +  E_Mass(i__),e_Mass_1(i__),rho(i__),E_rho(i__),e_rho_1(i__),
     +  Dist(i__),E_Dist(i__),e_Dist_1(i__),Av(i__),E_Av(i__),
     +  e_Av_1(i__),Mod(i__),Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Possible M dwarfs according to
*                             Gaidos+ 2016MNRAS.457.2877G

C  Format for file interpretation

    3 format(
     +  I8,1X,I4,1X,F4.2,1X,F4.2,1X,I4,1X,F4.2,1X,F6.2,1X,A3,I2,1X,A3,
     +  I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,53
        read(1,'(A52)')ar__2
        read(ar__2,3)
     +  KIC_2(i__),Teff_(i__),logg_(i__),Rad_(i__),Teff_1(i__),
     +  logg(i__),Rad_1(i__),n_Teff(i__),r_Teff(i__),n_logg(i__),
     +  r_logg(i__)
        if(ar__2(45:46) .EQ. '') r_Teff(i__) = iNULL__
        if(ar__2(51:52) .EQ. '') r_logg(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  KIC_2(i__),Teff_(i__),logg_(i__),Rad_(i__),Teff_1(i__),
     +  logg(i__),Rad_1(i__),n_Teff(i__),r_Teff(i__),n_logg(i__),
     +  r_logg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References (from table 5)

C  Format for file interpretation

    4 format(I2,1X,A19,1X,A23,1X,A60)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,91
        read(1,'(A107)')ar__3
        read(ar__3,4)Ref_1(i__),BibCode(i__),Auth(i__),Comm(i__)
c    ..............Just test output...........
        write(6,4)Ref_1(i__),BibCode(i__),Auth(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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