SPMST is a package to conduct 2-D surface wave tomography with topography by using shortest path method. In this method, surface wave is propagating along the (curved) surface.
The package currently support both spherical and cartesian coordinates to adapt surface wave tomography at different scales. The option could be set in the input file spmst.in. Some examples are included in the directory EXAMPLES.
set your topography in topo.dat
nlon nlat
lonmin lonmax latmin latmax
... (your topography here, in m)
...Where your topography will be started from (lonmin,latmin) to (lonmax,latmax) and is sorted like:
for i in range(nlat):
for j in range(nlon):
f.write("%f\n"%z[i,j])Set all parameters in spmst.in. It is a self-explained file.
In surfdata.txt, repeat the snippet:
# evlo evla nsta
stlo stla v0 (nsta lines )
The initial model and the True model are in the same format. It's just a single column (please see velocinit.in). The file will be started from (lonmin,latmin) to (lonmax,latmax) and is sorted like
for i in range(ny):
for j in range(nx):
f.write("%f\n"%veloc[i,j])Where nlat (ny) and nlon (nx) should be the same with that in spmst.in.
Download Eigen, and compile it by using:
mkdir build; cd build;
cmake .. -DCXX=g++ -DFC=gfortran -DEIGEN_INC=/path/to/your/eigen
make -j4to compile the package. After that three binaries will be in the bin directory: tomo2d, tomo3d and travel. The tomo2d/3d is to conduct 2-D/3-D tomography and the travel is to synthetic the whole time field for a specific event. Just feel free to go to EXAMPLES to run all the examples here.
