Direct surface wave radial anisotropy tomography package (DSurfRTomo)

please refer to:

Hu, S., H. Yao, and H. Huang (2020), Direct surface wave radial 
anisotropy tomography in the crust of the eastern Himalayan 
syntaxis. Journal of Geophysical Research: Solid Earth, 125, 
e2019JB018257. https://doi.org/10.1029/2019JB018257

for details of the algorithm

The code is based on previous studies, especially on the implementation of DSurfTomo:

Fang, H., H. Yao, H. Zhang, Y-C Huang, and R. D. van der Hilst (2015)
Direct inversion of surface wave dispersion for three-dimensional 
shallow crustal structure based on ray tracing methodology and 
application. Geophysical Journal International, 201, 1251-1263.

Please also refer to:

Rawlinson, N. and M. Sambridge (2004) Wave front evolution in 
strongly heterogeneous layered media using the fast marching method,
Geophysical Journal International, 156(3), 631-647

for implementation of the fast marching method, and

Herrmann, R. B. (2013) Computer programs in seismology: An evolving
tool for instruction and research. Seismological Research Letter,
84(6), 1081-1088

for implementation of the 1-D surface wave dispersion kernel.

The dispersion data (ALLR.dat for Rayleigh wave and ALLT.dat for Love wave), resulting model (DSurfRTomo.inMeasurement.dat) in the crust of the eastern Himalayan syntaxis is provided in example/

#############

2019/05/18 The code may still need some modification

##############

output (default DRadiSurfTomo.inMeasure.dat) is in the format : lon lat dep vsv gamma

To compute average shear wave velocity (vs) and radial anisotropy (xi), use the following equaitons:

1. vs=vsv*(1+gamma)/2.0

2. xi=2*(gamma-1)/(gamma+1)*100%

##############

2021/08/08

1. check parallel computation can be used
2. add roughness computation
3. output Rayleigh/Love raypath at the final iteration
4. add some useful scripts in utils/
5. check noiselevel can be used in the synthetic tests

##############