Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/709/168 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-Aug-15
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/709/168 Eccentric orbits in exoplanets (Anglada-Escude+, 2010)
*================================================================================
*How eccentric orbital solutions can hide planetary systems in 2:1 resonant
*orbits.
* Anglada-Escude G., Lopez-Morales M., Chambers J.E.
* <Astrophys. J., 709, 168-178 (2010)>
* =2010ApJ...709..168A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! List of eccentric planets and the amplitudes of
their correspondent eccentric harmonics
integer*4 nr__
parameter (nr__=252) ! Number of records
character*70 ar__ ! Full-size record
character*13 Name (nr__) ! Star name
character*1 m_Name (nr__) ! [b-f] Planet designation
real*8 msin_i (nr__) ! Planet sin(i) mass
real*8 P (nr__) ! (d) Planet orbital period
real*4 e (nr__) ! Planet eccentricity
real*8 K (nr__) ! (m/s) Planet radial velocity semi-amplitude
real*4 Ke (nr__) ! (m/s) Semi-amplitude of first eccentric harmonic
real*4 Ke2 (nr__) ! (m/s) Semi-amplitude of second order harmonic
real*8 mhsin_i (nr__) ! (Mgeo) Mass of the possible inner companion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Statistical comparison of circular, resonant,
and eccentric orbital solutions
integer*4 nr__1
parameter (nr__1=161) ! Number of records
character*75 ar__1 ! Full-size record
character*13 Name_1 (nr__1) ! Star name
character*1 f_Name (nr__1) ! [c] Individual notes (1)
real*4 chi_c (nr__1) ! ({chi}^2^)^0.5^ of circular orbital solution
real*4 chi_r (nr__1) ! ({chi}^2^)^0.5^ of resonant orbital solution
real*4 chi_e (nr__1) ! ({chi}^2^)^0.5^ of eccentric orbital solution
real*4 Conf (nr__1) ! (%) Confidence level (2)
integer*4 Nobs (nr__1) ! Number of observations
character*1 l_Nreq (nr__1) ! [+] Limit flag on Nreq
integer*4 Nreq (nr__1) ! ? Estimated number of required observations to
* detect the second harmonic (3)
character*1 f_Nreq (nr__1) ! [ab] Flag on Nreq (1)
character*9 Pref (nr__1) ! Preferred solution (Circular, Eccentric or
* Resonant)
character*1 u_FAP (nr__1) ! [~] Uncertainty flag on FAP
real*4 FAP (nr__1) ! (%) ? False alarm probability (FAP) of the favored
* solution
character*3 Qual (nr__1) ! [U* ] Significance of the solution (4)
*Note (1): Flag as follows:
* a = Other data sets and astrometry seem to confirm the eccentricity. Very
* noisy star.
* b = The obtained solution differs from the published one significantly.
* c = According to (Bean et al. 2007AJ....134..749B), astrometric observations
* indicate that the candidate is a star indeed with an orbital inclination
* close to 0.
*Note (2): To decide which of the orbital solutions is the best, we apply the
* confidence level test given by Lucy & Sweeney (1971AJ.....76..544L)
* and Lucy (2005A&A...439..663L). We only accept one of the non-circular
* solutions if the confidence level (c.l.) is better than 95%.
*Note (3): N_req_ is only given if the orbital solution is significantly
* non-circular.
*Note (4): Quality flag according to false alarm probability (FAP) as follows:
* *** = secure solution (FAP<0.1%)
* ** = FAP<1%
* * = FAP<5%
* U = undecided (the solutions are not significantly different: FAP>5%)
C=============================================================================
C Loading file 'table2.dat' ! List of eccentric planets and the amplitudes of
* their correspondent eccentric harmonics
C Format for file interpretation
1 format(
+ A13,1X,A1,5X,F7.2,1X,F7.2,1X,F4.2,1X,F7.2,1X,F6.2,1X,F6.2,1X,
+ F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,252
read(1,'(A70)')ar__
read(ar__,1)
+ Name(i__),m_Name(i__),msin_i(i__),P(i__),e(i__),K(i__),
+ Ke(i__),Ke2(i__),mhsin_i(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),m_Name(i__),msin_i(i__),P(i__),e(i__),K(i__),
+ Ke(i__),Ke2(i__),mhsin_i(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Statistical comparison of circular, resonant,
* and eccentric orbital solutions
C Format for file interpretation
2 format(
+ A13,4X,A1,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,A1,I4,A1,
+ 1X,A9,1X,A1,F5.2,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,161
read(1,'(A75)')ar__1
read(ar__1,2)
+ Name_1(i__),f_Name(i__),chi_c(i__),chi_r(i__),chi_e(i__),
+ Conf(i__),Nobs(i__),l_Nreq(i__),Nreq(i__),f_Nreq(i__),
+ Pref(i__),u_FAP(i__),FAP(i__),Qual(i__)
if(ar__1(50:53) .EQ. '') Nreq(i__) = iNULL__
if(ar__1(67:71) .EQ. '') FAP(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),f_Name(i__),chi_c(i__),chi_r(i__),chi_e(i__),
+ Conf(i__),Nobs(i__),l_Nreq(i__),Nreq(i__),f_Nreq(i__),
+ Pref(i__),u_FAP(i__),FAP(i__),Qual(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end