Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/605/A89 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-14
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/605/A89 Abundances of microlensed Bulge dwarf stars. VI. (Bensby+, 2017)
*================================================================================
*Chemical evolution of the Galactic bulge as traced by microlensed dwarfs and
*subgiant stars.
*VI. Age and abundance structure of the stellar populations in the central
*sub-kpc of the Milky Way.
* Bensby T., Feltzing S., Gould A., Yee, J.C., Johnson J.A., Asplund M.,
* Melendez J., Lucatello S., Howes L.M., McWilliam A., Udalski A.,
* Szymanski M.K., Soszynski I., Poleski R., Wyrzykowski L., Ulaczyk K.,
* Kozlowski S., Pietrukowicz P., Skowron J., Mroz P., Pawlak, M., Abe F.,
* Asakura Y., Bhattacharya A., Bond I.A., Bennett D.P., Hirao Y.,
* Nagakane M., Koshimoto N., Sumi T., Suzuki D., Tristram P.J.
* <Astron. Astrophys. 605, A89 (2017)>
* =2017A&A...605A..89B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Abundance ratios and uncertainties, stellar
parameters, radial velocities, stellar ages,
masses, photometry, and more
integer*4 nr__
parameter (nr__=91) ! Number of records
character*561 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)
integer*4 Num (nr__) ! Internal numbering of star
character*19 Name (nr__) ! Name of microlensing source star
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) [] Galactic longitude
real*4 GLAT (nr__) ! (deg) Galactic latitude
integer*4 Tevent (nr__) ! (d) ? Duration of microlensing event
real*8 Tmax (nr__) ! (d) Time of maximum magnification (HJD)
integer*4 Amax (nr__) ! ? Maximum magnification
real*8 Tobs (nr__) ! (d) Time when observed (MJD)
integer*4 ExpTime (nr__) ! (s) Exposure time
integer*4 snr (nr__) ! S/N per pixel at 640 nm
character*7 Spec (nr__) ! Spectrograph that was used
integer*4 R (nr__) ! Spectral resolving power
integer*4 Teff (nr__) ! (K) Effective temperature
integer*4 e_Teff (nr__) ! (K) Effective temperature uncertainty
real*4 logg (nr__) ! ([cm/s2]) Surface gravity
real*4 e_logg (nr__) ! ([cm/s2]) Surface gravity uncertainty
real*4 xit (nr__) ! (km/s) Microturbulence parameter
real*4 e_xit (nr__) ! (km/s) Microturbulence uncertainty
real*4 v_Fe_H_ (nr__) ! ([Sun]) Iron abundance normalised to Sun
real*4 e__Fe_H_ (nr__) ! ([Sun]) Iron abundance uncertainty
integer*4 o_FeI (nr__) ! Number of Fe I lines
integer*4 o_FeII (nr__) ! ? Number of Fe II lines
real*4 s_FeI (nr__) ! ([Sun]) Fe I line-to-line scatter
real*4 s_FeII (nr__) ! ([Sun]) ? Fe II line-to-line scatter
real*8 HRV (nr__) ! (km/s) Heliocentric radial velocity
real*8 vGC (nr__) ! (km/s) Galactocentric radial velocity
real*4 M (nr__) ! (Sun) Stellar mass
real*4 b_M (nr__) ! (Sun) ? Lower limit on stellar mass
real*4 B_M_1 (nr__) ! (Sun) ? Upper limit on stellar mass
real*4 Age (nr__) ! (Gyr) Stellar age
real*4 b_Age (nr__) ! (Gyr) ? Lower limit on stellar age
real*4 B_Age_1 (nr__) ! (Gyr) ? Upper limit on stellar age
real*4 v_V_I_0 (nr__) ! (mag) (V-I)_0 colour from spectroscopy
real*4 b__V_I_0 (nr__) ! (mag) Lower limit on (V-I)_0 colour
real*4 B__V_I_0_1 (nr__) ! (mag) Upper limit on (V-I)_0 colour
real*4 VMAG (nr__) ! (mag) Absolute V magnitude from spectroscopy
real*4 b_VMAG (nr__) ! (mag) Lower limit on absolute V magnitude
real*4 B_VMAG_1 (nr__) ! (mag) Upper limit on absolute V magnitude
real*4 IMAGmu (nr__) ! (mag) ? Absolute I magnitude from microlensing
real*4 v_V_I_0mu (nr__) ! (mag) ? Colour from microlensing
integer*4 Teffmu (nr__) ! (K) ? Effective temperature from microlensing
real*4 v_Na_Fe_ (nr__) ! ([Sun]) ? [Na/Fe] normalised abundance ratio
real*4 e__Na_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio
integer*4 o_Na (nr__) ! ? Number of Na lines used
real*4 s_Na (nr__) ! ([Sun]) ? Line-to-line scatter from Na lines
real*4 v_Mg_Fe_ (nr__) ! ([Sun]) ? [Mg/Fe] normalised abundance ratio
real*4 s__Mg_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio
integer*4 o_Mg (nr__) ! ? Number of Mg lines used
real*4 s_Mg (nr__) ! ([Sun]) ? Line-to-line scatter from Mg lines
real*4 v_Al_Fe_ (nr__) ! ([Sun]) ? [Al/Fe] normalised abundance ratio
real*4 e__Al_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Al/Fe]
integer*4 o_Al (nr__) ! ? Number of Al lines used
real*4 s_Al (nr__) ! ([Sun]) ? Line-to-line scatter from Al lines
real*4 v_Si_Fe_ (nr__) ! ([Sun]) ? [Si/Fe] normalised abundance ratio
real*4 e__Si_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Si/Fe]
integer*4 o_Si (nr__) ! ? Number of Si lines used
real*4 s_Si (nr__) ! ([Sun]) ? Line-to-line scatter from Si lines
real*4 v_Ca_Fe_ (nr__) ! ([Sun]) ? [Ca/Fe] normalised abundance ratio
real*4 e__Ca_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Ca/Fe]
integer*4 o_Ca (nr__) ! ? Number of Ca lines used
real*4 s_Ca (nr__) ! ([Sun]) ? Line-to-line scatter from Ca lines
real*4 v_Ti_Fe_ (nr__) ! ([Sun]) ? [Ti/Fe] normalised abundance ratio
real*4 e__Ti_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Ti/Fe]
integer*4 o_Ti (nr__) ! ? Number of Ti lines used
real*4 s_Ti (nr__) ! ([Sun]) ? Line-to-line scatter from Ti lines
real*4 v_Cr_Fe_ (nr__) ! ([Sun]) ? [Cr/Fe] normalised abundance ratio
real*4 e__Cr_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Cr/Fe]
integer*4 o_Cr (nr__) ! ? Number of Cr lines used
real*4 s_Cr (nr__) ! ([Sun]) ? Line-to-line scatter from Cr lines
real*4 v_Ni_Fe_ (nr__) ! ([Sun]) ? [Ni/Fe] normalised abundance ratio
real*4 e__Ni_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Ni/Fe]
integer*4 o_Ni (nr__) ! ? Number of Ni lines used
real*4 s_Ni (nr__) ! ([Sun]) ? Line-to-line scatter from Ni lines
real*4 v_Zn_Fe_ (nr__) ! ([Sun]) ? [Zn/Fe] normalised abundance ratio
real*4 e__Zn_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Zn/Fe]
integer*4 o_Zn (nr__) ! ? Number of Zn lines used
real*4 s_Zn (nr__) ! ([Sun]) ? Line-to-line scatter from Zn lines
real*4 v_Y_Fe_ (nr__) ! ([Sun]) ? [Y/Fe] normalised abundance ratio
real*4 e__Y_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Y/Fe]
integer*4 o_Y (nr__) ! ? Number of Y lines used
real*4 s_Y (nr__) ! ([Sun]) ? Line-to-line scatter from Y lines
real*4 v_Ba_Fe_ (nr__) ! ([Sun]) ? [Ba/Fe] normalised abundance ratio
real*4 e__Ba_Fe_ (nr__) ! ([Sun]) ? Error in abundance ratio [Ba/Fe]
integer*4 o_Ba (nr__) ! ? Number of Ba lines used
real*4 s_Ba (nr__) ! ([Sun]) ? Line-to-line scatter from Ba lines
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Equivalent widths and absolute abundances for
individual lines
integer*4 nr__1
parameter (nr__1=494) ! Number of records
character*1038 ar__1 ! Full-size record
character*2 Element (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 loggf (nr__1) ! Oscillator strength
real*4 EWSun (nr__1) ! (0.1pm) ? Equivalent width for Sun
real*4 ASun (nr__1) ! ([-]) ? Sun absolute abundance
* (log scale, A(H)=12)
real*4 EW1 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-174
real*4 A1 (nr__1) ! ? Absolute abundance MOA-2009-BLG-174
real*4 EW2 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-259
real*4 A2 (nr__1) ! ? Absolute abundance MOA-2009-BLG-259
real*4 EW3 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-167
real*4 A3 (nr__1) ! ? Absolute abundance MOA-2010-BLG-167
real*4 EW4 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-311
real*4 A4 (nr__1) ! ? Absolute abundance MOA-2010-BLG-311
real*4 EW5 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-446
real*4 A5 (nr__1) ! ? Absolute abundance MOA-2010-BLG-446
real*4 EW6 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-523
real*4 A6 (nr__1) ! ? Absolute abundance MOA-2010-BLG-523
real*4 EW7 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2011-BLG-0950
real*4 A7 (nr__1) ! ? Absolute abundance OGLE-2011-BLG-0950
real*4 EW8 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2011-BLG-0969
real*4 A8 (nr__1) ! ? Absolute abundance OGLE-2011-BLG-0969
real*4 EW9 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-034
real*4 A9 (nr__1) ! ? Absolute abundance MOA-2011-BLG-034
real*4 EW10 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-058
real*4 A10 (nr__1) ! ? Absolute abundance MOA-2011-BLG-058
real*4 EW11 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2011-BLG-1072
real*4 A11 (nr__1) ! ? Absolute abundance OGLE-2011-BLG-1072
real*4 EW12 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-090
real*4 A12 (nr__1) ! ? Absolute abundance MOA-2011-BLG-090
real*4 EW13 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-104
real*4 A13 (nr__1) ! ? Absolute abundance MOA-2011-BLG-104
real*4 EW14 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2011-BLG-1105
real*4 A14 (nr__1) ! ? Absolute abundance OGLE-2011-BLG-1105
real*4 EW15 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-174
real*4 A15 (nr__1) ! ? Absolute abundance MOA-2011-BLG-174
real*4 EW16 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-191
real*4 A16 (nr__1) ! ? Absolute abundance MOA-2011-BLG-191
real*4 EW17 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-234
real*4 A17 (nr__1) ! ? Absolute abundance MOA-2011-BLG-234
real*4 EW18 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-278
real*4 A18 (nr__1) ! ? Absolute abundance MOA-2011-BLG-278
real*4 EW19 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2011-BLG-1410
real*4 A19 (nr__1) ! ? Absolute abundance OGLE-2011-BLG-1410
real*4 EW20 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2011-BLG-445
real*4 A20 (nr__1) ! ? Absolute abundance MOA-2011-BLG-445
real*4 EW21 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-022
real*4 A21 (nr__1) ! ? Absolute abundance MOA-2012-BLG-022
real*4 EW22 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0026
real*4 A22 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0026
real*4 EW23 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-187
real*4 A23 (nr__1) ! ? Absolute abundance MOA-2012-BLG-187
real*4 EW24 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-202
real*4 A24 (nr__1) ! ? Absolute abundance MOA-2012-BLG-202
real*4 EW25 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0211
real*4 A25 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0211
real*4 EW26 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0270
real*4 A26 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0270
real*4 EW27 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-291
real*4 A27 (nr__1) ! ? Absolute abundance MOA-2012-BLG-291
real*4 EW28 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-391
real*4 A28 (nr__1) ! ? Absolute abundance MOA-2012-BLG-391
real*4 EW29 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-410
real*4 A29 (nr__1) ! ? Absolute abundance MOA-2012-BLG-410
real*4 EW30 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0521
real*4 A30 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0521
real*4 EW31 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2012-BLG-532
real*4 A31 (nr__1) ! ? Absolute abundance MOA-2012-BLG-532
real*4 EW32 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0563
real*4 A32 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0563
real*4 EW33 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0617
real*4 A33 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0617
real*4 EW34 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-0816
real*4 A34 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-0816
real*4 EW35 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1156
real*4 A35 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1156
real*4 EW36 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1217
real*4 A36 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1217
real*4 EW37 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1274
real*4 A37 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1274
real*4 EW38 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1279
real*4 A38 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1279
real*4 EW39 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1526
real*4 A39 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1526
real*4 EW40 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2012-BLG-1534
real*4 A40 (nr__1) ! ? Absolute abundance OGLE-2012-BLG-1534
real*4 EW41 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-063
real*4 A41 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-063
real*4 EW42 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-068
real*4 A42 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-068
real*4 EW43 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-299
real*4 A43 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-299
real*4 EW44 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-402
real*4 A44 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-402
real*4 EW45 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-0446
real*4 A45 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-0446
real*4 EW46 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-517
real*4 A46 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-517
real*4 EW47 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-524
real*4 A47 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-524
real*4 EW48 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2013-BLG-605
real*4 A48 (nr__1) ! ? Absolute abundance for MOA-2013-BLG-605
real*4 EW49 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-0692
real*4 A49 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-0692
real*4 EW50 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-0835
real*4 A50 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-0835
real*4 EW51 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-0911
real*4 A51 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-0911
real*4 EW52 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1015
real*4 A52 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1015
real*4 EW53 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1114
real*4 A53 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1114
real*4 EW54 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1125
real*4 A54 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1125
real*4 EW55 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1147
real*4 A55 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1147
real*4 EW56 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1259
real*4 A56 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1259
real*4 EW57 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1768
real*4 A57 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1768
real*4 EW58 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1793
real*4 A58 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1793
real*4 EW59 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1868
real*4 A59 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1868
real*4 EW60 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2013-BLG-1938
real*4 A60 (nr__1) ! ? Absolute abundance for OGLE-2013-BLG-1938
real*4 EW61 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2014-BLG-131
real*4 A61 (nr__1) ! ? Absolute abundance for MOA-2014-BLG-131
real*4 EW62 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-0157
real*4 A62 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-0157
real*4 EW63 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-0801
real*4 A63 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-0801
real*4 EW64 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-0953
real*4 A64 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-0953
real*4 EW65 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-0987
real*4 A65 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-0987
real*4 EW66 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-1122
real*4 A66 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-1122
real*4 EW67 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-1370
real*4 A67 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-1370
real*4 EW68 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-1418
real*4 A68 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-1418
real*4 EW69 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-1469
real*4 A69 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-1469
real*4 EW70 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2014-BLG-2040
real*4 A70 (nr__1) ! ? Absolute abundance for OGLE-2014-BLG-2040
real*4 EW71 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2015-BLG-0078
real*4 A71 (nr__1) ! ? Absolute abundance for OGLE-2015-BLG-0078
real*4 EW72 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2015-BLG-111
real*4 A72 (nr__1) ! ? Absolute abundance for MOA-2015-BLG-111
real*4 EW73 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2015-BLG-0159
real*4 A73 (nr__1) ! ? Absolute abundance for OGLE-2015-BLG-0159
real*4 EW74 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-037
real*4 A74 (nr__1) ! ? Absolute abundance MOA-2010-BLG-037
real*4 EW75 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-049
real*4 A75 (nr__1) ! ? Absolute abundance MOA-2010-BLG-049
real*4 EW76 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2009-BLG-076
real*4 A76 (nr__1) ! ? Absolute abundance OGLE-2009-BLG-076
real*4 EW77 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-078
real*4 A77 (nr__1) ! ? Absolute abundance MOA-2010-BLG-078
real*4 EW78 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2006-BLG-099
real*4 A78 (nr__1) ! ? Absolute abundance MOA-2006-BLG-099
real*4 EW79 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-133
real*4 A79 (nr__1) ! ? Absolute abundance MOA-2009-BLG-133
real*4 EW80 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2008-BLG-209
real*4 A80 (nr__1) ! ? Absolute abundance OGLE-2008-BLG-209
real*4 EW81 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2006-BLG-265
real*4 A81 (nr__1) ! ? Absolute abundance OGLE-2006-BLG-265
real*4 EW82 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2010-BLG-285
real*4 A82 (nr__1) ! ? Absolute abundance MOA-2010-BLG-285
real*4 EW83 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2008-BLG-310
real*4 A83 (nr__1) ! ? Absolute abundance MOA-2008-BLG-310
real*4 EW84 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2008-BLG-311
real*4 A84 (nr__1) ! ? Absolute abundance MOA-2008-BLG-311
real*4 EW85 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2007-BLG-349
real*4 A85 (nr__1) ! ? Absolute abundance OGLE-2007-BLG-349
real*4 EW86 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-456
real*4 A86 (nr__1) ! ? Absolute abundance MOA-2009-BLG-456
real*4 EW87 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-475
real*4 A87 (nr__1) ! ? Absolute abundance MOA-2009-BLG-475
real*4 EW88 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-489
real*4 A88 (nr__1) ! ? Absolute abundance MOA-2009-BLG-489
real*4 EW89 (nr__1) ! (0.1pm) ? Equivalent width for MOA-2009-BLG-493
real*4 A89 (nr__1) ! ? Absolute abundance MOA-2009-BLG-493
real*4 EW90 (nr__1) ! (0.1pm) ? Equivalent width for OGLE-2007-BLG-514
real*4 A90 (nr__1) ! ? Absolute abundance OGLE-2007-BLG-514
real*4 EW91 (nr__1) ! (0.1pm) ? Equivalent width for MACH-1999-BLG-022
real*4 A91 (nr__1) ! ? Absolute abundance MACH-1999-BLG-022
C=============================================================================
C Loading file 'tablea2.dat' ! Abundance ratios and uncertainties, stellar
* parameters, radial velocities, stellar ages,
* masses, photometry, and more
C Format for file interpretation
1 format(
+ I6,4X,A19,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,I4,1X,F7.2,1X,I4,1X,F8.3,1X,I5,3X,I3,1X,A7,1X,I5,
+ 2X,I4,2X,I4,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,1X,F5.2,2X,F4.2,
+ 3X,I3,4X,I2,2X,F4.2,2X,F4.2,1X,F7.1,1X,F7.1,2X,F4.2,2X,F4.2,
+ 2X,F4.2,2X,F4.1,2X,F4.1,2X,F4.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,2X,I4,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,I3,2X,F4.2,1X,F5.2,2X,F4.2,3X,
+ I3,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,91
read(1,'(A561)')ar__
read(ar__,1)
+ Num(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Tevent(i__),
+ Tmax(i__),Amax(i__),Tobs(i__),ExpTime(i__),snr(i__),Spec(i__),
+ R(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),xit(i__),
+ e_xit(i__),v_Fe_H_(i__),e__Fe_H_(i__),o_FeI(i__),o_FeII(i__),
+ s_FeI(i__),s_FeII(i__),HRV(i__),vGC(i__),M(i__),b_M(i__),
+ B_M_1(i__),Age(i__),b_Age(i__),B_Age_1(i__),v_V_I_0(i__),
+ b__V_I_0(i__),B__V_I_0_1(i__),VMAG(i__),b_VMAG(i__),
+ B_VMAG_1(i__),IMAGmu(i__),v_V_I_0mu(i__),Teffmu(i__),
+ v_Na_Fe_(i__),e__Na_Fe_(i__),o_Na(i__),s_Na(i__),
+ v_Mg_Fe_(i__),s__Mg_Fe_(i__),o_Mg(i__),s_Mg(i__),
+ v_Al_Fe_(i__),e__Al_Fe_(i__),o_Al(i__),s_Al(i__),
+ v_Si_Fe_(i__),e__Si_Fe_(i__),o_Si(i__),s_Si(i__),
+ v_Ca_Fe_(i__),e__Ca_Fe_(i__),o_Ca(i__),s_Ca(i__),
+ v_Ti_Fe_(i__),e__Ti_Fe_(i__),o_Ti(i__),s_Ti(i__),
+ v_Cr_Fe_(i__),e__Cr_Fe_(i__),o_Cr(i__),s_Cr(i__),
+ v_Ni_Fe_(i__),e__Ni_Fe_(i__),o_Ni(i__),s_Ni(i__),
+ v_Zn_Fe_(i__),e__Zn_Fe_(i__),o_Zn(i__),s_Zn(i__),v_Y_Fe_(i__),
+ e__Y_Fe_(i__),o_Y(i__),s_Y(i__),v_Ba_Fe_(i__),e__Ba_Fe_(i__),
+ o_Ba(i__),s_Ba(i__)
if(ar__(68:71) .EQ. '') Tevent(i__) = iNULL__
if(ar__(81:84) .EQ. '') Amax(i__) = iNULL__
if(ar__(178:179) .EQ. '') o_FeII(i__) = iNULL__
if(ar__(188:191) .EQ. '') s_FeII(i__) = rNULL__
if(ar__(216:219) .EQ. '') b_M(i__) = rNULL__
if(ar__(222:225) .EQ. '') B_M_1(i__) = rNULL__
if(ar__(234:237) .EQ. '') b_Age(i__) = rNULL__
if(ar__(240:243) .EQ. '') B_Age_1(i__) = rNULL__
if(ar__(281:285) .EQ. '') IMAGmu(i__) = rNULL__
if(ar__(287:291) .EQ. '') v_V_I_0mu(i__) = rNULL__
if(ar__(294:297) .EQ. '') Teffmu(i__) = iNULL__
if(ar__(299:303) .EQ. '') v_Na_Fe_(i__) = rNULL__
if(ar__(306:309) .EQ. '') e__Na_Fe_(i__) = rNULL__
if(ar__(313:315) .EQ. '') o_Na(i__) = iNULL__
if(ar__(318:321) .EQ. '') s_Na(i__) = rNULL__
if(ar__(323:327) .EQ. '') v_Mg_Fe_(i__) = rNULL__
if(ar__(330:333) .EQ. '') s__Mg_Fe_(i__) = rNULL__
if(ar__(337:339) .EQ. '') o_Mg(i__) = iNULL__
if(ar__(342:345) .EQ. '') s_Mg(i__) = rNULL__
if(ar__(347:351) .EQ. '') v_Al_Fe_(i__) = rNULL__
if(ar__(354:357) .EQ. '') e__Al_Fe_(i__) = rNULL__
if(ar__(361:363) .EQ. '') o_Al(i__) = iNULL__
if(ar__(366:369) .EQ. '') s_Al(i__) = rNULL__
if(ar__(371:375) .EQ. '') v_Si_Fe_(i__) = rNULL__
if(ar__(378:381) .EQ. '') e__Si_Fe_(i__) = rNULL__
if(ar__(385:387) .EQ. '') o_Si(i__) = iNULL__
if(ar__(390:393) .EQ. '') s_Si(i__) = rNULL__
if(ar__(395:399) .EQ. '') v_Ca_Fe_(i__) = rNULL__
if(ar__(402:405) .EQ. '') e__Ca_Fe_(i__) = rNULL__
if(ar__(409:411) .EQ. '') o_Ca(i__) = iNULL__
if(ar__(414:417) .EQ. '') s_Ca(i__) = rNULL__
if(ar__(419:423) .EQ. '') v_Ti_Fe_(i__) = rNULL__
if(ar__(426:429) .EQ. '') e__Ti_Fe_(i__) = rNULL__
if(ar__(433:435) .EQ. '') o_Ti(i__) = iNULL__
if(ar__(438:441) .EQ. '') s_Ti(i__) = rNULL__
if(ar__(443:447) .EQ. '') v_Cr_Fe_(i__) = rNULL__
if(ar__(450:453) .EQ. '') e__Cr_Fe_(i__) = rNULL__
if(ar__(457:459) .EQ. '') o_Cr(i__) = iNULL__
if(ar__(462:465) .EQ. '') s_Cr(i__) = rNULL__
if(ar__(467:471) .EQ. '') v_Ni_Fe_(i__) = rNULL__
if(ar__(474:477) .EQ. '') e__Ni_Fe_(i__) = rNULL__
if(ar__(481:483) .EQ. '') o_Ni(i__) = iNULL__
if(ar__(486:489) .EQ. '') s_Ni(i__) = rNULL__
if(ar__(491:495) .EQ. '') v_Zn_Fe_(i__) = rNULL__
if(ar__(498:501) .EQ. '') e__Zn_Fe_(i__) = rNULL__
if(ar__(505:507) .EQ. '') o_Zn(i__) = iNULL__
if(ar__(510:513) .EQ. '') s_Zn(i__) = rNULL__
if(ar__(515:519) .EQ. '') v_Y_Fe_(i__) = rNULL__
if(ar__(522:525) .EQ. '') e__Y_Fe_(i__) = rNULL__
if(ar__(529:531) .EQ. '') o_Y(i__) = iNULL__
if(ar__(534:537) .EQ. '') s_Y(i__) = rNULL__
if(ar__(539:543) .EQ. '') v_Ba_Fe_(i__) = rNULL__
if(ar__(546:549) .EQ. '') e__Ba_Fe_(i__) = rNULL__
if(ar__(553:555) .EQ. '') o_Ba(i__) = iNULL__
if(ar__(558:561) .EQ. '') s_Ba(i__) = rNULL__
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)
+ Num(i__),Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Tevent(i__),
+ Tmax(i__),Amax(i__),Tobs(i__),ExpTime(i__),snr(i__),Spec(i__),
+ R(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),xit(i__),
+ e_xit(i__),v_Fe_H_(i__),e__Fe_H_(i__),o_FeI(i__),o_FeII(i__),
+ s_FeI(i__),s_FeII(i__),HRV(i__),vGC(i__),M(i__),b_M(i__),
+ B_M_1(i__),Age(i__),b_Age(i__),B_Age_1(i__),v_V_I_0(i__),
+ b__V_I_0(i__),B__V_I_0_1(i__),VMAG(i__),b_VMAG(i__),
+ B_VMAG_1(i__),IMAGmu(i__),v_V_I_0mu(i__),Teffmu(i__),
+ v_Na_Fe_(i__),e__Na_Fe_(i__),o_Na(i__),s_Na(i__),
+ v_Mg_Fe_(i__),s__Mg_Fe_(i__),o_Mg(i__),s_Mg(i__),
+ v_Al_Fe_(i__),e__Al_Fe_(i__),o_Al(i__),s_Al(i__),
+ v_Si_Fe_(i__),e__Si_Fe_(i__),o_Si(i__),s_Si(i__),
+ v_Ca_Fe_(i__),e__Ca_Fe_(i__),o_Ca(i__),s_Ca(i__),
+ v_Ti_Fe_(i__),e__Ti_Fe_(i__),o_Ti(i__),s_Ti(i__),
+ v_Cr_Fe_(i__),e__Cr_Fe_(i__),o_Cr(i__),s_Cr(i__),
+ v_Ni_Fe_(i__),e__Ni_Fe_(i__),o_Ni(i__),s_Ni(i__),
+ v_Zn_Fe_(i__),e__Zn_Fe_(i__),o_Zn(i__),s_Zn(i__),v_Y_Fe_(i__),
+ e__Y_Fe_(i__),o_Y(i__),s_Y(i__),v_Ba_Fe_(i__),e__Ba_Fe_(i__),
+ o_Ba(i__),s_Ba(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 'tablea1.dat' ! Equivalent widths and absolute abundances for
* individual lines
C Format for file interpretation
2 format(
+ A2,1X,I1,1X,F8.3,1X,F5.3,1X,F6.3,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,
+ 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,
+ F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,
+ 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,494
read(1,'(A1038)')ar__1
read(ar__1,2)
+ Element(i__),Ion(i__),Lambda(i__),LEP(i__),loggf(i__),
+ EWSun(i__),ASun(i__),EW1(i__),A1(i__),EW2(i__),A2(i__),
+ EW3(i__),A3(i__),EW4(i__),A4(i__),EW5(i__),A5(i__),EW6(i__),
+ A6(i__),EW7(i__),A7(i__),EW8(i__),A8(i__),EW9(i__),A9(i__),
+ EW10(i__),A10(i__),EW11(i__),A11(i__),EW12(i__),A12(i__),
+ EW13(i__),A13(i__),EW14(i__),A14(i__),EW15(i__),A15(i__),
+ EW16(i__),A16(i__),EW17(i__),A17(i__),EW18(i__),A18(i__),
+ EW19(i__),A19(i__),EW20(i__),A20(i__),EW21(i__),A21(i__),
+ EW22(i__),A22(i__),EW23(i__),A23(i__),EW24(i__),A24(i__),
+ EW25(i__),A25(i__),EW26(i__),A26(i__),EW27(i__),A27(i__),
+ EW28(i__),A28(i__),EW29(i__),A29(i__),EW30(i__),A30(i__),
+ EW31(i__),A31(i__),EW32(i__),A32(i__),EW33(i__),A33(i__),
+ EW34(i__),A34(i__),EW35(i__),A35(i__),EW36(i__),A36(i__),
+ EW37(i__),A37(i__),EW38(i__),A38(i__),EW39(i__),A39(i__),
+ EW40(i__),A40(i__),EW41(i__),A41(i__),EW42(i__),A42(i__),
+ EW43(i__),A43(i__),EW44(i__),A44(i__),EW45(i__),A45(i__),
+ EW46(i__),A46(i__),EW47(i__),A47(i__),EW48(i__),A48(i__),
+ EW49(i__),A49(i__),EW50(i__),A50(i__),EW51(i__),A51(i__),
+ EW52(i__),A52(i__),EW53(i__),A53(i__),EW54(i__),A54(i__),
+ EW55(i__),A55(i__),EW56(i__),A56(i__),EW57(i__),A57(i__),
+ EW58(i__),A58(i__),EW59(i__),A59(i__),EW60(i__),A60(i__),
+ EW61(i__),A61(i__),EW62(i__),A62(i__),EW63(i__),A63(i__),
+ EW64(i__),A64(i__),EW65(i__),A65(i__),EW66(i__),A66(i__),
+ EW67(i__),A67(i__),EW68(i__),A68(i__),EW69(i__),A69(i__),
+ EW70(i__),A70(i__),EW71(i__),A71(i__),EW72(i__),A72(i__),
+ EW73(i__),A73(i__),EW74(i__),A74(i__),EW75(i__),A75(i__),
+ EW76(i__),A76(i__),EW77(i__),A77(i__),EW78(i__),A78(i__),
+ EW79(i__),A79(i__),EW80(i__),A80(i__),EW81(i__),A81(i__),
+ EW82(i__),A82(i__),EW83(i__),A83(i__),EW84(i__),A84(i__),
+ EW85(i__),A85(i__),EW86(i__),A86(i__),EW87(i__),A87(i__),
+ EW88(i__),A88(i__),EW89(i__),A89(i__),EW90(i__),A90(i__),
+ EW91(i__),A91(i__)
if(ar__1(28:32) .EQ. '') EWSun(i__) = rNULL__
if(ar__1(34:37) .EQ. '') ASun(i__) = rNULL__
if(ar__1(39:43) .EQ. '') EW1(i__) = rNULL__
if(ar__1(45:48) .EQ. '') A1(i__) = rNULL__
if(ar__1(50:54) .EQ. '') EW2(i__) = rNULL__
if(ar__1(56:59) .EQ. '') A2(i__) = rNULL__
if(ar__1(61:65) .EQ. '') EW3(i__) = rNULL__
if(ar__1(67:70) .EQ. '') A3(i__) = rNULL__
if(ar__1(72:76) .EQ. '') EW4(i__) = rNULL__
if(ar__1(78:81) .EQ. '') A4(i__) = rNULL__
if(ar__1(83:87) .EQ. '') EW5(i__) = rNULL__
if(ar__1(89:92) .EQ. '') A5(i__) = rNULL__
if(ar__1(94:98) .EQ. '') EW6(i__) = rNULL__
if(ar__1(100:103) .EQ. '') A6(i__) = rNULL__
if(ar__1(105:109) .EQ. '') EW7(i__) = rNULL__
if(ar__1(111:114) .EQ. '') A7(i__) = rNULL__
if(ar__1(116:120) .EQ. '') EW8(i__) = rNULL__
if(ar__1(122:125) .EQ. '') A8(i__) = rNULL__
if(ar__1(127:131) .EQ. '') EW9(i__) = rNULL__
if(ar__1(133:136) .EQ. '') A9(i__) = rNULL__
if(ar__1(138:142) .EQ. '') EW10(i__) = rNULL__
if(ar__1(144:147) .EQ. '') A10(i__) = rNULL__
if(ar__1(149:153) .EQ. '') EW11(i__) = rNULL__
if(ar__1(155:158) .EQ. '') A11(i__) = rNULL__
if(ar__1(160:164) .EQ. '') EW12(i__) = rNULL__
if(ar__1(166:169) .EQ. '') A12(i__) = rNULL__
if(ar__1(171:175) .EQ. '') EW13(i__) = rNULL__
if(ar__1(177:180) .EQ. '') A13(i__) = rNULL__
if(ar__1(182:186) .EQ. '') EW14(i__) = rNULL__
if(ar__1(188:191) .EQ. '') A14(i__) = rNULL__
if(ar__1(193:197) .EQ. '') EW15(i__) = rNULL__
if(ar__1(199:202) .EQ. '') A15(i__) = rNULL__
if(ar__1(204:208) .EQ. '') EW16(i__) = rNULL__
if(ar__1(210:213) .EQ. '') A16(i__) = rNULL__
if(ar__1(215:219) .EQ. '') EW17(i__) = rNULL__
if(ar__1(221:224) .EQ. '') A17(i__) = rNULL__
if(ar__1(226:230) .EQ. '') EW18(i__) = rNULL__
if(ar__1(232:235) .EQ. '') A18(i__) = rNULL__
if(ar__1(237:241) .EQ. '') EW19(i__) = rNULL__
if(ar__1(243:246) .EQ. '') A19(i__) = rNULL__
if(ar__1(248:252) .EQ. '') EW20(i__) = rNULL__
if(ar__1(254:257) .EQ. '') A20(i__) = rNULL__
if(ar__1(259:263) .EQ. '') EW21(i__) = rNULL__
if(ar__1(265:268) .EQ. '') A21(i__) = rNULL__
if(ar__1(270:274) .EQ. '') EW22(i__) = rNULL__
if(ar__1(276:279) .EQ. '') A22(i__) = rNULL__
if(ar__1(281:285) .EQ. '') EW23(i__) = rNULL__
if(ar__1(287:290) .EQ. '') A23(i__) = rNULL__
if(ar__1(292:296) .EQ. '') EW24(i__) = rNULL__
if(ar__1(298:301) .EQ. '') A24(i__) = rNULL__
if(ar__1(303:307) .EQ. '') EW25(i__) = rNULL__
if(ar__1(309:312) .EQ. '') A25(i__) = rNULL__
if(ar__1(314:318) .EQ. '') EW26(i__) = rNULL__
if(ar__1(320:323) .EQ. '') A26(i__) = rNULL__
if(ar__1(325:329) .EQ. '') EW27(i__) = rNULL__
if(ar__1(331:334) .EQ. '') A27(i__) = rNULL__
if(ar__1(336:340) .EQ. '') EW28(i__) = rNULL__
if(ar__1(342:345) .EQ. '') A28(i__) = rNULL__
if(ar__1(347:351) .EQ. '') EW29(i__) = rNULL__
if(ar__1(353:356) .EQ. '') A29(i__) = rNULL__
if(ar__1(358:362) .EQ. '') EW30(i__) = rNULL__
if(ar__1(364:367) .EQ. '') A30(i__) = rNULL__
if(ar__1(369:373) .EQ. '') EW31(i__) = rNULL__
if(ar__1(375:378) .EQ. '') A31(i__) = rNULL__
if(ar__1(380:384) .EQ. '') EW32(i__) = rNULL__
if(ar__1(386:389) .EQ. '') A32(i__) = rNULL__
if(ar__1(391:395) .EQ. '') EW33(i__) = rNULL__
if(ar__1(397:400) .EQ. '') A33(i__) = rNULL__
if(ar__1(402:406) .EQ. '') EW34(i__) = rNULL__
if(ar__1(408:411) .EQ. '') A34(i__) = rNULL__
if(ar__1(413:417) .EQ. '') EW35(i__) = rNULL__
if(ar__1(419:422) .EQ. '') A35(i__) = rNULL__
if(ar__1(424:428) .EQ. '') EW36(i__) = rNULL__
if(ar__1(430:433) .EQ. '') A36(i__) = rNULL__
if(ar__1(435:439) .EQ. '') EW37(i__) = rNULL__
if(ar__1(441:444) .EQ. '') A37(i__) = rNULL__
if(ar__1(446:450) .EQ. '') EW38(i__) = rNULL__
if(ar__1(452:455) .EQ. '') A38(i__) = rNULL__
if(ar__1(457:461) .EQ. '') EW39(i__) = rNULL__
if(ar__1(463:466) .EQ. '') A39(i__) = rNULL__
if(ar__1(468:472) .EQ. '') EW40(i__) = rNULL__
if(ar__1(474:477) .EQ. '') A40(i__) = rNULL__
if(ar__1(479:483) .EQ. '') EW41(i__) = rNULL__
if(ar__1(485:488) .EQ. '') A41(i__) = rNULL__
if(ar__1(490:494) .EQ. '') EW42(i__) = rNULL__
if(ar__1(496:499) .EQ. '') A42(i__) = rNULL__
if(ar__1(501:505) .EQ. '') EW43(i__) = rNULL__
if(ar__1(507:510) .EQ. '') A43(i__) = rNULL__
if(ar__1(512:516) .EQ. '') EW44(i__) = rNULL__
if(ar__1(518:521) .EQ. '') A44(i__) = rNULL__
if(ar__1(523:527) .EQ. '') EW45(i__) = rNULL__
if(ar__1(529:532) .EQ. '') A45(i__) = rNULL__
if(ar__1(534:538) .EQ. '') EW46(i__) = rNULL__
if(ar__1(540:543) .EQ. '') A46(i__) = rNULL__
if(ar__1(545:549) .EQ. '') EW47(i__) = rNULL__
if(ar__1(551:554) .EQ. '') A47(i__) = rNULL__
if(ar__1(556:560) .EQ. '') EW48(i__) = rNULL__
if(ar__1(562:565) .EQ. '') A48(i__) = rNULL__
if(ar__1(567:571) .EQ. '') EW49(i__) = rNULL__
if(ar__1(573:576) .EQ. '') A49(i__) = rNULL__
if(ar__1(578:582) .EQ. '') EW50(i__) = rNULL__
if(ar__1(584:587) .EQ. '') A50(i__) = rNULL__
if(ar__1(589:593) .EQ. '') EW51(i__) = rNULL__
if(ar__1(595:598) .EQ. '') A51(i__) = rNULL__
if(ar__1(600:604) .EQ. '') EW52(i__) = rNULL__
if(ar__1(606:609) .EQ. '') A52(i__) = rNULL__
if(ar__1(611:615) .EQ. '') EW53(i__) = rNULL__
if(ar__1(617:620) .EQ. '') A53(i__) = rNULL__
if(ar__1(622:626) .EQ. '') EW54(i__) = rNULL__
if(ar__1(628:631) .EQ. '') A54(i__) = rNULL__
if(ar__1(633:637) .EQ. '') EW55(i__) = rNULL__
if(ar__1(639:642) .EQ. '') A55(i__) = rNULL__
if(ar__1(644:648) .EQ. '') EW56(i__) = rNULL__
if(ar__1(650:653) .EQ. '') A56(i__) = rNULL__
if(ar__1(655:659) .EQ. '') EW57(i__) = rNULL__
if(ar__1(661:664) .EQ. '') A57(i__) = rNULL__
if(ar__1(666:670) .EQ. '') EW58(i__) = rNULL__
if(ar__1(672:675) .EQ. '') A58(i__) = rNULL__
if(ar__1(677:681) .EQ. '') EW59(i__) = rNULL__
if(ar__1(683:686) .EQ. '') A59(i__) = rNULL__
if(ar__1(688:692) .EQ. '') EW60(i__) = rNULL__
if(ar__1(694:697) .EQ. '') A60(i__) = rNULL__
if(ar__1(699:703) .EQ. '') EW61(i__) = rNULL__
if(ar__1(705:708) .EQ. '') A61(i__) = rNULL__
if(ar__1(710:714) .EQ. '') EW62(i__) = rNULL__
if(ar__1(716:719) .EQ. '') A62(i__) = rNULL__
if(ar__1(721:725) .EQ. '') EW63(i__) = rNULL__
if(ar__1(727:730) .EQ. '') A63(i__) = rNULL__
if(ar__1(732:736) .EQ. '') EW64(i__) = rNULL__
if(ar__1(738:741) .EQ. '') A64(i__) = rNULL__
if(ar__1(743:747) .EQ. '') EW65(i__) = rNULL__
if(ar__1(749:752) .EQ. '') A65(i__) = rNULL__
if(ar__1(754:758) .EQ. '') EW66(i__) = rNULL__
if(ar__1(760:763) .EQ. '') A66(i__) = rNULL__
if(ar__1(765:769) .EQ. '') EW67(i__) = rNULL__
if(ar__1(771:774) .EQ. '') A67(i__) = rNULL__
if(ar__1(776:780) .EQ. '') EW68(i__) = rNULL__
if(ar__1(782:785) .EQ. '') A68(i__) = rNULL__
if(ar__1(787:791) .EQ. '') EW69(i__) = rNULL__
if(ar__1(793:796) .EQ. '') A69(i__) = rNULL__
if(ar__1(798:802) .EQ. '') EW70(i__) = rNULL__
if(ar__1(804:807) .EQ. '') A70(i__) = rNULL__
if(ar__1(809:813) .EQ. '') EW71(i__) = rNULL__
if(ar__1(815:818) .EQ. '') A71(i__) = rNULL__
if(ar__1(820:824) .EQ. '') EW72(i__) = rNULL__
if(ar__1(826:829) .EQ. '') A72(i__) = rNULL__
if(ar__1(831:835) .EQ. '') EW73(i__) = rNULL__
if(ar__1(837:840) .EQ. '') A73(i__) = rNULL__
if(ar__1(842:846) .EQ. '') EW74(i__) = rNULL__
if(ar__1(848:851) .EQ. '') A74(i__) = rNULL__
if(ar__1(853:857) .EQ. '') EW75(i__) = rNULL__
if(ar__1(859:862) .EQ. '') A75(i__) = rNULL__
if(ar__1(864:868) .EQ. '') EW76(i__) = rNULL__
if(ar__1(870:873) .EQ. '') A76(i__) = rNULL__
if(ar__1(875:879) .EQ. '') EW77(i__) = rNULL__
if(ar__1(881:884) .EQ. '') A77(i__) = rNULL__
if(ar__1(886:890) .EQ. '') EW78(i__) = rNULL__
if(ar__1(892:895) .EQ. '') A78(i__) = rNULL__
if(ar__1(897:901) .EQ. '') EW79(i__) = rNULL__
if(ar__1(903:906) .EQ. '') A79(i__) = rNULL__
if(ar__1(908:912) .EQ. '') EW80(i__) = rNULL__
if(ar__1(914:917) .EQ. '') A80(i__) = rNULL__
if(ar__1(919:923) .EQ. '') EW81(i__) = rNULL__
if(ar__1(925:928) .EQ. '') A81(i__) = rNULL__
if(ar__1(930:934) .EQ. '') EW82(i__) = rNULL__
if(ar__1(936:939) .EQ. '') A82(i__) = rNULL__
if(ar__1(941:945) .EQ. '') EW83(i__) = rNULL__
if(ar__1(947:950) .EQ. '') A83(i__) = rNULL__
if(ar__1(952:956) .EQ. '') EW84(i__) = rNULL__
if(ar__1(958:961) .EQ. '') A84(i__) = rNULL__
if(ar__1(963:967) .EQ. '') EW85(i__) = rNULL__
if(ar__1(969:972) .EQ. '') A85(i__) = rNULL__
if(ar__1(974:978) .EQ. '') EW86(i__) = rNULL__
if(ar__1(980:983) .EQ. '') A86(i__) = rNULL__
if(ar__1(985:989) .EQ. '') EW87(i__) = rNULL__
if(ar__1(991:994) .EQ. '') A87(i__) = rNULL__
if(ar__1(996:1000) .EQ. '') EW88(i__) = rNULL__
if(ar__1(1002:1005) .EQ. '') A88(i__) = rNULL__
if(ar__1(1007:1011) .EQ. '') EW89(i__) = rNULL__
if(ar__1(1013:1016) .EQ. '') A89(i__) = rNULL__
if(ar__1(1018:1022) .EQ. '') EW90(i__) = rNULL__
if(ar__1(1024:1027) .EQ. '') A90(i__) = rNULL__
if(ar__1(1029:1033) .EQ. '') EW91(i__) = rNULL__
if(ar__1(1035:1038) .EQ. '') A91(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Element(i__),Ion(i__),Lambda(i__),LEP(i__),loggf(i__),
+ EWSun(i__),ASun(i__),EW1(i__),A1(i__),EW2(i__),A2(i__),
+ EW3(i__),A3(i__),EW4(i__),A4(i__),EW5(i__),A5(i__),EW6(i__),
+ A6(i__),EW7(i__),A7(i__),EW8(i__),A8(i__),EW9(i__),A9(i__),
+ EW10(i__),A10(i__),EW11(i__),A11(i__),EW12(i__),A12(i__),
+ EW13(i__),A13(i__),EW14(i__),A14(i__),EW15(i__),A15(i__),
+ EW16(i__),A16(i__),EW17(i__),A17(i__),EW18(i__),A18(i__),
+ EW19(i__),A19(i__),EW20(i__),A20(i__),EW21(i__),A21(i__),
+ EW22(i__),A22(i__),EW23(i__),A23(i__),EW24(i__),A24(i__),
+ EW25(i__),A25(i__),EW26(i__),A26(i__),EW27(i__),A27(i__),
+ EW28(i__),A28(i__),EW29(i__),A29(i__),EW30(i__),A30(i__),
+ EW31(i__),A31(i__),EW32(i__),A32(i__),EW33(i__),A33(i__),
+ EW34(i__),A34(i__),EW35(i__),A35(i__),EW36(i__),A36(i__),
+ EW37(i__),A37(i__),EW38(i__),A38(i__),EW39(i__),A39(i__),
+ EW40(i__),A40(i__),EW41(i__),A41(i__),EW42(i__),A42(i__),
+ EW43(i__),A43(i__),EW44(i__),A44(i__),EW45(i__),A45(i__),
+ EW46(i__),A46(i__),EW47(i__),A47(i__),EW48(i__),A48(i__),
+ EW49(i__),A49(i__),EW50(i__),A50(i__),EW51(i__),A51(i__),
+ EW52(i__),A52(i__),EW53(i__),A53(i__),EW54(i__),A54(i__),
+ EW55(i__),A55(i__),EW56(i__),A56(i__),EW57(i__),A57(i__),
+ EW58(i__),A58(i__),EW59(i__),A59(i__),EW60(i__),A60(i__),
+ EW61(i__),A61(i__),EW62(i__),A62(i__),EW63(i__),A63(i__),
+ EW64(i__),A64(i__),EW65(i__),A65(i__),EW66(i__),A66(i__),
+ EW67(i__),A67(i__),EW68(i__),A68(i__),EW69(i__),A69(i__),
+ EW70(i__),A70(i__),EW71(i__),A71(i__),EW72(i__),A72(i__),
+ EW73(i__),A73(i__),EW74(i__),A74(i__),EW75(i__),A75(i__),
+ EW76(i__),A76(i__),EW77(i__),A77(i__),EW78(i__),A78(i__),
+ EW79(i__),A79(i__),EW80(i__),A80(i__),EW81(i__),A81(i__),
+ EW82(i__),A82(i__),EW83(i__),A83(i__),EW84(i__),A84(i__),
+ EW85(i__),A85(i__),EW86(i__),A86(i__),EW87(i__),A87(i__),
+ EW88(i__),A88(i__),EW89(i__),A89(i__),EW90(i__),A90(i__),
+ EW91(i__),A91(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end