Fixed interface for flux calculation across interfaces.

The flux accessor now takes the (possibly user-defined) flux function (vector-valued) as an input,
which is more flexible than taking any of the equations out of the PDE system.