Bug in time indices in run.py

Found a bug in run.py: tind+1 was being used in the loop to run tracmass instead of tinds[j+1]. Have updated so that tout contains one extra time index and the run loop goes uses only the indices stored in tinds.
Also small tweak to the interpolation logic.