BELFAST 
Belfast Data Base on Atomic and Molecular Physics
Department of Computer Science
The Queen's University of Belfast
Belfast BT7 1NN, U.K.

ORNL-CFADC 
Controlled Fusion Atomic Data Center
Oak Ridge National Laboratory
M.S. 6373, P.O. Box 2008
Oak Ridge, TN 38731-6372
USA

PPPL 
Princeton Plasma Physics Laboratory
P.O. Box 451
Princeton, NJ 08544 USA

RIC-NAGOYA 
Research Information Center
Institute of Plasma Physics
Nagoya University
Chikusa-ku, Nagoya 464, Japan

IAEA-AMD-UNIT 
Atomic and Molecular Data Unit
Nuclear Data Section
International Atomic Energy Agency
P.O. Box 100
Wagramer Strasse 5
A-1400 Vienna, Austria

JAERI-NDC 
Nuclear Data Center
Japanesse Atomic Energy Research Institute (JAERI)
Tokai-Mura Naka-Gun
Ibaraki-Ken 319-11, Japan

NIST-AEL-DC
Atomic Energy Levels Data Center
National Institute of Standards and Technology
United States Department of Commerce
Physics Building, Room A167
Gaithersburg
Maryland 20899, USA

NIST-TP-DC
Transition Probabilities Data Center
National Institute of Standards and Technology
United States Department of Commerce
Gaithersburg
Maryland 20899, USA

CRAAMD 
Chinese Research Association for Atomic and Molecular Data
P.O. Box 8009
Telex: 222499 IAPCM CN
Beijing
The People's Republic of China

DOC=H-HE-PLASMA 
"Elementary Processes in Hydrogen-Helium Plasmas," R.K. Janev, W.D.
Langer, K. Evans and D.E.  Post, Jr., Springer Series on Atoms and
Plasmas, Springer-Verlag (Berlin, Heidelberg, New York) 1987.

DOC=IPPJ-AM-52 
"Energy Dependence of Ion-induced Sputtering Yields of Monatomic solids
in the Low Energy Region", N. Matsunami, Y. Yamamura, N. Itoh, H.
Tawara and T. K. Kawamura, Report IPPJ-AM-52, Institute of Plasma
Physics, Nagoya University, July, 1987.

DOC=IPPJ-AM-56
"Total and Partial Cross Sections for Electron Capture for C[q+]
(q=6-2) and O[q+] (q=8-2) Ions in Collisions with H, H[2] and He
Atoms," H.  Tawara, Report IPPJ-AM-56, Institute of Plasma Physics,
Nagoya University, October, 1987.

DOC=JPCRD-12/83 
"Recommended Cross Sections and Rates for Electron Ionisation of
Light Atoms and Ions," K.L. Bell, H.B. Gilbody, J.G.  Hughes,
A.E. Kingston and F.J. Smith, J. Phys. Chem. Ref. Data 12,
891-916 (1983).

DOC=JPCRD-17/88 
"Recommended Data on the Electron Impact Ionization of Atoms and
Ions:  Fluorine to Nickel," M.A. Lennon, K.L. Bell, H.B.
Gilbody, J.G. Hughes, A.E. Kingston, M.J. Murray and F.J. Smith,
J.  Phys. Chem. Ref. Data 17, 1285-1363 (1988).

DOC=NUC-FUS-SUPP/87 
"Recommended Data on Atomic Collision Processes Involving Iron
and its Ions," C. Bottcher, D.C.  Griffin, H.T. Hunter, R.K.
Janev, A.E. Kingston, M.A. Lennon, R.A. Phaneuf, M.S. Pindzola
and S.M. Younger Special Supplement to Nuclear Fusion, 1987.

DOC=ORNL-6090 
Atomic Data for Fusion, Volume 5 "Collisions of Carbon and Oxygen Ions
with Electrons, H, H[2] and He" R.A. Phaneuf, R.K. Janev and M.S.
Pindzola Oak Ridge National Laboratory Report ORNL-6090, February,
1987.

