Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/497/563 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-18
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/A+A/497/563 Chemical abundances of 451 stars (Neves+, 2009)
*================================================================================
*Chemical abundances of 451 stars from the HARPS GTO planet search program.
*Thin disc, thick disc, and planets.
* Neves V., Santos N.C., Sousa S.G., Correia A.C.M., Israelian G.
* <Astron. Astrophys., 497, 563-581 (2009)>
* =2009A&A...497..563N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Parameters to assign the stars to a Galactic
population
integer*4 nr__
parameter (nr__=451) ! Number of records
character*93 ar__ ! Full-size record
character*10 Star (nr__) ! HD or HIP name of the star
real*8 Ulsr (nr__) ! (km/s) U velocity of the star relative to the LSR (1)
real*8 Vlsr (nr__) ! (km/s) V velocity of the star relative to the LSR (1)
real*8 Wlsr (nr__) ! (km/s) W velocity of the star relative to the LSR (1)
real*4 fD (nr__) ! Distribution of the thin disc population (2)
real*4 fTD (nr__) ! Distribution of the thick disc population (2)
real*4 fH (nr__) ! Distribution of the halo population (2)
real*8 TDD (nr__) ! Probability of a star belonging to the thick
* disc relative to the thin disk
real*8 TDH (nr__) ! Probability of a star belonging to the thick
* disc relative to the halo
character*10 pop (nr__) ! probable population where the star belongs
*Note (1): The Galactic space velocities were calculated using the procedure
* from Johnson & Soderblom (1987AJ.....93..864J) and corrected for the
* solar motion relative to the Local Standard of Rest (LSR) using
* (U',V',W')=(+10.00,+5.25,+7.17)km/s from Dehnen & Binney
* (1998MNRAS.298..387D).
*Note (2): The selection of the different populations was completed using
* the method proposed by Bensby (2003, Cat. J/A+A/410/527, 2005, Cat.
* J/A+A/433/185).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Atomic parameters for lines selected in the paper
integer*4 nr__1
parameter (nr__1=180) ! Number of records
character*32 ar__1 ! Full-size record
character*5 El (nr__1) ! Element (1)
real*8 lambda (nr__1) ! (0.1nm) Central wavelength
real*4 EP (nr__1) ! Excitation potential
real*4 loggf (nr__1) ! Oscilator strength
real*4 EWsun (nr__1) ! (0.1pm) Solar equivalent width
*Note (1): Abundances for selected elements as follows:
* Na I : log(eps0)=6.33
* Mg I : log(eps0)=7.58
* Al I : log(eps0)=6.47
* Si I : log(eps0)=7.55
* Ca I : log(eps0)=6.36
* Sc I : log(eps0)=3.10
* Sc II : log(eps0)=3.10
* Ti I : log(eps0)=4.99
* Ti II : log(eps0)=4.99
* Mn I : log(eps0)=5.39
* Cr I : log(eps0)=5.67
* Cr II : log(eps0)=5.67
* V I : log(eps0)=4.00
* Co I : log(eps0)=4.92
* Ni I : log(eps0)=6.25
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Abundances, rms, and number of lines for each
star and element
integer*4 nr__2
parameter (nr__2=451) ! Number of records
character*225 ar__2 ! Full-size record
character*10 Star_1 (nr__2) ! HD or HIP name of the star
real*4 v_Si_H_ (nr__2) ! ([Sun]) Abundance [Si/H] (Z=14)
real*4 e__Si_H_ (nr__2) ! ([Sun]) rms uncertainty of [Si/H]
integer*4 o__Si_H_ (nr__2) ! Number of SiI lines used
real*4 v_Ca_H_ (nr__2) ! ([Sun]) Abundance [Ca/H] (Z=20)
real*4 e__Ca_H_ (nr__2) ! ([Sun]) rms uncertainty of [Ca/H]
integer*4 o__Ca_H_ (nr__2) ! Number of CaI lines used
real*4 v_ScI_H_ (nr__2) ! ([Sun]) Abundance [ScI/H] (Z=21.0)
real*4 e__ScI_H_ (nr__2) ! ([Sun]) rms uncertainty of [ScI/H]
integer*4 o__ScI_H_ (nr__2) ! Number of ScI lines used
real*4 v_ScII_H_ (nr__2) ! ([Sun]) Abundance [ScII/H] (Z=21.1)
real*4 e__ScII_H_ (nr__2) ! ([Sun]) rms uncertainty of [ScII/H]
integer*4 o__ScII_H_ (nr__2) ! Number of ScII lines used
real*4 v_TiI_H_ (nr__2) ! ([Sun]) Abundance [TiI/H] (Z=22.0)
real*4 e__TiI_H_ (nr__2) ! ([Sun]) rms uncertainty of [TiI/H]
integer*4 o__TiI_H_ (nr__2) ! Number of TiI lines used
real*4 v_TiII_H_ (nr__2) ! ([Sun]) Abundance [TiII/H] (Z=22.1)
real*4 e__TiII_H_ (nr__2) ! ([Sun]) rms uncertainty of [TiII/H]
integer*4 o__TiII_H_ (nr__2) ! Number of TiII lines used
real*4 v_V_H_ (nr__2) ! ([Sun]) Abundance [V/H] (Z=23)
real*4 e__V_H_ (nr__2) ! ([Sun]) rms uncertainty of [V/H]
integer*4 o__V_H_ (nr__2) ! Number of VI lines used
real*4 v_CrI_H_ (nr__2) ! ([Sun]) Abundance [CrI/H] (Z=24.0)
real*4 e__CrI_H_ (nr__2) ! ([Sun]) rms uncertainty of [CrI/H]
integer*4 o__CrI_H_ (nr__2) ! Number of CrI lines used
real*4 v_CrII_H_ (nr__2) ! ([Sun]) Abundance [CrII/H] (Z=24.1)
real*4 e__CrII_H_ (nr__2) ! ([Sun]) rms uncertainty of [CrII/H]
integer*4 o__CrII_H_ (nr__2) ! Number of CrII lines used
real*4 v_Mn_H_ (nr__2) ! ([Sun]) Abundance [Mn/H] (Z=25)
real*4 e__Mn_H_ (nr__2) ! ([Sun]) rms uncertainty of [Mn/H]
integer*4 o__Mn_H_ (nr__2) ! Number of MnI lines used
real*4 v_Co_H_ (nr__2) ! ([Sun]) Abundance [Co/H] (Z=27)
real*4 e__Co_H_ (nr__2) ! ([Sun]) rms uncertainty of [Co/H]
integer*4 o__Co_H_ (nr__2) ! Number of CoI lines used
real*4 v_Ni_H_ (nr__2) ! ([Sun]) Abundance [Ni/H] (Z=28)
real*4 e__Ni_H_ (nr__2) ! ([Sun]) rms uncertainty of [Ni/H]
integer*4 o__Ni_H_ (nr__2) ! Number of NiI lines used
real*4 v_Na_H_ (nr__2) ! ([Sun]) Abundance [Na/H] (Z=11)
real*4 e__Na_H_ (nr__2) ! ([Sun]) rms uncertainty of [Na/H]
integer*4 o__Na_H_ (nr__2) ! Number of Na lines used
real*4 v_Mg_H_ (nr__2) ! ([Sun]) Abundance [Mg/H] (Z=12)
real*4 e__Mg_H_ (nr__2) ! ([Sun]) rms uncertainty of [Mg/H]
integer*4 o__Mg_H_ (nr__2) ! Number of MgI lines used
real*4 v_Al_H_ (nr__2) ! ([Sun]) Abundance [Al/H] (Z=13)
real*4 e__Al_H_ (nr__2) ! ([Sun]) rms uncertainty of [Al/H]
integer*4 o__Al_H_ (nr__2) ! Number of AlI lines used
C=============================================================================
C Loading file 'table1.dat' ! Parameters to assign the stars to a Galactic
* population
C Format for file interpretation
1 format(
+ A10,1X,F8.3,1X,F8.3,1X,F8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,F8.3,
+ 1X,F8.3,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,451
read(1,'(A93)')ar__
read(ar__,1)
+ Star(i__),Ulsr(i__),Vlsr(i__),Wlsr(i__),fD(i__),fTD(i__),
+ fH(i__),TDD(i__),TDH(i__),pop(i__)
c ..............Just test output...........
write(6,1)
+ Star(i__),Ulsr(i__),Vlsr(i__),Wlsr(i__),fD(i__),fTD(i__),
+ fH(i__),TDD(i__),TDH(i__),pop(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Atomic parameters for lines selected in the paper
C Format for file interpretation
2 format(A5,2X,F7.2,1X,F4.2,1X,F6.3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,180
read(1,'(A32)')ar__1
read(ar__1,2)El(i__),lambda(i__),EP(i__),loggf(i__),EWsun(i__)
c ..............Just test output...........
write(6,2)El(i__),lambda(i__),EP(i__),loggf(i__),EWsun(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Abundances, rms, and number of lines for each
* star and element
C Format for file interpretation
3 format(
+ A10,1X,F6.3,1X,F4.2,1X,I2,1X,F6.3,1X,F4.2,1X,I2,1X,F6.3,1X,
+ F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I2,1X,
+ F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,
+ I2,1X,F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,
+ F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I2,1X,F6.3,1X,F4.2,1X,I1,1X,
+ F6.3,1X,F4.2,1X,I1,1X,F6.3,1X,F4.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,451
read(1,'(A225)')ar__2
read(ar__2,3)
+ Star_1(i__),v_Si_H_(i__),e__Si_H_(i__),o__Si_H_(i__),
+ v_Ca_H_(i__),e__Ca_H_(i__),o__Ca_H_(i__),v_ScI_H_(i__),
+ e__ScI_H_(i__),o__ScI_H_(i__),v_ScII_H_(i__),e__ScII_H_(i__),
+ o__ScII_H_(i__),v_TiI_H_(i__),e__TiI_H_(i__),o__TiI_H_(i__),
+ v_TiII_H_(i__),e__TiII_H_(i__),o__TiII_H_(i__),v_V_H_(i__),
+ e__V_H_(i__),o__V_H_(i__),v_CrI_H_(i__),e__CrI_H_(i__),
+ o__CrI_H_(i__),v_CrII_H_(i__),e__CrII_H_(i__),o__CrII_H_(i__),
+ v_Mn_H_(i__),e__Mn_H_(i__),o__Mn_H_(i__),v_Co_H_(i__),
+ e__Co_H_(i__),o__Co_H_(i__),v_Ni_H_(i__),e__Ni_H_(i__),
+ o__Ni_H_(i__),v_Na_H_(i__),e__Na_H_(i__),o__Na_H_(i__),
+ v_Mg_H_(i__),e__Mg_H_(i__),o__Mg_H_(i__),v_Al_H_(i__),
+ e__Al_H_(i__),o__Al_H_(i__)
c ..............Just test output...........
write(6,3)
+ Star_1(i__),v_Si_H_(i__),e__Si_H_(i__),o__Si_H_(i__),
+ v_Ca_H_(i__),e__Ca_H_(i__),o__Ca_H_(i__),v_ScI_H_(i__),
+ e__ScI_H_(i__),o__ScI_H_(i__),v_ScII_H_(i__),e__ScII_H_(i__),
+ o__ScII_H_(i__),v_TiI_H_(i__),e__TiI_H_(i__),o__TiI_H_(i__),
+ v_TiII_H_(i__),e__TiII_H_(i__),o__TiII_H_(i__),v_V_H_(i__),
+ e__V_H_(i__),o__V_H_(i__),v_CrI_H_(i__),e__CrI_H_(i__),
+ o__CrI_H_(i__),v_CrII_H_(i__),e__CrII_H_(i__),o__CrII_H_(i__),
+ v_Mn_H_(i__),e__Mn_H_(i__),o__Mn_H_(i__),v_Co_H_(i__),
+ e__Co_H_(i__),o__Co_H_(i__),v_Ni_H_(i__),e__Ni_H_(i__),
+ o__Ni_H_(i__),v_Na_H_(i__),e__Na_H_(i__),o__Na_H_(i__),
+ v_Mg_H_(i__),e__Mg_H_(i__),o__Mg_H_(i__),v_Al_H_(i__),
+ e__Al_H_(i__),o__Al_H_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end