Pylith noisy displacement field

PyLith
1

Hello,

I am running the demo of pylith-4.1.1 crustal-strikeslip-3d step02_varslip.cfg. When checking the displacement by reading step02_varslip-domain.h5 using MATLAB (left figure), the result seems to be contaminated by points with small values(green ones), especially in the near-fault field. The plotting code is shown below.

cor = h5read(‘step02_varslip-domain.h5’,‘/geometry/vertices’);
data=h5read(‘step02_varslip-domain.h5’,‘/vertex_fields/displacement’);
figure;fastscatter(cor(1,:),cor(2,:); ,data(1,:));

If I increase the mesh resolution, the points become clearer to see (right figure).

Screenshot 2025-01-24 at 12.55.00
Screenshot 2025-01-24 at 12.55.001920×636 165 KB

Could this be due to incorrect plotting?
Thanks.

2

You are plotting values from the entire domain projected onto the xy plane, so you see values from the bottom of the domain up to the top of the domain all together. Because the displacement field is not uniform with depth, this introduces variability. If you just plot the values on the ground surface (step02_varslip-groundsurf.h5), you will see a smooth displacement field.

Note: PyLith v4.1.1 contains a bug that introduces inconsistent fault orientations. This is fixed in v4.1.2 and v4.2.0 (but not v4.1.3). We strongly recommend updating to the most recent release (v4.2.0).

3

Thank you for your reply.