Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/512/A41 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 2024-Apr-20
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/512/A41 Abundances of microlensed stars in the Bulge (Bensby+, 2010)
*================================================================================
*Chemical evolution of the Galactic bulge as traced by microlensed dwarfs and
*subgiant stars.
*II. Ages, metallicities, detailed elemental abundances, and connections to the
*Galactic thick disc.
* Bensby T., Feltzing S., Johnson J.A., Gould A., Aden D., Asplund M.,
* Melendez J., Gal-Yam A., Lucatello S., Sana H., Sumi T., Miyake N.,
* Suzuki N., Han C., Bond I., Udalski A.
* <Astron. Astrophys. 512, A41 (2010)>
* =2010A&A...512A..41B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Summary of the so far 15 dwarf star microlensing
events in the Bulge that have been observed
with high-resolution spectrographs
integer*4 nr__
parameter (nr__=15) ! Number of records
character*111 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*19 Name (nr__) ! Star name
character*1 n_Name (nr__) ! [*] indicates the star MACHO 109.20893.3423
character*1 EW (nr__) ! [a] a for stars with data in table4
character*1 Ab (nr__) ! [+] + for stars with data in table5
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
real*4 RAs (nr__) ! (s) Right ascension (J2000)
character*1 DE_ (nr__) ! Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Declination (J2000)
real*4 GLON (nr__) ! (deg) [-90/90] Galactic longitude
real*4 GLAT (nr__) ! (deg) Galactic latitude
real*4 TE (nr__) ! (d) Duration of the event
real*8 Tmax (nr__) ! (d) Time when maximum magnification occured
* (HJD-2450000)
integer*4 Amax (nr__) ! Maximum magnification
real*8 Tobs (nr__) ! (d) Time when event was observed with
* high-resolution spectrograph (MJD)
integer*4 Texp (nr__) ! (s) Exposure time
integer*4 S_N (nr__) ! Measured signal-to-noise ratio per pixel
* at ~6400{AA}
character*5 Spec (nr__) ! Spectrograph used
integer*4 R (nr__) ! Spectral resolution
character*4 Ref (nr__) ! Reference where the star first appeared (2)
*Note (2): References as follows:
* TW = This work
* B09a = Bensby et al. (2009ApJ...699L.174B)
* B09b = Bensby et al. (2009, Cat. <J/A+A/499/737)
* J07 = Johnson et al. (2007, Cat. <J/ApJ/655/L33>)
* J08 = Johnson et al. (2008ApJ...685..508J)
* C08 = Cohen et al. (2008ApJ...682.1029C)
* C09 = Cohen et al. (2009ApJ...699...66C)
* C03 = Cavallo et al. (2003, SPIE Conf. Ser., 4834, 66)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Equivalent widths and absolute abundances for
individual lines (for 8 stars)
integer*4 nr__1
parameter (nr__1=391) ! Number of records
character*169 ar__1 ! Full-size record
character*2 El (nr__1) ! Element
integer*4 Ion (nr__1) ! Ionisation (1=neutral)
real*8 lambda (nr__1) ! (0.1nm) Wavelength
real*4 LEP (nr__1) ! (eV) Lower Excitation Potential
real*4 EWsun (nr__1) ! (0.1pm) Solar equivalent width of the element
real*4 A_Sun (nr__1) ! ([---]) Solar absolute abundance (log scale, A(H)=12)
integer*4 f_A_Sun (nr__1) ! [0,1] 0 if line measured and used,
* 1 if average abundance from other lines used
real*4 EW_076S (nr__1) ! (0.1pm) ? Equivalent width OGLE-2009-BLG-076S
real*4 A_076S (nr__1) ! ([---]) ? Absolute abundance OGLE-2009-BLG-076S
real*4 EW_493S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-493S
real*4 A_493S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-493S
real*4 EW_133S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-133S
real*4 A_133S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-133S
real*4 EW_475S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-475S
real*4 A_475S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-475S
real*4 EW_022S (nr__1) ! (0.1pm) ? Equivalent width MACHO-1999-BLG-022S
real*4 A_022S (nr__1) ! ([---]) ? Absolute abundance MACHO-1999-BLG-022S
real*4 EW_489S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-489S
real*4 A_489S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-489S
real*4 EW_456S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-456S
real*4 A_456S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-456S
real*4 EW_259S (nr__1) ! (0.1pm) ? Equivalent width MOA-2009-BLG-259S
real*4 A_259S (nr__1) ! ([---]) ? Absolute abundance MOA-2009-BLG-259S
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Abundances and uncertainties for each star
integer*4 nr__2
parameter (nr__2=13) ! Number of records
character*334 ar__2 ! Full-size record
character*19 Name_1 (nr__2) ! Star name
real*4 v_Fe_H_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_O_Fe_ (nr__2) ! ([Sun]) ? Normalised abundance (not NLTE corrected)
real*4 v_Na_Fe_ (nr__2) ! ([Sun]) ? Normalised abundance
real*4 v_Mg_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Al_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Si_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Ca_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Ti_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Cr_Fe_ (nr__2) ! ([Sun]) ? Normalised abundance
real*4 v_Ni_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 v_Zn_Fe_ (nr__2) ! ([Sun]) ? Normalised abundance
real*4 v_Y_Fe_ (nr__2) ! ([Sun]) ? Normalised abundance
real*4 v_Ba_Fe_ (nr__2) ! ([Sun]) Normalised abundance
real*4 e__Fe_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__O_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Na_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Mg_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Al_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Si_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ca_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ti_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Cr_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ni_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Zn_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Y_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ba_Fe_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Fe_H_r (nr__2) ! ([Sun]) ? Error in abundance ratio (repeated)
real*4 e__O_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Na_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Mg_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Al_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Si_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ca_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ti_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Cr_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ni_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Zn_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Y_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
real*4 e__Ba_H_ (nr__2) ! ([Sun]) ? Error in abundance ratio
integer*4 n_Fe (nr__2) ! ? Number of Fe lines used
integer*4 n_O (nr__2) ! ? Number of O lines used
integer*4 n_Na (nr__2) ! ? Number of Na lines used
integer*4 n_Mg (nr__2) ! ? Number of Mg lines used
integer*4 n_Al (nr__2) ! ? Number of Al lines used
integer*4 n_Si (nr__2) ! ? Number of Si lines used
integer*4 n_Ca (nr__2) ! ? Number of Ca lines used
integer*4 n_Ti (nr__2) ! ? Number of Ti lines used
integer*4 n_Cr (nr__2) ! ? Number of Cr lines used
integer*4 n_Ni (nr__2) ! ? Number of Ni lines used
integer*4 n_Zn (nr__2) ! ? Number of Zn lines used
integer*4 n_Y (nr__2) ! ? Number of Y lines used
integer*4 n_Ba (nr__2) ! ? Number of Ba lines used
C=============================================================================
C Loading file 'table1.dat' ! Summary of the so far 15 dwarf star microlensing
* events in the Bulge that have been observed
* with high-resolution spectrographs
C Format for file interpretation
1 format(
+ A19,A1,A1,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ F5.2,1X,F5.2,1X,F5.1,1X,F7.2,1X,I4,1X,F8.3,1X,I5,1X,I2,1X,A5,
+ 1X,I5,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,15
read(1,'(A111)')ar__
read(ar__,1)
+ Name(i__),n_Name(i__),EW(i__),Ab(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),GLON(i__),
+ GLAT(i__),TE(i__),Tmax(i__),Amax(i__),Tobs(i__),Texp(i__),
+ S_N(i__),Spec(i__),R(i__),Ref(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)
+ Name(i__),n_Name(i__),EW(i__),Ab(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),GLON(i__),
+ GLAT(i__),TE(i__),Tmax(i__),Amax(i__),Tobs(i__),Texp(i__),
+ S_N(i__),Spec(i__),R(i__),Ref(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 'table4.dat' ! Equivalent widths and absolute abundances for
* individual lines (for 8 stars)
C Format for file interpretation
2 format(
+ 1X,A2,1X,I1,5X,F8.3,3X,F4.2,3X,F5.1,3X,F4.2,8X,I1,3X,F5.1,3X,
+ F4.2,3X,F5.1,3X,F4.2,3X,F5.1,3X,F4.2,3X,F5.1,3X,F4.2,3X,F5.1,
+ 3X,F4.2,3X,F5.1,3X,F4.2,3X,F5.1,3X,F4.2,3X,F5.1,3X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,391
read(1,'(A169)')ar__1
read(ar__1,2)
+ El(i__),Ion(i__),lambda(i__),LEP(i__),EWsun(i__),A_Sun(i__),
+ f_A_Sun(i__),EW_076S(i__),A_076S(i__),EW_493S(i__),
+ A_493S(i__),EW_133S(i__),A_133S(i__),EW_475S(i__),A_475S(i__),
+ EW_022S(i__),A_022S(i__),EW_489S(i__),A_489S(i__),
+ EW_456S(i__),A_456S(i__),EW_259S(i__),A_259S(i__)
if(ar__1(53:57) .EQ. '') EW_076S(i__) = rNULL__
if(ar__1(61:64) .EQ. '') A_076S(i__) = rNULL__
if(ar__1(68:72) .EQ. '') EW_493S(i__) = rNULL__
if(ar__1(76:79) .EQ. '') A_493S(i__) = rNULL__
if(ar__1(83:87) .EQ. '') EW_133S(i__) = rNULL__
if(ar__1(91:94) .EQ. '') A_133S(i__) = rNULL__
if(ar__1(98:102) .EQ. '') EW_475S(i__) = rNULL__
if(ar__1(106:109) .EQ. '') A_475S(i__) = rNULL__
if(ar__1(113:117) .EQ. '') EW_022S(i__) = rNULL__
if(ar__1(121:124) .EQ. '') A_022S(i__) = rNULL__
if(ar__1(128:132) .EQ. '') EW_489S(i__) = rNULL__
if(ar__1(136:139) .EQ. '') A_489S(i__) = rNULL__
if(ar__1(143:147) .EQ. '') EW_456S(i__) = rNULL__
if(ar__1(151:154) .EQ. '') A_456S(i__) = rNULL__
if(ar__1(158:162) .EQ. '') EW_259S(i__) = rNULL__
if(ar__1(166:169) .EQ. '') A_259S(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ El(i__),Ion(i__),lambda(i__),LEP(i__),EWsun(i__),A_Sun(i__),
+ f_A_Sun(i__),EW_076S(i__),A_076S(i__),EW_493S(i__),
+ A_493S(i__),EW_133S(i__),A_133S(i__),EW_475S(i__),A_475S(i__),
+ EW_022S(i__),A_022S(i__),EW_489S(i__),A_489S(i__),
+ EW_456S(i__),A_456S(i__),EW_259S(i__),A_259S(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Abundances and uncertainties for each star
C Format for file interpretation
3 format(
+ A19,2X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,6X,F4.2,2X,
+ F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,
+ 2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,6X,F4.2,2X,F4.2,2X,F4.2,2X,
+ F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,
+ 2X,F4.2,2X,F4.2,4X,I3,2X,I3,2X,I3,2X,I3,2X,I3,3X,I3,3X,I3,3X,
+ I3,3X,I3,3X,I3,2X,I3,2X,I3,2X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,13
read(1,'(A334)')ar__2
read(ar__2,3)
+ Name_1(i__),v_Fe_H_(i__),v_O_Fe_(i__),v_Na_Fe_(i__),
+ v_Mg_Fe_(i__),v_Al_Fe_(i__),v_Si_Fe_(i__),v_Ca_Fe_(i__),
+ v_Ti_Fe_(i__),v_Cr_Fe_(i__),v_Ni_Fe_(i__),v_Zn_Fe_(i__),
+ v_Y_Fe_(i__),v_Ba_Fe_(i__),e__Fe_H_(i__),e__O_Fe_(i__),
+ e__Na_Fe_(i__),e__Mg_Fe_(i__),e__Al_Fe_(i__),e__Si_Fe_(i__),
+ e__Ca_Fe_(i__),e__Ti_Fe_(i__),e__Cr_Fe_(i__),e__Ni_Fe_(i__),
+ e__Zn_Fe_(i__),e__Y_Fe_(i__),e__Ba_Fe_(i__),e__Fe_H_r(i__),
+ e__O_H_(i__),e__Na_H_(i__),e__Mg_H_(i__),e__Al_H_(i__),
+ e__Si_H_(i__),e__Ca_H_(i__),e__Ti_H_(i__),e__Cr_H_(i__),
+ e__Ni_H_(i__),e__Zn_H_(i__),e__Y_H_(i__),e__Ba_H_(i__),
+ n_Fe(i__),n_O(i__),n_Na(i__),n_Mg(i__),n_Al(i__),n_Si(i__),
+ n_Ca(i__),n_Ti(i__),n_Cr(i__),n_Ni(i__),n_Zn(i__),n_Y(i__),
+ n_Ba(i__)
if(ar__2(28:32) .EQ. '') v_O_Fe_(i__) = rNULL__
if(ar__2(34:38) .EQ. '') v_Na_Fe_(i__) = rNULL__
if(ar__2(70:74) .EQ. '') v_Cr_Fe_(i__) = rNULL__
if(ar__2(82:86) .EQ. '') v_Zn_Fe_(i__) = rNULL__
if(ar__2(88:92) .EQ. '') v_Y_Fe_(i__) = rNULL__
if(ar__2(105:108) .EQ. '') e__Fe_H_(i__) = rNULL__
if(ar__2(111:114) .EQ. '') e__O_Fe_(i__) = rNULL__
if(ar__2(117:120) .EQ. '') e__Na_Fe_(i__) = rNULL__
if(ar__2(123:126) .EQ. '') e__Mg_Fe_(i__) = rNULL__
if(ar__2(129:132) .EQ. '') e__Al_Fe_(i__) = rNULL__
if(ar__2(135:138) .EQ. '') e__Si_Fe_(i__) = rNULL__
if(ar__2(141:144) .EQ. '') e__Ca_Fe_(i__) = rNULL__
if(ar__2(147:150) .EQ. '') e__Ti_Fe_(i__) = rNULL__
if(ar__2(153:156) .EQ. '') e__Cr_Fe_(i__) = rNULL__
if(ar__2(159:162) .EQ. '') e__Ni_Fe_(i__) = rNULL__
if(ar__2(165:168) .EQ. '') e__Zn_Fe_(i__) = rNULL__
if(ar__2(171:174) .EQ. '') e__Y_Fe_(i__) = rNULL__
if(ar__2(177:180) .EQ. '') e__Ba_Fe_(i__) = rNULL__
if(ar__2(187:190) .EQ. '') e__Fe_H_r(i__) = rNULL__
if(ar__2(193:196) .EQ. '') e__O_H_(i__) = rNULL__
if(ar__2(199:202) .EQ. '') e__Na_H_(i__) = rNULL__
if(ar__2(205:208) .EQ. '') e__Mg_H_(i__) = rNULL__
if(ar__2(211:214) .EQ. '') e__Al_H_(i__) = rNULL__
if(ar__2(217:220) .EQ. '') e__Si_H_(i__) = rNULL__
if(ar__2(223:226) .EQ. '') e__Ca_H_(i__) = rNULL__
if(ar__2(229:232) .EQ. '') e__Ti_H_(i__) = rNULL__
if(ar__2(235:238) .EQ. '') e__Cr_H_(i__) = rNULL__
if(ar__2(241:244) .EQ. '') e__Ni_H_(i__) = rNULL__
if(ar__2(247:250) .EQ. '') e__Zn_H_(i__) = rNULL__
if(ar__2(253:256) .EQ. '') e__Y_H_(i__) = rNULL__
if(ar__2(259:262) .EQ. '') e__Ba_H_(i__) = rNULL__
if(ar__2(267:269) .EQ. '') n_Fe(i__) = iNULL__
if(ar__2(272:274) .EQ. '') n_O(i__) = iNULL__
if(ar__2(277:279) .EQ. '') n_Na(i__) = iNULL__
if(ar__2(282:284) .EQ. '') n_Mg(i__) = iNULL__
if(ar__2(287:289) .EQ. '') n_Al(i__) = iNULL__
if(ar__2(293:295) .EQ. '') n_Si(i__) = iNULL__
if(ar__2(299:301) .EQ. '') n_Ca(i__) = iNULL__
if(ar__2(305:307) .EQ. '') n_Ti(i__) = iNULL__
if(ar__2(311:313) .EQ. '') n_Cr(i__) = iNULL__
if(ar__2(317:319) .EQ. '') n_Ni(i__) = iNULL__
if(ar__2(322:324) .EQ. '') n_Zn(i__) = iNULL__
if(ar__2(327:329) .EQ. '') n_Y(i__) = iNULL__
if(ar__2(332:334) .EQ. '') n_Ba(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Name_1(i__),v_Fe_H_(i__),v_O_Fe_(i__),v_Na_Fe_(i__),
+ v_Mg_Fe_(i__),v_Al_Fe_(i__),v_Si_Fe_(i__),v_Ca_Fe_(i__),
+ v_Ti_Fe_(i__),v_Cr_Fe_(i__),v_Ni_Fe_(i__),v_Zn_Fe_(i__),
+ v_Y_Fe_(i__),v_Ba_Fe_(i__),e__Fe_H_(i__),e__O_Fe_(i__),
+ e__Na_Fe_(i__),e__Mg_Fe_(i__),e__Al_Fe_(i__),e__Si_Fe_(i__),
+ e__Ca_Fe_(i__),e__Ti_Fe_(i__),e__Cr_Fe_(i__),e__Ni_Fe_(i__),
+ e__Zn_Fe_(i__),e__Y_Fe_(i__),e__Ba_Fe_(i__),e__Fe_H_r(i__),
+ e__O_H_(i__),e__Na_H_(i__),e__Mg_H_(i__),e__Al_H_(i__),
+ e__Si_H_(i__),e__Ca_H_(i__),e__Ti_H_(i__),e__Cr_H_(i__),
+ e__Ni_H_(i__),e__Zn_H_(i__),e__Y_H_(i__),e__Ba_H_(i__),
+ n_Fe(i__),n_O(i__),n_Na(i__),n_Mg(i__),n_Al(i__),n_Si(i__),
+ n_Ca(i__),n_Ti(i__),n_Cr(i__),n_Ni(i__),n_Zn(i__),n_Y(i__),
+ n_Ba(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end