Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/438/139 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-20
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/438/139 Abundances in Milky Way's disk (Soubiran+, 2005)
*================================================================================
*Abundance trends in kinematical groups of the Milky Way's disk.
* Soubiran C., Girard P.
* <Astron. Astrophys. 438, 139 (2005)>
* =2005A&A...438..139S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! Stellar data
integer*4 nr__
parameter (nr__=743) ! Number of records
character*135 ar__ ! Full-size record
character*11 Star (nr__) ! HD number (BD, CD, CPD, HIP, Giclas if missing)
real*4 v_Fe_H_ (nr__) ! (Sun) Metallicity
real*4 v_O_Fe_ (nr__) ! (Sun) ? Oxygen abundance relative to iron
real*4 v_Na_Fe_ (nr__) ! (Sun) ? Sodium abundance relative to iron
real*4 v_Mg_Fe_ (nr__) ! (Sun) ? Magnesium abundance relative to iron
real*4 v_Al_Fe_ (nr__) ! (Sun) ? Aluminium abundance relative to iron
real*4 v_Si_Fe_ (nr__) ! (Sun) Silicon abundance relative to iron
real*4 v_Ca_Fe_ (nr__) ! (Sun) ? Calcium abundance relative to iron
real*4 v_Ti_Fe_ (nr__) ! (Sun) ? Titanium abundance relative to iron
real*4 v_Ni_Fe_ (nr__) ! (Sun) ? Nickel abundance relative to iron
real*4 Uvel (nr__) ! (km/s) ? Galactic radial velocity (1)
real*4 Vvel (nr__) ! (km/s) ? Galactic rotational velocity (1)
real*4 Wvel (nr__) ! (km/s) ? Galactic vertical velocity (1)
real*4 Rmin (nr__) ! (kpc) ? Perigalactic distance
real*4 Rmax (nr__) ! (kpc) ? Apogalactic distance
real*4 Zmax (nr__) ! (kpc) ? Maximum distance from galactic plane
real*4 ecc (nr__) ! ? Eccentricity of galactic orbit
real*4 age (nr__) ! (Gyr) ? Age, in Gigayears
real*4 e_age (nr__) ! (%) ? Relative error on age
real*4 Pthin (nr__) ! (%) ? Thin disk membership probability
real*4 Pthick (nr__) ! (%) ? Thick disk membership probability
real*4 PHer (nr__) ! (%) ? Hercules stream membership probability
*Note (1): The three components of galactic velocity (U,V,W) are given with
* respect to the Sun, the U axis pointing towards the Galactic Center
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'param.dat' ! Stellar parameters
integer*4 nr__1
parameter (nr__1=743) ! Number of records
character*57 ar__1 ! Full-size record
character*11 Star_1 (nr__1) ! HD number (BD, CD, CPD, HIP, Giclas if missing)
integer*4 Teff (nr__1) ! (K) Effective temperature
real*4 logg (nr__1) ! ([cm/s2]) Gravity
real*4 v_Fe_H__1 (nr__1) ! (Sun) Metallicity
real*4 v_alpha_Fe_(nr__1) ! (Sun) Mean alpha abundance (1)
real*4 VMAG (nr__1) ! (mag) ? Absolute magnitude (2)
real*4 e_VMAG (nr__1) ! (mag) ? One sigma lower limit of Mv (2)
real*4 E_VMAG_1 (nr__1) ! (mag) ? One sigma upper limit of Mv (2)
*Note (1): [alpha/Fe]=0.25*([Mg/Fe]+[Si/Fe]+[Ca/Fe]+[Ti/Fe])
*Note (2): Absolute magnitudes are given for stars having an Hipparcos
* parallax greater than 2mas and a relative error lower than 50%.
* Apparent magnitudes have been deduced from Hipparcos magnitudes, Hp,
* transformed into V Jonhson according to the calibration by Harmanec
* (1998A&A...335..173H). One sigma lower and upper limits of Mv have
* been computed with {pi}+/-{sigma}_{pi}_.
C=============================================================================
C Loading file 'catalog.dat' ! Stellar data
C Format for file interpretation
1 format(
+ 1X,A11,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F6.1,1X,F6.1,1X,F6.1,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F4.2,1X,F4.1,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,743
read(1,'(A135)')ar__
read(ar__,1)
+ Star(i__),v_Fe_H_(i__),v_O_Fe_(i__),v_Na_Fe_(i__),
+ v_Mg_Fe_(i__),v_Al_Fe_(i__),v_Si_Fe_(i__),v_Ca_Fe_(i__),
+ v_Ti_Fe_(i__),v_Ni_Fe_(i__),Uvel(i__),Vvel(i__),Wvel(i__),
+ Rmin(i__),Rmax(i__),Zmax(i__),ecc(i__),age(i__),e_age(i__),
+ Pthin(i__),Pthick(i__),PHer(i__)
if(ar__(20:24) .EQ. '') v_O_Fe_(i__) = rNULL__
if(ar__(26:30) .EQ. '') v_Na_Fe_(i__) = rNULL__
if(ar__(32:36) .EQ. '') v_Mg_Fe_(i__) = rNULL__
if(ar__(38:42) .EQ. '') v_Al_Fe_(i__) = rNULL__
if(ar__(50:54) .EQ. '') v_Ca_Fe_(i__) = rNULL__
if(ar__(56:60) .EQ. '') v_Ti_Fe_(i__) = rNULL__
if(ar__(62:66) .EQ. '') v_Ni_Fe_(i__) = rNULL__
if(ar__(68:73) .EQ. '') Uvel(i__) = rNULL__
if(ar__(75:80) .EQ. '') Vvel(i__) = rNULL__
if(ar__(82:87) .EQ. '') Wvel(i__) = rNULL__
if(ar__(89:93) .EQ. '') Rmin(i__) = rNULL__
if(ar__(95:99) .EQ. '') Rmax(i__) = rNULL__
if(ar__(101:105) .EQ. '') Zmax(i__) = rNULL__
if(ar__(107:110) .EQ. '') ecc(i__) = rNULL__
if(ar__(112:115) .EQ. '') age(i__) = rNULL__
if(ar__(117:120) .EQ. '') e_age(i__) = rNULL__
if(ar__(122:125) .EQ. '') Pthin(i__) = rNULL__
if(ar__(127:130) .EQ. '') Pthick(i__) = rNULL__
if(ar__(132:135) .EQ. '') PHer(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Star(i__),v_Fe_H_(i__),v_O_Fe_(i__),v_Na_Fe_(i__),
+ v_Mg_Fe_(i__),v_Al_Fe_(i__),v_Si_Fe_(i__),v_Ca_Fe_(i__),
+ v_Ti_Fe_(i__),v_Ni_Fe_(i__),Uvel(i__),Vvel(i__),Wvel(i__),
+ Rmin(i__),Rmax(i__),Zmax(i__),ecc(i__),age(i__),e_age(i__),
+ Pthin(i__),Pthick(i__),PHer(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'param.dat' ! Stellar parameters
C Format for file interpretation
2 format(
+ 1X,A11,2X,I4,2X,F4.2,1X,F5.2,1X,F5.2,2X,F5.3,2X,F5.3,2X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'param.dat')
write(6,*) '....Loading file: param.dat'
do i__=1,743
read(1,'(A57)')ar__1
read(ar__1,2)
+ Star_1(i__),Teff(i__),logg(i__),v_Fe_H__1(i__),
+ v_alpha_Fe_(i__),VMAG(i__),e_VMAG(i__),E_VMAG_1(i__)
if(ar__1(39:43) .EQ. '') VMAG(i__) = rNULL__
if(ar__1(46:50) .EQ. '') e_VMAG(i__) = rNULL__
if(ar__1(53:57) .EQ. '') E_VMAG_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Star_1(i__),Teff(i__),logg(i__),v_Fe_H__1(i__),
+ v_alpha_Fe_(i__),VMAG(i__),e_VMAG(i__),E_VMAG_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end