Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/608/A145 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 2026-Oct-08
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/608/A145 What is the Milky Way outer halo made of? (Battaglia+, 2017)
*================================================================================
*What is the Milky Way outer halo made of?
*High resolution spectroscopy of distant red giants.
* Battaglia G., North P.L., Jablonka P., Shetrone M., Minniti D., Diaz M.,
* Starkenburg E., Savoy M.
* <Astron. Astrophys. 608, A145 (2017)>
* =2017A&A...608A.145B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Atmospheric stellar parameters and details of
observations for the sample stars
integer*4 nr__
parameter (nr__=31) ! Number of records
character*208 ar__ ! Full-size record
integer*4 Seq ! [1/28] Star sequential number (1)
integer*4 Teff ! (K) Spectroscopic effective temperature
real*4 v_Fe_H_ ! ([-]) Metallicity of the atmosphere model
real*4 logg ! ([cm/s2]) Spectroscopic surface gravity
real*4 vturb ! (km/s) Microturbulence velocity
real*4 loggph ! ([cm/s2]) Photometric surface gravity
real*4 e_loggph ! ([cm/s2]) rms uncertainty on loggph (2)
character*9 JDInt ! Days of the first and last exposures
integer*4 texp ! (s) Total exposure time
real*8 RVavg ! (km/s) Average radial velocity of all exposures
real*4 e_RVavg ! (km/s) rms uncertainty on RVavg
integer*4 S_N4000 ! ?=- Signal-to-noise ratios at ~4000{AA}
integer*4 S_N4900 ! Signal-to-noise ratios at ~4900{AA}
integer*4 S_N6700 ! Signal-to-noise ratios at ~6700{AA}
character*28 Inst ! Instrument
real*8 JD1 ! (d) Julian day for the middle of the exposure
real*8 RV1 ! (km/s) Radial velocity of JD1 exposure
real*4 e_RV1 ! (km/s) rms uncertainty on RV1
real*8 JD2 ! (d) ? Julian day for the middle of the exposure
real*8 RV2 ! (km/s) ? Radial velocity of JD2 exposure
real*4 e_RV2 ! (km/s) ? rms uncertainty on RV2
real*8 JD3 ! (d) ? Julian day for the middle of the exposure
real*8 RV3 ! (km/s) ? Radial velocity of JD3 exposure
real*4 e_RV3 ! (km/s) ? rms uncertainty on RV3
real*8 JD4 ! (d) ? Julian day for the middle of the exposure
real*8 RV4 ! (km/s) ? Radial velocity of JD4 exposure
real*4 e_RV4 ! (km/s) ? rms uncertainty on RV4
real*8 JD5 ! (d) ? Julian day for the middle of the exposure
real*8 RV5 ! (km/s) ? Radial velocity of JD5 exposure
real*4 e_RV5 ! (km/s) ? rms uncertainty on RV5
*Note (1): Sources: all stars observed with HRS and UVES, and stars #8-12 with
* MIKE, were selected from the Spaghetti survey, as listed in
* (Starkenburg, 2011, Ph.D. Thesis, Kapteyn Astronomical Institute,
* University of Groningen); star #11 comes from the APOGEE survey;
* stars from #13-18 were drawn from the Xue et al. (2014, Cat. J/ApJ/784/170)
* catalog.
*Note (2): The photometric logg of most HRS stars had asymmetric error bars;
* the error quoted here is the larger one of the upper and lower errors.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Coordinates, magnitude, distance and velocity
of the stars in the sample
integer*4 nr__1
parameter (nr__1=31) ! Number of records
character*181 ar__1 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2015.0)
integer*4 Seq_1 ! [1/28] Star sequential number
character*8 Name ! Star name
integer*8 Source ! Gaia source ID
real*8 RAdeg ! (deg) Gaia right ascension (ICRS) at epoch 2015.0
real*8 DEdeg ! (deg) Gaia declination (ICRS) at epoch 2015.0
real*8 GLON ! (deg) Gaia Galactic longitude at epoch 2015.0
real*8 GLAT ! (deg) Gaia Galactic latitude at epoch 2015.0
real*4 Gmag ! (mag) Gaia G-band mean magnitude
real*4 dH ! (kpc) ? Heliocentric distance
real*4 Rgal ! (kpc) ? Galactocentric distance
real*4 VGSR ! (km/s) ? Galactic Standard of Rest velocity
real*8 Lambda ! (deg) ? Longitude in a system aligned with the
* Sag stream
real*8 B ! (deg) ? Latitude in a system aligned with the
* Sag stream
character*7 Sub ! Substructure
character*32 Com ! Comments
character*12 File ! Name of the table with abundances
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea5.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #1, #2, #3, #4, and #21
observed with HRS
integer*4 nr__2
parameter (nr__2=215) ! Number of records
character*111 ar__2 ! Full-size record
character*2 El ! Element designation
integer*4 Ion ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_1 ! (0.1nm) Laboratory line wavelength
real*4 Chiex ! (eV) Excitation potential
real*8 log_gf ! Statistical weight times oscillator strength
real*4 EW01 ! (0.1pm) ? Observed equivalent width of star #01
real*4 e_EW01 ! (0.1pm) ? Error on EW01
real*4 Abund01 ! ? log of the abundance (H abundance = 12.00)
real*4 EW02 ! (0.1pm) ? Observed equivalent width of star #02
real*4 e_EW02 ! (0.1pm) ? Error on EW02
real*4 Abund02 ! ? log of the abundance (H abundance = 12.00)
real*4 EW03 ! (0.1pm) ? Observed equivalent width of star #03
real*4 e_EW03 ! (0.1pm) ? Error on EW03
real*4 Abund03 ! ? log of the abundance (H abundance = 12.00)
real*4 EW04 ! (0.1pm) ? Observed equivalent width of star #04
real*4 e_EW04 ! (0.1pm) ? Error on EW04
real*4 Abund04 ! ? log of the abundance (H abundance = 12.00)
real*4 EW21 ! (0.1pm) ? Observed equivalent width of star #21
real*4 e_EW21 ! (0.1pm) ? Error on EW21
real*4 Abund21 ! ? log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea6.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #26, #5, #6, #7, and #28
observed with HRS
integer*4 nr__3
parameter (nr__3=214) ! Number of records
character*111 ar__3 ! Full-size record
character*2 Element ! Element designation
integer*4 Ion_1 ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_2 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_1 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW26 ! (0.1pm) ? Observed equivalent width of star #26
real*4 e_EW26 ! (0.1pm) ? Error on EW26
real*4 Abund26 ! ? log of the abundance (H abundance = 12.00)
real*4 EW05 ! (0.1pm) ? Observed equivalent width of star #05
real*4 e_EW05 ! (0.1pm) ? Error on EW05
real*4 Abund05 ! ? log of the abundance (H abundance = 12.00)
real*4 EW06 ! (0.1pm) ? Observed equivalent width of star #06
real*4 e_EW06 ! (0.1pm) ? Error on EW06
real*4 Abund06 ! ? log of the abundance (H abundance = 12.00)
real*4 EW07 ! (0.1pm) ? Observed equivalent width of star #07
real*4 e_EW07 ! (0.1pm) ? Error on EW07
real*4 Abund07 ! ? log of the abundance (H abundance = 12.00)
real*4 EW28 ! (0.1pm) ? Observed equivalent width of star #28
real*4 e_EW28 ! (0.1pm) ? Error on EW28
real*4 Abund28 ! ? log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea7.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #8, #9, #10, #11, and #12
observed with MIKE
integer*4 nr__4
parameter (nr__4=374) ! Number of records
character*111 ar__4 ! Full-size record
character*2 El_1 ! Element designation
integer*4 Ion_2 ! Ionization state (1: neutral, 2: once ionized)
real*8 Lambda_3 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_2 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW08 ! (0.1pm) ? Observed equivalent width of star #08
real*4 e_EW08 ! (0.1pm) ? Error on EW08
real*4 Abund08 ! ? log of the abundance (H abundance = 12.00)
real*4 EW09 ! (0.1pm) ? Observed equivalent width of star #09
real*4 e_EW09 ! (0.1pm) ? Error on EW09
real*4 Abund09 ! ? log of the abundance (H abundance = 12.00)
real*4 EW10 ! (0.1pm) ? Observed equivalent width of star #10
real*4 e_EW10 ! (0.1pm) ? Error on EW10
character*1 l_Abund10 ! Limit flag on Abund10
real*4 Abund10 ! ? log of the abundance (H abundance = 12.00)
real*4 EW11 ! (0.1pm) ? Observed equivalent width of star #11
real*4 e_EW11 ! (0.1pm) ? Error on EW11
real*4 Abund11 ! ? log of the abundance (H abundance = 12.00)
real*4 EW12 ! (0.1pm) ? Observed equivalent width of star #12
real*4 e_EW12 ! (0.1pm) ? Error on EW12
character*1 l_Abund12 ! Limit flag on Abund12
real*4 Abund12 ! ? log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea8.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #13, #14, #15, #16, and #17
observed with MIKE
integer*4 nr__5
parameter (nr__5=283) ! Number of records
character*111 ar__5 ! Full-size record
character*2 El_2 ! Element designation
integer*4 Ion_3 ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_4 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_3 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW13 ! (0.1pm) ? Observed equivalent width of star #13
real*4 e_EW13 ! (0.1pm) ? Error on EW13
real*4 Abund13 ! ? log of the abundance (H abundance = 12.00)
real*4 EW14 ! (0.1pm) ? Observed equivalent width of star #14
real*4 e_EW14 ! (0.1pm) ? Error on EW14
real*4 Abund14 ! ? log of the abundance (H abundance = 12.00)
real*4 EW15 ! (0.1pm) ? Observed equivalent width of star #15
real*4 e_EW15 ! (0.1pm) ? Error on EW15
real*4 Abund15 ! ? log of the abundance (H abundance = 12.00)
real*4 EW16 ! (0.1pm) ? Observed equivalent width of star #16
real*4 e_EW16 ! (0.1pm) ? Error on EW16
real*4 Abund16 ! ? log of the abundance (H abundance = 12.00)
real*4 EW17 ! (0.1pm) ? Observed equivalent width of star #17
real*4 e_EW17 ! (0.1pm) ? Error on EW17
real*4 Abund17 ! ? log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea9.dat' ! Lines parameters, observed equivalent widths and
abundances of star #18 observed with MIKE
integer*4 nr__6
parameter (nr__6=104) ! Number of records
character*43 ar__6 ! Full-size record
character*2 El_3 ! Element designation
integer*4 Ion_4 ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_5 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_4 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW18 ! (0.1pm) Observed equivalent width of star #18
real*4 e_EW18 ! (0.1pm) Error on EW18
real*4 Abund18 ! log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea10.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #19, #20, #21, #22, and #23
observed with UVES
integer*4 nr__7
parameter (nr__7=312) ! Number of records
character*111 ar__7 ! Full-size record
character*2 El_4 ! Element designation
integer*4 Ion_5 ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_6 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_5 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW19 ! (0.1pm) ? Observed equivalent width of star #19
real*4 e_EW19 ! (0.1pm) ? Error on EW19
real*4 Abund19 ! ? log of the abundance (H abundance = 12.00)
real*4 EW20 ! (0.1pm) ? Observed equivalent width of star #20
real*4 e_EW20 ! (0.1pm) ? Error on EW20
real*4 Abund20 ! ? log of the abundance (H abundance = 12.00)
real*4 EW21_1 ! (0.1pm) ? Observed equivalent width of star #21
real*4 e_EW21_1 ! (0.1pm) ? Error on EW21
real*4 Abund21_1 ! ? log of the abundance (H abundance = 12.00)
real*4 EW22 ! (0.1pm) ? Observed equivalent width of star #22
real*4 e_EW22 ! (0.1pm) ? Error on EW22
real*4 Abund22 ! ? log of the abundance (H abundance = 12.00)
real*4 EW23 ! (0.1pm) ? Observed equivalent width of star #23
real*4 e_EW23 ! (0.1pm) ? Error on EW23
character*1 l_Abund23 ! Limit flag on Abund23
real*4 Abund23 ! ? log of the abundance (H abundance = 12.00)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea11.dat' ! Lines parameters, observed equivalent widths and
abundances of stars #24, #25, #26, #27, and #28
observed with UVES
integer*4 nr__8
parameter (nr__8=297) ! Number of records
character*111 ar__8 ! Full-size record
character*2 El_5 ! Element designation
integer*4 Ion_6 ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_7 ! (0.1nm) Laboratory line wavelength
real*4 Chiex_6 ! (eV) Excitation potential
real*8 log_gf_1 ! Statistical weight times oscillator strength
real*4 EW24 ! (0.1pm) ? Observed equivalent width of star #24
real*4 e_EW24 ! (0.1pm) ? Error on EW24
real*4 Abund24 ! ? log of the abundance (H abundance = 12.00)
real*4 EW25 ! (0.1pm) ? Observed equivalent width of star #25
real*4 e_EW25 ! (0.1pm) ? Error on EW25
real*4 Abund25 ! ? log of the abundance (H abundance = 12.00)
real*4 EW26_1 ! (0.1pm) ? Observed equivalent width of star #26
real*4 e_EW26_1 ! (0.1pm) ? Error on EW26
real*4 Abund26_1 ! ? log of the abundance (H abundance = 12.00)
real*4 EW27 ! (0.1pm) ? Observed equivalent width of star #27
real*4 e_EW27 ! (0.1pm) ? Error on EW27
real*4 Abund27 ! ? log of the abundance (H abundance = 12.00)
real*4 EW28_1 ! (0.1pm) ? Observed equivalent width of star #28
real*4 e_EW28_1 ! (0.1pm) ? Error on EW28
real*4 Abund28_1 ! ? log of the abundance (H abundance = 12.00)
C=============================================================================
C Loading file 'tablea1.dat' ! Atmospheric stellar parameters and details of
* observations for the sample stars
C Format for file interpretation
1 format(
+ I2,1X,I4,1X,F3.1,1X,F5.2,1X,F4.2,1X,F3.1,1X,F3.1,1X,A9,1X,I5,
+ 1X,F7.2,1X,F4.2,1X,I3,1X,I2,1X,I3,1X,A28,F8.3,1X,F7.2,1X,F4.2,
+ 1X,F8.3,1X,F7.2,1X,F4.2,1X,F8.3,1X,F7.2,1X,F4.2,1X,F8.3,1X,
+ F7.2,1X,F4.2,1X,F8.3,1X,F7.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,31
read(1,'(A208)')ar__
read(ar__,1)
+ Seq,Teff,v_Fe_H_,logg,vturb,loggph,e_loggph,JDInt,texp,RVavg,
+ e_RVavg,S_N4000,S_N4900,S_N6700,Inst,JD1,RV1,e_RV1,JD2,RV2,
+ e_RV2,JD3,RV3,e_RV3,JD4,RV4,e_RV4,JD5,RV5,e_RV5
if (S_N4000 .EQ. 45) S_N4000 = iNULL__
if(ar__(122:129) .EQ. '') JD2 = rNULL__
if(ar__(131:137) .EQ. '') RV2 = rNULL__
if(ar__(139:142) .EQ. '') e_RV2 = rNULL__
if(ar__(144:151) .EQ. '') JD3 = rNULL__
if(ar__(153:159) .EQ. '') RV3 = rNULL__
if(ar__(161:164) .EQ. '') e_RV3 = rNULL__
if(ar__(166:173) .EQ. '') JD4 = rNULL__
if(ar__(175:181) .EQ. '') RV4 = rNULL__
if(ar__(183:186) .EQ. '') e_RV4 = rNULL__
if(ar__(188:195) .EQ. '') JD5 = rNULL__
if(ar__(197:203) .EQ. '') RV5 = rNULL__
if(ar__(205:208) .EQ. '') e_RV5 = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq,Teff,v_Fe_H_,logg,vturb,loggph,e_loggph,JDInt,texp,RVavg,
+ e_RVavg,S_N4000,S_N4900,S_N6700,Inst,JD1,RV1,e_RV1,JD2,RV2,
+ e_RV2,JD3,RV3,e_RV3,JD4,RV4,e_RV4,JD5,RV5,e_RV5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Coordinates, magnitude, distance and velocity
* of the stars in the sample
C Format for file interpretation
2 format(
+ I2,1X,A8,1X,I19,1X,F13.9,F13.9,1X,F13.9,F13.9,1X,F6.3,1X,F4.1,
+ 1X,F4.1,1X,F6.1,1X,F8.4,1X,F8.4,1X,A7,1X,A32,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,31
read(1,'(A181)')ar__1
read(ar__1,2)
+ Seq_1,Name,Source,RAdeg,DEdeg,GLON,GLAT,Gmag,dH,Rgal,VGSR,
+ Lambda,B,Sub,Com,File
if(ar__1(94:97) .EQ. '') dH = rNULL__
if(ar__1(99:102) .EQ. '') Rgal = rNULL__
if(ar__1(104:109) .EQ. '') VGSR = rNULL__
if(ar__1(111:118) .EQ. '') Lambda = rNULL__
if(ar__1(120:127) .EQ. '') B = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1,Name,Source,RAdeg,DEdeg,GLON,GLAT,Gmag,dH,Rgal,VGSR,
+ Lambda,B,Sub,Com,File
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea5.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #1, #2, #3, #4, and #21
* observed with HRS
C Format for file interpretation
3 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea5.dat')
write(6,*) '....Loading file: tablea5.dat'
do i__=1,215
read(1,'(A111)')ar__2
read(ar__2,3)
+ El,Ion,lambda_1,Chiex,log_gf,EW01,e_EW01,Abund01,EW02,e_EW02,
+ Abund02,EW03,e_EW03,Abund03,EW04,e_EW04,Abund04,EW21,e_EW21,
+ Abund21
if(ar__2(28:32) .EQ. '') EW01 = rNULL__
if(ar__2(34:37) .EQ. '') e_EW01 = rNULL__
if(ar__2(39:43) .EQ. '') Abund01 = rNULL__
if(ar__2(45:49) .EQ. '') EW02 = rNULL__
if(ar__2(51:54) .EQ. '') e_EW02 = rNULL__
if(ar__2(56:60) .EQ. '') Abund02 = rNULL__
if(ar__2(62:66) .EQ. '') EW03 = rNULL__
if(ar__2(68:71) .EQ. '') e_EW03 = rNULL__
if(ar__2(73:77) .EQ. '') Abund03 = rNULL__
if(ar__2(79:83) .EQ. '') EW04 = rNULL__
if(ar__2(85:88) .EQ. '') e_EW04 = rNULL__
if(ar__2(90:94) .EQ. '') Abund04 = rNULL__
if(ar__2(96:100) .EQ. '') EW21 = rNULL__
if(ar__2(102:105) .EQ. '') e_EW21 = rNULL__
if(ar__2(107:111) .EQ. '') Abund21 = rNULL__
c ..............Just test output...........
write(6,3)
+ El,Ion,lambda_1,Chiex,log_gf,EW01,e_EW01,Abund01,EW02,e_EW02,
+ Abund02,EW03,e_EW03,Abund03,EW04,e_EW04,Abund04,EW21,e_EW21,
+ Abund21
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea6.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #26, #5, #6, #7, and #28
* observed with HRS
C Format for file interpretation
4 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea6.dat')
write(6,*) '....Loading file: tablea6.dat'
do i__=1,214
read(1,'(A111)')ar__3
read(ar__3,4)
+ Element,Ion_1,lambda_2,Chiex_1,log_gf_1,EW26,e_EW26,Abund26,
+ EW05,e_EW05,Abund05,EW06,e_EW06,Abund06,EW07,e_EW07,Abund07,
+ EW28,e_EW28,Abund28
if(ar__3(28:32) .EQ. '') EW26 = rNULL__
if(ar__3(34:37) .EQ. '') e_EW26 = rNULL__
if(ar__3(39:43) .EQ. '') Abund26 = rNULL__
if(ar__3(45:49) .EQ. '') EW05 = rNULL__
if(ar__3(51:54) .EQ. '') e_EW05 = rNULL__
if(ar__3(56:60) .EQ. '') Abund05 = rNULL__
if(ar__3(62:66) .EQ. '') EW06 = rNULL__
if(ar__3(68:71) .EQ. '') e_EW06 = rNULL__
if(ar__3(73:77) .EQ. '') Abund06 = rNULL__
if(ar__3(79:83) .EQ. '') EW07 = rNULL__
if(ar__3(85:88) .EQ. '') e_EW07 = rNULL__
if(ar__3(90:94) .EQ. '') Abund07 = rNULL__
if(ar__3(96:100) .EQ. '') EW28 = rNULL__
if(ar__3(102:105) .EQ. '') e_EW28 = rNULL__
if(ar__3(107:111) .EQ. '') Abund28 = rNULL__
c ..............Just test output...........
write(6,4)
+ Element,Ion_1,lambda_2,Chiex_1,log_gf_1,EW26,e_EW26,Abund26,
+ EW05,e_EW05,Abund05,EW06,e_EW06,Abund06,EW07,e_EW07,Abund07,
+ EW28,e_EW28,Abund28
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea7.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #8, #9, #10, #11, and #12
* observed with MIKE
C Format for file interpretation
5 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,A1,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,A1,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea7.dat')
write(6,*) '....Loading file: tablea7.dat'
do i__=1,374
read(1,'(A111)')ar__4
read(ar__4,5)
+ El_1,Ion_2,Lambda_3,Chiex_2,log_gf_1,EW08,e_EW08,Abund08,EW09,
+ e_EW09,Abund09,EW10,e_EW10,l_Abund10,Abund10,EW11,e_EW11,
+ Abund11,EW12,e_EW12,l_Abund12,Abund12
if(ar__4(28:32) .EQ. '') EW08 = rNULL__
if(ar__4(34:37) .EQ. '') e_EW08 = rNULL__
if(ar__4(39:43) .EQ. '') Abund08 = rNULL__
if(ar__4(45:49) .EQ. '') EW09 = rNULL__
if(ar__4(51:54) .EQ. '') e_EW09 = rNULL__
if(ar__4(56:60) .EQ. '') Abund09 = rNULL__
if(ar__4(62:66) .EQ. '') EW10 = rNULL__
if(ar__4(68:71) .EQ. '') e_EW10 = rNULL__
if(ar__4(73:77) .EQ. '') Abund10 = rNULL__
if(ar__4(79:83) .EQ. '') EW11 = rNULL__
if(ar__4(85:88) .EQ. '') e_EW11 = rNULL__
if(ar__4(90:94) .EQ. '') Abund11 = rNULL__
if(ar__4(96:100) .EQ. '') EW12 = rNULL__
if(ar__4(102:105) .EQ. '') e_EW12 = rNULL__
if(ar__4(107:111) .EQ. '') Abund12 = rNULL__
c ..............Just test output...........
write(6,5)
+ El_1,Ion_2,Lambda_3,Chiex_2,log_gf_1,EW08,e_EW08,Abund08,EW09,
+ e_EW09,Abund09,EW10,e_EW10,l_Abund10,Abund10,EW11,e_EW11,
+ Abund11,EW12,e_EW12,l_Abund12,Abund12
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea8.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #13, #14, #15, #16, and #17
* observed with MIKE
C Format for file interpretation
6 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea8.dat')
write(6,*) '....Loading file: tablea8.dat'
do i__=1,283
read(1,'(A111)')ar__5
read(ar__5,6)
+ El_2,Ion_3,lambda_4,Chiex_3,log_gf_1,EW13,e_EW13,Abund13,EW14,
+ e_EW14,Abund14,EW15,e_EW15,Abund15,EW16,e_EW16,Abund16,EW17,
+ e_EW17,Abund17
if(ar__5(28:32) .EQ. '') EW13 = rNULL__
if(ar__5(34:37) .EQ. '') e_EW13 = rNULL__
if(ar__5(39:43) .EQ. '') Abund13 = rNULL__
if(ar__5(45:49) .EQ. '') EW14 = rNULL__
if(ar__5(51:54) .EQ. '') e_EW14 = rNULL__
if(ar__5(56:60) .EQ. '') Abund14 = rNULL__
if(ar__5(62:66) .EQ. '') EW15 = rNULL__
if(ar__5(68:71) .EQ. '') e_EW15 = rNULL__
if(ar__5(73:77) .EQ. '') Abund15 = rNULL__
if(ar__5(79:83) .EQ. '') EW16 = rNULL__
if(ar__5(85:88) .EQ. '') e_EW16 = rNULL__
if(ar__5(90:94) .EQ. '') Abund16 = rNULL__
if(ar__5(96:100) .EQ. '') EW17 = rNULL__
if(ar__5(102:105) .EQ. '') e_EW17 = rNULL__
if(ar__5(107:111) .EQ. '') Abund17 = rNULL__
c ..............Just test output...........
write(6,6)
+ El_2,Ion_3,lambda_4,Chiex_3,log_gf_1,EW13,e_EW13,Abund13,EW14,
+ e_EW14,Abund14,EW15,e_EW15,Abund15,EW16,e_EW16,Abund16,EW17,
+ e_EW17,Abund17
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea9.dat' ! Lines parameters, observed equivalent widths and
* abundances of star #18 observed with MIKE
C Format for file interpretation
7 format(A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea9.dat')
write(6,*) '....Loading file: tablea9.dat'
do i__=1,104
read(1,'(A43)')ar__6
read(ar__6,7)
+ El_3,Ion_4,lambda_5,Chiex_4,log_gf_1,EW18,e_EW18,Abund18
c ..............Just test output...........
write(6,7)
+ El_3,Ion_4,lambda_5,Chiex_4,log_gf_1,EW18,e_EW18,Abund18
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea10.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #19, #20, #21, #22, and #23
* observed with UVES
C Format for file interpretation
8 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,A1,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea10.dat')
write(6,*) '....Loading file: tablea10.dat'
do i__=1,312
read(1,'(A111)')ar__7
read(ar__7,8)
+ El_4,Ion_5,lambda_6,Chiex_5,log_gf_1,EW19,e_EW19,Abund19,EW20,
+ e_EW20,Abund20,EW21_1,e_EW21_1,Abund21_1,EW22,e_EW22,Abund22,
+ EW23,e_EW23,l_Abund23,Abund23
if(ar__7(28:32) .EQ. '') EW19 = rNULL__
if(ar__7(34:37) .EQ. '') e_EW19 = rNULL__
if(ar__7(39:43) .EQ. '') Abund19 = rNULL__
if(ar__7(45:49) .EQ. '') EW20 = rNULL__
if(ar__7(51:54) .EQ. '') e_EW20 = rNULL__
if(ar__7(56:60) .EQ. '') Abund20 = rNULL__
if(ar__7(62:66) .EQ. '') EW21_1 = rNULL__
if(ar__7(68:71) .EQ. '') e_EW21_1 = rNULL__
if(ar__7(73:77) .EQ. '') Abund21_1 = rNULL__
if(ar__7(79:83) .EQ. '') EW22 = rNULL__
if(ar__7(85:88) .EQ. '') e_EW22 = rNULL__
if(ar__7(90:94) .EQ. '') Abund22 = rNULL__
if(ar__7(96:100) .EQ. '') EW23 = rNULL__
if(ar__7(102:105) .EQ. '') e_EW23 = rNULL__
if(ar__7(107:111) .EQ. '') Abund23 = rNULL__
c ..............Just test output...........
write(6,8)
+ El_4,Ion_5,lambda_6,Chiex_5,log_gf_1,EW19,e_EW19,Abund19,EW20,
+ e_EW20,Abund20,EW21_1,e_EW21_1,Abund21_1,EW22,e_EW22,Abund22,
+ EW23,e_EW23,l_Abund23,Abund23
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea11.dat' ! Lines parameters, observed equivalent widths and
* abundances of stars #24, #25, #26, #27, and #28
* observed with UVES
C Format for file interpretation
9 format(
+ A2,1X,I1,1X,F8.3,1X,F4.2,1X,F7.3,1X,F5.1,1X,F4.1,1X,F5.2,1X,
+ F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.2,1X,F5.1,1X,F4.1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea11.dat')
write(6,*) '....Loading file: tablea11.dat'
do i__=1,297
read(1,'(A111)')ar__8
read(ar__8,9)
+ El_5,Ion_6,lambda_7,Chiex_6,log_gf_1,EW24,e_EW24,Abund24,EW25,
+ e_EW25,Abund25,EW26_1,e_EW26_1,Abund26_1,EW27,e_EW27,Abund27,
+ EW28_1,e_EW28_1,Abund28_1
if(ar__8(28:32) .EQ. '') EW24 = rNULL__
if(ar__8(34:37) .EQ. '') e_EW24 = rNULL__
if(ar__8(39:43) .EQ. '') Abund24 = rNULL__
if(ar__8(45:49) .EQ. '') EW25 = rNULL__
if(ar__8(51:54) .EQ. '') e_EW25 = rNULL__
if(ar__8(56:60) .EQ. '') Abund25 = rNULL__
if(ar__8(62:66) .EQ. '') EW26_1 = rNULL__
if(ar__8(68:71) .EQ. '') e_EW26_1 = rNULL__
if(ar__8(73:77) .EQ. '') Abund26_1 = rNULL__
if(ar__8(79:83) .EQ. '') EW27 = rNULL__
if(ar__8(85:88) .EQ. '') e_EW27 = rNULL__
if(ar__8(90:94) .EQ. '') Abund27 = rNULL__
if(ar__8(96:100) .EQ. '') EW28_1 = rNULL__
if(ar__8(102:105) .EQ. '') e_EW28_1 = rNULL__
if(ar__8(107:111) .EQ. '') Abund28_1 = rNULL__
c ..............Just test output...........
write(6,9)
+ El_5,Ion_6,lambda_7,Chiex_6,log_gf_1,EW24,e_EW24,Abund24,EW25,
+ e_EW25,Abund25,EW26_1,e_EW26_1,Abund26_1,EW27,e_EW27,Abund27,
+ EW28_1,e_EW28_1,Abund28_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end