with entry tally.
collects and writes statistics on parameter adjustments
three levels of statistics can be compiled.  level 3 stats
will include several groups of level 2 stats which, in turn,
will include several groups of level 1 stats.  an example of
level 1 stats would be the s series of harmonic coefficients
which could be subsumed at level 2 with other harmonic coeff.
and at level 3 with initial conditions and other parameters
for a planet.

at the main entry to the routine there are three steps in
processing.  first, a summary is written for the level indi-
cated by plevel.  then the level indicated by ulevel is
updated by adding the totals of level plevel to those of
ulevel, where ulevel.gt.plevel.  finally, level plevel is
reinitialized (zeroed).  if ulevel.le.plevel then plevel
is summarized and reinitialized but no higher level is
updated.  the special case where plevel = 0 results in
(re)initialization of levels 1 through ulevel.  title is an
8 character string which can serve as a label on the written
summary.
to avoid redundant summarys, nu(ulevel) is incremented upon
updating and, therefore, nu(plevel) must be.gt.1 for
summarizing.

local variables
nu(level) = number of updates at level
m(level) = number of adjustments in group for which stats
included (sig.ne.0)

t(1) = average of fracts
t(2) = average of abs of fracts
t(3) = average of squares of fracts
t(4) = root mean square of fracts