Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/654/625 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.81 (2015-09-23), on 2026-Jun-08
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/654/625 Radial velocities of 9 long-period objects (Wittenmyer+, 2007)
*================================================================================
*Long-period objects in the extrasolar planetary systems 47 Ursae Majoris
*and 14 Herculis.
* Wittenmyer R.A., Endl M., Cochran W.D.
* <Astrophys. J., 654, 625-632 (2007)>
* =2007ApJ...654..625W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Stellar parameters
integer*4 nr__
parameter (nr__=9) ! Number of records
character*50 ar__ ! Full-size record
integer*4 HD ! HD number of star
real*4 Mass ! (solMass) Stellar mass
real*4 v_Fe_H_ ! ([Sun]) metallicity
real*4 e__Fe_H_ ! ([Sun]) Mean error on [Fe/H]
integer*4 Teff ! (K) Effective temperature of star
integer*4 e_Teff ! (K) Mean error on Teff
real*4 logR_HK ! Activity index (ratio of Ca H+K bands
* relative to the total stellar flux)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Planet parameters
integer*4 nr__1
parameter (nr__1=14) ! Number of records
character*104 ar__1 ! Full-size record
integer*4 HD_1 ! HD number of parent star
character*1 m_HD ! [bcd] Planet designation
real*8 P ! (d) Oribtal period
real*8 e_P ! (d) Mean error on P
real*8 T0_2400000 ! (d) Periastron epoch (JD-2400000)
real*4 e_T0_2400000 ! (d) Mean error on T0-2400000
real*4 e ! Eccentricity of orbit (1)
real*4 e_e ! ? Mean error on e (1)
real*4 omega ! (deg) Periastron angle (1)
real*4 e_omega ! (deg) ? Mean error on omega (1)
real*4 K ! (m/s) Velocity semi-amplitude
real*4 e_K ! (m/s) Mean error on K
real*4 Msini ! (jovMass) Projected mass, in units of Jupiter mass
real*4 e_Msini ! (jovMass) Mean error on Msini
real*4 a ! (AU) Semi-major axis of orbit
real*4 e_a ! (AU) Mean error on a
*Note (1): value of parameter is fixed when the error is blank.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'veloc.dat' ! Radial Velocities (tables 5 to 14 of paper)
integer*4 nr__2
parameter (nr__2=450) ! Number of records
character*39 ar__2 ! Full-size record
integer*4 HD_2 ! HD number
character*3 Tel ! Telescope
real*8 JD ! Julian Date of observation
real*4 RV ! (m/s) Radial velocity
real*4 e_RV ! (m/s) Uncertainty in RV
C=============================================================================
C Loading file 'table1.dat' ! Stellar parameters
C Format for file interpretation
1 format(3X,I6,10X,F4.2,2X,F4.2,1X,F4.2,1X,I5,1X,I2,2X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,9
read(1,'(A50)')ar__
read(ar__,1)HD,Mass,v_Fe_H_,e__Fe_H_,Teff,e_Teff,logR_HK
c ..............Just test output...........
write(6,1)HD,Mass,v_Fe_H_,e__Fe_H_,Teff,e_Teff,logR_HK
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Planet parameters
C Format for file interpretation
2 format(
+ 3X,I6,1X,A1,1X,F10.5,1X,F9.5,1X,F8.2,1X,F6.2,1X,F5.3,1X,F5.3,
+ 1X,F5.1,1X,F4.1,1X,F5.1,1X,F3.1,1X,F5.2,1X,F4.2,1X,F5.3,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,14
read(1,'(A104)')ar__1
read(ar__1,2)
+ HD_1,m_HD,P,e_P,T0_2400000,e_T0_2400000,e,e_e,omega,e_omega,K,
+ e_K,Msini,e_Msini,a,e_a
if(ar__1(56:60) .EQ. '') e_e = rNULL__
if(ar__1(68:71) .EQ. '') e_omega = rNULL__
c ..............Just test output...........
write(6,2)
+ HD_1,m_HD,P,e_P,T0_2400000,e_T0_2400000,e,e_e,omega,e_omega,K,
+ e_K,Msini,e_Msini,a,e_a
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'veloc.dat' ! Radial Velocities (tables 5 to 14 of paper)
C Format for file interpretation
3 format(3X,I6,1X,A3,1X,F13.5,1X,F6.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'veloc.dat')
write(6,*) '....Loading file: veloc.dat'
do i__=1,450
read(1,'(A39)')ar__2
read(ar__2,3)HD_2,Tel,JD,RV,e_RV
c ..............Just test output...........
write(6,3)HD_2,Tel,JD,RV,e_RV
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end