Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/131/2274 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 2024-Apr-16
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/131/2274 TrES-1 equivalent widths (Sozzetti+, 2006)
*================================================================================
*Chemical composition of the planet-harboring star TrES-1.
* Sozzetti A., Yong D., Carney B.W., Laird J.B., Latham D.W., Torres G.
* <Astron. J., 131, 2274-2289 (2006)>
* =2006AJ....131.2274S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! TrES-1 equivalent widths
integer*4 nr__
parameter (nr__=103) ! Number of records
character*42 ar__ ! Full-size record
real*8 lambda (nr__) ! (0.1nm) Wavelength in Angstroms
character*4 Ion (nr__) ! Species name
integer*4 Z (nr__) ! Atomic number
real*4 logEps (nr__) ! ([Sun]) Log of the nominal solar abundance
real*4 chi (nr__) ! (eV) Lower excitation potential
real*4 log_gf (nr__) ! ([-]) Log of the oscillator strength
real*4 EW (nr__) ! (10-13m) ? Equivalent width in milliAngstroms
character*1 f_EW (nr__) ! [a] EW from spectrum synthesis (1)
integer*4 r_log_gf (nr__) ! Reference for log(gf)
*Note (1): A "a" indicates that the abundance was determined from spectrum
* synthesis as no reliable EW measurement was available.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__1
parameter (nr__1=16) ! Number of records
character*98 ar__1 ! Full-size record
integer*4 Ref (nr__1) ! Reference number
character*19 Bibcode (nr__1) ! Bibcode
character*17 Aut (nr__1) ! Author's name
character*57 Com (nr__1) ! Comment
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Kinematic data of stars with planets and TrES-1
integer*4 nr__2
parameter (nr__2=121) ! Number of records
character*132 ar__2 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__2) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__2) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*9 Name (nr__2) ! Source name (1)
real*4 plx (nr__2) ! (mas) Parallax
real*8 HRV (nr__2) ! (km/s) Heliocentric Radial Velocity
real*8 pmRA (nr__2) ! (mas/yr) Proper motion in Right Ascension
real*8 pmDE (nr__2) ! (mas/yr) Proper motion in Declination
integer*4 RAh (nr__2) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm (nr__2) ! (min) Minute of Right Ascension (J2000)
real*4 RAs (nr__2) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__2) ! Sign of the Declination (J2000)
integer*4 DEd (nr__2) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__2) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__2) ! (arcsec) Arcsecond of Declination (J2000)
real*4 Uvel (nr__2) ! (km/s) Velocity U towards the Galactic anticenter
real*4 Vvel (nr__2) ! (km/s) Velocity V in the direction of galactic rotation
real*4 Wvel (nr__2) ! (km/s) Velocity W toward the North Galactic Pole
real*4 Pthin (nr__2) ! Mishenina et al. (2004, Cat. <J/A+A/418/551>)
* thin disk membership probability
real*4 Pthick (nr__2) ! Mishenina et al. (2004, Cat. <J/A+A/418/551>)
* thick disk membership probability
real*4 Ptd_d (nr__2) ! Bensby et al. (2003, Cat. <J/A+A/410/527>)
* thick-to-thin disk probability ratio
real*4 P1thin (nr__2) ! Venn et al. (2004, Cat. <J/AJ/128/1177>)
* thin disk membership probability (2)
real*4 P1thick (nr__2) ! Venn et al. (2004, Cat. <J/AJ/128/1177>)
* thick disk membership probability (2)
real*4 P2thin (nr__2) ! Venn et al. (2004, Cat. <J/AJ/128/1177>)
* thin disk membership probability (3)
real*4 P2thick (nr__2) ! Venn et al. (2004, Cat. <J/AJ/128/1177>)
* thick disk membership probability (3)
real*4 X (nr__2) ! ? Schuster et al. (1993, Cat. <J/A+AS/97/951>)
* stellar population criterion (X>=0 for halo
* stars)
real*4 a_Fe (nr__2) ! ([Sun]) ? Log of the {alpha}/Fe abundance ratio (4)
*Note (1): All values of positions, proper motions, and parallax are taken
* from Hipparcos (ESA 1997, Cat. <I/239>) and the SIMBAD database,
* except for TrES-1, for which RA, DE, pmRA and pmDE were taken from the
* USNO-B1.0 catalog (Monet et al. 2003, Cat. <I/284>), while the
* photometric parallax estimate was taken from Sozzetti et al.
* (2004MNRAS.354.1194S). The values of heliocentric radial-velocity are
* taken from SIMBAD and from the corresponding planet discovery papers.
*Note (2): Assuming uniform probability distributions.
*Note (3): Assuming Bensby's prior probability distributions.
*Note (4): From Bodaghee et al. (2003A&A...404..715B) and Beirao et al.
* (2005, Cat. <J/A+A/438/251>).
C=============================================================================
C Loading file 'table1.dat' ! TrES-1 equivalent widths
C Format for file interpretation
1 format(
+ F7.2,1X,A4,1X,I2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.1,1X,A1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,103
read(1,'(A42)')ar__
read(ar__,1)
+ lambda(i__),Ion(i__),Z(i__),logEps(i__),chi(i__),log_gf(i__),
+ EW(i__),f_EW(i__),r_log_gf(i__)
if(ar__(33:37) .EQ. '') EW(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ lambda(i__),Ion(i__),Z(i__),logEps(i__),chi(i__),log_gf(i__),
+ EW(i__),f_EW(i__),r_log_gf(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
2 format(I2,1X,A19,1X,A17,1X,A57)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,16
read(1,'(A98)')ar__1
read(ar__1,2)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,2)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Kinematic data of stars with planets and TrES-1
C Format for file interpretation
3 format(
+ A9,1X,F6.2,1X,F7.3,1X,F7.2,1X,F7.2,1X,I2,1X,I2,1X,F5.2,1X,A1,
+ I2,1X,I2,1X,F4.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F4.2,1X,F4.2,1X,
+ F6.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,121
read(1,'(A132)')ar__2
read(ar__2,3)
+ Name(i__),plx(i__),HRV(i__),pmRA(i__),pmDE(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),Pthin(i__),Pthick(i__),
+ Ptd_d(i__),P1thin(i__),P1thick(i__),P2thin(i__),P2thick(i__),
+ X(i__),a_Fe(i__)
if(ar__2(121:126) .EQ. '') X(i__) = rNULL__
if(ar__2(128:132) .EQ. '') a_Fe(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,3)
+ Name(i__),plx(i__),HRV(i__),pmRA(i__),pmDE(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ Uvel(i__),Vvel(i__),Wvel(i__),Pthin(i__),Pthick(i__),
+ Ptd_d(i__),P1thin(i__),P1thick(i__),P2thin(i__),P2thick(i__),
+ X(i__),a_Fe(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end