Improve discrete directions script, use polar plot

Use polar projection.
phi is now plotted in radians, the ticks still show degrees.
Correct theta axis limits.
Add isotime lines, where the time is constant for 2 detectors.