DOC=CLM-R294/1989 
"Atomic and Molecular Data for Fusion, Part 3 "Recommended Cross
Sections and Rates for Electron Impact Ionization of Atoms and Ions:
Copper to Uranium ", M.J. Higgins, M.A. Lennon, J.G. Hughes, K.L. Bell,
H.B.  Gilbody, A.E. Kingston and F.J. Smith, Culham Laboratory Report,
CLM-R294, UKAEA, Culham Laboratory, Abingdon, Oxon. OX14 3DB, England.

DOC=ADNDT,33,149,(1985) 
"Electron-Impact Cross Sections and Rate Coefficients for Excitations of
Carbon and Oxygen Ions," Y. Itikawa, S. Hara, T. Kato, S. Nakazaki,
M. S. Pindzola, D. H. Crandall, Atomic Data and Nuclear Data Tables,
33,149,(1985).

DOC=ADNDT,42,313,(1989) 
"Recommended Data for Excitation Rate Coefficients of Helium Atoms and
Helium-like Ions by Electron Impact," T. Kato and S. Nakazaki, Atomic
Data and Nuclear Data Tables, 42,313,(1989).

DOC=ORNL-6086 
Atomic Data for Fusion, Volume 1, "Collisions of H, H{2}, He and Li
Atoms and Ions with Atoms and Molecules," C.F. Barnett, Oak Ridge
National Laboratory Report ORNL-6086, July, 1990.

DOC=IAEA-AM-REP,249,(1991)
"Particle Reflection from Surfaces - A Recommended Data Base",
E.W. Thomas, R. K. Janev , J. J. Smith, International Atomic Energy
Agency, Atomic and Molecular Data Unit Report IAEA-AM-REP,249
June 1991.

#MEWE 
GENERALIZED GAUNT FACTOR LINE EXCITATION CALCULATION
REFERENCE: R. MEWE, ASTRON. AND ASTROPHYS. 20, 215 (1972)

INPUT: ELECTRON TEMPERATURE (keV)
RETURNS: SPECTRAL LINE EXCITATION RATE COEFFICIENT (CM3/S)

#MEWE DATA TYPE COEFFICIENTS ARE (IN ORDER):

1)  EXCITATION WAVELENGTH (ANGSTROMS)
2)  EXCITATION ENERGY (EV)
    (IF PEIJ=0, ASSUME A RESONANCE LINE AND CALCULATE EIJ FROM PWVL)
3)  OSCILLATOR STRENGTH
4)  BRANCHING RATIO
5-8) A, B, C, D COEFFICIENTS FOR MEWE GAUNT FACTOR FORMULA

#TAB1D 
General one-dimensional tabular data type

Arbitary  number of (X,Y) data pairs are supplied as coefficients

    X1, Y1, X2, Y2, ...

These (x,y) data pairs must be entered in order with X increasing.
Evaluation of #TAB1D data uses linear interpolation between the
(X,Y) data pairs to return a Y value for a specified X.  If the given
X is below the first X data value or larger than the last X data
value, an error return occurs.

#CHEB 
ORNL Chebyshev polynomial fitting subroutine.  The parameter list
consists of nine polynomial coefficients followed by two additional
parameters which define the lower and upper limits for the independent
variable over which the fit is valid.  For a reference, see page 40 of
the document DOC=NUC-FUS-SUPP/87.

#KING 
IAEA subroutine to evaluate the polynomial fit used to represent
electron-impact excitation rate coefficients from Kingston et al.
(See DOC=NUC-FUS-SUPP/87, page 47 ). The first three coefficient data
values are

1. The excitation energy for the transition in eV.
2. The lower limit of validity for the fit.
3. The upper limit of validity of the fit.

The remaining data values, of maximum number 9, are the fitting
coefficients for the polynomial.

#PHACX 
IAEA subroutine to evaluate the analytic function used to represent
the charge transfer cross sections by R. A. Phaneuf et al.
(See DOC=NUC-FUS-SUPP/87, page 14 ). The coefficient data values are (in
order);

