Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/130/2140 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 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/AJ/130/2140 Abundances and velocities in globular clusters (Pritzl+, 2005)
*================================================================================
*A comparison of elemental abundance ratios in globular clusters, field stars,
*and dwarf spheroidal galaxies.
* Pritzl B.J., Venn K.A., Irwin M.
* <Astron. J., 130, 2140-2165 (2005)>
* =2005AJ....130.2140P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Globular cluster mean abundance ratios
integer*4 nr__
parameter (nr__=106) ! Number of records
character*270 ar__ ! Full-size record
character*8 Cluster ! Globular cluster name
real*4 v_Fe_H_ ! ([Sun]) Log of the Fe-to-H ratio
character*1 f__Fe_H_ ! [ib] Flag on [Fe/H] (1)
real*4 v_Mg_Fe_ ! ([Sun]) ? Log of the Mg-to-Fe ratio
real*4 e__Mg_Fe_ ! ([Sun]) ? Standard error of the mean for [Mg/Fe]
character*1 f__Mg_Fe_ ! [ib] Flag on [Mg/Fe] (1)
real*4 v_Si_Fe_ ! ([Sun]) ? Log of the Si-to-Fe ratio
real*4 e__Si_Fe_ ! ([Sun]) ? Standard error of the mean for [Si/Fe]
character*1 f__Si_Fe_ ! [ib] Flag on [Si/Fe] (1)
real*4 v_Ca_Fe_ ! ([Sun]) ? Log of the Ca-to-Fe ratio
real*4 e__Ca_Fe_ ! ([Sun]) ? Standard error of the mean for [Ca/Fe]
character*1 f__Ca_Fe_ ! [ib] Flag on [Ca/Fe] (1)
real*4 v_TiI_Fe_ ! ([Sun]) ? Log of the TiI-to-Fe ratio
real*4 e__TiI_Fe_ ! ([Sun]) ? Standard error of the mean for [TiI/Fe]
character*1 f__TiI_Fe_ ! [ib] Flag on [TiI/Fe] (1)
real*4 v_TiII_Fe_ ! ([Sun]) ? Log of the TiII-to-Fe ratio
real*4 e__TiII_Fe_ ! ([Sun]) ? Standard error of the mean for [TiII/Fe]
character*1 f__TiII_Fe_ ! [ib] Flag on [TiII/Fe] (1)
real*4 v_Y_Fe_ ! ([Sun]) ? Log of the Y-to-Fe ratio
real*4 e__Y_Fe_ ! ([Sun]) ? Standard error of the mean for [Y/Fe]
character*1 f__Y_Fe_ ! [ib] Flag on [Y/Fe] (1)
real*4 v_Ba_Fe_ ! ([Sun]) ? Log of the Ba-to-Fe ratio
real*4 e__Ba_Fe_ ! ([Sun]) ? Standard error of the mean for [Ba/Fe]
character*1 f__Ba_Fe_ ! [ib] Flag on [Ba/Fe] (1)
real*4 v_La_Fe_ ! ([Sun]) ? Log of the La-to-Fe ratio
real*4 e__La_Fe_ ! ([Sun]) ? Standard error of the mean for [La/Fe]
character*1 f__La_Fe_ ! [ib] Flag on [La/Fe] (1)
real*4 v_Eu_Fe_ ! ([Sun]) ? Log of the Eu-to-Fe ratio
real*4 e__Eu_Fe_ ! ([Sun]) ? Standard error of the mean for [Eu/Fe]
character*1 f__Eu_Fe_ ! [ib] Flag on [Eu/Fe] (1)
real*4 v_Ti_Fe_ ! ([Sun]) ? Log of the Ti-to-Fe ratio
real*4 e__Ti_Fe_ ! ([Sun]) ? Standard error of the mean for [Ti/Fe]
character*1 f__Ti_Fe_ ! [ib] Flag on [Ti/Fe] (1)
real*4 v_a_Fe_ ! ([Sun]) ? [{alpha}/Fe] ratio ({alpha} elements)
character*1 f__a_Fe_ ! [ib] Flag on [alpha/Fe] (1)
real*4 v_Ba_Y_ ! ([Sun]) ? Log of the Ba-to-Y ratio
character*1 f__Ba_Y_ ! [ib] Flag on [Ba/Y] (1)
real*4 v_Y_Eu_ ! ([Sun]) ? Log of the Y-to-Eu ratio
character*1 f__Y_Eu_ ! [ib] Flag on [Y/Eu] (1)
real*4 v_Ba_Eu_ ! ([Sun]) ? Log of the Ba-to-Eu ratio
character*1 f__Ba_Eu_ ! [ib] Flag on [Ba/Eu] (1)
real*4 v_La_Eu_ ! ([Sun]) ? Log of the La-to-Eu ratio
character*1 f__La_Eu_ ! [ib] Flag on [La/Eu] (1)
character*83 Ref ! Reference
*Note (1): Flag defined as follows:
* i = ratio was not used in calculating the final weighted mean
* abundance ratios;
* b = final individual abundance ratio.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Globular cluster velocity and positional vector
components
integer*4 nr__1
parameter (nr__1=45) ! Number of records
character*71 ar__1 ! Full-size record
character*8 Cluster_1 ! Cluster name
real*4 Xgal ! (kpc) X component of the galactocentric
* positional vector
real*4 Ygal ! (kpc) Y component of the galactocentric
* positional vector
real*4 Zgal ! (kpc) Z component of the galactocentric
* positional vector
real*4 Pi ! (km/s) ? {Pi} component of the galactocentric
* velocity vector
real*4 e_Pi ! (km/s) ? Error in Pi
real*4 Theta ! (km/s) ? {Theta} component of the galactocentric
* velocity vector
real*4 e_Theta ! (km/s) ? Error in Theta
real*4 W ! (km/s) ? W component of the galactocentric
* velocity vector
real*4 e_W ! (km/s) ? Error in W
character*2 Cl1 ! Galactic component the globular cluster
* belongs to, according to our calculations (1)
character*2 Cl2 ! Galactic component the globular cluster
* belongs to, taken from Mackey & Gilmore
* (2004MNRAS.355..504M) (2)
*Note (1): Class Cl1 is defined as follows:
* B = bulge
* H = halo
* TK = thick disk
* TN = thin disk
*Note (2): Class Cl2 is defined as follows:
* YH = young halo
* OH = old halo
* BD = bulge/disk
* SG = Sgr cluster
C=============================================================================
C Loading file 'table2.dat' ! Globular cluster mean abundance ratios
C Format for file interpretation
1 format(
+ A8,1X,F5.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,
+ 1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,
+ 1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,
+ F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,
+ F5.2,1X,A1,1X,F5.2,1X,A1,1X,F5.2,1X,A1,1X,F5.2,1X,A1,1X,F5.2,
+ 1X,A1,1X,A83)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,106
read(1,'(A270)')ar__
read(ar__,1)
+ Cluster,v_Fe_H_,f__Fe_H_,v_Mg_Fe_,e__Mg_Fe_,f__Mg_Fe_,
+ v_Si_Fe_,e__Si_Fe_,f__Si_Fe_,v_Ca_Fe_,e__Ca_Fe_,f__Ca_Fe_,
+ v_TiI_Fe_,e__TiI_Fe_,f__TiI_Fe_,v_TiII_Fe_,e__TiII_Fe_,
+ f__TiII_Fe_,v_Y_Fe_,e__Y_Fe_,f__Y_Fe_,v_Ba_Fe_,e__Ba_Fe_,
+ f__Ba_Fe_,v_La_Fe_,e__La_Fe_,f__La_Fe_,v_Eu_Fe_,e__Eu_Fe_,
+ f__Eu_Fe_,v_Ti_Fe_,e__Ti_Fe_,f__Ti_Fe_,v_a_Fe_,f__a_Fe_,
+ v_Ba_Y_,f__Ba_Y_,v_Y_Eu_,f__Y_Eu_,v_Ba_Eu_,f__Ba_Eu_,v_La_Eu_,
+ f__La_Eu_,Ref
if(ar__(18:22) .EQ. '') v_Mg_Fe_ = rNULL__
if(ar__(24:27) .EQ. '') e__Mg_Fe_ = rNULL__
if(ar__(31:35) .EQ. '') v_Si_Fe_ = rNULL__
if(ar__(37:40) .EQ. '') e__Si_Fe_ = rNULL__
if(ar__(44:48) .EQ. '') v_Ca_Fe_ = rNULL__
if(ar__(50:53) .EQ. '') e__Ca_Fe_ = rNULL__
if(ar__(57:61) .EQ. '') v_TiI_Fe_ = rNULL__
if(ar__(63:66) .EQ. '') e__TiI_Fe_ = rNULL__
if(ar__(70:74) .EQ. '') v_TiII_Fe_ = rNULL__
if(ar__(76:79) .EQ. '') e__TiII_Fe_ = rNULL__
if(ar__(83:87) .EQ. '') v_Y_Fe_ = rNULL__
if(ar__(89:92) .EQ. '') e__Y_Fe_ = rNULL__
if(ar__(96:100) .EQ. '') v_Ba_Fe_ = rNULL__
if(ar__(102:105) .EQ. '') e__Ba_Fe_ = rNULL__
if(ar__(109:113) .EQ. '') v_La_Fe_ = rNULL__
if(ar__(115:118) .EQ. '') e__La_Fe_ = rNULL__
if(ar__(122:126) .EQ. '') v_Eu_Fe_ = rNULL__
if(ar__(128:131) .EQ. '') e__Eu_Fe_ = rNULL__
if(ar__(135:139) .EQ. '') v_Ti_Fe_ = rNULL__
if(ar__(141:144) .EQ. '') e__Ti_Fe_ = rNULL__
if(ar__(148:152) .EQ. '') v_a_Fe_ = rNULL__
if(ar__(156:160) .EQ. '') v_Ba_Y_ = rNULL__
if(ar__(164:168) .EQ. '') v_Y_Eu_ = rNULL__
if(ar__(172:176) .EQ. '') v_Ba_Eu_ = rNULL__
if(ar__(180:184) .EQ. '') v_La_Eu_ = rNULL__
c ..............Just test output...........
write(6,1)
+ Cluster,v_Fe_H_,f__Fe_H_,v_Mg_Fe_,e__Mg_Fe_,f__Mg_Fe_,
+ v_Si_Fe_,e__Si_Fe_,f__Si_Fe_,v_Ca_Fe_,e__Ca_Fe_,f__Ca_Fe_,
+ v_TiI_Fe_,e__TiI_Fe_,f__TiI_Fe_,v_TiII_Fe_,e__TiII_Fe_,
+ f__TiII_Fe_,v_Y_Fe_,e__Y_Fe_,f__Y_Fe_,v_Ba_Fe_,e__Ba_Fe_,
+ f__Ba_Fe_,v_La_Fe_,e__La_Fe_,f__La_Fe_,v_Eu_Fe_,e__Eu_Fe_,
+ f__Eu_Fe_,v_Ti_Fe_,e__Ti_Fe_,f__Ti_Fe_,v_a_Fe_,f__a_Fe_,
+ v_Ba_Y_,f__Ba_Y_,v_Y_Eu_,f__Y_Eu_,v_Ba_Eu_,f__Ba_Eu_,v_La_Eu_,
+ f__La_Eu_,Ref
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Globular cluster velocity and positional vector
* components
C Format for file interpretation
2 format(
+ A8,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,
+ F6.1,1X,F4.1,1X,A2,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,45
read(1,'(A71)')ar__1
read(ar__1,2)
+ Cluster_1,Xgal,Ygal,Zgal,Pi,e_Pi,Theta,e_Theta,W,e_W,Cl1,Cl2
if(ar__1(31:36) .EQ. '') Pi = rNULL__
if(ar__1(38:41) .EQ. '') e_Pi = rNULL__
if(ar__1(43:48) .EQ. '') Theta = rNULL__
if(ar__1(50:53) .EQ. '') e_Theta = rNULL__
if(ar__1(55:60) .EQ. '') W = rNULL__
if(ar__1(62:65) .EQ. '') e_W = rNULL__
c ..............Just test output...........
write(6,2)
+ Cluster_1,Xgal,Ygal,Zgal,Pi,e_Pi,Theta,e_Theta,W,e_W,Cl1,Cl2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end