lileipku00 / axisem3d

AxiSEM3D - A Quickstart Guide

I look ugly to you? I am a Markdown file. If you do not have a Markdown reader, read me on github, or paste me here.

Report issues to [email protected].

1 Get AxiSEM3D

git clone https://github.com/kuangdai/AxiSEM3D
export CURRENT_WORK_DIR=$PWD

2 Installing dependencies

AxiSEM3D has been built upon a few modern numerical packages for its performance and sustainability, as listed below (newer versions are acceptable):

name	version	build-from-source instructions
mpi	---	Common implementations such as open-mpi or mpich2 are recommended.
boost	1.60	Install Boost C++ Libraries with boost::mpi, following here.
eigen3	3.3-rc1	Simply download and unzip to your path of treasures.
fftw	3.3	Install both double- and single-precision (with --enable-float) versions to the same path. See here.
metis	5.1	Download and install following Install.txt. Skip step 3 in Install.txt to use the 32-bit build.
Exodus	---	Must be built from source, upon both HDF5 and NetCDF4. Refer to axisem3d_depends.sh.

Don't panic! Most of these popular packages (except Exodus) may be handily installed with free package management software, such as Homebrew for Mac OS X, Linuxbrew for common Linux distributions, and Conda for both. Here we introduce the wizard axisem3d_depends.sh.

• Get Brew (required)
  • Homebrew for Mac OS X
  • Linuxbrew for Linux (Finish all the three steps!)
• Get Conda (optional)
• Edit the first few lines in axisem3d_depends.sh and run it.
• Check your ~/.bash_profile (or ~/.bashrc) and ~/.axisem3d_roots. Grats, you are done!

Some notes:

• Depending on what you need and your network speed, axisem3d_depends.sh may take long time (especially for open-mpi and HDF5) to finish.
• Installing HDF5 and NetCDF4 might be problematic on some systems. If Brew fails or AxiSEM3D shows I/O runtime error later on, reinstall them with conda (see options in axisem3d_depends.sh), and re-run axisem3d_depends.sh. If the problem persists, you have to build them from source.

3 Building AxiSEM3D

• Edit SOLVER/CMakeLists.txt if needed (normally not), including

  to-be-edited	notes
  compiler suit	Changes are usually required only for HPC clusters, unless you have a bizarre MPI.
  dependency roots	No need if you have used axisem3d_depends.sh to install the dependencies.
  FFTW_WISDOM_DIR	Just specify any directory you like, or leave it as it.

• Build AxiSEM3D (go step by step to see what's happening)

  • Style 1: work under source

    cd $CURRENT_WORK_DIR/AxiSEM3D
    # make a simulation directory
    mkdir my_first_run
    cd my_first_run
    # cmake (*** see the notes below ***)
    cmake -DCMAKE_BUILD_TYPE=release ../SOLVER
    # compile and link
    make -j4
    # copy the input folder
    cp -R ../template/input ./
    # run it with any number of processors
    mpirun -np 4 ./axisem3d
    # check the outputs
    ls output

    Read the output of cmake carefully, which should end up with (If not, verify dependency installations as prompted)

    -- Configuring done
    -- Generating done
    -- Build files have been written to: ...

  • Style 2: keep source clean (suggested)

    To offer maximal flexibility for various infrastructures, AxiSEM3D is organized such that the directories of source, build and simulations are fully independent of one another. For exmaple:

    ########## build ##########
    cd $CURRENT_WORK_DIR
    # make a build directory
    mkdir my_axisem3d_build
    cd my_axisem3d_build
    # cmake
    cmake -DCMAKE_BUILD_TYPE=release ../AxiSEM3D/SOLVER
    # compile and link
    make -j4
    
    ########## run ##########
    # make a simulation directory
    cd $CURRENT_WORK_DIR
    mkdir my_second_run
    cd my_second_run
    # copy the executable 
    cp ../my_axisem3d_build/axisem3d ./
    # copy the input folder
    cp -R ../AxiSEM3D/template/input ./
    # Though optional, it is always a good practice 
    # to backup the source code for reproducibility:
    cp -r ../AxiSEM3D/SOLVER ./
    # run it with any number of processors
    mpirun -np 4 ./axisem3d
    # check the outputs
    ls output

4 The MESHER

In the above examples, we use the mesh file prem_ani_one_crust_spherical_2D_50s.e (anisotropic PREM model with one crustal layer and a 50 s period), located at SOLVER/template/input. To generate an AxiSEM3D mesh like this, you will need the salvus_mesher, a python-based command-line tool to generate several types of 2-D and 3-D finite element meshes.

Currently, the salvus_mesher is not published. Here we provide a stable version for AxiSEM3D:

• Install pyexodus (mainly by Lion Krischer at ETH Zurich)
• Install salvus_mesher (mainly by van Driel Martin at ETH Zurich)

These two packages have their own requirements and manual.