Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/608/A145 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 2026-Oct-06
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 (nr__) ! [1/28] Star sequential number (1)
integer*4 Teff (nr__) ! (K) Spectroscopic effective temperature
real*4 v_Fe_H_ (nr__) ! ([-]) Metallicity of the atmosphere model
real*4 logg (nr__) ! ([cm/s2]) Spectroscopic surface gravity
real*4 vturb (nr__) ! (km/s) Microturbulence velocity
real*4 loggph (nr__) ! ([cm/s2]) Photometric surface gravity
real*4 e_loggph (nr__) ! ([cm/s2]) rms uncertainty on loggph (2)
character*9 JDInt (nr__) ! Days of the first and last exposures
integer*4 texp (nr__) ! (s) Total exposure time
real*8 RVavg (nr__) ! (km/s) Average radial velocity of all exposures
real*4 e_RVavg (nr__) ! (km/s) rms uncertainty on RVavg
integer*4 S_N4000 (nr__) ! ?=- Signal-to-noise ratios at ~4000{AA}
integer*4 S_N4900 (nr__) ! Signal-to-noise ratios at ~4900{AA}
integer*4 S_N6700 (nr__) ! Signal-to-noise ratios at ~6700{AA}
character*28 Inst (nr__) ! Instrument
real*8 JD1 (nr__) ! (d) Julian day for the middle of the exposure
real*8 RV1 (nr__) ! (km/s) Radial velocity of JD1 exposure
real*4 e_RV1 (nr__) ! (km/s) rms uncertainty on RV1
real*8 JD2 (nr__) ! (d) ? Julian day for the middle of the exposure
real*8 RV2 (nr__) ! (km/s) ? Radial velocity of JD2 exposure
real*4 e_RV2 (nr__) ! (km/s) ? rms uncertainty on RV2
real*8 JD3 (nr__) ! (d) ? Julian day for the middle of the exposure
real*8 RV3 (nr__) ! (km/s) ? Radial velocity of JD3 exposure
real*4 e_RV3 (nr__) ! (km/s) ? rms uncertainty on RV3
real*8 JD4 (nr__) ! (d) ? Julian day for the middle of the exposure
real*8 RV4 (nr__) ! (km/s) ? Radial velocity of JD4 exposure
real*4 e_RV4 (nr__) ! (km/s) ? rms uncertainty on RV4
real*8 JD5 (nr__) ! (d) ? Julian day for the middle of the exposure
real*8 RV5 (nr__) ! (km/s) ? Radial velocity of JD5 exposure
real*4 e_RV5 (nr__) ! (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 (nr__1) ! [1/28] Star sequential number
character*8 Name (nr__1) ! Star name
integer*8 Source (nr__1) ! Gaia source ID
real*8 RAdeg (nr__1) ! (deg) Gaia right ascension (ICRS) at epoch 2015.0
real*8 DEdeg (nr__1) ! (deg) Gaia declination (ICRS) at epoch 2015.0
real*8 GLON (nr__1) ! (deg) Gaia Galactic longitude at epoch 2015.0
real*8 GLAT (nr__1) ! (deg) Gaia Galactic latitude at epoch 2015.0
real*4 Gmag (nr__1) ! (mag) Gaia G-band mean magnitude
real*4 dH (nr__1) ! (kpc) ? Heliocentric distance
real*4 Rgal (nr__1) ! (kpc) ? Galactocentric distance
real*4 VGSR (nr__1) ! (km/s) ? Galactic Standard of Rest velocity
real*8 Lambda (nr__1) ! (deg) ? Longitude in a system aligned with the
* Sag stream
real*8 B (nr__1) ! (deg) ? Latitude in a system aligned with the
* Sag stream
character*7 Sub (nr__1) ! Substructure
character*32 Com (nr__1) ! Comments
character*12 File (nr__1) ! 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 (nr__2) ! Element designation
integer*4 Ion (nr__2) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_1 (nr__2) ! (0.1nm) Laboratory line wavelength
real*4 Chiex (nr__2) ! (eV) Excitation potential
real*8 log_gf (nr__2) ! Statistical weight times oscillator strength
real*4 EW01 (nr__2) ! (0.1pm) ? Observed equivalent width of star #01
real*4 e_EW01 (nr__2) ! (0.1pm) ? Error on EW01
real*4 Abund01 (nr__2) ! ? log of the abundance (H abundance = 12.00)
real*4 EW02 (nr__2) ! (0.1pm) ? Observed equivalent width of star #02
real*4 e_EW02 (nr__2) ! (0.1pm) ? Error on EW02
real*4 Abund02 (nr__2) ! ? log of the abundance (H abundance = 12.00)
real*4 EW03 (nr__2) ! (0.1pm) ? Observed equivalent width of star #03
real*4 e_EW03 (nr__2) ! (0.1pm) ? Error on EW03
real*4 Abund03 (nr__2) ! ? log of the abundance (H abundance = 12.00)
real*4 EW04 (nr__2) ! (0.1pm) ? Observed equivalent width of star #04
real*4 e_EW04 (nr__2) ! (0.1pm) ? Error on EW04
real*4 Abund04 (nr__2) ! ? log of the abundance (H abundance = 12.00)
real*4 EW21 (nr__2) ! (0.1pm) ? Observed equivalent width of star #21
real*4 e_EW21 (nr__2) ! (0.1pm) ? Error on EW21
real*4 Abund21 (nr__2) ! ? 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 (nr__3) ! Element designation
integer*4 Ion_1 (nr__3) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_2 (nr__3) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_1 (nr__3) ! (eV) Excitation potential
real*8 log_gf_1 (nr__3) ! Statistical weight times oscillator strength
real*4 EW26 (nr__3) ! (0.1pm) ? Observed equivalent width of star #26
real*4 e_EW26 (nr__3) ! (0.1pm) ? Error on EW26
real*4 Abund26 (nr__3) ! ? log of the abundance (H abundance = 12.00)
real*4 EW05 (nr__3) ! (0.1pm) ? Observed equivalent width of star #05
real*4 e_EW05 (nr__3) ! (0.1pm) ? Error on EW05
real*4 Abund05 (nr__3) ! ? log of the abundance (H abundance = 12.00)
real*4 EW06 (nr__3) ! (0.1pm) ? Observed equivalent width of star #06
real*4 e_EW06 (nr__3) ! (0.1pm) ? Error on EW06
real*4 Abund06 (nr__3) ! ? log of the abundance (H abundance = 12.00)
real*4 EW07 (nr__3) ! (0.1pm) ? Observed equivalent width of star #07
real*4 e_EW07 (nr__3) ! (0.1pm) ? Error on EW07
real*4 Abund07 (nr__3) ! ? log of the abundance (H abundance = 12.00)
real*4 EW28 (nr__3) ! (0.1pm) ? Observed equivalent width of star #28
real*4 e_EW28 (nr__3) ! (0.1pm) ? Error on EW28
real*4 Abund28 (nr__3) ! ? 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 (nr__4) ! Element designation
integer*4 Ion_2 (nr__4) ! Ionization state (1: neutral, 2: once ionized)
real*8 Lambda_3 (nr__4) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_2 (nr__4) ! (eV) Excitation potential
real*8 log_gf_1 (nr__4) ! Statistical weight times oscillator strength
real*4 EW08 (nr__4) ! (0.1pm) ? Observed equivalent width of star #08
real*4 e_EW08 (nr__4) ! (0.1pm) ? Error on EW08
real*4 Abund08 (nr__4) ! ? log of the abundance (H abundance = 12.00)
real*4 EW09 (nr__4) ! (0.1pm) ? Observed equivalent width of star #09
real*4 e_EW09 (nr__4) ! (0.1pm) ? Error on EW09
real*4 Abund09 (nr__4) ! ? log of the abundance (H abundance = 12.00)
real*4 EW10 (nr__4) ! (0.1pm) ? Observed equivalent width of star #10
real*4 e_EW10 (nr__4) ! (0.1pm) ? Error on EW10
character*1 l_Abund10 (nr__4) ! Limit flag on Abund10
real*4 Abund10 (nr__4) ! ? log of the abundance (H abundance = 12.00)
real*4 EW11 (nr__4) ! (0.1pm) ? Observed equivalent width of star #11
real*4 e_EW11 (nr__4) ! (0.1pm) ? Error on EW11
real*4 Abund11 (nr__4) ! ? log of the abundance (H abundance = 12.00)
real*4 EW12 (nr__4) ! (0.1pm) ? Observed equivalent width of star #12
real*4 e_EW12 (nr__4) ! (0.1pm) ? Error on EW12
character*1 l_Abund12 (nr__4) ! Limit flag on Abund12
real*4 Abund12 (nr__4) ! ? 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 (nr__5) ! Element designation
integer*4 Ion_3 (nr__5) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_4 (nr__5) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_3 (nr__5) ! (eV) Excitation potential
real*8 log_gf_1 (nr__5) ! Statistical weight times oscillator strength
real*4 EW13 (nr__5) ! (0.1pm) ? Observed equivalent width of star #13
real*4 e_EW13 (nr__5) ! (0.1pm) ? Error on EW13
real*4 Abund13 (nr__5) ! ? log of the abundance (H abundance = 12.00)
real*4 EW14 (nr__5) ! (0.1pm) ? Observed equivalent width of star #14
real*4 e_EW14 (nr__5) ! (0.1pm) ? Error on EW14
real*4 Abund14 (nr__5) ! ? log of the abundance (H abundance = 12.00)
real*4 EW15 (nr__5) ! (0.1pm) ? Observed equivalent width of star #15
real*4 e_EW15 (nr__5) ! (0.1pm) ? Error on EW15
real*4 Abund15 (nr__5) ! ? log of the abundance (H abundance = 12.00)
real*4 EW16 (nr__5) ! (0.1pm) ? Observed equivalent width of star #16
real*4 e_EW16 (nr__5) ! (0.1pm) ? Error on EW16
real*4 Abund16 (nr__5) ! ? log of the abundance (H abundance = 12.00)
real*4 EW17 (nr__5) ! (0.1pm) ? Observed equivalent width of star #17
real*4 e_EW17 (nr__5) ! (0.1pm) ? Error on EW17
real*4 Abund17 (nr__5) ! ? 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 (nr__6) ! Element designation
integer*4 Ion_4 (nr__6) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_5 (nr__6) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_4 (nr__6) ! (eV) Excitation potential
real*8 log_gf_1 (nr__6) ! Statistical weight times oscillator strength
real*4 EW18 (nr__6) ! (0.1pm) Observed equivalent width of star #18
real*4 e_EW18 (nr__6) ! (0.1pm) Error on EW18
real*4 Abund18 (nr__6) ! 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 (nr__7) ! Element designation
integer*4 Ion_5 (nr__7) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_6 (nr__7) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_5 (nr__7) ! (eV) Excitation potential
real*8 log_gf_1 (nr__7) ! Statistical weight times oscillator strength
real*4 EW19 (nr__7) ! (0.1pm) ? Observed equivalent width of star #19
real*4 e_EW19 (nr__7) ! (0.1pm) ? Error on EW19
real*4 Abund19 (nr__7) ! ? log of the abundance (H abundance = 12.00)
real*4 EW20 (nr__7) ! (0.1pm) ? Observed equivalent width of star #20
real*4 e_EW20 (nr__7) ! (0.1pm) ? Error on EW20
real*4 Abund20 (nr__7) ! ? log of the abundance (H abundance = 12.00)
real*4 EW21_1 (nr__7) ! (0.1pm) ? Observed equivalent width of star #21
real*4 e_EW21_1 (nr__7) ! (0.1pm) ? Error on EW21
real*4 Abund21_1 (nr__7) ! ? log of the abundance (H abundance = 12.00)
real*4 EW22 (nr__7) ! (0.1pm) ? Observed equivalent width of star #22
real*4 e_EW22 (nr__7) ! (0.1pm) ? Error on EW22
real*4 Abund22 (nr__7) ! ? log of the abundance (H abundance = 12.00)
real*4 EW23 (nr__7) ! (0.1pm) ? Observed equivalent width of star #23
real*4 e_EW23 (nr__7) ! (0.1pm) ? Error on EW23
character*1 l_Abund23 (nr__7) ! Limit flag on Abund23
real*4 Abund23 (nr__7) ! ? 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 (nr__8) ! Element designation
integer*4 Ion_6 (nr__8) ! Ionization state (1: neutral, 2: once ionized)
real*8 lambda_7 (nr__8) ! (0.1nm) Laboratory line wavelength
real*4 Chiex_6 (nr__8) ! (eV) Excitation potential
real*8 log_gf_1 (nr__8) ! Statistical weight times oscillator strength
real*4 EW24 (nr__8) ! (0.1pm) ? Observed equivalent width of star #24
real*4 e_EW24 (nr__8) ! (0.1pm) ? Error on EW24
real*4 Abund24 (nr__8) ! ? log of the abundance (H abundance = 12.00)
real*4 EW25 (nr__8) ! (0.1pm) ? Observed equivalent width of star #25
real*4 e_EW25 (nr__8) ! (0.1pm) ? Error on EW25
real*4 Abund25 (nr__8) ! ? log of the abundance (H abundance = 12.00)
real*4 EW26_1 (nr__8) ! (0.1pm) ? Observed equivalent width of star #26
real*4 e_EW26_1 (nr__8) ! (0.1pm) ? Error on EW26
real*4 Abund26_1 (nr__8) ! ? log of the abundance (H abundance = 12.00)
real*4 EW27 (nr__8) ! (0.1pm) ? Observed equivalent width of star #27
real*4 e_EW27 (nr__8) ! (0.1pm) ? Error on EW27
real*4 Abund27 (nr__8) ! ? log of the abundance (H abundance = 12.00)
real*4 EW28_1 (nr__8) ! (0.1pm) ? Observed equivalent width of star #28
real*4 e_EW28_1 (nr__8) ! (0.1pm) ? Error on EW28
real*4 Abund28_1 (nr__8) ! ? 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(i__),Teff(i__),v_Fe_H_(i__),logg(i__),vturb(i__),
+ loggph(i__),e_loggph(i__),JDInt(i__),texp(i__),RVavg(i__),
+ e_RVavg(i__),S_N4000(i__),S_N4900(i__),S_N6700(i__),Inst(i__),
+ JD1(i__),RV1(i__),e_RV1(i__),JD2(i__),RV2(i__),e_RV2(i__),
+ JD3(i__),RV3(i__),e_RV3(i__),JD4(i__),RV4(i__),e_RV4(i__),
+ JD5(i__),RV5(i__),e_RV5(i__)
if (S_N4000(i__) .EQ. 45) S_N4000(i__) = iNULL__
if(ar__(122:129) .EQ. '') JD2(i__) = rNULL__
if(ar__(131:137) .EQ. '') RV2(i__) = rNULL__
if(ar__(139:142) .EQ. '') e_RV2(i__) = rNULL__
if(ar__(144:151) .EQ. '') JD3(i__) = rNULL__
if(ar__(153:159) .EQ. '') RV3(i__) = rNULL__
if(ar__(161:164) .EQ. '') e_RV3(i__) = rNULL__
if(ar__(166:173) .EQ. '') JD4(i__) = rNULL__
if(ar__(175:181) .EQ. '') RV4(i__) = rNULL__
if(ar__(183:186) .EQ. '') e_RV4(i__) = rNULL__
if(ar__(188:195) .EQ. '') JD5(i__) = rNULL__
if(ar__(197:203) .EQ. '') RV5(i__) = rNULL__
if(ar__(205:208) .EQ. '') e_RV5(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),Teff(i__),v_Fe_H_(i__),logg(i__),vturb(i__),
+ loggph(i__),e_loggph(i__),JDInt(i__),texp(i__),RVavg(i__),
+ e_RVavg(i__),S_N4000(i__),S_N4900(i__),S_N6700(i__),Inst(i__),
+ JD1(i__),RV1(i__),e_RV1(i__),JD2(i__),RV2(i__),e_RV2(i__),
+ JD3(i__),RV3(i__),e_RV3(i__),JD4(i__),RV4(i__),e_RV4(i__),
+ JD5(i__),RV5(i__),e_RV5(i__)
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(i__),Name(i__),Source(i__),RAdeg(i__),DEdeg(i__),
+ GLON(i__),GLAT(i__),Gmag(i__),dH(i__),Rgal(i__),VGSR(i__),
+ Lambda(i__),B(i__),Sub(i__),Com(i__),File(i__)
if(ar__1(94:97) .EQ. '') dH(i__) = rNULL__
if(ar__1(99:102) .EQ. '') Rgal(i__) = rNULL__
if(ar__1(104:109) .EQ. '') VGSR(i__) = rNULL__
if(ar__1(111:118) .EQ. '') Lambda(i__) = rNULL__
if(ar__1(120:127) .EQ. '') B(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),Name(i__),Source(i__),RAdeg(i__),DEdeg(i__),
+ GLON(i__),GLAT(i__),Gmag(i__),dH(i__),Rgal(i__),VGSR(i__),
+ Lambda(i__),B(i__),Sub(i__),Com(i__),File(i__)
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(i__),Ion(i__),lambda_1(i__),Chiex(i__),log_gf(i__),
+ EW01(i__),e_EW01(i__),Abund01(i__),EW02(i__),e_EW02(i__),
+ Abund02(i__),EW03(i__),e_EW03(i__),Abund03(i__),EW04(i__),
+ e_EW04(i__),Abund04(i__),EW21(i__),e_EW21(i__),Abund21(i__)
if(ar__2(28:32) .EQ. '') EW01(i__) = rNULL__
if(ar__2(34:37) .EQ. '') e_EW01(i__) = rNULL__
if(ar__2(39:43) .EQ. '') Abund01(i__) = rNULL__
if(ar__2(45:49) .EQ. '') EW02(i__) = rNULL__
if(ar__2(51:54) .EQ. '') e_EW02(i__) = rNULL__
if(ar__2(56:60) .EQ. '') Abund02(i__) = rNULL__
if(ar__2(62:66) .EQ. '') EW03(i__) = rNULL__
if(ar__2(68:71) .EQ. '') e_EW03(i__) = rNULL__
if(ar__2(73:77) .EQ. '') Abund03(i__) = rNULL__
if(ar__2(79:83) .EQ. '') EW04(i__) = rNULL__
if(ar__2(85:88) .EQ. '') e_EW04(i__) = rNULL__
if(ar__2(90:94) .EQ. '') Abund04(i__) = rNULL__
if(ar__2(96:100) .EQ. '') EW21(i__) = rNULL__
if(ar__2(102:105) .EQ. '') e_EW21(i__) = rNULL__
if(ar__2(107:111) .EQ. '') Abund21(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ El(i__),Ion(i__),lambda_1(i__),Chiex(i__),log_gf(i__),
+ EW01(i__),e_EW01(i__),Abund01(i__),EW02(i__),e_EW02(i__),
+ Abund02(i__),EW03(i__),e_EW03(i__),Abund03(i__),EW04(i__),
+ e_EW04(i__),Abund04(i__),EW21(i__),e_EW21(i__),Abund21(i__)
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(i__),Ion_1(i__),lambda_2(i__),Chiex_1(i__),
+ log_gf_1(i__),EW26(i__),e_EW26(i__),Abund26(i__),EW05(i__),
+ e_EW05(i__),Abund05(i__),EW06(i__),e_EW06(i__),Abund06(i__),
+ EW07(i__),e_EW07(i__),Abund07(i__),EW28(i__),e_EW28(i__),
+ Abund28(i__)
if(ar__3(28:32) .EQ. '') EW26(i__) = rNULL__
if(ar__3(34:37) .EQ. '') e_EW26(i__) = rNULL__
if(ar__3(39:43) .EQ. '') Abund26(i__) = rNULL__
if(ar__3(45:49) .EQ. '') EW05(i__) = rNULL__
if(ar__3(51:54) .EQ. '') e_EW05(i__) = rNULL__
if(ar__3(56:60) .EQ. '') Abund05(i__) = rNULL__
if(ar__3(62:66) .EQ. '') EW06(i__) = rNULL__
if(ar__3(68:71) .EQ. '') e_EW06(i__) = rNULL__
if(ar__3(73:77) .EQ. '') Abund06(i__) = rNULL__
if(ar__3(79:83) .EQ. '') EW07(i__) = rNULL__
if(ar__3(85:88) .EQ. '') e_EW07(i__) = rNULL__
if(ar__3(90:94) .EQ. '') Abund07(i__) = rNULL__
if(ar__3(96:100) .EQ. '') EW28(i__) = rNULL__
if(ar__3(102:105) .EQ. '') e_EW28(i__) = rNULL__
if(ar__3(107:111) .EQ. '') Abund28(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Element(i__),Ion_1(i__),lambda_2(i__),Chiex_1(i__),
+ log_gf_1(i__),EW26(i__),e_EW26(i__),Abund26(i__),EW05(i__),
+ e_EW05(i__),Abund05(i__),EW06(i__),e_EW06(i__),Abund06(i__),
+ EW07(i__),e_EW07(i__),Abund07(i__),EW28(i__),e_EW28(i__),
+ Abund28(i__)
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(i__),Ion_2(i__),Lambda_3(i__),Chiex_2(i__),log_gf_1(i__),
+ EW08(i__),e_EW08(i__),Abund08(i__),EW09(i__),e_EW09(i__),
+ Abund09(i__),EW10(i__),e_EW10(i__),l_Abund10(i__),
+ Abund10(i__),EW11(i__),e_EW11(i__),Abund11(i__),EW12(i__),
+ e_EW12(i__),l_Abund12(i__),Abund12(i__)
if(ar__4(28:32) .EQ. '') EW08(i__) = rNULL__
if(ar__4(34:37) .EQ. '') e_EW08(i__) = rNULL__
if(ar__4(39:43) .EQ. '') Abund08(i__) = rNULL__
if(ar__4(45:49) .EQ. '') EW09(i__) = rNULL__
if(ar__4(51:54) .EQ. '') e_EW09(i__) = rNULL__
if(ar__4(56:60) .EQ. '') Abund09(i__) = rNULL__
if(ar__4(62:66) .EQ. '') EW10(i__) = rNULL__
if(ar__4(68:71) .EQ. '') e_EW10(i__) = rNULL__
if(ar__4(73:77) .EQ. '') Abund10(i__) = rNULL__
if(ar__4(79:83) .EQ. '') EW11(i__) = rNULL__
if(ar__4(85:88) .EQ. '') e_EW11(i__) = rNULL__
if(ar__4(90:94) .EQ. '') Abund11(i__) = rNULL__
if(ar__4(96:100) .EQ. '') EW12(i__) = rNULL__
if(ar__4(102:105) .EQ. '') e_EW12(i__) = rNULL__
if(ar__4(107:111) .EQ. '') Abund12(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ El_1(i__),Ion_2(i__),Lambda_3(i__),Chiex_2(i__),log_gf_1(i__),
+ EW08(i__),e_EW08(i__),Abund08(i__),EW09(i__),e_EW09(i__),
+ Abund09(i__),EW10(i__),e_EW10(i__),l_Abund10(i__),
+ Abund10(i__),EW11(i__),e_EW11(i__),Abund11(i__),EW12(i__),
+ e_EW12(i__),l_Abund12(i__),Abund12(i__)
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(i__),Ion_3(i__),lambda_4(i__),Chiex_3(i__),log_gf_1(i__),
+ EW13(i__),e_EW13(i__),Abund13(i__),EW14(i__),e_EW14(i__),
+ Abund14(i__),EW15(i__),e_EW15(i__),Abund15(i__),EW16(i__),
+ e_EW16(i__),Abund16(i__),EW17(i__),e_EW17(i__),Abund17(i__)
if(ar__5(28:32) .EQ. '') EW13(i__) = rNULL__
if(ar__5(34:37) .EQ. '') e_EW13(i__) = rNULL__
if(ar__5(39:43) .EQ. '') Abund13(i__) = rNULL__
if(ar__5(45:49) .EQ. '') EW14(i__) = rNULL__
if(ar__5(51:54) .EQ. '') e_EW14(i__) = rNULL__
if(ar__5(56:60) .EQ. '') Abund14(i__) = rNULL__
if(ar__5(62:66) .EQ. '') EW15(i__) = rNULL__
if(ar__5(68:71) .EQ. '') e_EW15(i__) = rNULL__
if(ar__5(73:77) .EQ. '') Abund15(i__) = rNULL__
if(ar__5(79:83) .EQ. '') EW16(i__) = rNULL__
if(ar__5(85:88) .EQ. '') e_EW16(i__) = rNULL__
if(ar__5(90:94) .EQ. '') Abund16(i__) = rNULL__
if(ar__5(96:100) .EQ. '') EW17(i__) = rNULL__
if(ar__5(102:105) .EQ. '') e_EW17(i__) = rNULL__
if(ar__5(107:111) .EQ. '') Abund17(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ El_2(i__),Ion_3(i__),lambda_4(i__),Chiex_3(i__),log_gf_1(i__),
+ EW13(i__),e_EW13(i__),Abund13(i__),EW14(i__),e_EW14(i__),
+ Abund14(i__),EW15(i__),e_EW15(i__),Abund15(i__),EW16(i__),
+ e_EW16(i__),Abund16(i__),EW17(i__),e_EW17(i__),Abund17(i__)
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(i__),Ion_4(i__),lambda_5(i__),Chiex_4(i__),log_gf_1(i__),
+ EW18(i__),e_EW18(i__),Abund18(i__)
c ..............Just test output...........
write(6,7)
+ El_3(i__),Ion_4(i__),lambda_5(i__),Chiex_4(i__),log_gf_1(i__),
+ EW18(i__),e_EW18(i__),Abund18(i__)
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(i__),Ion_5(i__),lambda_6(i__),Chiex_5(i__),log_gf_1(i__),
+ EW19(i__),e_EW19(i__),Abund19(i__),EW20(i__),e_EW20(i__),
+ Abund20(i__),EW21_1(i__),e_EW21_1(i__),Abund21_1(i__),
+ EW22(i__),e_EW22(i__),Abund22(i__),EW23(i__),e_EW23(i__),
+ l_Abund23(i__),Abund23(i__)
if(ar__7(28:32) .EQ. '') EW19(i__) = rNULL__
if(ar__7(34:37) .EQ. '') e_EW19(i__) = rNULL__
if(ar__7(39:43) .EQ. '') Abund19(i__) = rNULL__
if(ar__7(45:49) .EQ. '') EW20(i__) = rNULL__
if(ar__7(51:54) .EQ. '') e_EW20(i__) = rNULL__
if(ar__7(56:60) .EQ. '') Abund20(i__) = rNULL__
if(ar__7(62:66) .EQ. '') EW21_1(i__) = rNULL__
if(ar__7(68:71) .EQ. '') e_EW21_1(i__) = rNULL__
if(ar__7(73:77) .EQ. '') Abund21_1(i__) = rNULL__
if(ar__7(79:83) .EQ. '') EW22(i__) = rNULL__
if(ar__7(85:88) .EQ. '') e_EW22(i__) = rNULL__
if(ar__7(90:94) .EQ. '') Abund22(i__) = rNULL__
if(ar__7(96:100) .EQ. '') EW23(i__) = rNULL__
if(ar__7(102:105) .EQ. '') e_EW23(i__) = rNULL__
if(ar__7(107:111) .EQ. '') Abund23(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ El_4(i__),Ion_5(i__),lambda_6(i__),Chiex_5(i__),log_gf_1(i__),
+ EW19(i__),e_EW19(i__),Abund19(i__),EW20(i__),e_EW20(i__),
+ Abund20(i__),EW21_1(i__),e_EW21_1(i__),Abund21_1(i__),
+ EW22(i__),e_EW22(i__),Abund22(i__),EW23(i__),e_EW23(i__),
+ l_Abund23(i__),Abund23(i__)
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(i__),Ion_6(i__),lambda_7(i__),Chiex_6(i__),log_gf_1(i__),
+ EW24(i__),e_EW24(i__),Abund24(i__),EW25(i__),e_EW25(i__),
+ Abund25(i__),EW26_1(i__),e_EW26_1(i__),Abund26_1(i__),
+ EW27(i__),e_EW27(i__),Abund27(i__),EW28_1(i__),e_EW28_1(i__),
+ Abund28_1(i__)
if(ar__8(28:32) .EQ. '') EW24(i__) = rNULL__
if(ar__8(34:37) .EQ. '') e_EW24(i__) = rNULL__
if(ar__8(39:43) .EQ. '') Abund24(i__) = rNULL__
if(ar__8(45:49) .EQ. '') EW25(i__) = rNULL__
if(ar__8(51:54) .EQ. '') e_EW25(i__) = rNULL__
if(ar__8(56:60) .EQ. '') Abund25(i__) = rNULL__
if(ar__8(62:66) .EQ. '') EW26_1(i__) = rNULL__
if(ar__8(68:71) .EQ. '') e_EW26_1(i__) = rNULL__
if(ar__8(73:77) .EQ. '') Abund26_1(i__) = rNULL__
if(ar__8(79:83) .EQ. '') EW27(i__) = rNULL__
if(ar__8(85:88) .EQ. '') e_EW27(i__) = rNULL__
if(ar__8(90:94) .EQ. '') Abund27(i__) = rNULL__
if(ar__8(96:100) .EQ. '') EW28_1(i__) = rNULL__
if(ar__8(102:105) .EQ. '') e_EW28_1(i__) = rNULL__
if(ar__8(107:111) .EQ. '') Abund28_1(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ El_5(i__),Ion_6(i__),lambda_7(i__),Chiex_6(i__),log_gf_1(i__),
+ EW24(i__),e_EW24(i__),Abund24(i__),EW25(i__),e_EW25(i__),
+ Abund25(i__),EW26_1(i__),e_EW26_1(i__),Abund26_1(i__),
+ EW27(i__),e_EW27(i__),Abund27(i__),EW28_1(i__),e_EW28_1(i__),
+ Abund28_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end