Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/154/105 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 2024-Apr-19
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/105 Parameters of 529 Kepler eclipsing binaries (Kjurkchieva+, 2017)
*================================================================================
*Orbital parameters of the eclipsing detached Kepler binaries with eccentric
*orbits.
* Kjurkchieva D., Vasileva D., Atanasova T.
* <Astron. J., 154, 105-105 (2017)>
* =2017AJ....154..105K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! *Information about the targets from the Eclipsing
Binary (EB) catalog
integer*4 nr__
parameter (nr__=529) ! Number of records
character*99 ar__ ! Full-size record
integer*4 KIC ! [1026032/12644769] Kepler Input Catalog identifier
real*8 Per ! (d) [1.2/1087.3] Period
real*8 T0 ! (d) [5487.3/55733.8] Initial epoch T_0_ (BJD-2400000)
real*8 Kpmag ! (mag) [7.3/18.4] Kepler magnitude (K_m_)
integer*4 Teff ! (K) [3717/9224]? Target effective temperature (T_eff_)
character*1 f_Teff ! [B] Flag on Teff (B=Binary with unknown T_eff_)
real*4 w1 ! [0/0.6] First eclipse width, in phase units (w_1_)
character*1 f_w1 ! [1] Flag on w1 (1=value derived first by us)
real*4 w2 ! [0/0.19] Second eclipse width, in phase units
* (w_2_)
character*1 f_w2 ! [1] Flag on w2 (1=value derived first by us)
real*4 d1 ! [0/1] First eclipse depth, in relative flux units
* (d_1_)
character*1 f_d1 ! [1] Flag on d1 (1=value derived first by us)
real*8 d2 ! [0/0.42] Second eclipse depth, in relative flux
* units (d_2_)
character*1 f_d2 ! [1] Flag on d2 (1=value derived first by us)
real*8 Phase2 ! [0.03/0.96] Secondary eclipse phase ({phi}_2_) (1)
character*1 f_Phase2 ! [1] Flag on Phase2 (1=value derived first by us)
*Note (1): It is assumed {phi}_1_=0.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *Stellar and orbital parameters of the targets
integer*4 nr__1
parameter (nr__1=529) ! Number of records
character*54 ar__1 ! Full-size record
integer*4 KIC_1 ! [1026032/12644769] Kepler Input Catalog identifier
integer*4 Teff1 ! (K) [3717/9224] Effective temperature of primary
* (T_1_) (2)
character*1 f_Teff1 ! [B] Flag on Teff1 (B=Binary with unknown T_eff_) (3)
integer*4 Teff2 ! (K) [2039/7805] Effective temperature of secondary
* (T_2_) (4)
real*4 q ! [0.18/1.4] Mass ratio (4)
real*4 r1 ! [0/0.96] Stellar radii of primary in relative
* units (r_1_) (4)
real*4 r2 ! [0/0.91] Stellar radii of secondary in relative
* units (r_2_) (4)
real*4 e ! [0.001/0.85] Orbital eccentricity (4)
real*4 omega ! (deg) Periastron angle {omega} (4)
*Note (2): We used the available target temperatures T_eff_ (Table1) from the
* Kepler Input Catalog (Kepler Mission Team 2009, Cat. V/133) and adopted
* T_1_=T_eff_.
*Note (3): We assumed the mean (solar) value of 5800K for the primary.
*Note (4):
* Based on the assumption that the stellar components are MS stars, we
* calculated approximate values of their parameters from empirical relations
* (Ivanov et al. 2010BASI...38...83I):
* Secondary temperature T_2_=T_1_(d_2_/d_1_)^1/4^,
* Mass ratio q=(T_2_/T_1_)^1.7^,
* Ratio of relative radii k=r_2_/r_1_=q^0.75^.
*
* From the last formula and the approximate expression for narrow eclipses
* r_1_+r_2_~{pi}<w>. Eq. (3).
* (<w> is the mean eclipse width from Table1) we derived approximate values
* of relative stellar radii:
* r_1_={pi}<w>/(1+k) r_2_=[{pi}<w>/(1+k)]k. Eq. (4).
*
* The calculated values of T_2_, q, r_1_, and r_2_ were used as fixed
* parameters for PHOEBE ((Prsa & Zwitter 2005ApJ...628..426P) while the
* calculated values e_0_ and {omega}_0_ were used as input parameters. The
* orbital inclination i was fixed to correspond approximately to the
* individual eclipse depths. By varying e and {omega} around e_0_ and
* {omega}_0_ as well as "phase shift" we searched for the perfect coincidence
* of the phases of eclipses.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Comparison with orbital parameters of Helminiak
et al. (2016MNRAS.461.2896H)
integer*4 nr__2
parameter (nr__2=3) ! Number of records
character*33 ar__2 ! Full-size record
integer*4 KIC_2 ! [3120320/10001167] Kepler Input Catalog identifier
real*4 eH ! [0.034/0.26] Orbital eccentricity, Helminiak
* et al. 2016MNRAS.461.2896H (e_H_)
real*4 omegaH ! (deg) [213/350.5] Periastron angle, Helminiak et al.
* 2016MNRAS.461.2896H ({omega}_H_)
real*4 e_1 ! [0.049/0.26] Orbital eccentricity, this paper
* (e_our_)
real*4 omega_1 ! (deg) [28.5/168.8] Periastron angle, this paper
* ({omega}_our_)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Comparison with orbital parameters of Borkovits
et al. (2015MNRAS.448..946B)
integer*4 nr__3
parameter (nr__3=21) ! Number of records
character*36 ar__3 ! Full-size record
integer*4 KIC_3 ! [4940201/12356914] Kepler Input Catalog
* identifier
real*4 eB ! [0.001/0.42] Orbital eccentricity, Borkovits
* et al. 2015MNRAS.448..946B (e_B_)
integer*4 omB0 ! (deg) [106/359] Lower value of periastron angle
* range, Borkovits et al. 2015MNRAS.448..946B
* ({omega}_B_)
integer*4 omB1 ! (deg) [108/457]? Upper value of periastron angle
* range, Borkovits et al. 2015MNRAS.448..946B
* ({omega}_B_)
real*4 e_2 ! [0.0025/0.44] Orbital eccentricity, this paper
* (e_our_)
real*4 omega_2 ! (deg) [14.1/346.6] Periastron angle, this paper
* ({omega}_our_)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Orbital parameters of eclipsing binary stars with
transiting circumbinary planet
integer*4 nr__4
parameter (nr__4=9) ! Number of records
character*92 ar__4 ! Full-size record
integer*4 KIC_4 ! [4862625/12644769] Kepler Input Catalog identifier
character*11 Kepler ! Kepler name, star with transiting circumbinary
* planet(s)
real*4 eP ! [0.02/0.53] Previous orbital eccentricity (e_pr_)
real*4 omegaP ! (deg) [89.1/300.9] Previous periastron angle
* ({omega}_pr_)
character*45 Ref ! Reference for previous eccentricity &
* periastron angle
real*4 e_3 ! [0.04/0.6] Orbital eccentricity, this paper
* (e_our_)
real*4 omega_3 ! (deg) [53/266.5] Periastron angle, this paper
* ({omega}_our_)
C=============================================================================
C Loading file 'table1.dat' ! *Information about the targets from the Eclipsing
* Binary (EB) catalog
C Format for file interpretation
1 format(
+ I8,1X,F12.7,1X,F12.6,1X,F10.4,1X,I4,1X,A1,1X,F6.4,1X,A1,1X,
+ F6.4,1X,A1,1X,F6.4,1X,A1,1X,F7.5,1X,A1,1X,F7.5,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,529
read(1,'(A99)')ar__
read(ar__,1)
+ KIC,Per,T0,Kpmag,Teff,f_Teff,w1,f_w1,w2,f_w2,d1,f_d1,d2,f_d2,
+ Phase2,f_Phase2
if(ar__(47:50) .EQ. '') Teff = iNULL__
c ..............Just test output...........
write(6,1)
+ KIC,Per,T0,Kpmag,Teff,f_Teff,w1,f_w1,w2,f_w2,d1,f_d1,d2,f_d2,
+ Phase2,f_Phase2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! *Stellar and orbital parameters of the targets
C Format for file interpretation
2 format(
+ I8,1X,I4,1X,A1,1X,I4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F5.3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,529
read(1,'(A54)')ar__1
read(ar__1,2)KIC_1,Teff1,f_Teff1,Teff2,q,r1,r2,e,omega
c ..............Just test output...........
write(6,2)KIC_1,Teff1,f_Teff1,Teff2,q,r1,r2,e,omega
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Comparison with orbital parameters of Helminiak
* et al. (2016MNRAS.461.2896H)
C Format for file interpretation
3 format(I8,1X,F5.3,1X,F5.1,1X,F5.3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,3
read(1,'(A33)')ar__2
read(ar__2,3)KIC_2,eH,omegaH,e_1,omega_1
c ..............Just test output...........
write(6,3)KIC_2,eH,omegaH,e_1,omega_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Comparison with orbital parameters of Borkovits
* et al. (2015MNRAS.448..946B)
C Format for file interpretation
4 format(I8,1X,F5.3,1X,I3,1X,I3,1X,F6.4,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,21
read(1,'(A36)')ar__3
read(ar__3,4)KIC_3,eB,omB0,omB1,e_2,omega_2
if(ar__3(20:22) .EQ. '') omB1 = iNULL__
c ..............Just test output...........
write(6,4)KIC_3,eB,omB0,omB1,e_2,omega_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Orbital parameters of eclipsing binary stars with
* transiting circumbinary planet
C Format for file interpretation
5 format(I8,1X,A11,1X,F5.3,1X,F6.2,1X,A45,1X,F5.3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,9
read(1,'(A92)')ar__4
read(ar__4,5)KIC_4,Kepler,eP,omegaP,Ref,e_3,omega_3
c ..............Just test output...........
write(6,5)KIC_4,Kepler,eP,omegaP,Ref,e_3,omega_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end