I am trying to understand the behavior of PyLith for a simple 2.5D model with gravity. In this set up, the model is 10 km x 2 km x 100m with an inclusion of an horizontal layer that is 50m thick. See the sketch below.
All sides, except the top boundary, have roller boundary conditions and I am using gravitational body forces. All materials are elastic.
I was expecting the model to show an uniform displacement in the Z direction, but instead what I find is that the thin layer shows a displacement that is much larger than the surrounding. See figure below with the results.
Could you comment on what could be causing this issue ? Do you think that this is correct ?
Below is a zip folder with the set up (.cfg file, .jason file and mesh journal file ) that can be used to reproduce the behavior.
PyLith_ProblemSetup.zip (347.3 KB)
PyLith version details:
Configured on 2017-08-03 14:01:44 -0500, GIT branch: branch-not-available, revision: v3.7.6-4826-gd686aaf
MPI standard: 3.0, implementation: MPICH, version: 3.1.3
HDF5 version: 1.8.11
NetCDF version: 184.108.40.206
Proj version: 480