TANDEM

Tsunami Adjoint simulator for Nonpoint sources using Dispersive wave equations on a self-gravitating Earth Model

Fig01

model_difference

Technical References

Takagawa, T., Allgeyer, S., & Cummins, P. (2024). Adjoint synthesis for trans‐oceanic tsunami waveforms and simultaneous inversion of fault geometry and slip distribution. Journal of Geophysical Research, [Solid Earth], 129(6).

Installation

TANDEM can be installed in the following 7 steps. This is an example for Ubuntu. Windows users are recommended to use WSL (Windows Subsystem for Linux).

Step 1. gfortran & pip

sudo apt update
sudo apt-get update
sudo apt install python3-pip
sudo apt install build-essential
sudo apt install gfortran

Step 2. Numpy

python -m pip install --user numpy

Step 3. FFTW

cd /tmp
wget https://www.fftw.org/fftw-3.3.10.tar.gz
tar xf fftw-3.3.10.tar.gz
cd fftw-3.3.10/
./configure --help # check options
./configure --enable-openmp --enable-shared --enable-avx
make -j8
sudo make install
# fftw3 is installed at /usr/local
sudo nano /etc/ld.so.conf
-----
include /etc/ld.so.conf.d/*.conf # default
/usr/local/lib                   # Add this line
-----
sudo ldconfig # update the library link
ldconfig -p | grep libfftw # check the library link

Step 4. Shtns

cd /tmp
wget https://gricad-gitlab.univ-grenoble-alpes.fr/schaeffn/shtns/-/archive/master/shtns-master.tar.gz
tar xf shtns-master.tar.gz
cd shtns-master/
./configure --help # check options
./configure --enable-openmp --enable-python
make -j8
sudo -E python setup.py install # -E option is needed for super user to find numpy module 

Step 5. OpenMPI & mpi4py

cd /tmp
wget "https://download.open-mpi.org/release/open-mpi/v5.0/openmpi-5.0.3.tar.gz"
tar xf openmpi-5.0.3.tar.gz
cd openmpi-5.0.3/
./configure --help # check options
./configure CC=gcc CXX=g++ F77=gfortran FC=gfortran # This may take a few minutes.
make -j8 # This may take a few minutes.
sudo make install # lib
sudo ldconfig # update the library link

# check installation
mpicc -v
mpiexec --version

# mpi4py
python -m pip install mpi4py

Step 6. petsc & petsc4py

# install requirements
sudo apt-get install libblas-dev liblapack-dev
python -m pip install --user cython --proxy=10.72.10.106:8080

cd /tmp
wget https://web.cels.anl.gov/projects/petsc/download/release-snapshots/petsc-3.21.0.tar.gz
tar xf petsc-3.21.0.tar.gz
cd petsc-3.21.0/
./configure --help
sudo mkdir /opt/petsc
sudo chown $USER /opt/petsc
./configure --prefix=/opt/petsc --download-hypre --download-petsc4py=1

# This may take several minutes.
make -j8 PETSC_DIR=/tmp/petsc-3.21.0 PETSC_ARCH=arch-linux-c-debug all
make PETSC_DIR=/tmp/petsc-3.21.0 PETSC_ARCH=arch-linux-c-debug install
# output
=====================================
To use petsc4py, add /opt/petsc/lib to PYTHONPATH
=====================================

nano ~/.bashrc
-----
export PYTHONPATH=/opt/petsc/lib
-----
source ~/.bashrc

Step 7. TANDEM

python -m pip install git+https://github.com/tomographyyy/tandem.git

How To Run

mpiexec -n 8 python tandem/run.py

How To View Results

view_results.ipynb

Acknowledgments

The sample topographic data was generated from GEBCO_2023 Grid.
- GEBCO Bathymetric Compilation Group 2023 (2023). The GEBCO_2023 Grid - a continuous terrain model of the global oceans and land. NERC EDS British Oceanographic Data Centre NOC. doi:10.5285/f98b053b-0cbc-6c23-e053-6c86abc0af7b

This site is open source. Improve this page.