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