Perform higher-order moment tensor inversions

The Neighbourhood Algorithm must be obtained from the original authors.

Python 3 is used for preprocessing and plotting.

This directory contains the problem-specific source code for higher-order moment tensor inversions.

In Inversion/src there is the main routine homti_na.f, some subroutines and a Makefile. To build the code, run

cd Inversion/src
make NA_DIR=</path/to/NA/package> all

The all target will build the necessary parts of the NA and link with the routines here. The default value of NA_DIR is ../NA.

This python package contains helpful scripts for data processing and plotting. Install in a python3 environment using pip

cd tempseis
python3 -m venv .venv
source .venv/bin/activate
pip install .

This will make sure all the correct packages and functions are available to the various scripts. Feel free to use a different environment name.

In pytempseis/dataprocessing, edit the parameters.yml as required then run

python process_data.py

to run through all the processing steps. Alternatively you can run the individual .py scripts if you wish to process step-by-step.

In pytempseis/plot you'll find helpful plotting scripts, most of which will take a result direction as input. TODO: use argparse

Once the data has been preprocessed, the easiest thing to do is

./total_proc_script

from the top level directory. This will copy data to the right place so that it can be found by the NA code, run the inversion, and plot the results. Everything should then be in the Results directory.