LLNL_ToFi is a small python program for tomographic filtering of hypothetical seismic mantle structure (vs or vp) using the resolution matrix of the LLNL-G3D-JPS model by Simmons et al. (2015). The routine LLNL_ToFi.py performs the matrix-vector multiplication R*m=m' to obtain the filtered version m' of the given seismic model m. To be able to perform this operation, m needs to be given in the parametrization of the LLNL-G3D-JPS model.
Main author: Bernhard Schuberth (Geophysics, LMU Munich, Germany, bernhard.schuberth@lmu.de)
Development is hosted on the GitLab server of the Leibniz Supercomputing Centre (LRZ) in the bschuberth/LLNL_ToFi repository.
Input data required by LLNL_ToFi.py are located on the LLNL server:
or by email request to simmons27@llnl.gov
-
Information on the resolution matrix and the parametrization of the LLNL-G3D-JPS tomographic model can be found in LLNL-G3D-JPS_R_Matrix_TomoFilter_README.pdf
-
The LLNL-G3D-JPS model is described in
Simmons, N. A., S. C. Myers, G. Johannesson, E. Matzel, and S. P. Grand (2015), Evidence for long-lived subduction of an ancient tectonic plate beneath the southern Indian Ocean, Geophys. Res. Lett., 42, doi:10.1002/2015GL066237. -
An example of applying the resolution matrix R to a geodynamic model is described in
N. A. Simmons, B. S. A. Schuberth, S. C. Myers, D. R. Knapp (2019), Resolution and Covariance of the LLNL-G3D-JPS Global Seismic Tomography Model: Applications to Travel Time Uncertainty and Tomographic Filtering of Geodynamic Models, Geophys. J. Int., in press -
Running the code:
-
To run the code, please first get the necessary input data (i.e., the resolution matrix files) from the source given above and put them into the directory ./DATA.
-
Edit the file model.py: Modify the code such that the dummy routine project_model_3D returns the values of your specific seismic velocity model at the given coordinates (radius, lat, lon).
-
Run the code in a python 2.7 environment using one of the following commands:
- python LLNL_ToFi.py
- ./LLNL_ToFi.py
for serial processing, or:
- mpirun -np {number of processes} LLNL_Tofi.py
for parallel processing.
Specify the option -n|--no-reparam if you run the code the several times and you do not want to perform the reparametrization again. This assumes that the reparametrized version of your seismic model (i.e., on the parametrization of the LLNL-G3D-JPS tomographic model) is already stored in the directory OUTPUT_FILES.
Output files will be stored in the directory ./OUTPUT_FiLES. Files containing the reparametrized model will be named according to the variable OUTFILE_PARM_PREFIX [default: 'LLNL_G3D_JPS_Parm_layer'] and the tomographically filtered model will be stored according to the variable OUTFILE_FILT_PREFIX [default: 'LLNL_G3D_JPS_ToFi_layer'].
-