Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/709/168 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.5, on 2013-May-22
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int 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 ! Star name
character*1 m_Name ! [b-f] Planet designation
real*8 msin_i ! Planet sin(i) mass
real*8 P ! (d) Planet orbital period
real*4 e ! Planet eccentricity
real*8 K ! (m/s) Planet radial velocity semi-amplitude
real*4 Ke ! (m/s) Semi-amplitude of first eccentric harmonic
real*4 Ke2 ! (m/s) Semi-amplitude of second order harmonic
real*8 mhsin_i ! (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 ! Star name
character*1 f_Name ! [c] Individual notes (1)
real*4 chi_c ! ({chi}^2^)^0.5^ of circular orbital solution
real*4 chi_r ! ({chi}^2^)^0.5^ of resonant orbital solution
real*4 chi_e ! ({chi}^2^)^0.5^ of eccentric orbital solution
real*4 Conf ! (%) Confidence level (2)
integer*4 Nobs ! Number of observations
character*1 l_Nreq ! [+] Limit flag on Nreq
integer*4 Nreq ! ? Estimated number of required observations to
* detect the second harmonic (3)
character*1 f_Nreq ! [ab] Flag on Nreq (1)
character*9 Pref ! Preferred solution (Circular, Eccentric or
* Resonant)
character*1 u_FAP ! [~] Uncertainty flag on FAP
real*4 FAP ! (%) ? False alarm probability (FAP) of the favored
* solution
character*3 Qual ! [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,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,252
read(1,'(A70)')ar__
read(ar__,1)Name,m_Name,msin_i,P,e,K,Ke,Ke2,mhsin_i
c ..............Just test output...........
write(6,1)Name,m_Name,msin_i,P,e,K,Ke,Ke2,mhsin_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,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,161
read(1,'(A75)')ar__1
read(ar__1,2)
+ Name_1,f_Name,chi_c,chi_r,chi_e,Conf,Nobs,l_Nreq,Nreq,f_Nreq,
+ Pref,u_FAP,FAP,Qual
if(ar__1(50:53) .EQ. '') Nreq = iNULL__
if(ar__1(67:71) .EQ. '') FAP = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,f_Name,chi_c,chi_r,chi_e,Conf,Nobs,l_Nreq,Nreq,f_Nreq,
+ Pref,u_FAP,FAP,Qual
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end