Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/167/103 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 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 (nr__) ! Source, A=Paper II; B=this work (1)
integer*4 Adopted (nr__) ! [0/1] 1=Adopted solution
integer*4 ID (nr__) ! [1/5] Run identifier
character*13 Name (nr__) ! Kepler Planet Name
real*8 KOI (nr__) ! Kepler object of interest
real*8 Per (nr__) ! (d) [0.7/500] Mean orbital period
real*8 Mp_M_ (nr__) ! (10-6) [0.73/1429] Planet-to-star mass ratio
real*4 E_Mp_M_ (nr__) ! (10-6) [0.1/58] Upper uncertainty on Mp/M*
real*4 e_Mp_M__1 (nr__) ! (10-6) [-67/0] Lower uncertainty on Mp/M*
real*8 Mp (nr__) ! (Mgeo) [0.2/597] Absolute planet mass (2)
real*8 E_Mp (nr__) ! (Mgeo) [0.04/37] Upper uncertainty on Mp
real*8 e_Mp_1 (nr__) ! (Mgeo) [-45/-0.04] Lower uncertainty on Mp
real*8 rhop (nr__) ! (g.cm-3) [0.04/49] Absolute planet density (3)
real*8 E_rhop (nr__) ! (g.cm-3) [0.004/20] Upper uncertainty on rhop
real*8 e_rhop_1 (nr__) ! (g.cm-3) [-45/-0.004] Lower uncertainty on rhop
real*4 Radp (nr__) ! (Rgeo) [0.96/12.1] Absolute planet radius (4)
real*4 E_Radp (nr__) ! (Rgeo) [0.01/2] Upper uncertainty on Radp
real*4 e_Radp_1 (nr__) ! (Rgeo) [-0.6/-0.02] Lower uncertainty on Radp
real*8 dex (nr__) ! [-0.3/0.26] Eccentricity difference to
* preceding planet in x (5)
real*8 E_dex (nr__) ! [0.0005/0.2] Upper uncertainty on dex
real*8 e_dex_1 (nr__) ! [-0.2/-0.0006] Lower uncertainty on dex
real*8 dey (nr__) ! [-0.47/0.58] Eccentricity difference to
* preceding planet in y (6)
real*8 E_dey (nr__) ! [0.0006/0.2] Upper uncertainty on dey
real*8 e_dey_1 (nr__) ! [-0.5/-0.0006] Lower uncertainty on dey
real*4 Ix (nr__) ! (deg) [-31.4/34.8] Inclination around the line of
* sight ("roll angle") (7)
real*4 E_Ix (nr__) ! (deg) [0/23] Upper uncertainty on Ix
real*4 e_Ix_1 (nr__) ! (deg) [-21/0] Lower uncertainty on Ix
real*8 Iy (nr__) ! (deg) [-4.39/7.1] Inclination relative to the line
* of sight (8)
real*8 E_Iy (nr__) ! (deg) [0.004/7] Upper uncertainty on Iy
real*8 e_Iy_1 (nr__) ! (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 (nr__1) ! [AB] Source, A=Paper II; B=this work (1)
integer*4 Adopted_1 (nr__1) ! [0/1] 1=Adopted solution
integer*4 ID_1 (nr__1) ! Run identifier
character*13 Name_1 (nr__1) ! Kepler Planet Name
real*8 KOI_1 (nr__1) ! Kepler object of interest
real*8 Per_1 (nr__1) ! (d) [0.7/1072] Mean orbital period
character*10 TTVStd (nr__1) ! (min) Best model TTV std
integer*4 J (nr__1) ! [0/5] J of the nearest first-order resonance
* J:J-1 (2)
real*8 Delta (nr__1) ! [-1/0.48] Normalized distance from the J:J-1
* resonance (3)
real*8 SPer (nr__1) ! (d) [0/6612] Super-period of the J:J-1 resonance (4)
real*8 T0 (nr__1) ! (h) [0/18] Transit duration at first data time stamp
real*8 E_T0 (nr__1) ! (h) [0.002/1.13] Upper uncertainty on T0
real*8 e_T0_1 (nr__1) ! (h) [-2/0] Lower uncertainty on T0
real*8 Tau0 (nr__1) ! (min) [0.44/161] ingress-egress time at first data
* time stamp
real*8 E_Tau0 (nr__1) ! (min) [0/25] Upper uncertainty on Tau0
real*8 e_Tau0_1 (nr__1) ! (min) [-21/0] Lower uncertainty on Tau0
real*8 b0 (nr__1) ! [-1.01/1.01] impact parameter at BKJD=0
real*8 E_b0 (nr__1) ! [0/2] Upper uncertainty on b0
real*8 e_b0_1 (nr__1) ! [-2/0] Lower uncertainty on b0
real*8 dbdt (nr__1) ! (yr-1) [-0.82/0.82] db/dt median
real*8 E_dbdt (nr__1) ! (yr-1) [0/0.1]Upper uncertainty on db/dt
real*8 e_dbdt_1 (nr__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 (nr__2) ! Source, A=Paper II; B=this work (1)
integer*4 KOI_2 (nr__2) ! Kepler object of interest
character*10 Name_2 (nr__2) ! Kepler Star Name
integer*4 KIC (nr__2) ! Kepler Input Catalog number
character*19 Ref (nr__2) ! Reference (ADS bibcode or NExScI) for Mstar, Rstar
real*4 u1 (nr__2) ! [0.23/0.67] first limb-darkening coeff, taken
* from NExScI
real*4 u2 (nr__2) ! [0.07/0.45] second limb-darkening coeff, taken
* from NExScI
real*4 Mass_ (nr__2) ! (Msun) [0.5/1.5] Absolute stellar mass, from Ref
real*4 E_Mass_ (nr__2) ! (Msun) [0.01/0.2] Upper uncertainty on Mass*
real*4 e_Mass__1 (nr__2) ! (Msun) [-0.2/-0.012] Lower uncertainty on Mass*
real*4 Rad_ (nr__2) ! (Rsun) [0.51/1.84] Absolute stellar radius from Ref
real*8 E_Rad_ (nr__2) ! (Rsun) [0.006/0.5] Upper uncertainty on Rad*
real*8 e_Rad__1 (nr__2) ! (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 (nr__3) ! Source, A=Paper II; B=this work (1)
integer*4 Adopted_2 (nr__3) ! [0/1] 1=Adopted solution
character*12 Name_3 (nr__3) ! Kepler Planet Name
real*8 KOI_3 (nr__3) ! Kepler object of interest
real*8 Per_2 (nr__3) ! (d) [0.7/500] mean orbital period
integer*4 Npl (nr__3) ! [2/6] model number of planets
integer*4 Pos (nr__3) ! [1/6] planet position within the system, from
* inside to outside
real*8 b0_1 (nr__3) ! [-0.94/0.92] impact parameter at BKJD=0
real*8 E_b0_2 (nr__3) ! [0.0004/0.2] Upper uncertainty on b0
real*8 e_b0_3 (nr__3) ! [-0.3/-0.0004] Lower uncertainty on b0
real*8 dbdt_1 (nr__3) ! (yr-1) [-0.77/0.09] db/dt median
real*8 E_dbdt_2 (nr__3) ! (yr-1) [3e-05/0.05] Upper uncertainty on db/dt
real*8 e_dbdt_3 (nr__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 (nr__4) ! (m) Instantaneous Cartesian x coordinate (1)
real*8 y (nr__4) ! (m) Instantaneous Cartesian y coordinate (1)
real*8 z (nr__4) ! (m) Instantaneous Cartesian z coordinate (1)
real*8 vx (nr__4) ! (m.s-1) Instantaneous Cartesian x velocity
real*8 vy (nr__4) ! (m.s-1) Instantaneous Cartesian y velocity
real*8 vz (nr__4) ! (m.s-1) Instantaneous Cartesian z velocity
character*20 Mass__1 (nr__4) ! (kg) Star mass
character*20 Massp (nr__4) ! (kg) Planet mass
integer*4 T0_1 (nr__4) ! (d) [120/353] Epoch, days since Barycentric Kepler
* Julian date
real*8 a (nr__4) ! (m) Orbital semi-major axis
real*8 e (nr__4) ! Eccentricity
real*8 Pomega (nr__4) ! (rad) Longitude of periastron
real*8 i (nr__4) ! (rad) Orbital inclination
real*8 Omega (nr__4) ! (rad) Longitude of the ascending node
real*8 Mean (nr__4) ! (rad) Mean anomaly
integer*4 Pos_1 (nr__4) ! [1/6] Planet position within the system,
* from inside to outside
real*8 KOI_4 (nr__4) ! Kepler object of interest
character*12 Name_4 (nr__4) ! Kepler Planet Name
character*20 sigma (nr__4) ! 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(i__),Adopted(i__),ID(i__),Name(i__),KOI(i__),Per(i__),
+ Mp_M_(i__),E_Mp_M_(i__),e_Mp_M__1(i__),Mp(i__),E_Mp(i__),
+ e_Mp_1(i__),rhop(i__),E_rhop(i__),e_rhop_1(i__),Radp(i__),
+ E_Radp(i__),e_Radp_1(i__),dex(i__),E_dex(i__),e_dex_1(i__),
+ dey(i__),E_dey(i__),e_dey_1(i__),Ix(i__),E_Ix(i__),
+ e_Ix_1(i__),Iy(i__),E_Iy(i__),e_Iy_1(i__)
c ..............Just test output...........
write(6,1)
+ Subset(i__),Adopted(i__),ID(i__),Name(i__),KOI(i__),Per(i__),
+ Mp_M_(i__),E_Mp_M_(i__),e_Mp_M__1(i__),Mp(i__),E_Mp(i__),
+ e_Mp_1(i__),rhop(i__),E_rhop(i__),e_rhop_1(i__),Radp(i__),
+ E_Radp(i__),e_Radp_1(i__),dex(i__),E_dex(i__),e_dex_1(i__),
+ dey(i__),E_dey(i__),e_dey_1(i__),Ix(i__),E_Ix(i__),
+ e_Ix_1(i__),Iy(i__),E_Iy(i__),e_Iy_1(i__)
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(i__),Adopted_1(i__),ID_1(i__),Name_1(i__),KOI_1(i__),
+ Per_1(i__),TTVStd(i__),J(i__),Delta(i__),SPer(i__),T0(i__),
+ E_T0(i__),e_T0_1(i__),Tau0(i__),E_Tau0(i__),e_Tau0_1(i__),
+ b0(i__),E_b0(i__),e_b0_1(i__),dbdt(i__),E_dbdt(i__),
+ e_dbdt_1(i__)
c ..............Just test output...........
write(6,2)
+ Subset_1(i__),Adopted_1(i__),ID_1(i__),Name_1(i__),KOI_1(i__),
+ Per_1(i__),TTVStd(i__),J(i__),Delta(i__),SPer(i__),T0(i__),
+ E_T0(i__),e_T0_1(i__),Tau0(i__),E_Tau0(i__),e_Tau0_1(i__),
+ b0(i__),E_b0(i__),e_b0_1(i__),dbdt(i__),E_dbdt(i__),
+ e_dbdt_1(i__)
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(i__),KOI_2(i__),Name_2(i__),KIC(i__),Ref(i__),
+ u1(i__),u2(i__),Mass_(i__),E_Mass_(i__),e_Mass__1(i__),
+ Rad_(i__),E_Rad_(i__),e_Rad__1(i__)
c ..............Just test output...........
write(6,3)
+ Subset_2(i__),KOI_2(i__),Name_2(i__),KIC(i__),Ref(i__),
+ u1(i__),u2(i__),Mass_(i__),E_Mass_(i__),e_Mass__1(i__),
+ Rad_(i__),E_Rad_(i__),e_Rad__1(i__)
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(i__),Adopted_2(i__),Name_3(i__),KOI_3(i__),
+ Per_2(i__),Npl(i__),Pos(i__),b0_1(i__),E_b0_2(i__),
+ e_b0_3(i__),dbdt_1(i__),E_dbdt_2(i__),e_dbdt_3(i__)
c ..............Just test output...........
write(6,4)
+ Subset_3(i__),Adopted_2(i__),Name_3(i__),KOI_3(i__),
+ Per_2(i__),Npl(i__),Pos(i__),b0_1(i__),E_b0_2(i__),
+ e_b0_3(i__),dbdt_1(i__),E_dbdt_2(i__),e_dbdt_3(i__)
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(i__),y(i__),z(i__),vx(i__),vy(i__),vz(i__),Mass__1(i__),
+ Massp(i__),T0_1(i__),a(i__),e(i__),Pomega(i__),i(i__),
+ Omega(i__),Mean(i__),Pos_1(i__),KOI_4(i__),Name_4(i__),
+ sigma(i__)
c ..............Just test output...........
write(6,5)
+ x(i__),y(i__),z(i__),vx(i__),vy(i__),vz(i__),Mass__1(i__),
+ Massp(i__),T0_1(i__),a(i__),e(i__),Pomega(i__),i(i__),
+ Omega(i__),Mean(i__),Pos_1(i__),KOI_4(i__),Name_4(i__),
+ sigma(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end