I\'m using FiPy to simulate a moderately large system of coupled PDEs (11 equations) that calculate the time-dependent evolution of a biochemical advection-diffusion-reactio