Converting Complex ODE Model from as_op to Native PyMC ODE Solvers

[aseyboldt]

I fixed the link in the above post.
Looks like a fun problem by the way! Is this some geochemical model? Reminds me a bit of reactive transport problems, but not quite