`RichardsJacobianEval` and `NlFunctionEval` both add in Dirichlet BC pressure values in the same manner. Should we have this behind a macro?
RichardsJacobianEvalandNlFunctionEvalboth add in Dirichlet BC pressure values in the same manner. Should we have this behind a macro?