Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/158/38 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-May-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/158/38 New substellar discoveries from Kepler and K2 (Carmichael+, 2019)
*================================================================================
*New substellar discoveries from Kepler and K2: is there a brown dwarf desert?
* Carmichael T.W., Latham D.W., Vanderburg A.M.
* <Astron. J., 158, 38 (2019)>
* =2019AJ....158...38C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! List of published transiting brown dwarfs
as of March 2019
integer*4 nr__
parameter (nr__=23) ! Number of records
character*136 ar__ ! Full-size record
character*14 Name (nr__) ! System identifier
character*1 m_Name (nr__) ! [abc] Multiplicity index on Name
character*1 Flag1 (nr__) ! [abcdefg] Flag on Name (G1)
real*8 Per (nr__) ! (d) [1.218/166.788] Orbital period
real*4 Mass (nr__) ! (Mjup) [11.79/95.1] Brown dwarf mass
real*4 e_Mass (nr__) ! (Mjup) [0.59/6.8] Uncertainty in Mass
real*4 R (nr__) ! (Rjup) [0.757/6.5] Brown dwarf radius
real*4 e_R (nr__) ! (Rjup) [0.015/0.53] Uncertainty in R
character*1 l_e (nr__) ! Limit flag on e
real*4 e (nr__) ! [0/0.689] Eccentricity
character*1 f_e (nr__) ! [a] Flag indicating an adopted value on e
real*4 e_e (nr__) ! [0.0002/0.057]? Uncertainty in e
real*4 Mp_ (nr__) ! (Msun) [0.188/1.78]? Primary stellar mass
real*4 Ms_ (nr__) ! (Msun) [0.3]? Secondary stellar mass
real*4 e_Mp_ (nr__) ! (Msun) [0.015/0.12]? Uncertainty in Mp*
real*4 Rp_ (nr__) ! (Rsun) [0.21/1.881]? Primary stellar radius
real*4 e_Rp_ (nr__) ! (Rsun) [0.005/0.31]? Uncertainty in Rp*
integer*4 Teff (nr__) ! (K) [3100/7930]? Stellar effective temperature
integer*4 e_Teff (nr__) ! (K) [40/820]? Uncertainty in Teff
real*4 v_Fe_H_ (nr__) ! ([Sun]) [-0.24/0.44]? Metallicity
real*4 e__Fe_H_ (nr__) ! ([Sun]) [0.053/0.25]? Uncertainty in [Fe/H]
character*19 Bibcode1 (nr__) ! Bibcode of reference
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Additional information on published transiting
brown dwarfs
integer*4 nr__1
parameter (nr__1=22) ! Number of records
character*97 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*14 Name_1 (nr__1) ! System identifier
character*1 Flag2 (nr__1) ! [bcd] Flag on Name (G1)
integer*4 RAh (nr__1) ! (h) ? Hour of Right Ascension (J2000)
integer*4 RAm (nr__1) ! (min) ? Minute of Right Ascension (J2000)
real*4 RAs (nr__1) ! (s) ? Second of Right Ascension (J2000)
character*1 DE_ (nr__1) ! Sign of the Declination (J2000)
integer*4 DEd (nr__1) ! (deg) ? Degree of Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) ? Arcminute of Declination (J2000)
real*4 DEs (nr__1) ! (arcsec) ? Arcsecond of Declination (J2000)
real*4 Vmag (nr__1) ! (mag) [9.8/18.73]? V band magnitude
real*4 Gmag (nr__1) ! (mag) [14.35/18.94]? G-band magnitude from the Gaia
* mission (Cat. I/345)
character*23 Ref (nr__1) ! Reference
character*19 Bibcode2 (nr__1) ! Bibcode of the reference
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Multi-order relative radial velocities
of CWW 89A from TRES
integer*4 nr__2
parameter (nr__2=18) ! Number of records
character*27 ar__2 ! Full-size record
real*8 BJD (nr__2) ! (d) [6799.92/8430.98] Barycentric Julian Date
* (BJD-2450000)
real*8 RV (nr__2) ! (m/s) [-8260.62/11493.66] Radial velocity
real*4 e_RV (nr__2) ! (m/s) [35.97/389.48] Uncertainty in RV
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Multi-order relative radial velocity
measurements of EPIC 212036875 from TRES
integer*4 nr__3
parameter (nr__3=14) ! Number of records
character*27 ar__3 ! Full-size record
real*8 BJD_1 (nr__3) ! (d) [6799.92/8430.98] Barycentric Julian Date
* (BJD-2450000)
real*8 RV_1 (nr__3) ! (m/s) [-8260.62/11493.66] Radial velocity
real*4 e_RV_1 (nr__3) ! (m/s) [35.97/389.48] Uncertainty in RV
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Multi-order relative radial velocity
measurements of KOI-607 from TRES
integer*4 nr__4
parameter (nr__4=14) ! Number of records
character*27 ar__4 ! Full-size record
real*8 BJD_2 (nr__4) ! (d) [6799.92/8430.98] Barycentric Julian Date
* (BJD-2450000)
real*8 RV_2 (nr__4) ! (m/s) [-8260.62/11493.66] Radial velocity
real*4 e_RV_2 (nr__4) ! (m/s) [35.97/389.48] Uncertainty in RV
C=============================================================================
C Loading file 'table7.dat' ! List of published transiting brown dwarfs
* as of March 2019
C Format for file interpretation
1 format(
+ A14,A1,1X,A1,1X,F7.3,1X,F5.2,1X,F4.2,1X,F5.3,1X,F5.3,1X,A1,
+ F6.4,1X,A1,1X,F6.4,1X,F5.3,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ I4,1X,I3,1X,F6.3,1X,F5.3,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,23
read(1,'(A136)')ar__
read(ar__,1)
+ Name(i__),m_Name(i__),Flag1(i__),Per(i__),Mass(i__),
+ e_Mass(i__),R(i__),e_R(i__),l_e(i__),e(i__),f_e(i__),e_e(i__),
+ Mp_(i__),Ms_(i__),e_Mp_(i__),Rp_(i__),e_Rp_(i__),Teff(i__),
+ e_Teff(i__),v_Fe_H_(i__),e__Fe_H_(i__),Bibcode1(i__)
if(ar__(60:65) .EQ. '') e_e(i__) = rNULL__
if(ar__(67:71) .EQ. '') Mp_(i__) = rNULL__
if(ar__(73:76) .EQ. '') Ms_(i__) = rNULL__
if(ar__(78:82) .EQ. '') e_Mp_(i__) = rNULL__
if(ar__(84:88) .EQ. '') Rp_(i__) = rNULL__
if(ar__(90:94) .EQ. '') e_Rp_(i__) = rNULL__
if(ar__(96:99) .EQ. '') Teff(i__) = iNULL__
if(ar__(101:103) .EQ. '') e_Teff(i__) = iNULL__
if(ar__(105:110) .EQ. '') v_Fe_H_(i__) = rNULL__
if(ar__(112:116) .EQ. '') e__Fe_H_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),m_Name(i__),Flag1(i__),Per(i__),Mass(i__),
+ e_Mass(i__),R(i__),e_R(i__),l_e(i__),e(i__),f_e(i__),e_e(i__),
+ Mp_(i__),Ms_(i__),e_Mp_(i__),Rp_(i__),e_Rp_(i__),Teff(i__),
+ e_Teff(i__),v_Fe_H_(i__),e__Fe_H_(i__),Bibcode1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Additional information on published transiting
* brown dwarfs
C Format for file interpretation
2 format(
+ A14,1X,A1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,A23,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,22
read(1,'(A97)')ar__1
read(ar__1,2)
+ Name_1(i__),Flag2(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),Gmag(i__),Ref(i__),
+ Bibcode2(i__)
if(ar__1(18:19) .EQ. '') RAh(i__) = iNULL__
if(ar__1(21:22) .EQ. '') RAm(i__) = iNULL__
if(ar__1(24:28) .EQ. '') RAs(i__) = rNULL__
if(ar__1(31:32) .EQ. '') DEd(i__) = iNULL__
if(ar__1(34:35) .EQ. '') DEm(i__) = iNULL__
if(ar__1(37:41) .EQ. '') DEs(i__) = rNULL__
if(ar__1(43:47) .EQ. '') Vmag(i__) = rNULL__
if(ar__1(49:53) .EQ. '') Gmag(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Flag2(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),Gmag(i__),Ref(i__),
+ Bibcode2(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Multi-order relative radial velocities
* of CWW 89A from TRES
C Format for file interpretation
3 format(F11.6,1X,F8.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,18
read(1,'(A27)')ar__2
read(ar__2,3)BJD(i__),RV(i__),e_RV(i__)
c ..............Just test output...........
write(6,3)BJD(i__),RV(i__),e_RV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Multi-order relative radial velocity
* measurements of EPIC 212036875 from TRES
C Format for file interpretation
4 format(F11.6,1X,F8.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,14
read(1,'(A27)')ar__3
read(ar__3,4)BJD_1(i__),RV_1(i__),e_RV_1(i__)
c ..............Just test output...........
write(6,4)BJD_1(i__),RV_1(i__),e_RV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Multi-order relative radial velocity
* measurements of KOI-607 from TRES
C Format for file interpretation
5 format(F11.6,1X,F8.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,14
read(1,'(A27)')ar__4
read(ar__4,5)BJD_2(i__),RV_2(i__),e_RV_2(i__)
c ..............Just test output...........
write(6,5)BJD_2(i__),RV_2(i__),e_RV_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end