Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/832/183 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-Apr-20
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/832/183 Binary stellar evolution data for Kepler systems (Kostov+, 2016)
*================================================================================
*Tatooine's future: the eccentric response of Kepler's circumbinary planets to
*common-envelope evolution of their host stars.
* Kostov V.B., Moore K., Tamayo D., Jayawardhana R., Rinehart S.A.
* <Astrophys. J., 832, 183-183 (2016)>
* =2016ApJ...832..183K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Binary Stellar Evolution (BSE) results for the
entire set of simulations for all Kepler systems
integer*4 nr__
parameter (nr__=760) ! Number of records
character*113 ar__ ! Full-size record
character*10 Name (nr__) ! Kepler system name
real*4 MA (nr__) ! (Msun) [0.6/1.6] Primary star, initial mass
real*4 MB (nr__) ! (Msun) [0.2/1.1] Secondary star, initial mass
real*4 Z (nr__) ! ([Sun]) [0.009/0.03] Metallicity
real*4 alpha (nr__) ! [0.5/10] Common Envelope (CE) efficiency parameter
character*4 TCP (nr__) ! Tides model applied (1)
character*11 deKool (nr__) ! de Kool CE Model: ON/OFF (always "DE KOOL OFF") (2)
real*8 Time (nr__) ! (Myr) [0/25000] Evolutionary time
real*4 M1 (nr__) ! (Msun) [0/2.5] Current mass, Primary star
real*4 M2 (nr__) ! (Msun) [0/1.1] Current mass, Secondary star
integer*4 K1 (nr__) ! [1/15] Current Stellar Type, Primary star (3)
integer*4 K2 (nr__) ! [0/15]? Current Stellar Type, Secondary star (3)
real*4 a (nr__) ! (Rsun) [0/49.3] Binary orbit, a_bin_
real*4 e (nr__) ! [0/0.6]?=-1 System eccentricity e_bin_
real*4 R1 (nr__) ! [0/8.8]?=-1 Fraction of Roche lobe filled,
* Primary
real*4 R2 (nr__) ! [0/2]?=-1 Fraction of Roche lobe filled,
* Secondary
character*8 Evol (nr__) ! Evolutionary Stage (4)
*Note (1): Tides model applied as follows:
* TCP = Tidal Circularization Path. Tides "ON" in BSE;
* NTCP = No Tidal Circularization Path. Tides "OFF" in BSE.
*Note (2): The BSE code also allows an alternate CE model to be used, the
* de Kool CE evolution model (de Kool 1990ApJ...358..189D), which first
* introduced the CE evolution binding energy factor {lambda}.
* See section 2.
*Note (3): Stellar Type (from Table 2) as follows:
* 0 = Deeply or fully convective low-mass MS star;
* 1 = Main Sequence star;
* 2 = Hertzsprung Gap (HG);
* 3 = First Giant Branch (GB);
* 4 = Core Helium Burning (CHeB);
* 5 = First/Early Asymptotic Giant Branch (EAGB);
* 6 = Second/Thermally Pulsing Asymptotic Giant Branch (TPAGB);
* 7 = Main Sequence Naked Helium star (HeMS);
* 8 = Hertzsprung Gap Naked Helium star (HeHG);
* 9 = Giant Branch Naked Helium Star (HeGB);
* 10 = Helium White Dwarf (HeWD);
* 11 = Carbon/Oxygen White Dwarf (COWD);
* 12 = Oxygen/Neon White Dwarf (ONeWD);
* 13 = Neutron Star (NS);
* 14 = Black Hole (BH);
* 15 = Massless Supernova/Massless Remnant.
*Note (4): Evolutionary Stage (from Table 3) as follows:
* INITIAL = Initial configuration;
* KW CHNGE = Stellar type change;
* BEG RCHE = Begin Roche lobe overflow;
* END RCHE = End Roche lobe overflow;
* CONTACT = Contact system;
* COELESCE = Coalescence of stars;
* COMENV = Common-envelope system;
* GNTAGE = New giant star from CE; appropriate age and initial mass set
* to match core-mass and stellar mass;
* NO REMNT = No remnant;
* MAX TIME = Max evolutionary time reached; end of program;
* DISRUPT = System is disrupted;
* BEG SYMB = Begin symbiotic system;
* END SYMB = End symbiotic system;
* BEG BSS = Begin blue stragglers.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Kepler system default initial parameters
integer*4 nr__1
parameter (nr__1=11) ! Number of records
character*69 ar__1 ! Full-size record
character*11 Name_1 (nr__1) ! Kepler system name
character*1 f_Name (nr__1) ! [ac] Flag on Name (1)
integer*4 KIC (nr__1) ! [4862625/12644769] Kepler ID number
real*4 M1_1 (nr__1) ! (Msun) [0.6/1.6] Primary mass
real*4 M2_1 (nr__1) ! (Msun) [0.2/1.1] Secondary mass
real*4 Pbin (nr__1) ! (d) [7.4/41.1] Binary period
real*4 e_1 (nr__1) ! [0.02/0.6] Binary eccentricity e_bin_
real*4 Z_1 (nr__1) ! ([Sun]) [0.009/0.4] Metallicity Z
character*1 f_Z (nr__1) ! [b] Flag on Z (1)
real*8 Pcbp (nr__1) ! (d) [49.5/1107.6] Circumbinary planet period,
* P_CBP_
real*4 ecbp (nr__1) ! [0.01/0.2] Eccentricity of the circumbinary
* planet, e_CBP_
character*2 n_ecbp (nr__1) ! Note on ecbp (2)
*Note (1): Flag as follows:
* a = J. Orosz (2016, private communication).
* b = BSE allows maximum metallicity of Z=0.03.
* c = Kepler-1647
*Note (2): P1/P2/P3 have 60%/15%/85% probability to become dynamically unstable
* within 1Myr after the Common Envelope phase.
C=============================================================================
C Loading file 'table4.dat' ! Binary Stellar Evolution (BSE) results for the
* entire set of simulations for all Kepler systems
C Format for file interpretation
1 format(
+ A10,1X,F4.2,1X,F4.2,1X,F5.3,1X,F4.1,1X,A4,1X,A11,1X,F10.4,1X,
+ F5.3,1X,F5.3,1X,I2,1X,I2,1X,F6.3,1X,F5.2,1X,F6.3,1X,F6.3,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,760
read(1,'(A113)')ar__
read(ar__,1)
+ Name(i__),MA(i__),MB(i__),Z(i__),alpha(i__),TCP(i__),
+ deKool(i__),Time(i__),M1(i__),M2(i__),K1(i__),K2(i__),a(i__),
+ e(i__),R1(i__),R2(i__),Evol(i__)
if(ar__(76:77) .EQ. '') K2(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),MA(i__),MB(i__),Z(i__),alpha(i__),TCP(i__),
+ deKool(i__),Time(i__),M1(i__),M2(i__),K1(i__),K2(i__),a(i__),
+ e(i__),R1(i__),R2(i__),Evol(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Kepler system default initial parameters
C Format for file interpretation
2 format(
+ A11,1X,A1,1X,I8,1X,F4.2,1X,F6.4,1X,F6.3,1X,F4.2,1X,F5.3,A1,1X,
+ F7.2,1X,F4.2,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,11
read(1,'(A69)')ar__1
read(ar__1,2)
+ Name_1(i__),f_Name(i__),KIC(i__),M1_1(i__),M2_1(i__),
+ Pbin(i__),e_1(i__),Z_1(i__),f_Z(i__),Pcbp(i__),ecbp(i__),
+ n_ecbp(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),f_Name(i__),KIC(i__),M1_1(i__),M2_1(i__),
+ Pbin(i__),e_1(i__),Z_1(i__),f_Z(i__),Pcbp(i__),ecbp(i__),
+ n_ecbp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end