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).