The following lines show how to create a simple travel time plot for a given distance range, selected phases and the iasp91 velocity model using the travelTimePlot() function of the module obspy.taup.
from obspy.taup.taup import travelTimePlot
travelTimePlot(min_degree=0, max_degree=50, phases=['P', 'S', 'PP'],
depth=120, model='iasp91')
[source code, hires.png, pdf]