Added finite differencing step to get cell fluxes