Made a change to Zhao equations, replacing 1.58e-4 with 1/r. This dramatically improved the accuracy of the code for mantle events.
Also fixed a bug that generated segmentation faults when source depth was exactly equal to Moho depth (to double precision).
Added a software disclaimer at the top of every source code file