Ok. So the error message was helpful in finding the error in the pde's. Concerning the numerical problems I learned in former times, that sometimes a rescaling of the variables (concentration, time and lengths), so that constants have about the same order of magnitude, turns out to make things running smoothly. But I can't make a proposal having not looked intensely at your system.
Regards Hans
