statphase2slowvec.py calculates for a given location on Earth, the horizontal slowness vectors associated with the outgoing seismic wavefield from an event at this particular location to a specified station and phase. It uses the tau-p routines in ObsPy.