Conversion of standardized ReadMe file for
file /./ftp/cats/J/AN/329/4 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.5, on 2013-May-20
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/AN/329/4 Elemental abundance analyses of nu Her & 41 Cyg (Adelman+, 2008)
*================================================================================
*Elemental abundance analyses with DAO spectrograms XXXI.
*The early F supergiants nu Her (F2 II) and 41 Cyg (F5 Ib-II).
* Adelman S.J., Cay I.H., Tektunali H.G., Gulliver A.F., Teker A.
* <Astron. Nachrichten, 329, 4 (2008)>
* =2008AN....329....4A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'abund.dat' ! Atomic abundances
integer*4 nr__
parameter (nr__=43) ! Number of records
character*34 ar__ ! Full-size record
character*6 Ion ! Ion
character*2 Element ! Element (identical to bytes 1-2)
real*4 log_Nel_Ntot_nu ! ?=- nu Her derived abundance of element (1)
real*4 e_log_Nel_Ntot_nu ! ? rms uncertainty on log(Nel/Ntot)
real*4 log_Nel_Ntot_41 ! ?=- 41 Cyg derived abundance of element (1)
real*4 e_log_Nel_Ntot_41 ! ? rms uncertainty on log(Nel/Ntot)
*Note (1): Ratio of number of atoms to the total number of atoms per unit volume
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Abundances for nu Her & 41 Cyg
integer*4 nr__1
parameter (nr__1=1910) ! Number of records
character*55 ar__1 ! Full-size record
character*5 Ion_1 ! Ion
character*7 Mult ! Multiplet number (1)
real*8 Lambda ! (0.1nm) Laboratory wavelength in air
real*4 loggf ! Logarithm of the oscillator strength
character*2 r_loggf ! References to the oscillator strength (2)
integer*4 EWnu ! (0.1pm) ?=- Equivalent width value for nu Her
real*4 log_Nel_Ntot_nu_1 ! ?=- Relative abundance for nu Her
integer*4 EW41 ! (0.1pm) ?=- Equivalent width value for 41 Cyg
real*4 log_Nel_Ntot_41_1 ! ?=- Relative abundance for 41 Cyg
*Note (1): Multiplet number from Moore (1945, A Multiplet Table of
* Astrophysical Interest, Princeton University Observatory, Princeton)
* on the first of the multiple. If not from Moore an initial(s)
* indicating the paper or the number from a later Multiple Table by
* Moore as follows:
* C = Catalan et al. (1964, J. Res. NBS, 68A, 9)
* D = Dworetsky (1971, Ph. D. thesis, University of California at
* Los Angeles)
* I = Iglesias et al. (1988, Publ. del Instituto de Optica de Madrid
* No.47) for V II and Iglesias & Velasco (1964, Publ. Del
* Instituto de Optica de Madrid No.23) for n II
* J = Johansson (1978, PhS, 18, 217),
* K = Kurucz & Bell (1995, Atomic Data for Opacity Calculations,
* Kurucz CD-Rom No. 23, Smithsonian Astrophysical Observatory,
* Cambridge, MA),
* MCS = Meggers, Corliss, & Scribner (1975, Tables of Spectral-Line
* Intensities Part I - Arranged by Elements, 2nd edition,
* (Washington, DC: Government Printing Office)
* N = Nave et al. (1994, Cat. <J/ApJS/94/221>) for Fe I and Nilsson et
* al. (1991, PhS, 44, 226) for Y II.
* For C I, O I, Si I, and Si II we used the more recent multiplet
* designations from Moore's studies (1965, NSRDS-NBS 3, Section 1
* (Washington, DC: US Government Printing Office), 1967, NSRDS-NBS 3,
* Section 2 (Washington, DC: US Government Printing Office), and 1993,
* Tables of Spectra of Hydrogen, Carbon, Nitrogen, and Oxygen Atoms
* and Ions, CRC Series in Evaluated Data in Atomic Physics, ed. J.W.
* Gallagher (Boca Raton, FL: CRC Press Inc.)
*Note (2): gf value references as follow:
* BG = Biemont et al. (1989A&A...209..391B) for V II and
* Biemont et al. (1981ApJ...248..867B) for Zr II
* BB = Blackwell-Whitehead & Bergemann (2007A&A...472L..43B)
* BK = Biemont et al. (1982A&A...107..166B)
* BL = Biemont et al. (2003, Eur. Phys. J. D, 27, 33)
* BR = Biemont et al. (1981ApJ...248..867B)
* DE = Den Hartog et al. (2003ApJS..148..543D)
* DL = Den Hartog et al. (2006ApJS..167..292D)
* FW = Fuhr & Wiese (1996, NIST Atomic Probability Tables, in Handbook
* of Chemistry and Physics, ed. David R. Lide, (Boca Raton: CRC
* Handbook of Chemistry and Physics, 10-128)
* GB = Grevesse et al. (1981, Upper Main Sequence Stars, 23rd Liege
* Astrophys. Colloq., Universite de Liege, p. 211)
* HL = Hannaford et al. (1982ApJ...261..736H)
* JK = Jonsson et al. (1984, Phys. Rev. A, 30, 2429)
* KF = Klose et al. (2002, JPCRD, 31, 217)
* KG = Kling & Griesmann (2000ApJ...531.1173K)
* KS = Kling et al. (2001ApJS..134..173K)
* KX = Kurucz & Bell (1995, Atomic Data for Opacity Calculations,
* Kurucz CD-Rom No. 23, Smithsonian Astrophysical Observatory,
* Cambridge, MA)
* LA = Lanz & Artru (1985, PhS, 32, 115)
* LB = Lawler et al. (2001ApJ...556..452L)
* LD = Lawler & Dakin (1989, JOSA B, 6, 1457)
* LN = Ljung et al. (2006A&A...456.1181L)
* LW = Lawler et al. (2001ApJ...563.1075L)
* F = Fuhr, Martin & Wiese (1988, JPCRD, 17, Suppl. 4) and
* Martin, Fuhr & Wiese (1988, JPCRD, 17, Suppl. 3)
* N4 = Fuhr & Wiese (2006, JPCRD, 35, 1669)
* PT = Pickering et al. (2001ApJS..132..403P)
* WF = Wiese, Fuhr & Deters (1996, JPCRD, Monograph 6)
* WL = Wickliffe et al. (2002, JQSRT, 66, 363)
* WM = Wiese & Martin (1980, NSRDS-NBS 68, Part 2, (Washington, DC: US
* Government Printing Office)
* WS = Wiese, Smith & Glennon (1966, NSRDS-NBS 4 (Washington, DC: US
* Government Printing Office) and
* Wiese, Smith & Miles (1969, NSRDS-NBS 22 (Washington, DC: US
* Government Printing Office)
C=============================================================================
C Loading file 'abund.dat' ! Atomic abundances
C Format for file interpretation
1 format(A6,1X,A2,1X,F6.2,1X,F4.2,2X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,file='abund.dat', status='old')
write(6,*) '....Loading file: abund.dat'
do i__=1,43
read(1,'(A34)')ar__
read(ar__,1)
+ Ion,Element,log_Nel_Ntot_nu,e_log_Nel_Ntot_nu,log_Nel_Ntot_41,
+ e_log_Nel_Ntot_41
if (log_Nel_Ntot_nu .EQ. 2.2233e-322) log_Nel_Ntot_nu = rNULL__
if(ar__(18:21) .EQ. '') e_log_Nel_Ntot_nu = rNULL__
if (log_Nel_Ntot_41 .EQ. 2.2233e-322) log_Nel_Ntot_41 = rNULL__
if(ar__(31:34) .EQ. '') e_log_Nel_Ntot_41 = rNULL__
c ..............Just test output...........
write(6,1)
+ Ion,Element,log_Nel_Ntot_nu,e_log_Nel_Ntot_nu,log_Nel_Ntot_41,
+ e_log_Nel_Ntot_41
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Abundances for nu Her & 41 Cyg
C Format for file interpretation
2 format(A5,A7,1X,F7.2,2X,F5.2,2X,A2,2X,I3,1X,F6.2,2X,I3,1X,F6.2)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,1910
read(1,'(A55)')ar__1
read(ar__1,2)
+ Ion_1,Mult,Lambda,loggf,r_loggf,EWnu,log_Nel_Ntot_nu_1,EW41,
+ log_Nel_Ntot_41_1
if (EWnu .EQ. 45) EWnu = iNULL__
if (log_Nel_Ntot_nu_1 .EQ. 2.2233e-322) log_Nel_Ntot_nu_1 = rNULL__
if (EW41 .EQ. 45) EW41 = iNULL__
if (log_Nel_Ntot_41_1 .EQ. 2.2233e-322) log_Nel_Ntot_41_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Ion_1,Mult,Lambda,loggf,r_loggf,EWnu,log_Nel_Ntot_nu_1,EW41,
+ log_Nel_Ntot_41_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end