Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/533/A69 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-19
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/533/A69 Spectroscopy of 124 RGB stars in NGC 1851 (Carretta+, 2011)
*================================================================================
*Multiple stellar populations in the globular cluster NGC 1851
* Carretta E., Lucatello S., Gratton R., Bragaglia A, D'Orazi V.
* <Astron. Astrophys. 533, A69 (2011)>
* =2011A&A...533A..69C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Identification, positions, B, V, K, RV, flag
integer*4 nr__
parameter (nr__=124) ! Number of records
character*68 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 ! 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) B magnitude
real*4 Vmag ! (mag) V magnitude
real*4 Kmag ! (mag) K magnitude (2MASS)
real*4 RV ! (km/s) Heliocentric radial velocity
character*9 Inst ! Intrument: UVES, Giraffe-HR11, or Giraffe-HR13
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Teff, logg, [A/H], Vt, [Fe/H]I, [Fe/H]II
integer*4 nr__1
parameter (nr__1=124) ! Number of records
character*58 ar__1 ! Full-size record
integer*4 ID_1 ! Identification number
integer*4 Teff ! (K) Effective temperature
real*4 logg ! ([cm/s2]) Gravity
real*4 v_A_H_ ! ([Sun]) Metallicity
real*4 vt ! (km/s) Microturbulent velocity {xi}
integer*4 o__Fe_H_1 ! Number of FeI lines
real*4 v_Fe_H_1 ! ([Sun]) Iron I abundance
real*4 e__Fe_H_1 ! ([Sun]) Error on [Fe/H]1
integer*4 o__Fe_H_2 ! ? Number of FeII lines
real*4 v_Fe_H_2 ! ([Sun]) ? Iron II abundance
real*4 e__Fe_H_2 ! ([Sun]) ? Error on [Fe/H]2
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe]
integer*4 nr__2
parameter (nr__2=124) ! Number of records
character*61 ar__2 ! Full-size record
integer*4 ID_2 ! Identification number
integer*4 o__O_Fe_ ! ? Number of O lines
real*4 v_O_Fe_ ! ([Sun]) ? Oxygen abundance
real*4 e__O_Fe_ ! ([Sun]) ? rms uncertainty on [O/Fe]
integer*4 o__Na_Fe_ ! ? 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 o__Mg_Fe_ ! ? 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 o__Al_Fe_ ! ? Number of Al lines
real*4 v_Al_Fe_ ! ([Sun]) ? Aluminum abundance
real*4 e__Al_Fe_ ! ([Sun]) ? rms uncertainty on [Al/Fe]
character*1 l__O_Fe_ ! [0,1] 0 = upper limit in [O/Fe]
character*1 l__Al_Fe_ ! [0,1]? 0 = upper limit in [Al/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! [Si/Fe], [Ca/Fe], [Ti/Fe] I & II
integer*4 nr__3
parameter (nr__3=124) ! Number of records
character*55 ar__3 ! Full-size record
integer*4 ID_3 ! Identification number
integer*4 o__Si_Fe_ ! ? 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 o__Ca_Fe_ ! ? 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 o__Ti_Fe_1 ! ? Number of Ti I lines
real*4 v_Ti_Fe_1 ! ([Sun]) ? Titanium I abundance
real*4 e__Ti_Fe_1 ! ([Sun]) ? rms uncertainty on [Si/Fe]
integer*4 o__Ti_Fe_2 ! ? Number of Ti II lines
real*4 v_Ti_Fe_2 ! ([Sun]) ? Titanium II abundance
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! [Sc/Fe], [V/Fe], [Cr/Fe], [Mn/Fe], [Co/Fe], [Ni/Fe],
[Cu/Fe]
integer*4 nr__4
parameter (nr__4=124) ! Number of records
character*97 ar__4 ! Full-size record
integer*4 ID_4 ! Identification number
integer*4 o__Sc_Fe_ ! 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 o__V_Fe_ ! ? 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 o__Cr_Fe_ ! ? Number of Cr lines
real*4 v_Cr_Fe_ ! ([Sun]) ? Cromium abundance
real*4 e__Cr_Fe_ ! ([Sun]) ? rms uncertainty on [Cr/Fe]
integer*4 o__Mn_Fe_ ! ? 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 o__Co_Fe_ ! ? 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 o__Ni_Fe_ ! ? 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 o__Cu_Fe_ ! ? Number of Cu lines
real*4 v_Cu_Fe_ ! ([Sun]) ? Copper abundance
real*4 e__Cu_Fe_ ! ([Sun]) ? rms uncertainty on [Cu/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! [Y/Fe], [Zr/Fe], [La/Fe], [Ce/Fe]. [Nd/Fe], [Eu/Fe],
[Dy/Fe]
integer*4 nr__5
parameter (nr__5=13) ! Number of records
character*96 ar__5 ! Full-size record
integer*4 ID_5 ! Identification number
integer*4 o__Y_Fe_ ! ? Number of Y II lines
real*4 v_Y_Fe_ ! ([Sun]) ? Yttrium abundance
real*4 e__Y_Fe_ ! ([Sun]) ? rms uncertainty on [Y/Fe]
integer*4 o__Zr_Fe_ ! ? Number of Zr II lines
real*4 v_Zr_Fe_ ! ([Sun]) ? Zirconium abundance
real*4 e__Zr_Fe_ ! ([Sun]) ? rms uncertainty on [Zr/Fe]
integer*4 o__La_Fe_ ! ? Number of La II lines
real*4 v_La_Fe_ ! ([Sun]) ? Lanthanum abundance
real*4 e__La_Fe_ ! ([Sun]) ? rms uncertainty on [La/Fe]
integer*4 o__Ce_Fe_ ! ? Number of Ce II lines
real*4 v_Ce_Fe_ ! ([Sun]) ? Cerium abundance
real*4 e__Ce_Fe_ ! ([Sun]) ? rms uncertainty on [Ce/Fe]
integer*4 o__Nd_Fe_ ! ? Number of Nd II lines
real*4 v_Nd_Fe_ ! ([Sun]) ? Neodimium abundance
real*4 e__Nd_Fe_ ! ([Sun]) ? rms uncertainty on [Nd/Fe]
integer*4 o__Eu_Fe_ ! ? Number of Eu II lines
real*4 v_Eu_Fe_ ! ([Sun]) ? Europium abundance
real*4 e__Eu_Fe_ ! ([Sun]) ? rms uncertainty on [Eu/Fe]
integer*4 o__Dy_Fe_ ! ? Number of Dy II lines
real*4 v_Dy_Fe_ ! ([Sun]) ? Dysprosium abundance
real*4 e__Dy_Fe_ ! ([Sun]) ? rms uncertainty on [Dy/Fe]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! [Ba/Fe]
integer*4 nr__6
parameter (nr__6=124) ! Number of records
character*18 ar__6 ! Full-size record
integer*4 ID_6 ! Identification number
integer*4 o__Ba_Fe_ ! ? Number of Ba II 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' ! Identification, positions, B, V, K, RV, flag
C Format for file interpretation
1 format(
+ I5,1X,I1,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F6.3,1X,F6.3,
+ 1X,F6.3,1X,F6.2,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,124
read(1,'(A68)')ar__
read(ar__,1)
+ ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Bmag,Vmag,Kmag,RV,Inst
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,Kmag,RV,Inst
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, [Fe/H]I, [Fe/H]II
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,124
read(1,'(A58)')ar__1
read(ar__1,2)
+ ID_1,Teff,logg,v_A_H_,vt,o__Fe_H_1,v_Fe_H_1,e__Fe_H_1,
+ o__Fe_H_2,v_Fe_H_2,e__Fe_H_2
if(ar__1(44:45) .EQ. '') o__Fe_H_2 = iNULL__
if(ar__1(47:52) .EQ. '') v_Fe_H_2 = rNULL__
if(ar__1(54:58) .EQ. '') e__Fe_H_2 = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1,Teff,logg,v_A_H_,vt,o__Fe_H_1,v_Fe_H_1,e__Fe_H_1,
+ o__Fe_H_2,v_Fe_H_2,e__Fe_H_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe]
C Format for file interpretation
3 format(
+ I5,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,A1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,124
read(1,'(A61)')ar__2
read(ar__2,3)
+ ID_2,o__O_Fe_,v_O_Fe_,e__O_Fe_,o__Na_Fe_,v_Na_Fe_,e__Na_Fe_,
+ o__Mg_Fe_,v_Mg_Fe_,e__Mg_Fe_,o__Al_Fe_,v_Al_Fe_,e__Al_Fe_,
+ l__O_Fe_,l__Al_Fe_
if(ar__2(7:7) .EQ. '') o__O_Fe_ = iNULL__
if(ar__2(9:13) .EQ. '') v_O_Fe_ = rNULL__
if(ar__2(15:18) .EQ. '') e__O_Fe_ = rNULL__
if(ar__2(20:20) .EQ. '') o__Na_Fe_ = iNULL__
if(ar__2(22:26) .EQ. '') v_Na_Fe_ = rNULL__
if(ar__2(28:31) .EQ. '') e__Na_Fe_ = rNULL__
if(ar__2(33:33) .EQ. '') o__Mg_Fe_ = iNULL__
if(ar__2(35:39) .EQ. '') v_Mg_Fe_ = rNULL__
if(ar__2(41:44) .EQ. '') e__Mg_Fe_ = rNULL__
if(ar__2(46:46) .EQ. '') o__Al_Fe_ = iNULL__
if(ar__2(48:52) .EQ. '') v_Al_Fe_ = rNULL__
if(ar__2(54:57) .EQ. '') e__Al_Fe_ = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_2,o__O_Fe_,v_O_Fe_,e__O_Fe_,o__Na_Fe_,v_Na_Fe_,e__Na_Fe_,
+ o__Mg_Fe_,v_Mg_Fe_,e__Mg_Fe_,o__Al_Fe_,v_Al_Fe_,e__Al_Fe_,
+ l__O_Fe_,l__Al_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.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,F4.2,1X,I1,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,124
read(1,'(A55)')ar__3
read(ar__3,4)
+ ID_3,o__Si_Fe_,v_Si_Fe_,e__Si_Fe_,o__Ca_Fe_,v_Ca_Fe_,
+ e__Ca_Fe_,o__Ti_Fe_1,v_Ti_Fe_1,e__Ti_Fe_1,o__Ti_Fe_2,v_Ti_Fe_2
if(ar__3(7:8) .EQ. '') o__Si_Fe_ = iNULL__
if(ar__3(10:14) .EQ. '') v_Si_Fe_ = rNULL__
if(ar__3(16:19) .EQ. '') e__Si_Fe_ = rNULL__
if(ar__3(21:22) .EQ. '') o__Ca_Fe_ = iNULL__
if(ar__3(24:28) .EQ. '') v_Ca_Fe_ = rNULL__
if(ar__3(30:33) .EQ. '') e__Ca_Fe_ = rNULL__
if(ar__3(35:36) .EQ. '') o__Ti_Fe_1 = iNULL__
if(ar__3(38:42) .EQ. '') v_Ti_Fe_1 = rNULL__
if(ar__3(44:47) .EQ. '') e__Ti_Fe_1 = rNULL__
if(ar__3(49:49) .EQ. '') o__Ti_Fe_2 = iNULL__
if(ar__3(51:55) .EQ. '') v_Ti_Fe_2 = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_3,o__Si_Fe_,v_Si_Fe_,e__Si_Fe_,o__Ca_Fe_,v_Ca_Fe_,
+ e__Ca_Fe_,o__Ti_Fe_1,v_Ti_Fe_1,e__Ti_Fe_1,o__Ti_Fe_2,v_Ti_Fe_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! [Sc/Fe], [V/Fe], [Cr/Fe], [Mn/Fe], [Co/Fe], [Ni/Fe],
* [Cu/Fe]
C Format for file interpretation
5 format(
+ I5,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)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,124
read(1,'(A97)')ar__4
read(ar__4,5)
+ ID_4,o__Sc_Fe_,v_Sc_Fe_,e__Sc_Fe_,o__V_Fe_,v_V_Fe_,e__V_Fe_,
+ o__Cr_Fe_,v_Cr_Fe_,e__Cr_Fe_,o__Mn_Fe_,v_Mn_Fe_,e__Mn_Fe_,
+ o__Co_Fe_,v_Co_Fe_,e__Co_Fe_,o__Ni_Fe_,v_Ni_Fe_,e__Ni_Fe_,
+ o__Cu_Fe_,v_Cu_Fe_,e__Cu_Fe_
if(ar__4(15:18) .EQ. '') e__Sc_Fe_ = rNULL__
if(ar__4(20:20) .EQ. '') o__V_Fe_ = iNULL__
if(ar__4(22:26) .EQ. '') v_V_Fe_ = rNULL__
if(ar__4(28:31) .EQ. '') e__V_Fe_ = rNULL__
if(ar__4(33:33) .EQ. '') o__Cr_Fe_ = iNULL__
if(ar__4(35:39) .EQ. '') v_Cr_Fe_ = rNULL__
if(ar__4(41:44) .EQ. '') e__Cr_Fe_ = rNULL__
if(ar__4(46:46) .EQ. '') o__Mn_Fe_ = iNULL__
if(ar__4(48:52) .EQ. '') v_Mn_Fe_ = rNULL__
if(ar__4(54:57) .EQ. '') e__Mn_Fe_ = rNULL__
if(ar__4(59:59) .EQ. '') o__Co_Fe_ = iNULL__
if(ar__4(61:65) .EQ. '') v_Co_Fe_ = rNULL__
if(ar__4(67:70) .EQ. '') e__Co_Fe_ = rNULL__
if(ar__4(72:73) .EQ. '') o__Ni_Fe_ = iNULL__
if(ar__4(75:79) .EQ. '') v_Ni_Fe_ = rNULL__
if(ar__4(81:84) .EQ. '') e__Ni_Fe_ = rNULL__
if(ar__4(86:86) .EQ. '') o__Cu_Fe_ = iNULL__
if(ar__4(88:92) .EQ. '') v_Cu_Fe_ = rNULL__
if(ar__4(94:97) .EQ. '') e__Cu_Fe_ = rNULL__
c ..............Just test output...........
write(6,5)
+ ID_4,o__Sc_Fe_,v_Sc_Fe_,e__Sc_Fe_,o__V_Fe_,v_V_Fe_,e__V_Fe_,
+ o__Cr_Fe_,v_Cr_Fe_,e__Cr_Fe_,o__Mn_Fe_,v_Mn_Fe_,e__Mn_Fe_,
+ o__Co_Fe_,v_Co_Fe_,e__Co_Fe_,o__Ni_Fe_,v_Ni_Fe_,e__Ni_Fe_,
+ o__Cu_Fe_,v_Cu_Fe_,e__Cu_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! [Y/Fe], [Zr/Fe], [La/Fe], [Ce/Fe]. [Nd/Fe], [Eu/Fe],
* [Dy/Fe]
C Format for file interpretation
6 format(
+ I5,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,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=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,13
read(1,'(A96)')ar__5
read(ar__5,6)
+ ID_5,o__Y_Fe_,v_Y_Fe_,e__Y_Fe_,o__Zr_Fe_,v_Zr_Fe_,e__Zr_Fe_,
+ o__La_Fe_,v_La_Fe_,e__La_Fe_,o__Ce_Fe_,v_Ce_Fe_,e__Ce_Fe_,
+ o__Nd_Fe_,v_Nd_Fe_,e__Nd_Fe_,o__Eu_Fe_,v_Eu_Fe_,e__Eu_Fe_,
+ o__Dy_Fe_,v_Dy_Fe_,e__Dy_Fe_
if(ar__5(7:7) .EQ. '') o__Y_Fe_ = iNULL__
if(ar__5(9:13) .EQ. '') v_Y_Fe_ = rNULL__
if(ar__5(15:18) .EQ. '') e__Y_Fe_ = rNULL__
if(ar__5(20:20) .EQ. '') o__Zr_Fe_ = iNULL__
if(ar__5(22:26) .EQ. '') v_Zr_Fe_ = rNULL__
if(ar__5(28:31) .EQ. '') e__Zr_Fe_ = rNULL__
if(ar__5(33:33) .EQ. '') o__La_Fe_ = iNULL__
if(ar__5(35:39) .EQ. '') v_La_Fe_ = rNULL__
if(ar__5(41:44) .EQ. '') e__La_Fe_ = rNULL__
if(ar__5(46:46) .EQ. '') o__Ce_Fe_ = iNULL__
if(ar__5(48:52) .EQ. '') v_Ce_Fe_ = rNULL__
if(ar__5(54:57) .EQ. '') e__Ce_Fe_ = rNULL__
if(ar__5(59:59) .EQ. '') o__Nd_Fe_ = iNULL__
if(ar__5(61:65) .EQ. '') v_Nd_Fe_ = rNULL__
if(ar__5(67:70) .EQ. '') e__Nd_Fe_ = rNULL__
if(ar__5(72:72) .EQ. '') o__Eu_Fe_ = iNULL__
if(ar__5(74:78) .EQ. '') v_Eu_Fe_ = rNULL__
if(ar__5(80:83) .EQ. '') e__Eu_Fe_ = rNULL__
if(ar__5(85:85) .EQ. '') o__Dy_Fe_ = iNULL__
if(ar__5(87:91) .EQ. '') v_Dy_Fe_ = rNULL__
if(ar__5(93:96) .EQ. '') e__Dy_Fe_ = rNULL__
c ..............Just test output...........
write(6,6)
+ ID_5,o__Y_Fe_,v_Y_Fe_,e__Y_Fe_,o__Zr_Fe_,v_Zr_Fe_,e__Zr_Fe_,
+ o__La_Fe_,v_La_Fe_,e__La_Fe_,o__Ce_Fe_,v_Ce_Fe_,e__Ce_Fe_,
+ o__Nd_Fe_,v_Nd_Fe_,e__Nd_Fe_,o__Eu_Fe_,v_Eu_Fe_,e__Eu_Fe_,
+ o__Dy_Fe_,v_Dy_Fe_,e__Dy_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.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=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,124
read(1,'(A18)')ar__6
read(ar__6,7)ID_6,o__Ba_Fe_,v_Ba_Fe_,e__Ba_Fe_
if(ar__6(7:7) .EQ. '') o__Ba_Fe_ = iNULL__
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,o__Ba_Fe_,v_Ba_Fe_,e__Ba_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end