Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/578/A116 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-15
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 (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 ID (nr__) ! [208/22172] Identification number
integer*4 RAh (nr__) ! (h) Right ascension (J2000.0)
integer*4 RAm (nr__) ! (min) Right ascension (J2000.0)
real*4 RAs (nr__) ! (s) Right ascension (J2000.0)
character*1 DE_ (nr__) ! Declination sign (J2000.0)
integer*4 DEd (nr__) ! (deg) Declination (J2000.0)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000.0)
real*4 DEs (nr__) ! (arcsec) Declination (J2000.0)
real*4 Bmag (nr__) ! (mag) [14.7/17.9] B magnitude
real*4 Vmag (nr__) ! (mag) [12.8/16.8] V magnitude
character*9 setup (nr__) ! Giraffe setup (HR11,HR13) and/or UVES
real*4 HRV (nr__) ! (km/s) [-94/260] Heliocentric radial velocity
real*4 e_HRV (nr__) ! (km/s) [0/6.6] error on HRV
real*4 Dist (nr__) ! (arcmin) [0.99/12.9] Distance from centre of NGC 6093
character*6 flag (nr__) ! [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 (nr__1) ! [208/22172] Identification number
integer*4 Teff (nr__1) ! (K) [4046/5138] Effective temperature
real*4 logg (nr__1) ! ([cm/s2]) [0.42/2.63] Surface gravity
real*4 v_A_H_ (nr__1) ! ([Sun]) [-1.9/-1.7] Metallicity
real*4 vt (nr__1) ! (km/s) [0.65/3.57] Microturbulent velocity
integer*4 nFeI (nr__1) ! [3/76] Number of FeI lines
real*4 v_FeI_H_ (nr__1) ! ([Sun]) [-1.9/-1.7] Iron I abundance
real*4 e__FeI_H_ (nr__1) ! ([Sun]) [0/0.21] Error on [FeI/H]
integer*4 nFeII (nr__1) ! [1/10]? Number of FeII lines
real*4 v_FeII_H_ (nr__1) ! ([Sun]) [-1.9/-1.6]? Iron II abundance
real*4 e__FeII_H_ (nr__1) ! ([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 (nr__2) ! [208/22172] Identification number
integer*4 nO (nr__2) ! number of O lines
character*1 l__O_Fe_ (nr__2) ! Limit flag in O (2)
real*4 v_O_Fe_ (nr__2) ! ([Sun]) ? Oxygen abundance
real*4 e__O_Fe_ (nr__2) ! ([Sun]) ? rms uncertainty on [O/Fe]
integer*4 nNa (nr__2) ! ? Number of Na lines
real*4 v_Na_Fe_ (nr__2) ! ([Sun]) ? Sodium abundance
real*4 e__Na_Fe_ (nr__2) ! ([Sun]) ? rms uncertainty on [Na/Fe]
integer*4 nMg (nr__2) ! Number of Mg lines
real*4 v_Mg_Fe_ (nr__2) ! ([Sun]) ? Magnesium abundance
real*4 e__Mg_Fe_ (nr__2) ! ([Sun]) ? rms uncertainty on [Mg/Fe]
integer*4 nAl (nr__2) ! Number of Al lines
character*1 l__Al_Fe_ (nr__2) ! Limit flag for Al (2)
real*4 v_Al_Fe_ (nr__2) ! ([Sun]) ? Aluminum abundance
real*4 e__Al_Fe_ (nr__2) ! ([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 (nr__3) ! [208/22172] Identification number
integer*4 nSi (nr__3) ! Number of Si lines
real*4 v_Si_Fe_ (nr__3) ! ([Sun]) ? Silicon abundance
real*4 e__Si_Fe_ (nr__3) ! ([Sun]) ? rms uncertainty on [Si/Fe]
integer*4 nCa (nr__3) ! Number of Ca lines
real*4 v_Ca_Fe_ (nr__3) ! ([Sun]) Calcium abundance
real*4 e__Ca_Fe_ (nr__3) ! ([Sun]) ? rms uncertainty on [ca/Fe]
integer*4 nTi1 (nr__3) ! Number of Ti I lines
real*4 v_TiI_Fe_ (nr__3) ! ([Sun]) ? Titanium I abundance
real*4 e__TiI_Fe_ (nr__3) ! ([Sun]) ? rms uncertainty on [TiI/Fe]
integer*4 nTi2 (nr__3) ! Number of Ti II lines
real*4 v_TiII_Fe_ (nr__3) ! ([Sun]) ? Titanium II abundance
real*4 e__TiII_Fe_(nr__3) ! ([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 (nr__4) ! [208/22172] Identification number
integer*4 nSc (nr__4) ! Number of Sc lines
real*4 v_Sc_Fe_ (nr__4) ! ([Sun]) Scandium abundance
real*4 e__Sc_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [Sc/Fe]
integer*4 nV (nr__4) ! Number of V lines
real*4 v_V_Fe_ (nr__4) ! ([Sun]) ? Vanadium abundance
real*4 e__V_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [V/Fe]
integer*4 nCrI (nr__4) ! Number of CrI lines
real*4 v_CrI_Fe_ (nr__4) ! ([Sun]) ? Cromium abundance
real*4 e__CrI_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [CrI/Fe]
integer*4 nCrII (nr__4) ! Number of CrII lines
real*4 v_CrII_Fe_ (nr__4) ! ([Sun]) ? Cromium abundance
real*4 e__CrII_Fe_(nr__4) ! ([Sun]) ? rms uncertainty on [CrII/Fe]
integer*4 nMn (nr__4) ! Number of Mn lines
real*4 v_Mn_Fe_ (nr__4) ! ([Sun]) ? Manganese abundance
real*4 e__Mn_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [Mn/Fe]
integer*4 nCo (nr__4) ! Number of Co lines
real*4 v_Co_Fe_ (nr__4) ! ([Sun]) ? Cobalt abundance
real*4 e__Co_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [Co/Fe]
integer*4 nNi (nr__4) ! Number of Ni lines
real*4 v_Ni_Fe_ (nr__4) ! ([Sun]) ? Nichel abundance
real*4 e__Ni_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [Ni/Fe]
integer*4 nCu (nr__4) ! Number of Cu lines
real*4 v_Cu_Fe_ (nr__4) ! ([Sun]) ? Copper abundance
real*4 e__Cu_Fe_ (nr__4) ! ([Sun]) ? rms uncertainty on [Cu/Fe]
integer*4 nZn (nr__4) ! Number of Zn lines
real*4 v_Zn_Fe_ (nr__4) ! ([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 (nr__5) ! [208/22172] Identification number
integer*4 nY (nr__5) ! Number of Y lines
real*4 v_Y_Fe_ (nr__5) ! ([Sun]) Yttrium abundance
real*4 e__Y_Fe_ (nr__5) ! ([Sun]) rms uncertainty on [Y/Fe]
integer*4 nZrI (nr__5) ! Number of ZrI lines
real*4 v_ZrI_Fe_ (nr__5) ! ([Sun]) ? Zirconium I abundance
real*4 e__ZrI_Fe_ (nr__5) ! ([Sun]) ? rms uncertainty on [ZrI/Fe]
integer*4 nZrII (nr__5) ! Number of ZrII lines
real*4 v_ZrII_Fe_ (nr__5) ! ([Sun]) Zirconium II abundance
real*4 e__ZrII_Fe_(nr__5) ! ([Sun]) ? rms uncertainty on [ZrII/Fe]
integer*4 nLa (nr__5) ! Number of La lines
real*4 v_La_Fe_ (nr__5) ! ([Sun]) Lanthanum abundance
real*4 e__La_Fe_ (nr__5) ! ([Sun]) ? rms uncertainty on [La/Fe]
integer*4 nCe (nr__5) ! Number of Ce lines
real*4 v_Ce_Fe_ (nr__5) ! ([Sun]) Cerium abundance
real*4 e__Ce_Fe_ (nr__5) ! ([Sun]) ? rms uncertainty on [Ce/Fe]
integer*4 nPr (nr__5) ! Number of Pr lines
real*4 v_Pr_Fe_ (nr__5) ! ([Sun]) Praseodymium abundance
real*4 e__Pr_Fe_ (nr__5) ! ([Sun]) rms uncertainty on [Pr/Fe]
integer*4 nNd (nr__5) ! Number of Nd lines
real*4 v_Nd_Fe_ (nr__5) ! ([Sun]) Neodymium abundance
real*4 e__Nd_Fe_ (nr__5) ! ([Sun]) rms uncertainty on [Nd/Fe]
integer*4 nSm (nr__5) ! Number of Sm lines
real*4 v_Sm_Fe_ (nr__5) ! ([Sun]) Samarium abundance
real*4 e__Sm_Fe_ (nr__5) ! ([Sun]) ? rms uncertainty on [Sm/Fe]
integer*4 nEu (nr__5) ! Number of Eu lines
real*4 v_Eu_Fe_ (nr__5) ! ([Sun]) Europium abundance
real*4 e__Eu_Fe_ (nr__5) ! ([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 (nr__6) ! [208/22172] Identification number
integer*4 nBa (nr__6) ! Number of Ba lines
real*4 v_Ba_Fe_ (nr__6) ! ([Sun]) ? barium abundance
real*4 e__Ba_Fe_ (nr__6) ! ([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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),Bmag(i__),Vmag(i__),setup(i__),HRV(i__),e_HRV(i__),
+ Dist(i__),flag(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ ID(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),Bmag(i__),Vmag(i__),setup(i__),HRV(i__),e_HRV(i__),
+ Dist(i__),flag(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),Teff(i__),logg(i__),v_A_H_(i__),vt(i__),nFeI(i__),
+ v_FeI_H_(i__),e__FeI_H_(i__),nFeII(i__),v_FeII_H_(i__),
+ e__FeII_H_(i__)
if(ar__1(44:45) .EQ. '') nFeII(i__) = iNULL__
if(ar__1(47:52) .EQ. '') v_FeII_H_(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),Teff(i__),logg(i__),v_A_H_(i__),vt(i__),nFeI(i__),
+ v_FeI_H_(i__),e__FeI_H_(i__),nFeII(i__),v_FeII_H_(i__),
+ e__FeII_H_(i__)
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(i__),nO(i__),l__O_Fe_(i__),v_O_Fe_(i__),e__O_Fe_(i__),
+ nNa(i__),v_Na_Fe_(i__),e__Na_Fe_(i__),nMg(i__),v_Mg_Fe_(i__),
+ e__Mg_Fe_(i__),nAl(i__),l__Al_Fe_(i__),v_Al_Fe_(i__),
+ e__Al_Fe_(i__)
if(ar__2(10:14) .EQ. '') v_O_Fe_(i__) = rNULL__
if(ar__2(16:19) .EQ. '') e__O_Fe_(i__) = rNULL__
if(ar__2(21:21) .EQ. '') nNa(i__) = iNULL__
if(ar__2(23:27) .EQ. '') v_Na_Fe_(i__) = rNULL__
if(ar__2(29:32) .EQ. '') e__Na_Fe_(i__) = rNULL__
if(ar__2(36:40) .EQ. '') v_Mg_Fe_(i__) = rNULL__
if(ar__2(42:45) .EQ. '') e__Mg_Fe_(i__) = rNULL__
if(ar__2(50:54) .EQ. '') v_Al_Fe_(i__) = rNULL__
if(ar__2(56:59) .EQ. '') e__Al_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_2(i__),nO(i__),l__O_Fe_(i__),v_O_Fe_(i__),e__O_Fe_(i__),
+ nNa(i__),v_Na_Fe_(i__),e__Na_Fe_(i__),nMg(i__),v_Mg_Fe_(i__),
+ e__Mg_Fe_(i__),nAl(i__),l__Al_Fe_(i__),v_Al_Fe_(i__),
+ e__Al_Fe_(i__)
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(i__),nSi(i__),v_Si_Fe_(i__),e__Si_Fe_(i__),nCa(i__),
+ v_Ca_Fe_(i__),e__Ca_Fe_(i__),nTi1(i__),v_TiI_Fe_(i__),
+ e__TiI_Fe_(i__),nTi2(i__),v_TiII_Fe_(i__),e__TiII_Fe_(i__)
if(ar__3(10:14) .EQ. '') v_Si_Fe_(i__) = rNULL__
if(ar__3(16:19) .EQ. '') e__Si_Fe_(i__) = rNULL__
if(ar__3(30:33) .EQ. '') e__Ca_Fe_(i__) = rNULL__
if(ar__3(38:42) .EQ. '') v_TiI_Fe_(i__) = rNULL__
if(ar__3(44:48) .EQ. '') e__TiI_Fe_(i__) = rNULL__
if(ar__3(53:57) .EQ. '') v_TiII_Fe_(i__) = rNULL__
if(ar__3(59:62) .EQ. '') e__TiII_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_3(i__),nSi(i__),v_Si_Fe_(i__),e__Si_Fe_(i__),nCa(i__),
+ v_Ca_Fe_(i__),e__Ca_Fe_(i__),nTi1(i__),v_TiI_Fe_(i__),
+ e__TiI_Fe_(i__),nTi2(i__),v_TiII_Fe_(i__),e__TiII_Fe_(i__)
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(i__),nSc(i__),v_Sc_Fe_(i__),e__Sc_Fe_(i__),nV(i__),
+ v_V_Fe_(i__),e__V_Fe_(i__),nCrI(i__),v_CrI_Fe_(i__),
+ e__CrI_Fe_(i__),nCrII(i__),v_CrII_Fe_(i__),e__CrII_Fe_(i__),
+ nMn(i__),v_Mn_Fe_(i__),e__Mn_Fe_(i__),nCo(i__),v_Co_Fe_(i__),
+ e__Co_Fe_(i__),nNi(i__),v_Ni_Fe_(i__),e__Ni_Fe_(i__),nCu(i__),
+ v_Cu_Fe_(i__),e__Cu_Fe_(i__),nZn(i__),v_Zn_Fe_(i__)
if(ar__4(16:19) .EQ. '') e__Sc_Fe_(i__) = rNULL__
if(ar__4(24:28) .EQ. '') v_V_Fe_(i__) = rNULL__
if(ar__4(30:33) .EQ. '') e__V_Fe_(i__) = rNULL__
if(ar__4(38:42) .EQ. '') v_CrI_Fe_(i__) = rNULL__
if(ar__4(44:47) .EQ. '') e__CrI_Fe_(i__) = rNULL__
if(ar__4(52:56) .EQ. '') v_CrII_Fe_(i__) = rNULL__
if(ar__4(58:61) .EQ. '') e__CrII_Fe_(i__) = rNULL__
if(ar__4(66:70) .EQ. '') v_Mn_Fe_(i__) = rNULL__
if(ar__4(72:75) .EQ. '') e__Mn_Fe_(i__) = rNULL__
if(ar__4(80:84) .EQ. '') v_Co_Fe_(i__) = rNULL__
if(ar__4(86:89) .EQ. '') e__Co_Fe_(i__) = rNULL__
if(ar__4(94:98) .EQ. '') v_Ni_Fe_(i__) = rNULL__
if(ar__4(100:103) .EQ. '') e__Ni_Fe_(i__) = rNULL__
if(ar__4(107:111) .EQ. '') v_Cu_Fe_(i__) = rNULL__
if(ar__4(113:116) .EQ. '') e__Cu_Fe_(i__) = rNULL__
if(ar__4(120:124) .EQ. '') v_Zn_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ ID_4(i__),nSc(i__),v_Sc_Fe_(i__),e__Sc_Fe_(i__),nV(i__),
+ v_V_Fe_(i__),e__V_Fe_(i__),nCrI(i__),v_CrI_Fe_(i__),
+ e__CrI_Fe_(i__),nCrII(i__),v_CrII_Fe_(i__),e__CrII_Fe_(i__),
+ nMn(i__),v_Mn_Fe_(i__),e__Mn_Fe_(i__),nCo(i__),v_Co_Fe_(i__),
+ e__Co_Fe_(i__),nNi(i__),v_Ni_Fe_(i__),e__Ni_Fe_(i__),nCu(i__),
+ v_Cu_Fe_(i__),e__Cu_Fe_(i__),nZn(i__),v_Zn_Fe_(i__)
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(i__),nY(i__),v_Y_Fe_(i__),e__Y_Fe_(i__),nZrI(i__),
+ v_ZrI_Fe_(i__),e__ZrI_Fe_(i__),nZrII(i__),v_ZrII_Fe_(i__),
+ e__ZrII_Fe_(i__),nLa(i__),v_La_Fe_(i__),e__La_Fe_(i__),
+ nCe(i__),v_Ce_Fe_(i__),e__Ce_Fe_(i__),nPr(i__),v_Pr_Fe_(i__),
+ e__Pr_Fe_(i__),nNd(i__),v_Nd_Fe_(i__),e__Nd_Fe_(i__),nSm(i__),
+ v_Sm_Fe_(i__),e__Sm_Fe_(i__),nEu(i__),v_Eu_Fe_(i__),
+ e__Eu_Fe_(i__)
if(ar__5(23:27) .EQ. '') v_ZrI_Fe_(i__) = rNULL__
if(ar__5(29:32) .EQ. '') e__ZrI_Fe_(i__) = rNULL__
if(ar__5(42:45) .EQ. '') e__ZrII_Fe_(i__) = rNULL__
if(ar__5(55:58) .EQ. '') e__La_Fe_(i__) = rNULL__
if(ar__5(68:71) .EQ. '') e__Ce_Fe_(i__) = rNULL__
if(ar__5(108:111) .EQ. '') e__Sm_Fe_(i__) = rNULL__
if(ar__5(121:124) .EQ. '') e__Eu_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ ID_5(i__),nY(i__),v_Y_Fe_(i__),e__Y_Fe_(i__),nZrI(i__),
+ v_ZrI_Fe_(i__),e__ZrI_Fe_(i__),nZrII(i__),v_ZrII_Fe_(i__),
+ e__ZrII_Fe_(i__),nLa(i__),v_La_Fe_(i__),e__La_Fe_(i__),
+ nCe(i__),v_Ce_Fe_(i__),e__Ce_Fe_(i__),nPr(i__),v_Pr_Fe_(i__),
+ e__Pr_Fe_(i__),nNd(i__),v_Nd_Fe_(i__),e__Nd_Fe_(i__),nSm(i__),
+ v_Sm_Fe_(i__),e__Sm_Fe_(i__),nEu(i__),v_Eu_Fe_(i__),
+ e__Eu_Fe_(i__)
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(i__),nBa(i__),v_Ba_Fe_(i__),e__Ba_Fe_(i__)
if(ar__6(9:13) .EQ. '') v_Ba_Fe_(i__) = rNULL__
if(ar__6(15:18) .EQ. '') e__Ba_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,7)ID_6(i__),nBa(i__),v_Ba_Fe_(i__),e__Ba_Fe_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end