Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/809/25 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-16
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/809/25 Stellar and planet properties for K2 candidates (Montet+, 2015)
*================================================================================
*Stellar and planetary properties of K2 campaign 1 candidates and validation
*of 17 planets, including a planet receiving earth-like insolation.
* Montet B.T., Morton T.D., Foreman-Mackey D., Johnson J.A., Hogg D.W.,
* Bowler B.P., Latham D.W., Bieryla A., Mann A.W.
* <Astrophys. J., 809, 25 (2015)>
* =2015ApJ...809...25M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! *Photometry for All Objects of Interest
integer*4 nr__
parameter (nr__=31) ! Number of records
character*130 ar__ ! Full-size record
integer*4 Star ! Star's name (EPIC NNNNNNNNN in Simbad)
real*4 Bmag ! (mag) ? APASS DR6 B-band magnitude (1)
real*4 e_Bmag ! (mag) ? Uncertainty in Bmag (1)
real*4 Vmag ! (mag) ? APASS DR6 V-band magnitude (1)
real*4 e_Vmag ! (mag) ? Uncertainty in Vmag (1)
real*4 gmag ! (mag) APASS DR6 g-band magnitude (1)
real*4 e_gmag ! (mag) Uncertainty in gmag (1)
real*4 rmag ! (mag) APASS DR6 r-band magnitude (1)
real*4 e_rmag ! (mag) Uncertainty in rmag (1)
real*4 imag ! (mag) APASS DR6 i-band magnitude (1)
real*4 e_imag ! (mag) Uncertainty in imag (1)
real*4 Jmag ! (mag) 2MASS J-band magnitude (2)
real*4 e_Jmag ! (mag) Uncertainty in Jmag (2)
real*4 Hmag ! (mag) 2MASS H-band magnitude (2)
real*4 e_Hmag ! (mag) Uncertainty in Hmag (2)
real*4 Kmag ! (mag) 2MASS K-band magnitude (2)
real*4 e_Kmag ! (mag) Uncertainty in Kmag (2)
real*4 W1mag ! (mag) AllWISE W1-band magnitude (3)
real*4 e_W1mag ! (mag) Uncertainty in W1mag (3)
real*4 W2mag ! (mag) AllWISE W2-band magnitude (3)
real*4 e_W2mag ! (mag) Uncertainty in W2mag (3)
real*4 W3mag ! (mag) ? AllWISE W3-band magnitude (3)
real*4 e_W3mag ! (mag) ? Uncertainty in W3mag (3)
*Note (1): Magnitude from the AAVSO Photometric All-sky Survey (APASS) DR6
* (Henden & Munari 2014CoSka..43..518H) as reported in the UCAC4 Catalog
* (Zacharias et al. 2012yCat.1322....0Z, Cat. I/322).
*Note (2): Magnitude from the 2MASS All-sky Catalog of Point Sources
* (Cutri et al. 2003, Cat. II/246).
*Note (3): Magnitude from the ALLWise Data Release
* (Cutri et al. 2013, Cat. II/328).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *Stellar Properties for All Objects of Interest
integer*4 nr__1
parameter (nr__1=31) ! Number of records
character*118 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Star_1 ! Star's name (EPIC NNNNNNNNN in Simbad)
character*5 Other ! Other name (K2-NN)
character*1 n_Other ! [a] Note on Other (G1)
real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000) (1)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000) (1)
real*4 Mass ! (Msun) Stellar mass
real*4 E_Mass ! (Msun) Upper limit uncertainty in Mass
real*4 e_Mass_1 ! (Msun) Lower limit uncertainty in Mass
real*4 Rad ! (Rsun) Stellar radius
real*4 E_Rad ! (Rsun) Upper limit uncertainty in Rad
real*4 e_Rad_1 ! (Rsun) Lower limit uncertainty in Rad
integer*4 Teff ! (K) Effective temperature
integer*4 E_Teff ! (K) Upper limit uncertainty in Teff
integer*4 e_Teff_1 ! (K) Lower limit uncertainty in Teff
real*4 v_Fe_H_ ! ([-]) Stellar metallicity (in dex)
real*4 E__Fe_H_ ! ([-]) Upper limit uncertainty in [Fe/H]
real*4 e__Fe_H__1 ! ([-]) Lower limit uncertainty in [Fe/H]
integer*4 Dist ! (pc) Distance
integer*4 E_Dist ! (pc) Upper limit uncertainty in Dist
integer*4 e_Dist_1 ! (pc) Lower limit uncertainty in Dist
*Note (1): The coordinates are retrieved directly from the EPIC (Huber et al.
* 2016, J/ApJS/224/2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! *Planet Properties for All Objects of Interest
integer*4 nr__2
parameter (nr__2=36) ! Number of records
character*112 ar__2 ! Full-size record
real*8 Planet ! Planet candidate name (EPIC NNNNNNNNN.NN in
* Simbad)
character*6 Other_1 ! Other name (K2-NNa)
character*1 n_Other_1 ! [a] Note on Other (G1)
real*8 Per ! (d) ? Period
real*8 e_Per ! (d) ? Uncertainty in Per
real*8 Epoch ! (d) ? Epoch (BJD-2456808)
real*4 e_Epoch ! (d) ? Uncertainty in Epoch
real*4 Rad_1 ! (Rgeo) ? Planet radius
real*4 e_Rad_2 ! (Rgeo) ? Uncertainty in Rad
real*4 a_R_ ! ? Semi-major axis to stellar radius ratio
real*4 e_a_R_ ! ? Uncertainty in a/R*
real*4 a ! (AU) ? Semi-major axis in AU
real*4 e_a ! (AU) ? Uncertainty in a
integer*4 Teq ! (K) ? Equilibrium temperature
integer*4 e_Teq ! (K) ? Uncertainty in Teq
character*9 Disp ! Disposition (Planet,Candidate,FP=False Positive)
character*1 n_Disp ! [b] Note on Disp (1)
*Note (1): Note as follows:
* b = Declared a false positive due to noise modeling systematics
* (see Section 5.2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Detected Companions to Candidate Host Stars
integer*4 nr__3
parameter (nr__3=31) ! Number of records
character*79 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Star_2 ! Primary star name (EPIC NNNNNNNNN in Simbad)
real*4 Aper ! (arcsec) Aperture used to create the K2 stellar light
* curve
real*8 RAdeg_1 ! (deg) ? Right Ascension in decimal degrees (J2000) (1)
real*8 DEdeg_1 ! (deg) ? Declination in decimal degrees (J2000) (1)
character*7 Detec ! Dataset used to detect the imaged companion
real*4 Sep ! (arcsec) ? Separation (2)
real*4 e_Sep ! (arcsec) ? Uncertainty in Sep
character*1 n_Sep ! [hi] Note on Sep (3)
real*4 drmag ! (mag) ? Difference in r-band magnitude between the
* primary K2 target star and the companion
real*4 e_drmag ! (mag) ? Uncertainty in drmag
character*1 n_drmag ! [j] Note on drmag (4)
real*4 Depth1 ! (10-3) ? Observed "transit" depth (in parts per
* thousand, ppt) (5)
real*4 Depth ! (10-3) Observed transit depth in the K2 dataset (in
* parts per thousand, ppt) (6)
*Note (1): Position of imaged companion.
*Note (2): Distance between the primary K2 target star and companion, in the
* dataset in which the companion is detected.
*Note (3): Note as follows:
* h = Separation from AO imaging;
* i = Separation from SDSS photometry.
*Note (4): Note as follows:
* j = {Delta}r inferred from JHK relative photometry.
*Note (5): If the imaged companion's flux were fully contained in the aperture
* and if it were an equal-mass eclipsing binary, leading to an eclipse
* depth of 50%. This is the maximum possible false positive eclipse
* depth, as described in Section 4.3.
*Note (6): If larger than the "max depth", this transit event cannot be caused by
* eclipses of the background star.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! False Positive Probability Calculation Results
integer*4 nr__4
parameter (nr__4=36) ! Number of records
character*75 ar__4 ! Full-size record
real*8 Planet_1 ! Planet candidate name (EPIC NNNNNNNNN.NN in
* Simbad)
character*6 Other_2 ! Other name (K2-NNa)
real*4 delta ! (10-3) Maximum depth of potential secondary eclipse
* signal (in parts per thousand, ppt)
character*1 AO ! [Y-] Adaptive optics observation presented in
* this paper ?
character*1 l_PrEB ! [<] Limit flag on PrEB
real*4 PrEB ! ? Relative probability for this false positive
* scenario : undiluted eclipsing binary
character*1 l_PrBEB ! [<] Limit flag on PrBEB
real*4 PrBEB ! ? Relative probability for this false positive
* scenario : chance-aligned background
* (/foreground) eclipsing binary
character*1 l_PrHEB ! [<] Limit flag on PrHEB
real*4 PrHEB ! ? Relative probability for this false positive
* scenario : hierarchical triple eclipsing binary
real*4 fp ! ? Integrated planet occurrence rate (1)
character*1 l_FPP ! [<] Limit flag on FPP
real*4 FPP ! ? False positive probability (2)
character*9 Disp_1 ! Disposition (Planet,Candidate,FP=False Positive)
character*1 n_Disp_1 ! [defg] Note on Disp (3)
*Note (1): Assumed between 0.7x and 1.3x the candidate radius.
*Note (2): Candidates are declared to be validated planets if FPP<0.01.
*Note (3): Note as follows:
* d = Despite low FPP, returned to candidate status out of abundance of caution
* due to secondary star detection within or near photometric aperture.
* e = Declared a false positive due to noise modeling systematics
* (see Section 5.2).
* f = Identified as planets by Crossfield et al. (2015ApJ...804...10C).
* g = Identified as planets by Armstrong et al. (2015, J/A+A/582/A33).
C=============================================================================
C Loading file 'table1.dat' ! *Photometry for All Objects of Interest
C Format for file interpretation
1 format(
+ I9,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.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,31
read(1,'(A130)')ar__
read(ar__,1)
+ Star,Bmag,e_Bmag,Vmag,e_Vmag,gmag,e_gmag,rmag,e_rmag,imag,
+ e_imag,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,W1mag,e_W1mag,
+ W2mag,e_W2mag,W3mag,e_W3mag
if(ar__(11:15) .EQ. '') Bmag = rNULL__
if(ar__(17:20) .EQ. '') e_Bmag = rNULL__
if(ar__(22:26) .EQ. '') Vmag = rNULL__
if(ar__(28:31) .EQ. '') e_Vmag = rNULL__
if(ar__(121:125) .EQ. '') W3mag = rNULL__
if(ar__(127:130) .EQ. '') e_W3mag = rNULL__
c ..............Just test output...........
write(6,1)
+ Star,Bmag,e_Bmag,Vmag,e_Vmag,gmag,e_gmag,rmag,e_rmag,imag,
+ e_imag,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,W1mag,e_W1mag,
+ W2mag,e_W2mag,W3mag,e_W3mag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! *Stellar Properties for All Objects of Interest
C Format for file interpretation
2 format(
+ I9,1X,A5,A1,1X,F10.6,1X,F9.6,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F5.3,1X,I4,1X,I3,1X,I3,1X,F5.2,1X,F4.2,1X,F4.2,1X,
+ I4,1X,I4,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,31
read(1,'(A118)')ar__1
read(ar__1,2)
+ Star_1,Other,n_Other,RAdeg,DEdeg,Mass,E_Mass,e_Mass_1,Rad,
+ E_Rad,e_Rad_1,Teff,E_Teff,e_Teff_1,v_Fe_H_,E__Fe_H_,
+ e__Fe_H__1,Dist,E_Dist,e_Dist_1
c ..............Just test output...........
write(6,2)
+ Star_1,Other,n_Other,RAdeg,DEdeg,Mass,E_Mass,e_Mass_1,Rad,
+ E_Rad,e_Rad_1,Teff,E_Teff,e_Teff_1,v_Fe_H_,E__Fe_H_,
+ e__Fe_H__1,Dist,E_Dist,e_Dist_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! *Planet Properties for All Objects of Interest
C Format for file interpretation
3 format(
+ F12.2,1X,A6,1X,A1,1X,F8.5,1X,F7.5,1X,F7.4,1X,F6.4,1X,F6.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F6.4,1X,F6.4,1X,I4,1X,I3,1X,A9,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,36
read(1,'(A112)')ar__2
read(ar__2,3)
+ Planet,Other_1,n_Other_1,Per,e_Per,Epoch,e_Epoch,Rad_1,
+ e_Rad_2,a_R_,e_a_R_,a,e_a,Teq,e_Teq,Disp,n_Disp
if(ar__2(23:30) .EQ. '') Per = rNULL__
if(ar__2(32:38) .EQ. '') e_Per = rNULL__
if(ar__2(40:46) .EQ. '') Epoch = rNULL__
if(ar__2(48:53) .EQ. '') e_Epoch = rNULL__
if(ar__2(55:60) .EQ. '') Rad_1 = rNULL__
if(ar__2(62:66) .EQ. '') e_Rad_2 = rNULL__
if(ar__2(68:72) .EQ. '') a_R_ = rNULL__
if(ar__2(74:78) .EQ. '') e_a_R_ = rNULL__
if(ar__2(80:85) .EQ. '') a = rNULL__
if(ar__2(87:92) .EQ. '') e_a = rNULL__
if(ar__2(94:97) .EQ. '') Teq = iNULL__
if(ar__2(99:101) .EQ. '') e_Teq = iNULL__
c ..............Just test output...........
write(6,3)
+ Planet,Other_1,n_Other_1,Per,e_Per,Epoch,e_Epoch,Rad_1,
+ e_Rad_2,a_R_,e_a_R_,a,e_a,Teq,e_Teq,Disp,n_Disp
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Detected Companions to Candidate Host Stars
C Format for file interpretation
4 format(
+ I9,1X,F4.1,1X,F10.6,1X,F9.6,1X,A7,1X,F5.2,1X,F4.2,A1,1X,F4.2,
+ 1X,F4.2,A1,1X,F5.1,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,31
read(1,'(A79)')ar__3
read(ar__3,4)
+ Star_2,Aper,RAdeg_1,DEdeg_1,Detec,Sep,e_Sep,n_Sep,drmag,
+ e_drmag,n_drmag,Depth1,Depth
if(ar__3(16:25) .EQ. '') RAdeg_1 = rNULL__
if(ar__3(27:35) .EQ. '') DEdeg_1 = rNULL__
if(ar__3(45:49) .EQ. '') Sep = rNULL__
if(ar__3(51:54) .EQ. '') e_Sep = rNULL__
if(ar__3(57:60) .EQ. '') drmag = rNULL__
if(ar__3(62:65) .EQ. '') e_drmag = rNULL__
if(ar__3(68:72) .EQ. '') Depth1 = rNULL__
c ..............Just test output...........
write(6,4)
+ Star_2,Aper,RAdeg_1,DEdeg_1,Detec,Sep,e_Sep,n_Sep,drmag,
+ e_drmag,n_drmag,Depth1,Depth
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! False Positive Probability Calculation Results
C Format for file interpretation
5 format(
+ F12.2,1X,A6,1X,F4.2,1X,A1,1X,A1,E6.2,1X,A1,E6.2,1X,A1,E6.2,1X,
+ F4.2,1X,A1,E6.2,1X,A9,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,36
read(1,'(A75)')ar__4
read(ar__4,5)
+ Planet_1,Other_2,delta,AO,l_PrEB,PrEB,l_PrBEB,PrBEB,l_PrHEB,
+ PrHEB,fp,l_FPP,FPP,Disp_1,n_Disp_1
if(ar__4(29:34) .EQ. '') PrEB = rNULL__
if(ar__4(37:42) .EQ. '') PrBEB = rNULL__
if(ar__4(45:50) .EQ. '') PrHEB = rNULL__
if(ar__4(52:55) .EQ. '') fp = rNULL__
if(ar__4(58:63) .EQ. '') FPP = rNULL__
c ..............Just test output...........
write(6,5)
+ Planet_1,Other_2,delta,AO,l_PrEB,PrEB,l_PrBEB,PrBEB,l_PrHEB,
+ PrHEB,fp,l_FPP,FPP,Disp_1,n_Disp_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end