Changes to maintain the existing functionality (and input arguments) whilst incorporating the increased flexibility in the times and input data you added.
Add support for writing out the data properly when we're bracketing a value. This now works for single scalars across non-scalar outputs as well as for non-scalar inputs on non-scalar outputs.
Remove the additional loading of the data when it is not strictly necessary. This follows the improvment Karen made in speed.
Merge the updates I did whilst on the temporary branch.