Your code is not running with the Monod2d benchmark

After setting a bit more strict tolerance criteria, the code fails to solve the Newton iteration even in the first time step. Therefore, it is not appropriate to say that the implementation of kinetic reactions are fully verified. Still need to work on it.