Output format for IRC and DRC

The IRC and DRC can produce several different forms of output. Because of the large size of these outputs, users are recommended to use search functions to extract information. To facilitate this, specific lines have specific characters. Thus, a search for the "%" symbol will summarize the energy profile while a search for "AA" will yield the coordinates of atom 1, whenever it is printed. The main flags to use in searches are:

`%`	Energies for all points calculated, excluding extrema
`%M`	Energies for all turning points
`%MAX`	Energies for all maxima
`%MIN`	Energies for all minima
`%`	Energies for all points calculated
*`AA`**	Internal coordinates for atom 1 for every point
*`AE`**	Internal coordinates for atom 5 for every point
*`123AB`**	Internal coordinates for atom 2 for point 123

As the keywords for the IRC/DRC are interdependent, the following list of keywords illustrates various options.

`DRC`	The Dynamic Reaction Coordinate is calculated. Energy is conserved, and no initial impetus.
`DRC=0.5`	In the DRC kinetic energy is lost with a half-life of femtoseconds.
`DRC=1.0`	Energy is put into a DRC with an half-life of -1.0 femtoseconds, i.e., the system gains energy.
`IRC`	The Intrinsic Reaction Coordinate is calculated. No initial impetus is given. Energy not conserved.
`IRC=4`	The IRC is run starting with an impetus in the negative of the 4th normal mode direction. The impetus is one quantum of vibrational energy.
`IRC=1 KINETIC=1`	The first normal mode is used in an IRC, with the initial impetus being 1.0 kcal/mol.
`DRC KINETIC=5`	In a DRC, after the velocity is defined, 5 kcal of kinetic energy is added in the direction of the initial velocity.
`IRC=1`	`DRC KINETIC=4` After starting with a 4 kcal impetus in the direction of the first normal mode, energy is conserved.
`DRC VELOCITY KINETIC=10`	Follow a DRC trajectory which starts with an initial velocity read in, normalized to a kinetic energy of 10 kcal/mol.

Instead of every point being printed, the option exists to print specific points determined by the keywords T-PRIORITY, X-PRIORITY and H-PRIORITY. If any one of these words is specified, then the calculated points are used to define quadratics in time for all variables normally printed. In addition, if the flag for the first atom is set to "T" then all kinetic energy turning points are printed. If the flag for any other internal coordinate is set to "T" then, when that coordinate passes through an extremum, that point will be printed. As with the PRIORITY's, the point will be calculated via a quadratic to minimize non-linear errors.

N.B.: Quadratics are unstable in the regions of inflection points; in these circumstances linear interpolation will be used. A result of this is that points printed in the region of an inflection may not correspond exactly to those requested. This is not an error and should not affect the quality of the results.