Fixed wrong coefficient in flow conservation.
Added new function to read solution.