Off fault receivers


Ascii receivers are enabled using the namelist output. Here is a commented example:

pickdt = 0.01 ! Pickpoint Sampling
pickDtType = 1 ! Pickpoint Type
RFileName = 'receivers.dat' ! Record Points in extra file

If pickDtType = 2, the output is generated every N time steps, where N is set by pickdt. If pickDtType = 1, output is generated every pickdt second.

receivers.dat is an ASCII file describing the coordinates of the receivers in the form:

x1 y1 z1
x2 y2 z2
xn yn zn

The receivers files contain the time-histories of the stress tensor (6 variables) and the particle velocities (3). Currently, there is no way to write only a subset of these variables.

The variable ReceiverOutputInterval (in the section Output of the Parameter File) controls the frequency of flushing receiver time-histories. If not specified, they are written at the end of the simulation.

Rotational Output

You can additionally choose to write the rotation of the velocity field by setting ReceiverComputeRotation=1 in the parameter file. The rotation of the vector field is defined as \(\text{rot} v = \begin{pmatrix} \partial_2 v_3 - \partial_3 v_2 \\ \partial_3 v_1 - \partial_1 v_3 \\ \partial_1 v_2 - \partial_2 v_1 \\ \end{pmatrix}\).

Placing free-surface receivers

Placing receivers on the free-surface requires special care when a realistic topography is used. The procedure to move receivers exactly to the surface is described here.

Compiling place_receivers on SuperMUC

Load the relevant modules.

git submodule update --init
mkdir build && cd build
cmake ..
make -j