this function cause run-time error while simulating monod2d problem. disbleing for now.
fall back to the backward Euler scheme for ODE.
a bit lose constrain on local problem.