Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/154/108 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-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/108 California-Kepler Survey (CKS). II. Properties (Johnson+, 2017)
*================================================================================
*The California-Kepler Survey. II. Precise physical properties of 2025 Kepler
*planets and their host stars.
* Johnson J.A., Petigura E.A., Fulton B.J., Marcy G.W., Howard A.W.,
* Isaacson H., Hebb L., Cargile P.A., Morton T.D., Weiss L.M., Winn J.N.,
* Rogers L.A., Sinukoff E., Hirsch L.A.
* <Astron. J., 154, 108-108 (2017)>
* =2017AJ....154..108J (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Stellar properties
integer*4 nr__
parameter (nr__=1305) ! Number of records
character*156 ar__ ! Full-size record
character*12 Name (nr__) ! Star name: mostly Kepler Object Identifier
* (KOI), or Kepler Input Calalog (KIC) number
character*11 TYC2 (nr__) ! Tycho-2 catalog identifier
real*4 Ksmag (nr__) ! [6.5/14.7] Apparent 2MASS Ks band magnitude
integer*4 Teff (nr__) ! (K) [4619/6651] Isochrones derived effective
* temperature (T_eff,iso_) (1)
integer*4 e_Teff (nr__) ! (K) [49/192] Lower uncertainty in Teff
integer*4 E_Teff_1 (nr__) ! (K) [48/182] Upper uncertainty in Teff
real*4 logg (nr__) ! ([cm/s2]) [2.7/4.7] Isochrones derived log surface
* gravity (logg_iso_) (1)
real*4 e_logg (nr__) ! ([cm/s2]) [0.02/0.33] Lower uncertainty in logg
real*4 E_logg_1 (nr__) ! ([cm/s2]) [0.02/0.15] Upper uncertainty in logg
real*4 v_Fe_H_ (nr__) ! ([Sun]) [-0.6/0.5] Isochrones derived metallicity
* ([Fe/H]_iso_) (1)
real*4 e__Fe_H_ (nr__) ! ([Sun]) [0.01/0.1] Lower uncertainty in [Fe/H]
real*4 E__Fe_H__1 (nr__) ! ([Sun]) [0.01/0.19] Upper uncertainty in [Fe/H]
real*4 Mass (nr__) ! (Msun) [0.6/2.2] Isochrones derived stellar mass
* (M_*_) (1)
real*4 e_Mass (nr__) ! (Msun) [0.02/0.23] Lower uncertainty in Mass
real*4 E_Mass_1 (nr__) ! (Msun) [0.02/0.2] Upper uncertainty in Mass
real*4 Rad (nr__) ! (Rsun) [0.6/10.4] Isochrones derived stellar
* radius (R_*_) (1)
real*4 e_Rad (nr__) ! (Rsun) [0.02/1.6] Lower uncertainty in Rad
real*4 E_Rad_1 (nr__) ! (Rsun) [0.02/4.6] Upper uncertainty in Rad
real*4 logA (nr__) ! ([yr]) [8.9/10.2] Isochrones derived log age
* (log_10_(age)) (1)
real*4 e_logA (nr__) ! ([yr]) [0.05/0.56] Lower uncertainty in logA
real*4 E_logA_1 (nr__) ! ([yr]) [0.04/0.35] Upper uncertainty in logA
real*4 plx1 (nr__) ! (mas) [0.1/35] Isochrones derived parallax (1)
real*4 e_plx1 (nr__) ! (mas) [0.01/1.4] Lower uncertainty in plx1
real*4 E_plx1_1 (nr__) ! (mas) [0.01/1.6] Upper uncertainty in plx1
real*4 plx2 (nr__) ! (mas) [1/15.3]? Tycho-Gaia Astrometric Solution
* trigonometric parallax
real*4 e_plx2 (nr__) ! (mas) [0.2/1]? Lower uncertainty in plx2
real*4 E_plx2_1 (nr__) ! (mas) [0.2/1]? Upper uncertainty in plx2
*Note (1):
* We used the isochrones Python package to derive the following physical
* parameters:
* Teff, logg, [Fe/H], Mass, Rad, logA, and plx1.
* Isochrones returns posterior distributions on effective temperature,
* surface gravity, and metallicity, which are distinguished from the purely
* spectroscopic measurements.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! California-Kepler Survey (CKS) planet parameters
integer*4 nr__1
parameter (nr__1=2025) ! Number of records
character*106 ar__1 ! Full-size record
real*8 KOI (nr__1) ! Kepler Object Identifier (KOI) of the planet
* candidate
real*8 Per (nr__1) ! (d) Orbital period
real*8 Rp_R_ (nr__1) ! ? Planetary to stellar radii ratio (1)
real*8 e_Rp_R_ (nr__1) ! ? Lower uncertainty in Rp/R*
real*8 E_Rp_R__1 (nr__1) ! ? Upper uncertainty in Rp/R*
real*8 Rad_1 (nr__1) ! (Rgeo) ? Planetary radius R_p_
real*8 e_Rad_2 (nr__1) ! (Rgeo) ? Lower uncertainty in Rad
real*8 E_Rad_3 (nr__1) ! (Rgeo) ? Upper uncertainty in Rad
integer*4 Sinc (nr__1) ! (Earth) ? Stellar irradiance (incident stellar flux)
* S_inc_ (2)
integer*4 e_Sinc (nr__1) ! (Earth) ? Lower uncertainty in Sinc
integer*4 E_Sinc_1 (nr__1) ! (Earth) ? Upper uncertainty in Sinc
integer*4 Teq (nr__1) ! (K) ? Equilibrium temperature T_eq_ (3)
integer*4 e_Teq (nr__1) ! (K) ? Lower uncertainty in Teq
integer*4 E_Teq_1 (nr__1) ! (K) ? Upper uncertainty in Teq
*Note (1): From the NASA's Exoplanet Archive Q1-Q16 KOI catalog (Mullally et al.
* 2015, Cat. J/ApJS/217/31).
*Note (2):
* Stellar irradiance received at the planet relative to the Earth.
* We compute the incident flux as:
* S_inc_/S_{Earth}_=(T_eff_/5778K)^4^(R_*_/R_{Sun}_)^2^(a/au)^-2^, Eq. (3).
*Note (3):
* Defined according to
* (T_eq_/280K)=(S_inc_/S_{Earth}_)^1/4^(1-{alpha}/4)^1/4^, Eq. (4).
* Assuming a Bond albedo {alpha} of 0.3, typical for super-Earth-size planets
* (Demory 2014ApJ...789L..20D).
C=============================================================================
C Loading file 'table1.dat' ! Stellar properties
C Format for file interpretation
1 format(
+ A12,5X,A11,1X,F4.1,1X,I4,1X,I3,1X,I3,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,1305
read(1,'(A156)')ar__
read(ar__,1)
+ Name(i__),TYC2(i__),Ksmag(i__),Teff(i__),e_Teff(i__),
+ E_Teff_1(i__),logg(i__),e_logg(i__),E_logg_1(i__),
+ v_Fe_H_(i__),e__Fe_H_(i__),E__Fe_H__1(i__),Mass(i__),
+ e_Mass(i__),E_Mass_1(i__),Rad(i__),e_Rad(i__),E_Rad_1(i__),
+ logA(i__),e_logA(i__),E_logA_1(i__),plx1(i__),e_plx1(i__),
+ E_plx1_1(i__),plx2(i__),e_plx2(i__),E_plx2_1(i__)
if(ar__(142:146) .EQ. '') plx2(i__) = rNULL__
if(ar__(148:151) .EQ. '') e_plx2(i__) = rNULL__
if(ar__(153:156) .EQ. '') E_plx2_1(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),TYC2(i__),Ksmag(i__),Teff(i__),e_Teff(i__),
+ E_Teff_1(i__),logg(i__),e_logg(i__),E_logg_1(i__),
+ v_Fe_H_(i__),e__Fe_H_(i__),E__Fe_H__1(i__),Mass(i__),
+ e_Mass(i__),E_Mass_1(i__),Rad(i__),e_Rad(i__),E_Rad_1(i__),
+ logA(i__),e_logA(i__),E_logA_1(i__),plx1(i__),e_plx1(i__),
+ E_plx1_1(i__),plx2(i__),e_plx2(i__),E_plx2_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! California-Kepler Survey (CKS) planet parameters
C Format for file interpretation
2 format(
+ 4X,F8.2,1X,F7.2,1X,F9.6,1X,F9.6,1X,F9.6,1X,F8.2,1X,F8.2,1X,
+ F8.2,1X,I5,1X,I4,1X,I4,1X,I4,1X,I3,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,2025
read(1,'(A106)')ar__1
read(ar__1,2)
+ KOI(i__),Per(i__),Rp_R_(i__),e_Rp_R_(i__),E_Rp_R__1(i__),
+ Rad_1(i__),e_Rad_2(i__),E_Rad_3(i__),Sinc(i__),e_Sinc(i__),
+ E_Sinc_1(i__),Teq(i__),e_Teq(i__),E_Teq_1(i__)
if(ar__1(22:30) .EQ. '') Rp_R_(i__) = rNULL__
if(ar__1(32:40) .EQ. '') e_Rp_R_(i__) = rNULL__
if(ar__1(42:50) .EQ. '') E_Rp_R__1(i__) = rNULL__
if(ar__1(52:59) .EQ. '') Rad_1(i__) = rNULL__
if(ar__1(61:68) .EQ. '') e_Rad_2(i__) = rNULL__
if(ar__1(70:77) .EQ. '') E_Rad_3(i__) = rNULL__
if(ar__1(79:83) .EQ. '') Sinc(i__) = iNULL__
if(ar__1(85:88) .EQ. '') e_Sinc(i__) = iNULL__
if(ar__1(90:93) .EQ. '') E_Sinc_1(i__) = iNULL__
if(ar__1(95:98) .EQ. '') Teq(i__) = iNULL__
if(ar__1(100:102) .EQ. '') e_Teq(i__) = iNULL__
if(ar__1(104:106) .EQ. '') E_Teq_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ KOI(i__),Per(i__),Rp_R_(i__),e_Rp_R_(i__),E_Rp_R__1(i__),
+ Rad_1(i__),e_Rad_2(i__),E_Rad_3(i__),Sinc(i__),e_Sinc(i__),
+ E_Sinc_1(i__),Teq(i__),e_Teq(i__),E_Teq_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end