Adding a new utility (range.hh) which provides adapters for pairs of

iterators that let us use range based loops more easily in Spheral.
Also adds a C++ range based version of Python's "enumerate" method to
iterate over iterators and provide an size_t index simultaneously,
which is a pattern we have a lot in Spheral.  Converted the loops in
Integrator to this new pattern as a test.