Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/153/208 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/153/208 LCES HIRES/Keck radial velocity Exoplanet Survey (Butler+, 2017)
*================================================================================
*The LCES HIRES/Keck precision radial velocity Exoplanet Survey.
* Butler R.P., Vogt S.S., Laughlin G., Burt J.A., Rivera E.J., Tuomi M.,
* Teske J., Arriagada P., Diaz M., Holden B., Keiser S.
* <Astron. J., 153, 208-208 (2017)>
* =2017AJ....153..208B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'targets.dat' ! Stars observed (corrected velocities in Tal-Or
et al, 2019, Cat. J/MNRAS/484/L8)
integer*4 nr__
parameter (nr__=1680) ! Number of records
character*14 ar__ ! Full-size record
character*14 Name (nr__) ! Target identifier
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Radial velocities and chromospheric activity
indices of all program stars
(corrected velocities in Tal-Or
et al, 2019, Cat. J/MNRAS/484/L8)
integer*4 nr__1
parameter (nr__1=62913) ! Number of records
character*74 ar__1 ! Full-size record
character*14 Name_1 (nr__1) ! Target identifier
real*8 BJD (nr__1) ! (d) Barycentric Julian Date
real*8 RV (nr__1) ! (m/s) [-39874.3/42363.6] Radial velocity
real*4 e_RV (nr__1) ! (m/s) [0.01/252.9] 1{sigma} uncertainty in RV
real*8 S (nr__1) ! [-5.5/68.9] Chromosperic activity index S-index (1)
real*8 H (nr__1) ! [-1/0.11]? Chromosperic activity index H-index (2)
integer*4 Count (nr__1) ! [67/616935] Median photons per pixel
integer*4 Exp (nr__1) ! (s) [-1/2551] Exposure time
*Note (1): Obtained from measurement of the emission reversal at the cores of the
* Fraunhofer H and K lines of Ca II at 3968{AA} and 3934{AA} respectively.
*Note (2): A measure of the chromospheric emission component at the H{alpha}
* Balmer line of hydrogen at 6563{AA} (Gomes da Silva et al.
* 2011A&A...534A..30G). Note that H-index values H>=0.062 are to be
* considered with caution (see Section 4 in the paper).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *Signals in the Keck data sets
integer*4 nr__2
parameter (nr__2=357) ! Number of records
character*73 ar__2 ! Full-size record
integer*4 Num (nr__2) ! [4/1672] Target number
character*10 Name_2 (nr__2) ! Target identifier
integer*4 Ns (nr__2) ! [1/4] Number of signals (n_signals_)
real*4 LR (nr__2) ! ([-]) [16.5/673.3] Likelihood ratio ({Delta}lnL_k_)
real*8 Per (nr__2) ! (d) [1/6264.7] Period
real*8 e_Per (nr__2) ! (d) [0/573.77] Uncertainty in Per
real*4 K (nr__2) ! (m/s) [1.25/840.71] Radial velocity amplitude
real*4 e_K (nr__2) ! (m/s) [0.17/165.1] Uncertainty in K
character*9 Signal (nr__2) ! Interpretation (either Activity, Candidate,
* Published, or SRC=Signal Requiring
* Confirmation) (1)
character*3 Note (nr__2) ! Additional note (1-7) (2)
*Note (1): The signals are classified as follows:
* Published = Signal that have already been published as planet claim;
* Candidate = Significant unpublished planet candidate that await
* photometric follow-up to rule out activity-related cause;
* SRC = Signal Requiring Confirmation by additional data before rising
* to classification as planet candidate;
* Activity = Keplerian-like signal that we classify as likely arising from
* stellar activity.
*Note (2): Additional notes are defined as follows:
* 1 = Signal coincides with a photometric rotation period based on ASAS
* V-band photometry;
* 2 = Corresponding signal in the S-index data at a very nearby period likely
* indicative of the rotation period of the star;
* 3 = Signal not constrained from above in the parameter space as the
* analyses were limited to 1 day<P<T_obs_;
* 4 = Coincides with the stellar magnetic cycle apparent in the S-index but
* strong enough as a signal to be interpreted as a candidate planet;
* 5 = Likely artifact signal caused by orbital evolution of the system;
* 6 = Daily alias of the known transiting planet orbiting the star with an
* orbital period of 0.73 days (Winn et al. 2011ApJ...737L..18W);
* 7 = The strongest S-index counterpart in our sample indicates this signal
* is likely caused by activity unless one argues that S-index variability
* cannot have counterparts in radial velocities. Signal reported as a
* planet candidate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Maximum likelihood parameter estimates and
standard errors for the "nuisance" parameters
integer*4 nr__3
parameter (nr__3=1010) ! Number of records
character*84 ar__3 ! Full-size record
integer*4 Num_1 (nr__3) ! [3/1673] Target number
character*14 Name_3 (nr__3) ! Target identifier
real*8 dV_dt (nr__3) ! (m/s/yr) [-900.31/622.51] Linear acceleration {gamma}
real*8 e_dV_dt (nr__3) ! (m/s/yr) [0.035/332.111] Standard error in dV/dt
real*4 sigmaJ (nr__3) ! (m/s) [0.32/98.7] Excess white noise or "jitter"
* ({sigma}_J_)
real*4 e_sigmaJ (nr__3) ! (m/s) [0.07/20.31] Standard error in sigmaJ
real*8 cS (nr__3) ! (m/s) [-4866.9/8437.7] Coefficient quantifying
* linear dependence of the radial velocities
* on the S-index (c_S_)
real*8 e_cS (nr__3) ! (m/s) [0.18/4622.1] Standard error in c_S_
real*8 LRS (nr__3) ! ([-]) [7.883/159.22]?=0 Likelihood ratio statistic
* {Delta}lnL_S_ (1)
real*8 Per_1 (nr__3) ! (d) [2.065/6308.8]?=0 Period of the signal in the
* S-index time-series (P_S_)
*Note (1): {Delta}lnL_S_=lnL_S_(k=1)-lnL_S_(k=0) represents the likelihood ratio
* statistic of model with a signal (k=1) in the S-index with respect to a
* model without a signal (k=0). Only signals exceeding a ratio of 7.82
* corresponding to 5% False Alarm Probability (FAP) are shown.
C=============================================================================
C Loading file 'targets.dat' ! Stars observed (corrected velocities in Tal-Or
* et al, 2019, Cat. J/MNRAS/484/L8)
C Format for file interpretation
1 format(A14)
C Effective file loading
open(unit=1,status='old',file=
+'targets.dat')
write(6,*) '....Loading file: targets.dat'
do i__=1,1680
read(1,'(A14)')ar__
read(ar__,1)Name(i__)
c ..............Just test output...........
write(6,1)Name(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! Radial velocities and chromospheric activity
* indices of all program stars
* (corrected velocities in Tal-Or
* et al, 2019, Cat. J/MNRAS/484/L8)
C Format for file interpretation
2 format(A14,1X,F13.5,1X,F9.2,1X,F6.2,1X,F7.4,1X,F8.5,1X,I6,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,62913
read(1,'(A74)')ar__1
read(ar__1,2)
+ Name_1(i__),BJD(i__),RV(i__),e_RV(i__),S(i__),H(i__),
+ Count(i__),Exp(i__)
if(ar__1(55:62) .EQ. '') H(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),BJD(i__),RV(i__),e_RV(i__),S(i__),H(i__),
+ Count(i__),Exp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! *Signals in the Keck data sets
C Format for file interpretation
3 format(
+ I4,1X,A10,1X,I1,1X,F6.2,1X,F10.5,1X,F9.5,1X,F6.2,1X,F6.2,1X,
+ A9,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,357
read(1,'(A73)')ar__2
read(ar__2,3)
+ Num(i__),Name_2(i__),Ns(i__),LR(i__),Per(i__),e_Per(i__),
+ K(i__),e_K(i__),Signal(i__),Note(i__)
c ..............Just test output...........
write(6,3)
+ Num(i__),Name_2(i__),Ns(i__),LR(i__),Per(i__),e_Per(i__),
+ K(i__),e_K(i__),Signal(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Maximum likelihood parameter estimates and
* standard errors for the "nuisance" parameters
C Format for file interpretation
4 format(
+ I4,1X,A14,1X,F8.3,1X,F7.3,1X,F5.2,1X,F5.2,1X,F9.3,1X,F8.3,1X,
+ F7.3,1X,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,1010
read(1,'(A84)')ar__3
read(ar__3,4)
+ Num_1(i__),Name_3(i__),dV_dt(i__),e_dV_dt(i__),sigmaJ(i__),
+ e_sigmaJ(i__),cS(i__),e_cS(i__),LRS(i__),Per_1(i__)
c ..............Just test output...........
write(6,4)
+ Num_1(i__),Name_3(i__),dV_dt(i__),e_dV_dt(i__),sigmaJ(i__),
+ e_sigmaJ(i__),cS(i__),e_cS(i__),LRS(i__),Per_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end