FORTRAN Generation
(/./ftp/cats/J/AN/329/4)

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.8, on 2014-Oct-22
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