1. The lower limit of the range of validity for the scaled projectile
   energy (in eV/(u*sqrt(q), where q is the charge state of the
   projectile ion).
2. The upper limit of the range of validity for the scaled projectile
   energy (in eV/(u*sqrt(q)).
3-6. The four parameters for the analytic fit.

As the seventh parameter in the coefficient array passed to #PHACX
the charge state of the projectile ion must be passed. For the data
referenced in DOC=NUC-FUS-SUPP/87 the limits for q are +5 <= q <= +26.

#BELI
IAEA subroutine to evaluate the analytic fit for the electron-impact
ionization cross sections from Bell et al. (See DOC=JPCRD-12/83).  The
number of fitting parameters specified in the coefficient data varies
depending on the number of terms taken in the numerical fit and on
whether allowance for excitation-autoionization (E-A) has been made.
For cross sections which include E-A contributions, two seperate fits
are defined. One from the ionisation threshold to the threshold for E-A
and a second fit for energies above the E-A threshold.  (The number of
parameters (in any entry) is given in the variable KNCF).  The
coefficient data values are

For the first region

1. The ionisation potential (eV) of the target.
2-7. The fitting parameters for the fit. (A, B1, B2 etc of Bell et al)

For the second region (for cross sections with E-A only)

8. The E-A threshold energy (eV).
9. An ionisation potential (eV) for the region.
10-15. The fitting parameters for the fit. (A, B1, B2 etc of Bell et al)

#JAN1 
IAEA subroutine to evaluate the 9-term polynomial fit used in the book
of Janev et al. (see DOC=H-HE-PLASMA). The first and second parameters
of the coefficient data are the lower and upper range of validity of the
fit, for cross sections Emin to Emax and for rate coefficients Tmin to
Tmax.  This range data is followed by the nine coefficients for the
polynomial.

The Boolean labels included with this evaluation function are

a. For cross sections, values at Emin in label XS(EMIN) while for
   reaction rate coefficients the value at Tmin in label RC(Tmin).
b. For cross sections, the maximum valure of the cross section in
   label XSMAX, while for rates the maximunm rate in label RCMAX.
c. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).

#ETABSQ 
IAEA subroutine to calaculate the yields (atoms/ions) of ion-induced
sputtering of monatomic solids of N. Matsunami et al. (see
DOC=IPPJ-AM-52). The data coefficients that must be passed to #ETABSQ
are:

1. The atomic number of the incident ion
2. The atomic mass of the incident ion
3. The atomic number of the target
4. The atomic mass of the target
5. Sublimation energy of the target (in eV)
6. The Q value, defined by in IPPJ-AM- 52, for the target

The coefficient data for this evaluation function are read by the
evaluation function #RTABSQ.

#RTABSQ 
IAEA subroutine to read and interpret the entry data for the yields
(atoms/ions) of ion-induced sputtering of monatomic solids of N.
Matsunami et al. (see DOC=IPPJ-AM-52). This routine requires as input
(in order) the specification of three atomic species of the incident ion
and of the target. Ther subroutine returns the coefficient data required
by #ETABSQ to determine the sputtering yield.

#JANRD 
IAEA subroutine to evaluate the double 9-term polynomial fit used in the
book of Janev et al.  (see DOC=H-HE-PLASMA) for reaction rate
coefficients (cm[3]/s) as a function of the particle impact energy (eV)
and the plasma temperature (eV) for heavy particle reactions.  The first
and second parameters of the coefficient data are the lower and upper
range of validity for the plasma temperature (in eV) of the fit, while
the third and fourth give the range of validity for the energy of the
incident particle (in eV).  This range data is followed by the 81
coefficients for the polynomial.  The order of apperance of the
coefficiennts in the coefficient array is all 9 coefficients for each
incident particle energy, ie A(E=0,T=0-9), A(E=1,T=0-9),...  ,
A(E=9,T=0-9).  The coefficient array is temporarily updated whren
calling JANRD.  PCF(86) = E, the energy of the incident particle.  PCF
(87-95) stores temperature independent data (ln (E)) ** i , where is 0
to 9.  See DOC=H-HE-PLASMA and evaluation function for details.

The Boolean labels which may be included with this evaluation function
are

1.   The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).
2.   The error class for the fit. (see DOC=H-HE-PLASMA for meaning).

