Dear PyLith team,
I want to simulate post-seismic deformation, and the fault geometry is shown in the figure. My PyLith version is 2.2.2. Based on my current understanding, the fault planes should be internal interfaces of the blocks. Could you please provide me with some advice on how to divide the blocks?
Thank you！
Is there a reason why you are using PyLith v2.2.2 instead of a recent release? We just released version 5.0.1, and as long as you don’t need fault friction, it will give much better performance.
For the fault geometry shown in the figure, using the Gmsh mesh generator (compatible with recent PyLith releases, you can embed the faults in the volume. Refer to examples/crustal-strikeslip-3 for a demonstration of how to do this.
Thank you for your reply!
