Added support for calculating fluxes at the boundaries
Added units tests for Neumann boundary conditions