Feature/pass fluxcache to neumann

This is still WIP because we have to check if that is all consistent with staggered or rather if staggered could do something similar. Also, we still need to build and run all tests!

Fixes #660 (closed)