Species names included automatically, kinetic reactants not working

The loop inside change_components now scan for the name of the
species, removes the trailing concentration names, and use it
as the input for the accumulatelines. total_O, total_H and cb
are done outside the loop.

There must be a match on the name of the species in the
SELECTED_OUTPUT and inside biogeochem code. Also the order is
important: the conc_array must match the order in the output
from Phreeqc.