Can I graph finite faults in Gempy, so that they don't cut off my entire strata series?

[shairethomg]

I'm trying to make this fault only cut me up to the second stratum, and I've tried placing the series in the order of the succession of events but it doesn't work. I also tried to modify the correlation matrix. Can someone please tell me how to graph the finite faults?
thanks.

[javoha]

Hi,
thanks for your question. Just to make sure: You want the fault to only displace the mioceno_medio element, right? mioceno_superior and plioceno_inferior are supposed to be unaffected?
Actually you did everything correct and the result is exactly what you want then. The only probelm is that gempy atm always displays the fautls over the whole thing. This is just a visualization artefact. In your model the top two units are not actually displaced by the fault.
There are ways to change the visualization using some pyvista magic I guess if that is really important to you.
Let me know if that helps or if I got the question wrong.
Cheers,
Jan

[shairethomg]

Hi, that's right, that's my question, thank you very much.
Do you know how you could use Pyvista to improve this?

[NilsChudalla]

Hi @shairethomg ,

if it comes to Pyvista usage, you might stumble upon boolean operations (Example: https://docs.pyvista.org/examples/01-filter/boolean-operations.html). Those might be not sufficient for you, as they require manifold/watertight geometries, which these surface are not. If it is a onetime time thing, you could export the scene of your pyvista plotter object with
plotter = gpv.plot_3d(geo_model)
plotter.p.export_obj("name_of_scene.obj").
and then remove the vertices in some modeling software such as Blender. If you have to do it multiple times, I cannot think of any other way than creating a manifold object (e.g. pv.Box(model xmin, model xmax, model ymin, model ymax, depth of intersection, model zmax)) and doing the boolean operation with that.

[shairethomg]

Thank you very much!