Conversion of standardized ReadMe file for
file /./ftp/cats/III/193 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. III/193 Chemical Abundances in Late-Type Stars (Thevenin, 1998)
*================================================================================
*Chemical Abundances in Late-Type Stars
* Thevenin F.
* <Bull. Inf. CDS 49 (unpublished) >
* =1998yCat.3193....0T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! Abundances from Li to Eu for individual stars
integer*4 nr__
parameter (nr__=1108) ! Number of records
character*184 ar__ ! Full-size record
character*8 Cluster (nr__) ! Name of cluster
character*11 Star (nr__) ! Star name
real*4 theta (nr__) ! (K-1) Effective temperature (5040/T) (Sun=0.87)
real*4 logg (nr__) ! ([cm/s2]) Gravity (Sun=4.4)
real*4 zeta (nr__) ! (km/s) Microturbulence velocity (Sun=1.1)
real*4 v_Li_ (nr__) ! ([Sun]) ? Lithium (Z= 3) abundance
real*4 v_O_ (nr__) ! ([Sun]) ? Oxygen (Z= 8) abundance (Sun=8.92)
real*4 v_Na_ (nr__) ! ([Sun]) ? Natrium (Z=11) abundance (Sun=6.28)
real*4 v_Mg_ (nr__) ! ([Sun]) ? Magnesium (Z=12) abundance (Sun=7.53)
real*4 v_Al_ (nr__) ! ([Sun]) ? Aluminium (Z=13) abundance (Sun=6.43)
real*4 v_Si_ (nr__) ! ([Sun]) ? Silicium (Z=14) abundance (Sun=7.50)
real*4 v_Ca_ (nr__) ! ([Sun]) ? Calcium (Z=20) abundance (Sun=6.36)
real*4 v_Sc_ (nr__) ! ([Sun]) ? Scandium (Z=21) abundance (Sun=2.99)
real*4 v_Ti_ (nr__) ! ([Sun]) ? Titanium (Z=22) abundance (Sun=4.88)
real*4 v_V_ (nr__) ! ([Sun]) ? Vanadium (Z=23) abundance (Sun=3.91)
real*4 v_Cr_ (nr__) ! ([Sun]) ? Chromium (Z=24) abundance (Sun=5.61)
real*4 v_Mn_ (nr__) ! ([Sun]) ? Maganese (Z=25) abundance (Sun=5.47)
real*4 v_Fe_ (nr__) ! ([Sun]) ? Iron (Z=26) abundance (Sun=7.46)
real*4 v_Co_ (nr__) ! ([Sun]) ? Cobalt (Z=27) abundance (Sun=4.85)
real*4 v_Ni_ (nr__) ! ([Sun]) ? Nickel (Z=28) abundance (Sun=6.18)
real*4 v_Sr_ (nr__) ! ([Sun]) ? Strontium (Z=38) abundance (Sun=2.93)
real*4 v_Y_ (nr__) ! ([Sun]) ? Yttrium (Z=39) abundance (Sun=2.18)
real*4 v_Zr_ (nr__) ! ([Sun]) ? Zirconium (Z=40) abundance (Sun=2.46)
real*4 v_Mo_ (nr__) ! ([Sun]) ? Molybdene (Z=42) abundance (Sun=2.10)
real*4 v_Ba_ (nr__) ! ([Sun]) ? Barium (Z=56) abundance (Sun=2.18)
real*4 v_La_ (nr__) ! ([Sun]) ? Lanthane (Z=57) abundance (Sun=1.07)
real*4 v_Ce_ (nr__) ! ([Sun]) ? Cerium (Z=58) abundance (Sun=1.58)
real*4 v_Nd_ (nr__) ! ([Sun]) ? Neodyme (Z=60) abundance (Sun=1.40)
real*4 v_Sm_ (nr__) ! ([Sun]) ? Samarium (Z=62) abundance (Sun=0.88)
real*4 v_Eu_ (nr__) ! ([Sun]) ? Europium (Z=63) abundance (Sun=0.48)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cluster.dat' ! [Fe/H] Global values for Clusters
integer*4 nr__1
parameter (nr__1=45) ! Number of records
character*45 ar__1 ! Full-size record
character*19 Cluster_1 (nr__1) ! Cluster Name
real*4 v_Fe_H_ (nr__1) ! ([Sun]) Iron abundance of the cluster
real*4 e__Fe_H_ (nr__1) ! ([Sun]) Mean error on [Fe/H]
C=============================================================================
C Loading file 'catalog.dat' ! Abundances from Li to Eu for individual stars
C Format for file interpretation
1 format(
+ A8,A11,1X,F4.2,1X,F3.1,1X,F4.1,2X,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,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,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,1108
read(1,'(A184)')ar__
read(ar__,1)
+ Cluster(i__),Star(i__),theta(i__),logg(i__),zeta(i__),
+ v_Li_(i__),v_O_(i__),v_Na_(i__),v_Mg_(i__),v_Al_(i__),
+ v_Si_(i__),v_Ca_(i__),v_Sc_(i__),v_Ti_(i__),v_V_(i__),
+ v_Cr_(i__),v_Mn_(i__),v_Fe_(i__),v_Co_(i__),v_Ni_(i__),
+ v_Sr_(i__),v_Y_(i__),v_Zr_(i__),v_Mo_(i__),v_Ba_(i__),
+ v_La_(i__),v_Ce_(i__),v_Nd_(i__),v_Sm_(i__),v_Eu_(i__)
if(ar__(36:40) .EQ. '') v_Li_(i__) = rNULL__
if(ar__(42:46) .EQ. '') v_O_(i__) = rNULL__
if(ar__(48:52) .EQ. '') v_Na_(i__) = rNULL__
if(ar__(54:58) .EQ. '') v_Mg_(i__) = rNULL__
if(ar__(60:64) .EQ. '') v_Al_(i__) = rNULL__
if(ar__(66:70) .EQ. '') v_Si_(i__) = rNULL__
if(ar__(72:76) .EQ. '') v_Ca_(i__) = rNULL__
if(ar__(78:82) .EQ. '') v_Sc_(i__) = rNULL__
if(ar__(84:88) .EQ. '') v_Ti_(i__) = rNULL__
if(ar__(90:94) .EQ. '') v_V_(i__) = rNULL__
if(ar__(96:100) .EQ. '') v_Cr_(i__) = rNULL__
if(ar__(102:106) .EQ. '') v_Mn_(i__) = rNULL__
if(ar__(108:112) .EQ. '') v_Fe_(i__) = rNULL__
if(ar__(114:118) .EQ. '') v_Co_(i__) = rNULL__
if(ar__(120:124) .EQ. '') v_Ni_(i__) = rNULL__
if(ar__(126:130) .EQ. '') v_Sr_(i__) = rNULL__
if(ar__(132:136) .EQ. '') v_Y_(i__) = rNULL__
if(ar__(138:142) .EQ. '') v_Zr_(i__) = rNULL__
if(ar__(144:148) .EQ. '') v_Mo_(i__) = rNULL__
if(ar__(150:154) .EQ. '') v_Ba_(i__) = rNULL__
if(ar__(156:160) .EQ. '') v_La_(i__) = rNULL__
if(ar__(162:166) .EQ. '') v_Ce_(i__) = rNULL__
if(ar__(168:172) .EQ. '') v_Nd_(i__) = rNULL__
if(ar__(174:178) .EQ. '') v_Sm_(i__) = rNULL__
if(ar__(180:184) .EQ. '') v_Eu_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Cluster(i__),Star(i__),theta(i__),logg(i__),zeta(i__),
+ v_Li_(i__),v_O_(i__),v_Na_(i__),v_Mg_(i__),v_Al_(i__),
+ v_Si_(i__),v_Ca_(i__),v_Sc_(i__),v_Ti_(i__),v_V_(i__),
+ v_Cr_(i__),v_Mn_(i__),v_Fe_(i__),v_Co_(i__),v_Ni_(i__),
+ v_Sr_(i__),v_Y_(i__),v_Zr_(i__),v_Mo_(i__),v_Ba_(i__),
+ v_La_(i__),v_Ce_(i__),v_Nd_(i__),v_Sm_(i__),v_Eu_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cluster.dat' ! [Fe/H] Global values for Clusters
C Format for file interpretation
2 format(A19,14X,F5.2,3X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'cluster.dat')
write(6,*) '....Loading file: cluster.dat'
do i__=1,45
read(1,'(A45)')ar__1
read(ar__1,2)Cluster_1(i__),v_Fe_H_(i__),e__Fe_H_(i__)
c ..............Just test output...........
write(6,2)Cluster_1(i__),v_Fe_H_(i__),e__Fe_H_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end