optimised t_ar construction
for GeoChem
Currently working when velocity is
ignored, posibility time scale is
out