Fault output
For information on specifying faults in PyLith v3 and later, I recommend the PyLith tutorials on YoutTube (for example, https://www.youtube.com/watch?v=fcHKL8PqFDk&list=PLNTb1ey7FFk7aaP7meQlevDqahYa-7S1h&index=3) and Horizontal Cross-Section of Strike-Slip Fault (2D) — PyLith 4.2.0dev documentation.
In PyLith v3 and later, in most cases you do not need to do as much for output. In examples/strikeslip-2d, we specify the fields we want output (we can select from slip, traction change, and any solution fields) in pylithapp.cfg:
# Output `slip` and the change in tractions on the fault.
observers.observer.data_fields = [slip, traction_change]
Poroelasticity solver error
Changing the material properties may be creating an ill-conditioned system due to the inconsistency between the material properties and the scales used for nondimensionalization. For a well-conditioned system for poroelasticity, the ratio of the nondimensional permeability to viscosity should be about 1. Assuming the stress and length scales are appropriate for your problem, then you would adjust the time scale so that the nondimensional permeability to viscosity ratio is about 1.