Major changes
- fixed NaN problem when source and station are on the same longitude
- improved precision and efficiency of computation for shallow events
- supported conditions with fluid or vacuum at center of planet
- converted from .f to .f90 format and organized code, leading to speed-up of run time