Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/394/187 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-19
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/394/187 {gamma} Ser & 101 Her elemental abundances (Caliskan+, 2002)
*================================================================================
*Elemental abundance analyses with DAO spectrograms. XXVI.
*The superficially normal stars {gamma} Serpentic (F6 V) and 101 Herculis (A7 V).
* Caliskan H., Adelman S.J., Cay M.T., Gulliver A.F., Tektunali H.G., Teker A.
* <Astron. Astrophys. 394, 187 (2002)>
* =2002A&A...394..187C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3a.dat' ! Elemental abundances
integer*4 nr__
parameter (nr__=850) ! Number of records
character*50 ar__ ! Full-size record
character*5 Ion (nr__) ! Ion
character*6 No (nr__) ! Multiplet number from Moore (1945) if first
* line of multiplet, if not from Moore an
* initial indicating the paper (1)
real*8 Lambda (nr__) ! (0.1nm) Laboratory wavelength (air)
real*4 log_gf (nr__) ! Oscillator strength
character*2 r_log_gf (nr__) ! Reference to oscillator strength (2)
integer*4 EW1 (nr__) ! (0.1pm) ? Equivalent width for {gamma} Ser
real*4 log_N_Ntot_1(nr__) ! ? Elemental abundance for {gamma} Ser (3)
integer*4 EW2 (nr__) ! (0.1pm) ? Equivalent width for 101 Her
real*4 log_N_Ntot_2(nr__) ! ? Elemental abundance for 101 Her (3)
*Note (1): the acronyms are:
* C: Catalan et al., 1964, J. Res. NBS, 68A, 9
* J: Johansson S., 1978, Phys. Scripta, 18, 217
* K: Kurucz & Bell 1995, Kurucz CD-Rom No. 23
* MCS: Meggers et al., 1975, Tables of Spectral-Line Intensities, US
* Goverment Printing Office, Washington, DC
* N: Nave et al., 1994, Cat. <J/ApJS/94/221>
*Note (2): Sources of gf-values:
* BG: Biemont et al. (1989A&A...209..391B) for V II;
* Biemont et al. (1981ApJ...248..867B) for Zr II
* BK: Biemont et al. (1982A&A...107..166B)
* FW: Fuhr & Wiese (1990, in Lide D.R. ed., CRC Handbook of Chemistry
* and Physics, CRC Press, Cleveland, OH)
* HL: Hannaford et al. (1982ApJ...261..736H)
* JK: Jonsson et al. (1984, Phys. Rev. A, 30, 2429)
* KX: Kurucz & Bell (1995, Atomic Data for Opacity Calculations, Kurucz
* CD-Rom No. 23, (Smithsonian Astrophysical Observatory, Cambridge, MA)
* LA: Lanz & Artru (1985PhyS...32..115L)
* LD: Lawler & Dakin (1989, JOSA B, 6, 1457)
* MF: Fuhr, Martin & Wiese (1988, Cat. <VI/72>) &
* Martin, Fuhr & Wiese (1988, Cat. <VI/72>)
* MC: Magazzu & Cowley (1986ApJ...308..254M)
* MW: Miles & Wiese (1969, NBS Technical Note No. 474)
* YF: Younger et al. (1978, J. Phys. Chem. Ref. Data, 7, 495)
* WF: Wiese, Fuhr & Deters (1996, J. Phys. Chem. Ref. Data, Monograph 6)
* WM: Wiese & Martin (1980, NSRDS-NBS 68, Part.2, US Governement
* Printing Office, Washington)
* WS Wiese, Smith & Glennon (1966, NSRDS-NBS 4, US Governement
* Printing Office, Washington) &
* Wiese, Smith & Miles (1969, NSRDS-NBS, D.C.: US Department of
* Commerce, National Bureau of Standards)
* WV Ward (1985MNRAS.213...71W)
*Note (3): the designations are:
* N = number of atoms/volume of particular element
* Ntot = number of atoms/volume of all elements
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3b.dat' ! Relative abundances
integer*4 nr__1
parameter (nr__1=35) ! Number of records
character*33 ar__1 ! Full-size record
character*5 Ion_1 (nr__1) ! Ion
character*2 El (nr__1) ! Element
real*4 log_El_Ntot_1(nr__1) ! ? Elemental abundance for {gamma} Ser (1)
real*4 e_log_El_Ntot_1(nr__1) ! ? rms uncertainty on log(El/Ntot)1
real*4 log_El_Ntot_2(nr__1) ! ? Elemental abundance for 101 Her (1)
real*4 e_log_El_Ntot_2(nr__1) ! ? rms uncertainty on log(El/Ntot)2
*Note (1): El = Element
* Ntot = number of atoms/volume of all elements
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Comparison of derived and solar abundances
integer*4 nr__2
parameter (nr__2=35) ! Number of records
character*58 ar__2 ! Full-size record
character*5 Ion_2 (nr__2) ! Ion
real*4 logN_H1 (nr__2) ! ([Sun]) ? Elemental abundance for 99 Her (1)
real*4 logN_H2 (nr__2) ! ([Sun]) ? Elemental abundance for {gamma} Ser (1)
real*4 e_logN_H2 (nr__2) ! ([Sun]) ? rms uncertainty on logN/H2
integer*4 o_logN_H1 (nr__2) ! Number of lines used to derive
* logN/H1 and logN/H2
real*4 logN_H3 (nr__2) ! ([Sun]) ? Elemental abundance for 28 And (1)
real*4 logN_H4 (nr__2) ! ([Sun]) ? Elemental abundance for 101 Her (1)
real*4 e_logN_H4 (nr__2) ! ([Sun]) ? rms uncertainty on logN/H4
integer*4 o_logN_H3 (nr__2) ! ? Number of lines used to derive
* logN/H3 and logN/H4
real*4 Sun (nr__2) ! ([Sun]) Solar elemental abundance
*Note (1): Effective temperatures:
* 99 Her: Teff=6100K
* 28 And: Teff=7350K
* 101 Her: Teff=8061K
* {gamma} Ser: Teff=6300K
C=============================================================================
C Loading file 'table3a.dat' ! Elemental abundances
C Format for file interpretation
1 format(A5,A6,1X,F7.2,1X,F5.2,1X,A2,1X,I3,1X,F6.2,1X,I3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3a.dat')
write(6,*) '....Loading file: table3a.dat'
do i__=1,850
read(1,'(A50)')ar__
read(ar__,1)
+ Ion(i__),No(i__),Lambda(i__),log_gf(i__),r_log_gf(i__),
+ EW1(i__),log_N_Ntot_1(i__),EW2(i__),log_N_Ntot_2(i__)
if(ar__(30:32) .EQ. '') EW1(i__) = iNULL__
if(ar__(34:39) .EQ. '') log_N_Ntot_1(i__) = rNULL__
if(ar__(41:43) .EQ. '') EW2(i__) = iNULL__
if(ar__(45:50) .EQ. '') log_N_Ntot_2(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Ion(i__),No(i__),Lambda(i__),log_gf(i__),r_log_gf(i__),
+ EW1(i__),log_N_Ntot_1(i__),EW2(i__),log_N_Ntot_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3b.dat' ! Relative abundances
C Format for file interpretation
2 format(A5,2X,A2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3b.dat')
write(6,*) '....Loading file: table3b.dat'
do i__=1,35
read(1,'(A33)')ar__1
read(ar__1,2)
+ Ion_1(i__),El(i__),log_El_Ntot_1(i__),e_log_El_Ntot_1(i__),
+ log_El_Ntot_2(i__),e_log_El_Ntot_2(i__)
if(ar__1(11:16) .EQ. '') log_El_Ntot_1(i__) = rNULL__
if(ar__1(18:21) .EQ. '') e_log_El_Ntot_1(i__) = rNULL__
if(ar__1(23:28) .EQ. '') log_El_Ntot_2(i__) = rNULL__
if(ar__1(30:33) .EQ. '') e_log_El_Ntot_2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Ion_1(i__),El(i__),log_El_Ntot_1(i__),e_log_El_Ntot_1(i__),
+ log_El_Ntot_2(i__),e_log_El_Ntot_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Comparison of derived and solar abundances
C Format for file interpretation
3 format(
+ A5,1X,F6.2,1X,F6.2,1X,F4.2,1X,I3,1X,F6.2,1X,F6.2,1X,F4.2,1X,
+ I3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,35
read(1,'(A58)')ar__2
read(ar__2,3)
+ Ion_2(i__),logN_H1(i__),logN_H2(i__),e_logN_H2(i__),
+ o_logN_H1(i__),logN_H3(i__),logN_H4(i__),e_logN_H4(i__),
+ o_logN_H3(i__),Sun(i__)
if(ar__2(7:12) .EQ. '') logN_H1(i__) = rNULL__
if(ar__2(14:19) .EQ. '') logN_H2(i__) = rNULL__
if(ar__2(21:24) .EQ. '') e_logN_H2(i__) = rNULL__
if(ar__2(30:35) .EQ. '') logN_H3(i__) = rNULL__
if(ar__2(37:42) .EQ. '') logN_H4(i__) = rNULL__
if(ar__2(44:47) .EQ. '') e_logN_H4(i__) = rNULL__
if(ar__2(49:51) .EQ. '') o_logN_H3(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Ion_2(i__),logN_H1(i__),logN_H2(i__),e_logN_H2(i__),
+ o_logN_H1(i__),logN_H3(i__),logN_H4(i__),e_logN_H4(i__),
+ o_logN_H3(i__),Sun(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end