Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/578/A116 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-07
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/578/A116 Abundances of red giant stars in NGC6093 (M80) (Carretta+, 2015)
*================================================================================
*Detailed abundances of a large sample of giant stars in NGC 6093 (M80).
* Carretta E. Bragaglia A., Gratton R., D'Orazi V., Lucatello S., Sollima A.,
* Momany Y., Catanzaro G., Leone F.
* <Astron. Astrophys. 578, A116 (2015)>
* =2015A&A...578A.116C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Positions, BV magnitudes and radial velocities
integer*4 nr__
parameter (nr__=143) ! Number of records
character*80 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000.0
real*8 DEdeg ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID ! [208/22172] Identification number
integer*4 RAh ! (h) Right ascension (J2000.0)
integer*4 RAm ! (min) Right ascension (J2000.0)
real*4 RAs ! (s) Right ascension (J2000.0)
character*1 DE_ ! Declination sign (J2000.0)
integer*4 DEd ! (deg) Declination (J2000.0)
integer*4 DEm ! (arcmin) Declination (J2000.0)
real*4 DEs ! (arcsec) Declination (J2000.0)
real*4 Bmag ! (mag) [14.7/17.9] B magnitude
real*4 Vmag ! (mag) [12.8/16.8] V magnitude
character*9 setup ! Giraffe setup (HR11,HR13) and/or UVES
real*4 HRV ! (km/s) [-94/260] Heliocentric radial velocity
real*4 e_HRV ! (km/s) [0/6.6] error on HRV
real*4 Dist ! (arcmin) [0.99/12.9] Distance from centre of NGC 6093
character*6 flag ! [M NM,MET RV] Membership flag (1)
*Note (1): Membership flag as follows:
* M = member
* NM,MET = not-a-member, from metallicity
* NM,RV = not-a-member, from radial velocity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Teff, logg, [A/H], Vt, [FeI/H], [FeII/H]
integer*4 nr__1
parameter (nr__1=82) ! Number of records
character*58 ar__1 ! Full-size record
integer*4 ID_1 ! [208/22172] Identification number
integer*4 Teff ! (K) [4046/5138] Effective temperature
real*4 logg ! ([cm/s2]) [0.42/2.63] Surface gravity
real*4 v_A_H_ ! ([Sun]) [-1.9/-1.7] Metallicity
real*4 vt ! (km/s) [0.65/3.57] Microturbulent velocity
integer*4 nFeI ! [3/76] Number of FeI lines
real*4 v_FeI_H_ ! ([Sun]) [-1.9/-1.7] Iron I abundance
real*4 e__FeI_H_ ! ([Sun]) [0/0.21] Error on [FeI/H]
integer*4 nFeII ! [1/10]? Number of FeII lines
real*4 v_FeII_H_ ! ([Sun]) [-1.9/-1.6]? Iron II abundance
real*4 e__FeII_H_ ! ([Sun]) [0/0.22]?=9.999 Error on [FeII/H]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe]
integer*4 nr__2
parameter (nr__2=82) ! Number of records
character*59 ar__2 ! Full-size record
integer*4 ID_2 ! [208/22172] Identification number
integer*4 nO ! number of O lines
character*1 l__O_Fe_ ! Limit flag in O (2)
real*4 v_O_Fe_ ! ([Sun]) ? Oxygen abundance
real*4 e__O_Fe_ ! ([Sun]) ? rms uncertainty on [O/Fe]
integer*4 nNa ! ? Number of Na lines
real*4 v_Na_Fe_ ! ([Sun]) ? Sodium abundance
real*4 e__Na_Fe_ ! ([Sun]) ? rms uncertainty on [Na/Fe]
integer*4 nMg ! Number of Mg lines
real*4 v_Mg_Fe_ ! ([Sun]) ? Magnesium abundance
real*4 e__Mg_Fe_ ! ([Sun]) ? rms uncertainty on [Mg/Fe]
integer*4 nAl ! Number of Al lines
character*1 l__Al_Fe_ ! Limit flag for Al (2)
real*4 v_Al_Fe_ ! ([Sun]) ? Aluminum abundance
real*4 e__Al_Fe_ ! ([Sun]) ? rms uncertainty on [Al/Fe]
*Note (2): < for limit, else measure.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! [Si/Fe], [Ca/Fe], [Ti/Fe] I & II
integer*4 nr__3
parameter (nr__3=82) ! Number of records
character*62 ar__3 ! Full-size record
integer*4 ID_3 ! [208/22172] Identification number
integer*4 nSi ! Number of Si lines
real*4 v_Si_Fe_ ! ([Sun]) ? Silicon abundance
real*4 e__Si_Fe_ ! ([Sun]) ? rms uncertainty on [Si/Fe]
integer*4 nCa ! Number of Ca lines
real*4 v_Ca_Fe_ ! ([Sun]) Calcium abundance
real*4 e__Ca_Fe_ ! ([Sun]) ? rms uncertainty on [ca/Fe]
integer*4 nTi1 ! Number of Ti I lines
real*4 v_TiI_Fe_ ! ([Sun]) ? Titanium I abundance
real*4 e__TiI_Fe_ ! ([Sun]) ? rms uncertainty on [TiI/Fe]
integer*4 nTi2 ! Number of Ti II lines
real*4 v_TiII_Fe_ ! ([Sun]) ? Titanium II abundance
real*4 e__TiII_Fe_ ! ([Sun]) ? rms uncertainty on [TiII/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! [Sc/Fe], [V/Fe], [Cr/Fe], [Mn/Fe], [Co/Fe],
[Ni/Fe], [Cu/Fe], [Zn.Fe]
integer*4 nr__4
parameter (nr__4=82) ! Number of records
character*124 ar__4 ! Full-size record
integer*4 ID_4 ! [208/22172] Identification number
integer*4 nSc ! Number of Sc lines
real*4 v_Sc_Fe_ ! ([Sun]) Scandium abundance
real*4 e__Sc_Fe_ ! ([Sun]) ? rms uncertainty on [Sc/Fe]
integer*4 nV ! Number of V lines
real*4 v_V_Fe_ ! ([Sun]) ? Vanadium abundance
real*4 e__V_Fe_ ! ([Sun]) ? rms uncertainty on [V/Fe]
integer*4 nCrI ! Number of CrI lines
real*4 v_CrI_Fe_ ! ([Sun]) ? Cromium abundance
real*4 e__CrI_Fe_ ! ([Sun]) ? rms uncertainty on [CrI/Fe]
integer*4 nCrII ! Number of CrII lines
real*4 v_CrII_Fe_ ! ([Sun]) ? Cromium abundance
real*4 e__CrII_Fe_ ! ([Sun]) ? rms uncertainty on [CrII/Fe]
integer*4 nMn ! Number of Mn lines
real*4 v_Mn_Fe_ ! ([Sun]) ? Manganese abundance
real*4 e__Mn_Fe_ ! ([Sun]) ? rms uncertainty on [Mn/Fe]
integer*4 nCo ! Number of Co lines
real*4 v_Co_Fe_ ! ([Sun]) ? Cobalt abundance
real*4 e__Co_Fe_ ! ([Sun]) ? rms uncertainty on [Co/Fe]
integer*4 nNi ! Number of Ni lines
real*4 v_Ni_Fe_ ! ([Sun]) ? Nichel abundance
real*4 e__Ni_Fe_ ! ([Sun]) ? rms uncertainty on [Ni/Fe]
integer*4 nCu ! Number of Cu lines
real*4 v_Cu_Fe_ ! ([Sun]) ? Copper abundance
real*4 e__Cu_Fe_ ! ([Sun]) ? rms uncertainty on [Cu/Fe]
integer*4 nZn ! Number of Zn lines
real*4 v_Zn_Fe_ ! ([Sun]) ? Zinc abundance
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! [Y/Fe], [Zr/Fe]I & II, [La/Fe], [Ce/Fe], [Pr/Fe],
[Nd/fe], [Sm/Fe], [Eu/Fe]
integer*4 nr__5
parameter (nr__5=14) ! Number of records
character*124 ar__5 ! Full-size record
integer*4 ID_5 ! [208/22172] Identification number
integer*4 nY ! Number of Y lines
real*4 v_Y_Fe_ ! ([Sun]) Yttrium abundance
real*4 e__Y_Fe_ ! ([Sun]) rms uncertainty on [Y/Fe]
integer*4 nZrI ! Number of ZrI lines
real*4 v_ZrI_Fe_ ! ([Sun]) ? Zirconium I abundance
real*4 e__ZrI_Fe_ ! ([Sun]) ? rms uncertainty on [ZrI/Fe]
integer*4 nZrII ! Number of ZrII lines
real*4 v_ZrII_Fe_ ! ([Sun]) Zirconium II abundance
real*4 e__ZrII_Fe_ ! ([Sun]) ? rms uncertainty on [ZrII/Fe]
integer*4 nLa ! Number of La lines
real*4 v_La_Fe_ ! ([Sun]) Lanthanum abundance
real*4 e__La_Fe_ ! ([Sun]) ? rms uncertainty on [La/Fe]
integer*4 nCe ! Number of Ce lines
real*4 v_Ce_Fe_ ! ([Sun]) Cerium abundance
real*4 e__Ce_Fe_ ! ([Sun]) ? rms uncertainty on [Ce/Fe]
integer*4 nPr ! Number of Pr lines
real*4 v_Pr_Fe_ ! ([Sun]) Praseodymium abundance
real*4 e__Pr_Fe_ ! ([Sun]) rms uncertainty on [Pr/Fe]
integer*4 nNd ! Number of Nd lines
real*4 v_Nd_Fe_ ! ([Sun]) Neodymium abundance
real*4 e__Nd_Fe_ ! ([Sun]) rms uncertainty on [Nd/Fe]
integer*4 nSm ! Number of Sm lines
real*4 v_Sm_Fe_ ! ([Sun]) Samarium abundance
real*4 e__Sm_Fe_ ! ([Sun]) ? rms uncertainty on [Sm/Fe]
integer*4 nEu ! Number of Eu lines
real*4 v_Eu_Fe_ ! ([Sun]) Europium abundance
real*4 e__Eu_Fe_ ! ([Sun]) ? rms uncertainty on [Eu/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! [Ba/Fe]
integer*4 nr__6
parameter (nr__6=82) ! Number of records
character*18 ar__6 ! Full-size record
integer*4 ID_6 ! [208/22172] Identification number
integer*4 nBa ! Number of Ba lines
real*4 v_Ba_Fe_ ! ([Sun]) ? barium abundance
real*4 e__Ba_Fe_ ! ([Sun]) ? rms uncertainty on [Ba/Fe]
C=============================================================================
C Loading file 'table2.dat' ! Positions, BV magnitudes and radial velocities
C Format for file interpretation
1 format(
+ I5,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F6.3,1X,F6.3,
+ 1X,A9,1X,F6.2,1X,F4.2,1X,F5.2,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,143
read(1,'(A80)')ar__
read(ar__,1)
+ ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,Vmag,setup,HRV,e_HRV,Dist,
+ flag
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,Vmag,setup,HRV,e_HRV,Dist,
+ flag
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Teff, logg, [A/H], Vt, [FeI/H], [FeII/H]
C Format for file interpretation
2 format(
+ I5,1X,I4,1X,F4.2,1X,F5.2,1X,F4.2,1X,I2,1X,F6.3,1X,F5.3,1X,I2,
+ 1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,82
read(1,'(A58)')ar__1
read(ar__1,2)
+ ID_1,Teff,logg,v_A_H_,vt,nFeI,v_FeI_H_,e__FeI_H_,nFeII,
+ v_FeII_H_,e__FeII_H_
if(ar__1(44:45) .EQ. '') nFeII = iNULL__
if(ar__1(47:52) .EQ. '') v_FeII_H_ = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1,Teff,logg,v_A_H_,vt,nFeI,v_FeI_H_,e__FeI_H_,nFeII,
+ v_FeII_H_,e__FeII_H_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe]
C Format for file interpretation
3 format(
+ I5,1X,I1,1X,A1,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,I1,1X,
+ F5.2,1X,F4.2,1X,I1,1X,A1,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,82
read(1,'(A59)')ar__2
read(ar__2,3)
+ ID_2,nO,l__O_Fe_,v_O_Fe_,e__O_Fe_,nNa,v_Na_Fe_,e__Na_Fe_,nMg,
+ v_Mg_Fe_,e__Mg_Fe_,nAl,l__Al_Fe_,v_Al_Fe_,e__Al_Fe_
if(ar__2(10:14) .EQ. '') v_O_Fe_ = rNULL__
if(ar__2(16:19) .EQ. '') e__O_Fe_ = rNULL__
if(ar__2(21:21) .EQ. '') nNa = iNULL__
if(ar__2(23:27) .EQ. '') v_Na_Fe_ = rNULL__
if(ar__2(29:32) .EQ. '') e__Na_Fe_ = rNULL__
if(ar__2(36:40) .EQ. '') v_Mg_Fe_ = rNULL__
if(ar__2(42:45) .EQ. '') e__Mg_Fe_ = rNULL__
if(ar__2(50:54) .EQ. '') v_Al_Fe_ = rNULL__
if(ar__2(56:59) .EQ. '') e__Al_Fe_ = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_2,nO,l__O_Fe_,v_O_Fe_,e__O_Fe_,nNa,v_Na_Fe_,e__Na_Fe_,nMg,
+ v_Mg_Fe_,e__Mg_Fe_,nAl,l__Al_Fe_,v_Al_Fe_,e__Al_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! [Si/Fe], [Ca/Fe], [Ti/Fe] I & II
C Format for file interpretation
4 format(
+ I5,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,F5.2,
+ 1X,F5.2,1X,I2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,82
read(1,'(A62)')ar__3
read(ar__3,4)
+ ID_3,nSi,v_Si_Fe_,e__Si_Fe_,nCa,v_Ca_Fe_,e__Ca_Fe_,nTi1,
+ v_TiI_Fe_,e__TiI_Fe_,nTi2,v_TiII_Fe_,e__TiII_Fe_
if(ar__3(10:14) .EQ. '') v_Si_Fe_ = rNULL__
if(ar__3(16:19) .EQ. '') e__Si_Fe_ = rNULL__
if(ar__3(30:33) .EQ. '') e__Ca_Fe_ = rNULL__
if(ar__3(38:42) .EQ. '') v_TiI_Fe_ = rNULL__
if(ar__3(44:48) .EQ. '') e__TiI_Fe_ = rNULL__
if(ar__3(53:57) .EQ. '') v_TiII_Fe_ = rNULL__
if(ar__3(59:62) .EQ. '') e__TiII_Fe_ = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_3,nSi,v_Si_Fe_,e__Si_Fe_,nCa,v_Ca_Fe_,e__Ca_Fe_,nTi1,
+ v_TiI_Fe_,e__TiI_Fe_,nTi2,v_TiII_Fe_,e__TiII_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! [Sc/Fe], [V/Fe], [Cr/Fe], [Mn/Fe], [Co/Fe],
* [Ni/Fe], [Cu/Fe], [Zn.Fe]
C Format for file interpretation
5 format(
+ I5,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,F5.2,
+ 1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I2,1X,
+ F5.2,1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,
+ I1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,82
read(1,'(A124)')ar__4
read(ar__4,5)
+ ID_4,nSc,v_Sc_Fe_,e__Sc_Fe_,nV,v_V_Fe_,e__V_Fe_,nCrI,
+ v_CrI_Fe_,e__CrI_Fe_,nCrII,v_CrII_Fe_,e__CrII_Fe_,nMn,
+ v_Mn_Fe_,e__Mn_Fe_,nCo,v_Co_Fe_,e__Co_Fe_,nNi,v_Ni_Fe_,
+ e__Ni_Fe_,nCu,v_Cu_Fe_,e__Cu_Fe_,nZn,v_Zn_Fe_
if(ar__4(16:19) .EQ. '') e__Sc_Fe_ = rNULL__
if(ar__4(24:28) .EQ. '') v_V_Fe_ = rNULL__
if(ar__4(30:33) .EQ. '') e__V_Fe_ = rNULL__
if(ar__4(38:42) .EQ. '') v_CrI_Fe_ = rNULL__
if(ar__4(44:47) .EQ. '') e__CrI_Fe_ = rNULL__
if(ar__4(52:56) .EQ. '') v_CrII_Fe_ = rNULL__
if(ar__4(58:61) .EQ. '') e__CrII_Fe_ = rNULL__
if(ar__4(66:70) .EQ. '') v_Mn_Fe_ = rNULL__
if(ar__4(72:75) .EQ. '') e__Mn_Fe_ = rNULL__
if(ar__4(80:84) .EQ. '') v_Co_Fe_ = rNULL__
if(ar__4(86:89) .EQ. '') e__Co_Fe_ = rNULL__
if(ar__4(94:98) .EQ. '') v_Ni_Fe_ = rNULL__
if(ar__4(100:103) .EQ. '') e__Ni_Fe_ = rNULL__
if(ar__4(107:111) .EQ. '') v_Cu_Fe_ = rNULL__
if(ar__4(113:116) .EQ. '') e__Cu_Fe_ = rNULL__
if(ar__4(120:124) .EQ. '') v_Zn_Fe_ = rNULL__
c ..............Just test output...........
write(6,5)
+ ID_4,nSc,v_Sc_Fe_,e__Sc_Fe_,nV,v_V_Fe_,e__V_Fe_,nCrI,
+ v_CrI_Fe_,e__CrI_Fe_,nCrII,v_CrII_Fe_,e__CrII_Fe_,nMn,
+ v_Mn_Fe_,e__Mn_Fe_,nCo,v_Co_Fe_,e__Co_Fe_,nNi,v_Ni_Fe_,
+ e__Ni_Fe_,nCu,v_Cu_Fe_,e__Cu_Fe_,nZn,v_Zn_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! [Y/Fe], [Zr/Fe]I & II, [La/Fe], [Ce/Fe], [Pr/Fe],
* [Nd/fe], [Sm/Fe], [Eu/Fe]
C Format for file interpretation
6 format(
+ I5,1X,I2,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,
+ 1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,I1,1X,
+ F5.2,1X,F4.2,1X,I2,1X,F5.2,1X,F4.2,1X,I1,1X,F5.2,1X,F4.2,1X,
+ I1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,14
read(1,'(A124)')ar__5
read(ar__5,6)
+ ID_5,nY,v_Y_Fe_,e__Y_Fe_,nZrI,v_ZrI_Fe_,e__ZrI_Fe_,nZrII,
+ v_ZrII_Fe_,e__ZrII_Fe_,nLa,v_La_Fe_,e__La_Fe_,nCe,v_Ce_Fe_,
+ e__Ce_Fe_,nPr,v_Pr_Fe_,e__Pr_Fe_,nNd,v_Nd_Fe_,e__Nd_Fe_,nSm,
+ v_Sm_Fe_,e__Sm_Fe_,nEu,v_Eu_Fe_,e__Eu_Fe_
if(ar__5(23:27) .EQ. '') v_ZrI_Fe_ = rNULL__
if(ar__5(29:32) .EQ. '') e__ZrI_Fe_ = rNULL__
if(ar__5(42:45) .EQ. '') e__ZrII_Fe_ = rNULL__
if(ar__5(55:58) .EQ. '') e__La_Fe_ = rNULL__
if(ar__5(68:71) .EQ. '') e__Ce_Fe_ = rNULL__
if(ar__5(108:111) .EQ. '') e__Sm_Fe_ = rNULL__
if(ar__5(121:124) .EQ. '') e__Eu_Fe_ = rNULL__
c ..............Just test output...........
write(6,6)
+ ID_5,nY,v_Y_Fe_,e__Y_Fe_,nZrI,v_ZrI_Fe_,e__ZrI_Fe_,nZrII,
+ v_ZrII_Fe_,e__ZrII_Fe_,nLa,v_La_Fe_,e__La_Fe_,nCe,v_Ce_Fe_,
+ e__Ce_Fe_,nPr,v_Pr_Fe_,e__Pr_Fe_,nNd,v_Nd_Fe_,e__Nd_Fe_,nSm,
+ v_Sm_Fe_,e__Sm_Fe_,nEu,v_Eu_Fe_,e__Eu_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! [Ba/Fe]
C Format for file interpretation
7 format(I5,1X,I1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,82
read(1,'(A18)')ar__6
read(ar__6,7)ID_6,nBa,v_Ba_Fe_,e__Ba_Fe_
if(ar__6(9:13) .EQ. '') v_Ba_Fe_ = rNULL__
if(ar__6(15:18) .EQ. '') e__Ba_Fe_ = rNULL__
c ..............Just test output...........
write(6,7)ID_6,nBa,v_Ba_Fe_,e__Ba_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end