Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/471/625 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-16
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/471/625 VLT-FLAMES survey of massive stars (Trundle+, 2007)
*================================================================================
*The VLT-FLAMES survey of massive stars: Evolution of surface N abundances and
*effective temperature scales in the Galaxy and Magellanic Clouds.
* Trundle C., Dufton P.L., Hunter I., Evans C.J., Lennon D.J., Smartt S.J.,
* Ryans R.S.I.
* <Astron. Astrophys., 471, 625-643 (2007)>
* =2007A&A...471..625T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'stars.dat' ! Star positions, from Simbad
integer*4 nr__
parameter (nr__=77) ! Number of records
character*79 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*11 Name (nr__) ! Star identification
integer*4 RAh (nr__) ! (h) Simbad Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Simbad Right ascension (J2000)
real*4 RAs (nr__) ! (s) Simbad Right ascension (J2000)
character*1 DE_ (nr__) ! Simbad Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Simbad Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Simbad Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Simbad Declination (J2000)
real*4 Bmag (nr__) ! (mag) Simbad B magnitude
real*4 Vmag (nr__) ! (mag) Simbad V magnitude
character*10 SpType (nr__) ! Simbad MK spectral type
character*20 SName (nr__) ! Simbad name
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Mean absolute abundances using initial atmospheric
parameters
integer*4 nr__1
parameter (nr__1=61) ! Number of records
character*105 ar__1 ! Full-size record
character*11 Name_1 (nr__1) ! Star Identification (G1)
character*11 SpType_1 (nr__1) ! MK Spectral Classification
real*4 CII (nr__1) ! mean carbon abundance (G2)
real*4 NII (nr__1) ! mean nitrogen abundance (G2)
real*4 e_NII (nr__1) ! Uncertainty in nitrogen abundance (G3)
integer*4 o_NII (nr__1) ! Number of lines used for mean N abundance
real*4 OII (nr__1) ! mean oxygen abundance (G2)
real*4 e_OII (nr__1) ! Uncertainty in oxygen abundance (G3)
integer*4 o_OII (nr__1) ! Number of lines used for mean O abundance
real*4 MgII (nr__1) ! mean magnesium abundance (G2)
real*4 e_MgII (nr__1) ! Uncertainty in magnesium abundance (G3)
real*4 SiII (nr__1) ! ? mean Si II abundance (G2)
real*4 e_SiII (nr__1) ! ? Uncertainty in Si II abundance (G3)
real*4 SiIII (nr__1) ! mean Si III abundance (G2)
real*4 e_SiIII (nr__1) ! Uncertainty in Si III abundance (G3)
real*4 SiIV (nr__1) ! ? mean Si IV abundance (G2)
real*4 e_SiIV (nr__1) ! ? Uncertainty in Si IV abundance (G3)
real*4 FeIII (nr__1) ! ? mean iron abundance (G2)
real*4 e_FeIII (nr__1) ! ? Uncertainty in iron abundance (G3)
integer*4 o_FeIII (nr__1) ! ? Number of lines used for mean iron abundance
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Mean absolute abundances adjusting microturbulence
integer*4 nr__2
parameter (nr__2=61) ! Number of records
character*106 ar__2 ! Full-size record
character*11 Name_2 (nr__2) ! Star Identification (G1)
character*11 SpType_2 (nr__2) ! MK Spectral Classification
character*1 l_Vtc (nr__2) ! Limit flag on Vtc
integer*4 Vtc (nr__2) ! (km/s) Corrected microturbulence to give mean Si
character*1 l_Vturb (nr__2) ! Limit flag on Vturb
integer*4 Vturb (nr__2) ! (km/s) Initial microturbulence
real*4 CIIa (nr__2) ! mean carbon abundance (G2)
real*4 NIIa (nr__2) ! mean nitrogen abundance (G2)
real*4 e_NIIa (nr__2) ! Uncertainty in nitrogen abundance (G3)
real*4 OIIa (nr__2) ! mean oxygen abundance (G2)
real*4 e_OIIa (nr__2) ! Uncertainty in oxygen abundance (G3)
real*4 MgIIa (nr__2) ! mean magnesium abundance (G2)
real*4 e_MgIIa (nr__2) ! Uncertainty in magnesium abundance (G3)
real*4 SiIIa (nr__2) ! ? mean Si II abundance (G2)
real*4 e_SiIIa (nr__2) ! ? Uncertainty in Si II abundance (G3)
real*4 SiIIIa (nr__2) ! mean Si III abundance (G2)
real*4 e_SiIIIa (nr__2) ! Uncertainty in Si III abundance (G3)
real*4 SiIVa (nr__2) ! ? mean Si IV abundance (G2)
real*4 e_SiIVa (nr__2) ! ? Uncertainty in Si IV abundance (G3)
real*4 FeIIIa (nr__2) ! ? mean iron abundance (G2)
real*4 e_FeIIIa (nr__2) ! ? Uncertainty in iron abundance (G3)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Fe III Equivalent Widths and abundances in the young
clusters from the VLT-FLAMES survey of massive stars
integer*4 nr__3
parameter (nr__3=16) ! Number of records
character*37 ar__3 ! Full-size record
character*11 Name_3 (nr__3) ! Star Identification (G1)
integer*4 EWFe4419 (nr__3) ! (0.1pm) Equivalent width of FeIII 4419 line
real*4 AFe4419 (nr__3) ! Absolute abundance of FeIII 4419 line (G2)
integer*4 EWFe4430 (nr__3) ! (0.1pm) ? Equivalent width of FeIII 4430 line
real*4 AFe4430 (nr__3) ! ? Absolute abundance of FeIII 4430 line (G2)
real*4 AFeIII (nr__3) ! Mean FeIII abundance (G2)
real*4 e_AFeIII (nr__3) ! Uncertainty in FeIII abundance (G3)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec1.dat' ! NGC 3293 line by line equivalent widths & abundances
integer*4 nr__4
parameter (nr__4=424) ! Number of records
character*53 ar__4 ! Full-size record
character*11 Name_4 (nr__4) ! Star Identification (G1)
character*8 Ion (nr__4) ! Element/Ion (1)
real*8 lambda (nr__4) ! (0.1nm) Wavelength of line (1)
integer*4 EW (nr__4) ! (0.1pm) ? Equivalent width of line
real*4 Ab1 (nr__4) ! ? Abundance using initial parameters (G2)(2)
real*4 Ab2 (nr__4) ! ? Abundance with corrected microturbulence (G2)(2)
real*4 Ab3 (nr__4) ! ? Abundance using corrected parameters (G2)(2)
*Note (1): The energy levels and oscillator strengths relating to the
* transitions for the metal lines are available online at
* http://star.pst.qub.ac.uk/~pld/line\_identifications.html.
*Note (2): The abundances are estimated from the lines and they show the results
* from the 3 steps in determining the abundances as described in
* Sect. 4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec2.dat' ! NGC 4755 line by line equivalent widths & abundances
integer*4 nr__5
parameter (nr__5=530) ! Number of records
character*53 ar__5 ! Full-size record
character*11 Name_5 (nr__5) ! Star Identification (G1)
character*8 Ion_1 (nr__5) ! Element/Ion (1)
real*8 lambda_1 (nr__5) ! (0.1nm) Wavelength of line (1)
integer*4 EW_1 (nr__5) ! (0.1pm) ? Equivalent width of line
real*4 Ab1_1 (nr__5) ! ? Abundance using initial parameters (G2)(2)
real*4 Ab2_1 (nr__5) ! ? Abundance with corrected microturbulence (G2)(2)
real*4 Ab3_1 (nr__5) ! ? Abundance using corrected parameters (G2)(2)
*Note (1): The energy levels and oscillator strengths relating to the
* transitions for the metal lines are available online at
* http://star.pst.qub.ac.uk/~pld/line\_identifications.html.
*Note (2): The abundances are estimated from the lines and they show the results
* from the 3 steps in determining the abundances as described in
* Sect. 4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec3.dat' ! NGC 2004 line by line equivalent widths & abundances
integer*4 nr__6
parameter (nr__6=1219) ! Number of records
character*53 ar__6 ! Full-size record
character*11 Name_6 (nr__6) ! Star Identification (G1)
character*8 Ion_2 (nr__6) ! Element/Ion (1)
real*8 lambda_2 (nr__6) ! (0.1nm) Wavelength of line (1)
integer*4 EW_2 (nr__6) ! (0.1pm) ? Equivalent width of line
real*4 Ab1_2 (nr__6) ! ? Abundance using initial parameters (G2)(2)
real*4 Ab2_2 (nr__6) ! ? Abundance with corrected microturbulence (G2)(2)
real*4 Ab3_2 (nr__6) ! ? Abundance using corrected parameters (G2)(2)
*Note (1): The energy levels and oscillator strengths relating to the
* transitions for the metal lines are available online at
* http://star.pst.qub.ac.uk/~pld/line\_identifications.html.
*Note (2): The abundances are estimated from the lines and they show the results
* from the 3 steps in determining the abundances as described in
* Sect. 4 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec4.dat' ! NGC 330 line by line equivalent widths & abundances
integer*4 nr__7
parameter (nr__7=1060) ! Number of records
character*53 ar__7 ! Full-size record
character*11 Name_7 (nr__7) ! Star Identification (G1)
character*8 Ion_3 (nr__7) ! Element/Ion (1)
real*8 lambda_3 (nr__7) ! (0.1nm) Wavelength of line (1)
integer*4 EW_3 (nr__7) ! (0.1pm) ? Equivalent width of line
real*4 Ab1_3 (nr__7) ! ? Abundance using initial parameters (G2)(2)
real*4 Ab2_3 (nr__7) ! ? Abundance with corrected microturbulence (G2)(2)
real*4 Ab3_3 (nr__7) ! ? Abundance using corrected parameters (G2)(2)
*Note (1): The energy levels and oscillator strengths relating to the
* transitions for the metal lines are available online at
* http://star.pst.qub.ac.uk/~pld/line\_identifications.html.
*Note (2): The abundances are estimated from the lines and they show the results
* from the 3 steps in determining the abundances as described in
* Sect. 4 of the paper.
C=============================================================================
C Loading file 'stars.dat' ! Star positions, from Simbad
C Format for file interpretation
1 format(
+ A11,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,
+ F5.2,1X,A10,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'stars.dat')
write(6,*) '....Loading file: stars.dat'
do i__=1,77
read(1,'(A79)')ar__
read(ar__,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Bmag(i__),Vmag(i__),SpType(i__),SName(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__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Bmag(i__),Vmag(i__),SpType(i__),SName(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' ! Mean absolute abundances using initial atmospheric
* parameters
C Format for file interpretation
2 format(
+ A11,1X,A11,1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,1X,F4.2,1X,F4.2,1X,
+ I2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,61
read(1,'(A105)')ar__1
read(ar__1,2)
+ Name_1(i__),SpType_1(i__),CII(i__),NII(i__),e_NII(i__),
+ o_NII(i__),OII(i__),e_OII(i__),o_OII(i__),MgII(i__),
+ e_MgII(i__),SiII(i__),e_SiII(i__),SiIII(i__),e_SiIII(i__),
+ SiIV(i__),e_SiIV(i__),FeIII(i__),e_FeIII(i__),o_FeIII(i__)
if(ar__1(65:68) .EQ. '') SiII(i__) = rNULL__
if(ar__1(70:73) .EQ. '') e_SiII(i__) = rNULL__
if(ar__1(85:88) .EQ. '') SiIV(i__) = rNULL__
if(ar__1(90:93) .EQ. '') e_SiIV(i__) = rNULL__
if(ar__1(95:98) .EQ. '') FeIII(i__) = rNULL__
if(ar__1(100:103) .EQ. '') e_FeIII(i__) = rNULL__
if(ar__1(105:105) .EQ. '') o_FeIII(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),SpType_1(i__),CII(i__),NII(i__),e_NII(i__),
+ o_NII(i__),OII(i__),e_OII(i__),o_OII(i__),MgII(i__),
+ e_MgII(i__),SiII(i__),e_SiII(i__),SiIII(i__),e_SiIII(i__),
+ SiIV(i__),e_SiIV(i__),FeIII(i__),e_FeIII(i__),o_FeIII(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Mean absolute abundances adjusting microturbulence
C Format for file interpretation
3 format(
+ A11,1X,A11,1X,A1,I2,1X,A1,I2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,61
read(1,'(A106)')ar__2
read(ar__2,3)
+ Name_2(i__),SpType_2(i__),l_Vtc(i__),Vtc(i__),l_Vturb(i__),
+ Vturb(i__),CIIa(i__),NIIa(i__),e_NIIa(i__),OIIa(i__),
+ e_OIIa(i__),MgIIa(i__),e_MgIIa(i__),SiIIa(i__),e_SiIIa(i__),
+ SiIIIa(i__),e_SiIIIa(i__),SiIVa(i__),e_SiIVa(i__),FeIIIa(i__),
+ e_FeIIIa(i__)
if(ar__2(68:71) .EQ. '') SiIIa(i__) = rNULL__
if(ar__2(73:76) .EQ. '') e_SiIIa(i__) = rNULL__
if(ar__2(88:91) .EQ. '') SiIVa(i__) = rNULL__
if(ar__2(93:96) .EQ. '') e_SiIVa(i__) = rNULL__
if(ar__2(98:101) .EQ. '') FeIIIa(i__) = rNULL__
if(ar__2(103:106) .EQ. '') e_FeIIIa(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),SpType_2(i__),l_Vtc(i__),Vtc(i__),l_Vturb(i__),
+ Vturb(i__),CIIa(i__),NIIa(i__),e_NIIa(i__),OIIa(i__),
+ e_OIIa(i__),MgIIa(i__),e_MgIIa(i__),SiIIa(i__),e_SiIIa(i__),
+ SiIIIa(i__),e_SiIIIa(i__),SiIVa(i__),e_SiIVa(i__),FeIIIa(i__),
+ e_FeIIIa(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Fe III Equivalent Widths and abundances in the young
* clusters from the VLT-FLAMES survey of massive stars
C Format for file interpretation
4 format(A11,1X,I2,1X,F4.2,1X,I2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,16
read(1,'(A37)')ar__3
read(ar__3,4)
+ Name_3(i__),EWFe4419(i__),AFe4419(i__),EWFe4430(i__),
+ AFe4430(i__),AFeIII(i__),e_AFeIII(i__)
if(ar__3(21:22) .EQ. '') EWFe4430(i__) = iNULL__
if(ar__3(24:27) .EQ. '') AFe4430(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),EWFe4419(i__),AFe4419(i__),EWFe4430(i__),
+ AFe4430(i__),AFeIII(i__),e_AFeIII(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec1.dat' ! NGC 3293 line by line equivalent widths & abundances
C Format for file interpretation
5 format(A11,2X,A8,2X,F7.2,2X,I3,2X,F4.2,1X,F5.2,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec1.dat')
write(6,*) '....Loading file: tablec1.dat'
do i__=1,424
read(1,'(A53)')ar__4
read(ar__4,5)
+ Name_4(i__),Ion(i__),lambda(i__),EW(i__),Ab1(i__),Ab2(i__),
+ Ab3(i__)
if(ar__4(33:35) .EQ. '') EW(i__) = iNULL__
if(ar__4(38:41) .EQ. '') Ab1(i__) = rNULL__
if(ar__4(43:47) .EQ. '') Ab2(i__) = rNULL__
if(ar__4(50:53) .EQ. '') Ab3(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),Ion(i__),lambda(i__),EW(i__),Ab1(i__),Ab2(i__),
+ Ab3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec2.dat' ! NGC 4755 line by line equivalent widths & abundances
C Format for file interpretation
6 format(A11,2X,A8,2X,F7.2,2X,I3,2X,F4.2,1X,F5.2,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec2.dat')
write(6,*) '....Loading file: tablec2.dat'
do i__=1,530
read(1,'(A53)')ar__5
read(ar__5,6)
+ Name_5(i__),Ion_1(i__),lambda_1(i__),EW_1(i__),Ab1_1(i__),
+ Ab2_1(i__),Ab3_1(i__)
if(ar__5(33:35) .EQ. '') EW_1(i__) = iNULL__
if(ar__5(38:41) .EQ. '') Ab1_1(i__) = rNULL__
if(ar__5(43:47) .EQ. '') Ab2_1(i__) = rNULL__
if(ar__5(50:53) .EQ. '') Ab3_1(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5(i__),Ion_1(i__),lambda_1(i__),EW_1(i__),Ab1_1(i__),
+ Ab2_1(i__),Ab3_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec3.dat' ! NGC 2004 line by line equivalent widths & abundances
C Format for file interpretation
7 format(A11,2X,A8,2X,F7.2,2X,I3,2X,F4.2,1X,F5.2,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec3.dat')
write(6,*) '....Loading file: tablec3.dat'
do i__=1,1219
read(1,'(A53)')ar__6
read(ar__6,7)
+ Name_6(i__),Ion_2(i__),lambda_2(i__),EW_2(i__),Ab1_2(i__),
+ Ab2_2(i__),Ab3_2(i__)
if(ar__6(33:35) .EQ. '') EW_2(i__) = iNULL__
if(ar__6(38:41) .EQ. '') Ab1_2(i__) = rNULL__
if(ar__6(43:47) .EQ. '') Ab2_2(i__) = rNULL__
if(ar__6(50:53) .EQ. '') Ab3_2(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Name_6(i__),Ion_2(i__),lambda_2(i__),EW_2(i__),Ab1_2(i__),
+ Ab2_2(i__),Ab3_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablec4.dat' ! NGC 330 line by line equivalent widths & abundances
C Format for file interpretation
8 format(A11,2X,A8,2X,F7.2,2X,I3,2X,F4.2,1X,F5.2,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablec4.dat')
write(6,*) '....Loading file: tablec4.dat'
do i__=1,1060
read(1,'(A53)')ar__7
read(ar__7,8)
+ Name_7(i__),Ion_3(i__),lambda_3(i__),EW_3(i__),Ab1_3(i__),
+ Ab2_3(i__),Ab3_3(i__)
if(ar__7(33:35) .EQ. '') EW_3(i__) = iNULL__
if(ar__7(38:41) .EQ. '') Ab1_3(i__) = rNULL__
if(ar__7(43:47) .EQ. '') Ab2_3(i__) = rNULL__
if(ar__7(50:53) .EQ. '') Ab3_3(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ Name_7(i__),Ion_3(i__),lambda_3(i__),EW_3(i__),Ab1_3(i__),
+ Ab2_3(i__),Ab3_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end