More description incoming soon. I had the idea to create surfaces for faults using GemPy, using the multiple faults per scalarfield. This means that fewer scalarfields have to be computed, but also, that faults from one scalarfield will not cross cut each other.