Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/167/103 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-Sep-09
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/167/103 System properties & impact parameter variations (Judkovsky+, 2024)
*================================================================================
*Kepler Multitransiting System Physical Properties and Impact Parameter
*Variations.
* Judkovsky Y., Ofir A., Aharonson O.
* <Astron. J., 167, 103 (2024)>
* =2024AJ....167..103J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Physical and orbital elements of our solutions that
pass all tests
integer*4 nr__
parameter (nr__=470) ! Number of records
character*231 ar__ ! Full-size record
character*1 Subset ! Source, A=Paper II; B=this work (1)
integer*4 Adopted ! [0/1] 1=Adopted solution
integer*4 ID ! [1/5] Run identifier
character*13 Name ! Kepler Planet Name
real*8 KOI ! Kepler object of interest
real*8 Per ! (d) [0.7/500] Mean orbital period
real*8 Mp_M_ ! (10-6) [0.73/1429] Planet-to-star mass ratio
real*4 E_Mp_M_ ! (10-6) [0.1/58] Upper uncertainty on Mp/M*
real*4 e_Mp_M__1 ! (10-6) [-67/0] Lower uncertainty on Mp/M*
real*8 Mp ! (Mgeo) [0.2/597] Absolute planet mass (2)
real*8 E_Mp ! (Mgeo) [0.04/37] Upper uncertainty on Mp
real*8 e_Mp_1 ! (Mgeo) [-45/-0.04] Lower uncertainty on Mp
real*8 rhop ! (g.cm-3) [0.04/49] Absolute planet density (3)
real*8 E_rhop ! (g.cm-3) [0.004/20] Upper uncertainty on rhop
real*8 e_rhop_1 ! (g.cm-3) [-45/-0.004] Lower uncertainty on rhop
real*4 Radp ! (Rgeo) [0.96/12.1] Absolute planet radius (4)
real*4 E_Radp ! (Rgeo) [0.01/2] Upper uncertainty on Radp
real*4 e_Radp_1 ! (Rgeo) [-0.6/-0.02] Lower uncertainty on Radp
real*8 dex ! [-0.3/0.26] Eccentricity difference to
* preceding planet in x (5)
real*8 E_dex ! [0.0005/0.2] Upper uncertainty on dex
real*8 e_dex_1 ! [-0.2/-0.0006] Lower uncertainty on dex
real*8 dey ! [-0.47/0.58] Eccentricity difference to
* preceding planet in y (6)
real*8 E_dey ! [0.0006/0.2] Upper uncertainty on dey
real*8 e_dey_1 ! [-0.5/-0.0006] Lower uncertainty on dey
real*4 Ix ! (deg) [-31.4/34.8] Inclination around the line of
* sight ("roll angle") (7)
real*4 E_Ix ! (deg) [0/23] Upper uncertainty on Ix
real*4 e_Ix_1 ! (deg) [-21/0] Lower uncertainty on Ix
real*8 Iy ! (deg) [-4.39/7.1] Inclination relative to the line
* of sight (8)
real*8 E_Iy ! (deg) [0.004/7] Upper uncertainty on Iy
real*8 e_Iy_1 ! (deg) [-3/-0.004] Lower uncertainty on Iy
*Note (1): This table combines results from (Subset=A) Table A1 of
* Paper II (Judkovsky+, 2022AJ....163...91J) for Keplers two- and
* three-transiting planets systems and results from this work (Subset=B)
* for Kepler systems containing four or more transiting planets. The
* results for Subset=A are reproduced exactly from Paper II.
*Note (2): Absolute planet mass derived from Mp/M*, the planet-to-star mass
* ratio, and literature stellar mass.
*Note (3): Absolute planet density derived from RhoRel, the
* planet-to-star density ratio and literature stellar density.
*Note (4): Absolute planet radius derived from ror, the planet-to-star
* radius ratio and literature stellar radius.
*Note (5): Eccentricity difference to preceding planet; x is pointing at
* the observer.
*Note (6): Eccentricity difference to preceding planet; y is the
* direction perpendicular to the line of sight pointing with the
* planets motion along transit.
*Note (7): Ix=I*cos(Omega) is the inclination component corresponding to
* the inclination around to the line of sight ("roll angle"), where
* I is the inclination and Omega is the longitude of ascending node.
*Note (8): Iy=I*sin(Omega) is the inclination component corresponding to
* the inclination relative to the line of sight, where
* I is the inclination and Omega is the longitude of ascending node.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Light curve transit and MMR-proximity parameters of
all analyzed planets
integer*4 nr__1
parameter (nr__1=2324) ! Number of records
character*194 ar__1 ! Full-size record
character*1 Subset_1 ! [AB] Source, A=Paper II; B=this work (1)
integer*4 Adopted_1 ! [0/1] 1=Adopted solution
integer*4 ID_1 ! Run identifier
character*13 Name_1 ! Kepler Planet Name
real*8 KOI_1 ! Kepler object of interest
real*8 Per_1 ! (d) [0.7/1072] Mean orbital period
character*10 TTVStd ! (min) Best model TTV std
integer*4 J ! [0/5] J of the nearest first-order resonance
* J:J-1 (2)
real*8 Delta ! [-1/0.48] Normalized distance from the J:J-1
* resonance (3)
real*8 SPer ! (d) [0/6612] Super-period of the J:J-1 resonance (4)
real*8 T0 ! (h) [0/18] Transit duration at first data time stamp
real*8 E_T0 ! (h) [0.002/1.13] Upper uncertainty on T0
real*8 e_T0_1 ! (h) [-2/0] Lower uncertainty on T0
real*8 Tau0 ! (min) [0.44/161] ingress-egress time at first data
* time stamp
real*8 E_Tau0 ! (min) [0/25] Upper uncertainty on Tau0
real*8 e_Tau0_1 ! (min) [-21/0] Lower uncertainty on Tau0
real*8 b0 ! [-1.01/1.01] impact parameter at BKJD=0
real*8 E_b0 ! [0/2] Upper uncertainty on b0
real*8 e_b0_1 ! [-2/0] Lower uncertainty on b0
real*8 dbdt ! (yr-1) [-0.82/0.82] db/dt median
real*8 E_dbdt ! (yr-1) [0/0.1]Upper uncertainty on db/dt
real*8 e_dbdt_1 ! (yr-1) [-0.2/0] Lower uncertainty on db/dt
*Note (1): This table combines results from (Subset=A) Table A2 of
* Paper II (Judkovsky+, 2022AJ....163...91J) for Keplers two- and
* three-transiting planets systems and results from this work (Subset=B)
* for Kepler systems containing four or more transiting planets. The
* results for Subset=A are a slight revision of the earlier results,
* including a slight correction to the values of T0, Tau0, and b0.
*Note (2): J of the nearest first-order resonance J:J-1 with the planet
* from the inside. If the planet is innermost, the value is 0.
*Note (3): Normalized distance from the J:J-1 resonance with the planet
* from the inside. If the planet is innermost, the value is 0.
*Note (4): Super-period of the J:J-1 resonance with the planet from the
* inside. If the planet is innermost, the value is 0.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Stellar parameters of the systems for which valid
dynamical solutions were found
integer*4 nr__2
parameter (nr__2=77) ! Number of records
character*103 ar__2 ! Full-size record
character*1 Subset_2 ! Source, A=Paper II; B=this work (1)
integer*4 KOI_2 ! Kepler object of interest
character*10 Name_2 ! Kepler Star Name
integer*4 KIC ! Kepler Input Catalog number
character*19 Ref ! Reference (ADS bibcode or NExScI) for Mstar, Rstar
real*4 u1 ! [0.23/0.67] first limb-darkening coeff, taken
* from NExScI
real*4 u2 ! [0.07/0.45] second limb-darkening coeff, taken
* from NExScI
real*4 Mass_ ! (Msun) [0.5/1.5] Absolute stellar mass, from Ref
real*4 E_Mass_ ! (Msun) [0.01/0.2] Upper uncertainty on Mass*
real*4 e_Mass__1 ! (Msun) [-0.2/-0.012] Lower uncertainty on Mass*
real*4 Rad_ ! (Rsun) [0.51/1.84] Absolute stellar radius from Ref
real*8 E_Rad_ ! (Rsun) [0.006/0.5] Upper uncertainty on Rad*
real*8 e_Rad__1 ! (Rsun) [-0.3/-0.009] Lower uncertainty on Rad*
*Note (1): This table combines results from (Subset=A) Table A3 of
* Paper II (Judkovsky+, 2022AJ....163...91J) for Keplers two- and
* three-transiting planets systems and results from this work (Subset=B)
* for Kepler systems containing four or more transiting planets. The
* results for Subset=A are reproduced exactly from Paper II.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Planets with significant impact parameter
variations
integer*4 nr__3
parameter (nr__3=131) ! Number of records
character*98 ar__3 ! Full-size record
character*1 Subset_3 ! Source, A=Paper II; B=this work (1)
integer*4 Adopted_2 ! [0/1] 1=Adopted solution
character*12 Name_3 ! Kepler Planet Name
real*8 KOI_3 ! Kepler object of interest
real*8 Per_2 ! (d) [0.7/500] mean orbital period
integer*4 Npl ! [2/6] model number of planets
integer*4 Pos ! [1/6] planet position within the system, from
* inside to outside
real*8 b0_1 ! [-0.94/0.92] impact parameter at BKJD=0
real*8 E_b0_2 ! [0.0004/0.2] Upper uncertainty on b0
real*8 e_b0_3 ! [-0.3/-0.0004] Lower uncertainty on b0
real*8 dbdt_1 ! (yr-1) [-0.77/0.09] db/dt median
real*8 E_dbdt_2 ! (yr-1) [3e-05/0.05] Upper uncertainty on db/dt
real*8 e_dbdt_3 ! (yr-1) [-0.04/-3e-05] Lower uncertainty on db/dt
*Note (1): This table combines results from (Subset=A) Table A4 of
* Paper II (Judkovsky+, 2022AJ....163...91J) for Keplers two- and
* three-transiting planets systems and results from this work (Subset=B)
* for Kepler systems containing four or more transiting planets. The
* results for Subset=A are a slight revision of the earlier results,
* including a slight correction to the value of b0.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Instantaneous coordinates, velocities, and orbital
elements of the best-fitting parameters for the
adopted solutions in this work and in paper II
integer*4 nr__4
parameter (nr__4=241) ! Number of records
character*333 ar__4 ! Full-size record
real*8 x ! (m) Instantaneous Cartesian x coordinate (1)
real*8 y ! (m) Instantaneous Cartesian y coordinate (1)
real*8 z ! (m) Instantaneous Cartesian z coordinate (1)
real*8 vx ! (m.s-1) Instantaneous Cartesian x velocity
real*8 vy ! (m.s-1) Instantaneous Cartesian y velocity
real*8 vz ! (m.s-1) Instantaneous Cartesian z velocity
character*20 Mass__1 ! (kg) Star mass
character*20 Massp ! (kg) Planet mass
integer*4 T0_1 ! (d) [120/353] Epoch, days since Barycentric Kepler
* Julian date
real*8 a ! (m) Orbital semi-major axis
real*8 e ! Eccentricity
real*8 Pomega ! (rad) Longitude of periastron
real*8 i ! (rad) Orbital inclination
real*8 Omega ! (rad) Longitude of the ascending node
real*8 Mean ! (rad) Mean anomaly
integer*4 Pos_1 ! [1/6] Planet position within the system,
* from inside to outside
real*8 KOI_4 ! Kepler object of interest
character*12 Name_4 ! Kepler Planet Name
character*20 sigma ! N-body {sigma}
*Note (1): The axis system is such that yz forms the plane of the sky,
* and x points towards the observer (transit occurs when x>0).
C=============================================================================
C Loading file 'table1.dat' ! Physical and orbital elements of our solutions that
* pass all tests
C Format for file interpretation
1 format(
+ A1,1X,I1,1X,I1,1X,A13,F7.2,1X,F12.8,1X,F7.2,1X,F6.2,1X,F6.2,
+ 1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.4,1X,F8.4,1X,F8.4,1X,F6.3,1X,
+ F6.3,1X,F6.3,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,
+ 1X,F6.2,1X,F6.2,1X,F6.2,1X,F7.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,470
read(1,'(A231)')ar__
read(ar__,1)
+ Subset,Adopted,ID,Name,KOI,Per,Mp_M_,E_Mp_M_,e_Mp_M__1,Mp,
+ E_Mp,e_Mp_1,rhop,E_rhop,e_rhop_1,Radp,E_Radp,e_Radp_1,dex,
+ E_dex,e_dex_1,dey,E_dey,e_dey_1,Ix,E_Ix,e_Ix_1,Iy,E_Iy,e_Iy_1
c ..............Just test output...........
write(6,1)
+ Subset,Adopted,ID,Name,KOI,Per,Mp_M_,E_Mp_M_,e_Mp_M__1,Mp,
+ E_Mp,e_Mp_1,rhop,E_rhop,e_rhop_1,Radp,E_Radp,e_Radp_1,dex,
+ E_dex,e_dex_1,dey,E_dey,e_dey_1,Ix,E_Ix,e_Ix_1,Iy,E_Iy,e_Iy_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Light curve transit and MMR-proximity parameters of
* all analyzed planets
C Format for file interpretation
2 format(
+ A1,2X,I1,1X,I1,1X,A13,1X,F7.2,1X,F13.8,1X,A10,1X,I1,1X,F11.8,
+ 1X,F12.7,1X,F7.4,1X,F7.4,1X,F7.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.5,1X,F8.5,1X,F8.5,1X,F11.8,1X,F11.8,1X,F11.8)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,2324
read(1,'(A194)')ar__1
read(ar__1,2)
+ Subset_1,Adopted_1,ID_1,Name_1,KOI_1,Per_1,TTVStd,J,Delta,
+ SPer,T0,E_T0,e_T0_1,Tau0,E_Tau0,e_Tau0_1,b0,E_b0,e_b0_1,dbdt,
+ E_dbdt,e_dbdt_1
c ..............Just test output...........
write(6,2)
+ Subset_1,Adopted_1,ID_1,Name_1,KOI_1,Per_1,TTVStd,J,Delta,
+ SPer,T0,E_T0,e_T0_1,Tau0,E_Tau0,e_Tau0_1,b0,E_b0,e_b0_1,dbdt,
+ E_dbdt,e_dbdt_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Stellar parameters of the systems for which valid
* dynamical solutions were found
C Format for file interpretation
3 format(
+ A1,1X,I4,1X,A10,1X,I8,1X,A19,1X,F6.4,1X,F6.4,1X,F5.3,1X,F6.3,
+ 1X,F6.3,1X,F6.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,77
read(1,'(A103)')ar__2
read(ar__2,3)
+ Subset_2,KOI_2,Name_2,KIC,Ref,u1,u2,Mass_,E_Mass_,e_Mass__1,
+ Rad_,E_Rad_,e_Rad__1
c ..............Just test output...........
write(6,3)
+ Subset_2,KOI_2,Name_2,KIC,Ref,u1,u2,Mass_,E_Mass_,e_Mass__1,
+ Rad_,E_Rad_,e_Rad__1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Planets with significant impact parameter
* variations
C Format for file interpretation
4 format(
+ A1,1X,I1,1X,A12,1X,F7.2,1X,F12.8,1X,I1,1X,I1,1X,F8.5,1X,F8.5,
+ 1X,F8.5,1X,F9.6,1X,F9.6,1X,F9.6)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,131
read(1,'(A98)')ar__3
read(ar__3,4)
+ Subset_3,Adopted_2,Name_3,KOI_3,Per_2,Npl,Pos,b0_1,E_b0_2,
+ e_b0_3,dbdt_1,E_dbdt_2,e_dbdt_3
c ..............Just test output...........
write(6,4)
+ Subset_3,Adopted_2,Name_3,KOI_3,Per_2,Npl,Pos,b0_1,E_b0_2,
+ e_b0_3,dbdt_1,E_dbdt_2,e_dbdt_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Instantaneous coordinates, velocities, and orbital
* elements of the best-fitting parameters for the
* adopted solutions in this work and in paper II
C Format for file interpretation
5 format(
+ F19.6,1X,F20.7,1X,F21.8,1X,F20.12,1X,F20.12,1X,F20.13,1X,A20,
+ 1X,A20,1X,I3,1X,F18.5,1X,F19.17,1X,F18.16,1X,F20.18,1X,F19.17,
+ 1X,F18.16,1X,I1,1X,F7.2,1X,A12,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,241
read(1,'(A333)')ar__4
read(ar__4,5)
+ x,y,z,vx,vy,vz,Mass__1,Massp,T0_1,a,e,Pomega,i,Omega,Mean,
+ Pos_1,KOI_4,Name_4,sigma
c ..............Just test output...........
write(6,5)
+ x,y,z,vx,vy,vz,Mass__1,Massp,T0_1,a,e,Pomega,i,Omega,Mean,
+ Pos_1,KOI_4,Name_4,sigma
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end