(1) Use obspy and download data for a shallow & strong earthquake
which generated a lot of surface waves.

(2) Run polfre_s1.66el on this data and plot the BAZ as function of
time and frequency. Verify the BAZ using station and event position. 
 
(3) Try to detect also the R2,R3,.. wave trains.