#JBORN1 
IAEA subroutine to evaluate the empirical Born-Bethe type fit Born,1
used for electron impact cross sections (cm[2]) as a function of
electron energy in the book of Janev et al., (see DOC=H-HE-PLASMA). This
fit is derived using uniform weighting in the least-squares procedure.
The coefficient data values are

1. The threshold energy for the reaction in eV.
2. The minimum energy , Emin, for application of the fitting formula.
3. The coefficient aBorn,1  used in DOC=H-HE-PLASMA.
4. The coefficient nBorn,1  used in DOC=H-HE-PLASMA.

The Boolean labels include

a. The cross section at Emin in label XS(EMIN)
b. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).

#JBORN2 
IAEA subroutine to evaluate the empirical Born-Bethe type fit Born,2
used for electron impact cross sections (cm[2]) as a function of
electron energy in the book of Janev et al., (see DOC=H-HE-PLASMA). This
fit is derived using a weighting which does not take into account the
points at energies beloew the peak in the cross section in the
least-squares procedure.  The coefficient data values are

1. The threshold energy for the reaction in eV.
2. The minimum energy , Emin, for application of the fitting formula.
3. The coefficient aBorn,2 used in DOC=H-HE-PLASMA.
4. The coefficient nBorn,2 used in DOC=H-HE-PLASMA.

The Boolean labels include

a. The cross section at Emin in label XS(EMIN)
b. The fitting error, ERROR. (see DOC=H-HE-PLASMA for meaning).

#EEXCH2 
IAEA subroutine to evaluate the analytic function used to represent the
electron impact excitation cross section (cm[2]) as a function of
electron energy (eV) of L. Vriens and H. M.  Smeets. (See
DOC=VRIENS,PRA,22,1980). The coefficient data values are

1. The principal quantum number (n) of the initial state.
2. The principal quantum number (m) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array ,

3. The threshold energy (eV) for the transition.
4. Parameter A(n,m) as quoted in the formulae of DOC=H-HE-PLASMA.
5. Parameter B(n,m) as quoted in the formulae of DOC=H-HE-PLASMA.
6. Parameter delta(n,m) as quoted in the formulae of DOC=H-HE-PLASMA.
7. Parameter gamma(n,m) as quoted in the formulae of DOC=H-HE-PLASMA.

#EIONH2 
IAEA subroutine to evaluate the analytic function used to represent the
electron impact ionization cross section (cm[2]) as a function of
electron energy (eV) of L. Vriens and H. M.  Smeets. (See
DOC=VRIENS,PRA,22,1980). The coefficient data values are

1. The principal quantum number (n) of the initial state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array ,

2. The threshold energy (eV) for the transition.

#EIONHR 
IAEA subroutine to evaluate the analytic function used to represent the
electron impact ionization rate coefficient (cm[3]/s) as a function of
plasma temperature (eV) of L. Vriens and H. M.  Smeets. (See
DOC=VRIENS,PRA,22,1980). The coefficient data values are

1. The principal quantum number (n) of the initial state.

#JRREC1 
IAEA subroutine to evaluate the analytic function used to represent the
radiative recombination rate coefficient (cm[3]/s), as a function of
temperature (eV), into the n=1 and n=2 levels as used in reaction 2.1.8
of Janev et al. (See DOC=H-HE-PLASMA). The coefficient data values are

1. The threshold energy (eV) for the transition.
2. Parameter A(n,l) as quoted in the formulae of DOC=H-HE-PLASMA.
3. Parameter chi (n,l) as quoted in the formlulae of DOC=H-HE-PLASMA.

#JRREC3 
IAEA subroutine to evaluate the analytic function used to represent the
radiative recombination rate coefficient (cm[3]), as a function of
temperature (eV), into the levels with principal quantum number, n>=3 as
used in reaction 2.1.8 of Janev et al. (See DOC=H-HE-PLASMA). The
evaluation function is only defined for the reactions covered in
DOC=H-HE-PLAMSA. It may need to be modified to cover reactions not
contained in the above database. The coefficient data values are

