add description about time grid shape mismatch
update infer_potential to fix grid shape mismatch
switch back to using arrays, but store units globally
in integrator, integrate up to and including t2