Miroslav Hallo  Software
Bayesian earthquake hypocenter location from direct P and Swaves arrival times
 Open source Matlab functions for assessment of the earthquake hypocenter location by grid search method in the homogeneous or 1D layered medium.
The posterior probability density function (PDF) of the earthquake hypocenter location is following theory by Tarantola (2005, Chapter 7.1).
Hence, the solution contains also information about hypocenter location uncertainty.
Codes are published under the GNU General Public License (GNU GPL) for noncommercial use.
Methodology described in:
 Hallo, M., Opršal, I., Asano, K., Gallovič, F. (2019): Seismotectonics of the 2018 Northern Osaka M6.1 earthquake and its aftershocks: joint movements on strikeslip and reverse faults in inland Japan, Earth, Planets and Space, 71:34.

List of the programs:
 File crustal.dat contains 1D velocity model of Earth's crust.
 The code a1_Pert_Model.m perturb velocity model and perform MonteCarlo simulations to get uncertainties of P and S wave arrival times.
 The code a2_LatLon2xyz.m converts coordinates of stations from WGS84 Lat/Lon into local coordinates Easting/Northing
 The code a3_hypo_homo.m compute earthquake hypocenter location by grid search method in the homogeneous medium
 The code a4_hypo_layer.m compute earthquake hypocenter location by grid search method in the 1D layered medium
 The a3_hypo_homo.m is much faster than a4_hypo_layer.m, but less precise as the homogeneous medium is used. We recommend to use a3_hypo_homo.m to rough search for hypocentre and the a4_hypo_layer.m for precise adjustments
Download
MATLAB (revision 12/2018):

Requirements
MATLAB:

The codes do not require any additional MatLab Toolboxes
The codes were tested on various versions of MatLab software, and it should work without any problems even on older versions. Tested versions: Matlab 2012b 32bit (Windows 7), Matlab 2018b 64bit (Windows 7).