Included effects of arrival time delay in the formulae for the waveforms (assuming a round earth)

Expressed output X sin(phi) + Y cos(phi) in the alternative form r sin(phi+offset)

Currently offset is printed in degrees.  Printing code assumes UTF-8 on the terminal. This may be fragile, I'm not sure.