Added interpolation function Flux_get_dF. Added checks for out of boundary values in it. Added test for reading (to do: make it a real test with boolean result)