Lagrangian simulation of drift trajectories
===========================================