SPMST: Surface Wave Tomography with Topography

I'm still organizing manual and examples ...

Introduction

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.

USAGE:

1. Set your Topography File

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])

2. Set Study Area and Inverse Problem Parameters

Set all parameters in spmst.in. It is a self-explained file.

3. Set dispersion data

In surfdata.txt, repeat the snippet:

# evlo evla nsta
stlo stla v0 (nsta lines )

4. Set Your Initial and True Model

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.

5. Compile and Run the Code

Download Eigen, and compile it by using:

mkdir build; cd build;
cmake .. -DCXX=g++ -DFC=gfortran -DEIGEN_INC=/path/to/your/eigen
make -j4

to 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.

Gallery

Topography Adaptive Mesh

Travel time Field with Topography

Checkerboard Recovered

An Image in Sichuan Basin (T=6s)