Allow BodyForce to use variables as arguments
This MR adds the ability to use variables in
BodyForce function definitions, so that certain simulations are now enabled with the use of IMEX time stepping -- in particular, buoyancy terms in
IncNavierStokesSolver for convective flows and simple reaction-diffusion equations in
ADRSolver (through the use of a new equation system that this MR introduces).