normalize now uses "RLOOP", to use openmp

I think.

"impose model constraints" should only be called *after* all Runge Kutta
substeps have been computed.
Otherwise it messes up the numerical method.