1. Z, the atomic number of the target.
2. Q, the charge state of the target ion.  (for Hydrogen Q=1,
   for Helium Q= 1 or 2).
3. The principal quantum number (n) of the final state. (n>=3)

#JEEXC1 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]),  as a function of the
electron energy (eV), for dipole allowed transitions The analytic
expression is defined in reaction 2.3.7 in DOC=H-HE-PLASMA.
The coefficient data values are

1. Parameter beta(n) as quoted in the formulae of DOC=H-HE-PLASMA.
2. Parameter gamma(n) as quoted in the formlulae of DOC=H-HE-PLASMA.
3. Parameter delta as quoted in the formlulae of DOC=H-HE-PLASMA.
4. Principal quantum number (n) of the initial state.
5. Orbital angular momentum (l) of the initial state.
6. Orbital angular momentum (l') of the final state.
7. Spin multiplicity (2S+1) of the initial state.
8. The total angular momentum of the initial state, only required
   for fine structure transitions.
9. Principal quantum number (n') of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

10. The threshold energy (eV) for the transition.
11. The oscillator strength for the transition.
12. Coefficient xi of the formulae in DOC=H-HE-PLASMA.

#JEEXC2 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]) as a function of the
electron energy (eV), for dipole forbidden transitions without a spin
change.  The analytic expression is defined in reaction 2.3.2 in
DOC=H-HE-PLASMA. The coefficient data values are

1. Parameter beta(n,l) as quoted in the formulae of DOC=H-HE-PLASMA.
2. Parameter gamma(n,l) quoted in the formlulae of DOC=H-HE-PLASMA.
3. Parameter delta(n,l) quoted in the formlulae of DOC=H-HE-PLASMA.
4. Principal quantum number (n) of the final state.
5. Orbital angular momentum (l) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

6. The threshold energy (eV) for the transition.
7. The energy of ionization of the (n,l) state.
8. Coefficient t of the formulae in DOC=H-HE-PLASMA.

#JEEXC3 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]) as a function of
electron energy (eV), for dipole forbidden transitions with a spin
change.  The analytic expression is defined in reaction 2.3.3 in
DOC=H-HE-PLASMA. The coefficient data values are

1. Parameter a(n,l) as quoted in the formulae of DOC=H-HE-PLASMA.
2. Parameter b quoted in the formlulae of DOC=H-HE-PLASMA.
3. Principal quantum number (n) of the final state.
4. Orbital angular momentum (l) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

5. The threshold energy (eV) for the transition.
6. The energy of ionization of the (n,l) state.

#JEEXC4 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]) as a function of
electron energy (eV), for dipole forbidden transitions with a spin
change. The analytic expression is defined in reaction 2.3.4 in
DOC=H-HE-PLASMA. The coefficient data values are

1. Parameter Q as quoted in the formulae of DOC=H-HE-PLASMA.
2. Orbital angular momentum (l) of the initial state.
3. Spin multiplicity (2S+1) of the initial state.
4. Orbital angular momentum (l') of the final state.
5. Spin multiplicity (2S'+1) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

6. The threshold energy (eV) for the transition.

#JEEXC6 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections as a function of electron
energy (eV), for transitions between states defined in terms of the
initial and final principal quantum numbers alone. The analytic
expression is defined in reaction 2.3.8 in DOC=H-HE-PLASMA.
The coefficient data values are

1. Principal quantum number (n) of the initial state.
2. Principal quantum number (m) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

3. The threshold energy (eV) for the transition.
4. Coefficient xi of the formulae in DOC=H-HE-PLASMA.
5. The oscillator strength for the transition, averaged over the
   orbital angular momentum and spin quantum numbers.

#JIONHE 
IAEA subroutine to evaluate the analytic function used to represent the
electron impact ionization cross sections (cm[2]), as a function of
electron energy (eV), as defined in reaction 2.3.10 in DOC=H-HE-PLASMA.
The coefficient data values are

1. An integer, ITRANS,  that defines the type of initial excited
   target state. For (n,n') transitions ITRANS = 1 while for
   (nl,n'l') transitions ITRANS = 2.
2. Coefficient beta of the formulae in DOC=H-HE-PLASMA.
3. Principal quantum number (n) of the initial state.

For (nl,n'l') transitions the following coefficients are also
required,

4. Orbital angular momentum (l) of the initial state.
5. Spin multiplicity (2S+1) of the initial state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

6. The threshold energy (eV) for the transition.

#JEION5 
IAEA subroutine to evaluate the analytic function used to represent the
electron impact ionization cross sections (cm[2]), as a function of
electron energy (eV), as defined in reaction 2.3.10 in DOC=H-HE-PLASMA.
The coefficient data values are

1. Coefficient A of the formulae in DOC=H-HE-PLASMA.
2. Principal quantum number (n) of the initial state.

#JRREC2 
IAEA subroutine to evaluate the analytic function used to represent the
radiative recombination rate coefficient (cm[3]) as a function of
temperature (eV), into n=2 levels as used in reaction 2.3.13 of Janev et
al. (See DOC=H-HE-PLASMA). The coefficient data values are

1. The threshold energy (eV) for the transition.
2. Parameter A(n,l) as quoted in the formulae of DOC=H-HE-PLASMA.

#JEEXC5 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]) as a function of
electron energy (eV), for a target of singly charged Helium. The basic
form of the analytic expression is given in reactions 2.3.17 in
DOC=H-HE-PLASMA. The coefficient data values are

1. An integer, TTYPE. If TTYPE = 1 the expression contains a
   natural logarithm involving the reduced energy u, if TTYPE=2
   this term is replaced by unity.
2. Normalisation coefficient (leading number in formula of
   DOC=H-HE-PLASMA)
3. Parameter C(n) as quoted in the formulae of DOC=H-HE-PLASMA.
4. Parameter phi(n) as quoted in the formulae of DOC=H-HE-PLASMA.
5. Principal quantum number (n) of the initial state.
6. Principal quantum number (m) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

7. The threshold energy (eV) for the transition.

#JEEXC7 
IAEA subroutine to evaluate the analytic function used to represent
electron impact excitation cross sections (cm[2]), as a function of
electron energy (eV), for a target of singly charged Helium for (n,n')
transitions with n>=3.  The analytic expression is given in reactions
2.3.18 in DOC=H-HE-PLASMA. The coefficient data values are

1. Principal quantum number (n) of the initial state.
2. Principal quantum number (m) of the final state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

3. The threshold energy (eV) for the transition.
4. The oscillator strength for the transition, averaged over the
   orbital angular momentum and spin quantum numbers.

#JDREC1 
IAEA subroutine to evaluate the analytic function used to represent the
reaction rate coefficient for dielectronic recombination.  (cm[3]) as a
function of temperature (eV), The analytic expression is given in
reactions 2.3.22 in DOC=H-HE-PLASMA. The coefficient data values are

1. Principal quantum number (n) of the initial state.
2. Principal quantum number (m) of the final state.

#JAN3 
IAEA subroutine which is a general form for a cross section with
simple logarithmic and inverse energy dependence, which is of the form

cross section = A * LN ( B*E + C) / E

where A, B and C are coefficients and E is the energy. For instance,
as used in reaction 2.1.3 in DOC=H-HE-PLASMA. The coefficient data values 
are

1. Threshold energy for the transition (eV).
2. Coefficient A.
3. Coefficient B.
4. Coefficient C.

#PEXCH4 
IAEA subroutine to evaluate the analytic function used to represent the
heavy particle impact excitation cross sections (cm[2]), as a function
of impact energy (eV), from the dipole close-coupling approximation.
The analytic expression is defined in reaction 3.1.5 in
DOC=H-HE-PLASMA. The coefficient data values are

1. The atomic number of the target.
2. The atomic number of the projectile.
3. An integer, ITRANS,  that defines the type of initial excited
   target state. For (n,n') transitions ITRANS = 1 while for
   (nl,n'l') transitions ITRANS = 2.
4. Principal quantum number (n) of the initial state.
5. Principal quantum number (n') of the final state.

For (nl,n'l') transitions the following coefficients are also
required,

6. An integer, SUMEN,  that defines the choice of ionization energy
   to be used. If SUMEN = 0 the ionization energy for the
   specific (n,l,s) state is taken, while if SUMEN = 1 the ionization
   energy averaged over possible (l,s) values is taken.
7. Orbital angular momentum (l) of the initial state.
8. Orbital angular momentum (l') of the final state.
9. Spin multiplicity (2S+1) of the initial state.
10. The total angular momentum of the initial state, only required
    for fine structure transitions.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

11. The threshold energy (eV) for the transition.
12. The energy difference between the initial and final states.
13. Coefficient lambda(n,n').

#IONBEA 
IAEA subroutine to evaluate the analytic function used to represent the
heavy particle impact ionization cross sections (cm[2]), as a function
of impact energy (eV), from the binary encounter approximation.  The
analytic expression is defined in reaction 3.1.7 in DOC=H-HE-PLASMA.
The coefficient data values are

1. The atomic number of the target.
2. The charge of the core of the target.
3. The binding energy of the target atom.
4. The charge of the projectile ion.
5. An integer, ITRANS,  that defines the type of initial excited
   target state. For (n,n') transitions ITRANS = 1 while for
   (nl,n'l') transitions ITRANS = 2.
6. An integer, SUMEN,  that defines the choice of ionization energy
   to be used. If SUMEN = 0 the ionization energy for the
   specific (n,l,s) state is taken, while if SUMEN = 1 the ionization
   energy averaged over possible (l,s) values is taken.
7. Principal quantum number (n) of the initial state.
8. Orbital angular momentum (l) of the initial state.
9. Spin multiplicity (2S+1) of the initial state.

On the first call to the evaluation function the coefficient array is
updated with the following data, where the integer refers to the index
for the coefficient array,

11.   The threshold energy (eV) for the transition.
12.   The "effective" value of the principal quantum number used in
      evaluating the cross section.

#JCX1 
IAEA subroutine to evaluate the analytic function used to represent the
charge transfer cross sections (cm[2]), as a function of impact energy
(eV), from the analytic expression defined in reaction 3.1.12 in
DOC=H-HE-PLASMA. The coefficient data values are

1. Principal quantum number (n) of the initial state.

#JCX2 
IAEA subroutine to evaluate the analytic function used to represent the
charge transfer cross sections (cm[2]), as a function of impact energy
(eV), from the analytic expression defined in reaction 6.1.7 in
DOC=H-HE-PLASMA. The coefficient data values are

1. Principal quantum number (n) of the initial state.
2. Principal quantum number (n') of the final state.

#DBETA 
IAEA function to evaluate the D(beta) function defined in Janev R.K.
and Presnyakov L.P., J. Phys. B., 13, 4233, (1980) and also as defined
on page 319 in DOC=H-HE-PLASMA. The argument passed to DBETA is the
value of beta.

#OSCSTH1 
IAEA subroutine to return the oscillator stength, F(nl,ml') for
hydrogen. The transitions included are those contained in table A.1 on
page 314 in DOC=H-HE-PLASMA.  The argument call of OSCSTH1 is, (in
order):

1. Principal quantum number (n) of the initial state. (Integer input)
2. Orbital angular momentum (l) of the initial state. (Integer input)
3. Principal quantum number (m) of the initial state. (Integer input)
4. Orbital angular momentum (l') of the final state. (Integer input)
5. Oscillator strength, F(nl,ml'), (Real output)

#OSCSTH2 
IAEA subroutine to evaluate the oscillator stength, F(n,m) for hydrogen
from the formula of Johnson L. C., Astrophys. J. 174, 227, (1972).
The argument call of OSCSTH2 is:

1. Principal quantum number (n) of the initial state. (Integer input)
2. Principal quantum number (m) of the final state. (Integer input)
3. Oscillator strength, F(n,m). (Real output)

#OSCSTHE 
IAEA subroutine to return the oscillator stength, F(nl,ml') for
helium. The transitions included are those contained in table A.6 on
page 317 in DOC=H-HE-PLASMA.  The argument call of OSCSTHE is, (in
order):

1. Principal quantum number (n) of the initial state. (Integer input)
2. Orbital angular momentum (l) of the initial state. (Integer input)
3. Principal quantum number (m) of the initial state. (Integer input)
4. Orbital angular momentum (l') of the final state. (Integer input)
5. Spin multiplicity (2S+1) of the state. (Integer input)
6. The total angular momentum of the initial state, only required
   for fine structure transitions (Integer  input).
7. Oscillator strength, F(nl,ml'), (Real output)
8. Error Indicator. (Character output field)

#EXINT 
IAEA function to evaluate the exponential integral function defined in
"Handbook of Mathematical Functions", Abramowitz M. and Stegun I. A,
National Bureau of standard, Washington DC, page 231, 1964.

#HEIONEN
IAEA subroutine to return the ionization energy for excited states of
helium taken from tables A.4 and A.5  on page 316 in DOC=H-HE-PLASMA.
The argument call of HEIONEN is

1. Principal quantum number (n) of the excited electron. (Integerinput)
2. Orbital angular momentum (l) of the excited electron. (Integer input)
3. Spin multiplicity (2S+1) of the state. (Integer input)
4. An input integer, SUMEN,  that defines the choice of ionization
   energy to be returned. If SUMEN = 0 the ionization energy for the
   specific (n,l,s) state is taken, while if SUMEN = 1 the ionization
   energy averaged over possible (l,s) values is returned.
5. The ionization energy (eV). (Real output field).
6. Error Indicator. (Character output field)

#HEEXCEN 
IAEA subroutine to return the excitation energy for excited states of
helium taken from tables A.4 and A.5  on page 316 in DOC=H-HE-PLASMA.
The argument call of HEEXCEN is

1. Principal quantum number (n) of the excited electron. (Integer input)
2. Orbital angular momentum (l) of the excited electron. (Integer input)
3. Spin multiplicity (2S+1) of the state. (Integer input)
4. An input integer, SUMEN,  that defines the choice of excitation
   energy to be returned. If SUMEN = 0 the excitation energy for the
   specific (n,l,s) state is taken, while if SUMEN = 1 the excitation
   energy averaged over possible (l,s) values is returned.
5. The ionization energy (eV). (Real output field).
6. Error Indicator. (Character output field)

#FNAGEX 
IAEA subroutine to return the electron impact excitation reaction rate
coefficients for a specified electron temperature in eV, using the
analytic expressions of Y. Itikawa et al, See DOC=ADNDT,33,149,(1985).
In this reference two type of analytic fits are employed to represent
the reaction rate coefficients. The basic forms are a power-log type
expression and an exponential type expression. See
DOC=ADNDT,33,149,(1985) for details. The argument call of FNAGEX is 

1. Index for the type of fit, either 1 or 2, Variable ITYPE
2. The excitation energy (eV)
3. The lower limit of the fit to the reduced energy, X
4. The upper limit of the fit to the reduced energy, X
5. statistical weight of the initial state (2s+1)*(2l+1)
6. Parameter A, 7.  Parameter B,  8. Parameter C, 9 Parameter D
10. Parameter E

If ITYPE = 1 this can be followed by three more parameters for
the region which contains resonaces represented by a linear term

11.  Parameter P, 12. Parameter Q, 13. Parameter X1

If ITYPE = 2 this can be followed by

11. Parameter F and possibly three more parameters for the region 
    which contains resonances represented by a linear term.

12. Parameter P, 13. Parameter Q, 14. Parameter X1

#REFL1 
IAEA subroutine to return the general form for the particle or energy
reflection coefficient as a function of projectile energy in keV.
DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL1 is the
six fitting coefficients needed to determine the reflection
coefficient.

#REFL2
IAEA subroutine to return the particle or energy reflection coefficient
as a function of projectile energy in keV for the extended energy fit.
DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL2 is the
eight fitting coefficients needed to determine the reflection
coefficient.

#REFL3 
IAEA subroutine to return the self particle or energy reflection
coefficient as a function of projectile energy in keV.
DOC=IAEA-AM-REP-249,(1991) The coefficient array passed to REFL3 is the
nine fitting coefficients needed to determine the reflection
coefficient.