RATES
=====

Rates is a module for integrating a cross section with an energy distribution,
to obtain rate coefficients between 100 and 10,000K. A Maxwellian velocity
distribution is assumed. The 'alternative extended Simpson's rule' is used
for the integration (from Numerical recipes).  
          
This program assumes that energies and cross sections are 
tabulated respectively EITHER in Rydbergs and Bohrs**2 (Default) 
OR in eV and Angstroms^2. Rates are converted from atomic units
to cm**3 sec**-1 by multiplying by 6.126163*10**(-9).

A sample input file for rates looks like the following::

    &RATESIN
     name="Rate Coefficient for xsec_1_doublet ",
     ieunit=2,
     infile = "xsec_1_doublet",
     outfile =" xsec_1_doublet_rates.dat",
     deffile = "xsec_1_doublet_rates.def",
    /

where ``xsec_1_doublet`` is the input file and the filename with the cross
section data. This file is created in each channel run from automatically
created **fort.402** during ``XSEC`` module run for example. This name
pattern is then used to create the output and the definition files for
the ``rates`` binary. The ``ieunit=2`` has the same meaning as in
module ``RSOLV`` and implies that the cross section is in eV units.

This module was updated by Keir Little in 2016 from the original rates.f
