Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/876/23 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/876/23 Multiple populations of extrasolar gas giants (Goda+, 2019)
*================================================================================
*Multiple populations of extrasolar gas giants.
* Goda S., Matsuo T.
* <Astrophys. J., 876, 23 (2019)>
* =2019ApJ...876...23G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Planet parameters
integer*4 nr__
parameter (nr__=569) ! Number of records
character*112 ar__ ! Full-size record
character*18 ID (nr__) ! Planet identifier
real*4 MajAxis (nr__) ! (AU) [0.02/22] Semi-major axis (1)
real*4 e_MajAxis (nr__) ! (AU) [0/3.4] Uncertainty in MajAxis
character*1 l_Msini (nr__) ! [l] l=lower limit of companion mass was
* applied as that listed in Ref
real*4 Msini (nr__) ! (Mjup) [0.3/103.1] Planet mass times sine
* of inclination
real*4 e_Msini (nr__) ! (Mjup) [0/28.7] Uncertainty in Msini
real*8 K (nr__) ! (m/s) [3.1/6114]? Stellar velocity semi-amplitude
real*8 e_K (nr__) ! (m/s) [0/1361]? Uncertainty in K
real*8 P (nr__) ! (d) [1.3/39785] Orbital period
real*8 e_P (nr__) ! (d) [0/6205] Uncertainty in P
real*4 e (nr__) ! [0/0.97] Eccentricity
real*4 e_e (nr__) ! [0/0.3] Uncertainty in e
character*19 BibCode (nr__) ! Bibcode of reference
*Note (1): Calculated based on Equations 1 and 2:
*
* Equation (1): a~9.996x10^-1^P^2/3^M_*_^1/3^ (au)
* where P is the orbital period in years.
*
* Equation (2): M_p_sini~4.919x10^-3^P^1/3^(1-e^2^)^1/2^KM_*_^2/3^ (M_J_)
* where P is in days, e is the orbital eccentricity, and K is the velocity
* semi-amplitude of the star in m/s.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Parameters of stars hosting 569 original samples
and parameters of radial-velocity observations
used in this paper
integer*4 nr__1
parameter (nr__1=485) ! Number of records
character*73 ar__1 ! Full-size record
character*16 ID_1 (nr__1) ! Star identifier
character*1 f__Fe_H_ (nr__1) ! Flag on [Fe/H] (1)
real*4 v_Fe_H_ (nr__1) ! ([Sun]) [-0.8/0.5] Metallicity
real*4 e__Fe_H_ (nr__1) ! ([Sun]) [0.01/0.2] Uncertainty in [Fe/H]
character*1 f_Mass (nr__1) ! Flag on Mass (1)
real*4 Mass (nr__1) ! (Msun) [0.3/4.2] Stellar mass
real*4 e_Mass (nr__1) ! (Msun) [0.02/1.5] Uncertainty in Mass
real*4 sigma (nr__1) ! (m/s) [0.9/720]? Companion detection accuracy (2)
real*4 tau (nr__1) ! (yr) [0.2/33]? Companion detection duration (2)
character*19 BibCode_1 (nr__1) ! Bibcode reference
*Note (1): Flag as follows:
* a = The host-star metallicity was calibrated based on linear regression
* derived from the host-star metallicity correlations between the
* SWEET-Cat and Geneva-Copenhagen catalogs.
*Note (2): A blank indicates that the accuracy and duration of the radial
* velocity measurement for the sample were not used for constructing
* the detection probabilities for a companion as derived from the
* radial-velocity measurements.
C=============================================================================
C Loading file 'table1.dat' ! Planet parameters
C Format for file interpretation
1 format(
+ A18,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F5.2,1X,F10.5,1X,F7.2,1X,
+ F8.2,1X,F7.2,1X,F5.3,1X,F5.3,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,569
read(1,'(A112)')ar__
read(ar__,1)
+ ID(i__),MajAxis(i__),e_MajAxis(i__),l_Msini(i__),Msini(i__),
+ e_Msini(i__),K(i__),e_K(i__),P(i__),e_P(i__),e(i__),e_e(i__),
+ BibCode(i__)
if(ar__(46:55) .EQ. '') K(i__) = rNULL__
if(ar__(57:63) .EQ. '') e_K(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),MajAxis(i__),e_MajAxis(i__),l_Msini(i__),Msini(i__),
+ e_Msini(i__),K(i__),e_K(i__),P(i__),e_P(i__),e(i__),e_e(i__),
+ BibCode(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Parameters of stars hosting 569 original samples
* and parameters of radial-velocity observations
* used in this paper
C Format for file interpretation
2 format(
+ A16,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,F6.2,1X,
+ F4.1,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,485
read(1,'(A73)')ar__1
read(ar__1,2)
+ ID_1(i__),f__Fe_H_(i__),v_Fe_H_(i__),e__Fe_H_(i__),
+ f_Mass(i__),Mass(i__),e_Mass(i__),sigma(i__),tau(i__),
+ BibCode_1(i__)
if(ar__1(43:48) .EQ. '') sigma(i__) = rNULL__
if(ar__1(50:53) .EQ. '') tau(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),f__Fe_H_(i__),v_Fe_H_(i__),e__Fe_H_(i__),
+ f_Mass(i__),Mass(i__),e_Mass(i__),sigma(i__),tau(i__),
+ BibCode_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end