you are right, we don't have (yet?) a tracer solver. We thought about it, but it can always be done as a post-processing step (using e.g. Forward Euler based on the computed velocity time series), in MATLAB/python. If you would like to implement one in C++ in ISSM, let us know and we can help you.
All the best