Use new interpolators in some output functions
Add get dim and ts_pt functions to caching interpolator
Use new interpolators in hms calculations
Fix stupid bug in magnetostatic field calculation

* I was copying more values than were initialised (h_ms is of length 3
  but dphi/dx is only of length dim!), again...