3d. Updating fault geometry

import numpy as np
import pandas as pd
import LoopStructural as LS

# Define a dataset for a fault

origin = [0, 0, 0]
extent = [10, 10, 10]

data = pd.DataFrame(
    [
        [5, 5, 5, 0, 0.70710678, 0.0, 0.70710678, 0, "fault"],
        [5, 5, 5, 0, -0.70710678, 0.0, 0.70710678, 1, "fault"],
        [8, 5, 5, 0, 0, 0, 1, np.nan, "strati"],
    ],
    columns=["X", "Y", "Z", "val", "nx", "ny", "nz", "coord", "feature_name"],
)

data

	X	Y	Z	val	nx	ny	nz	coord	feature_name
0	5	5	5	0	0.707107	0.0	0.707107	0.0	fault
1	5	5	5	0	-0.707107	0.0	0.707107	1.0	fault
2	8	5	5	0	0.000000	0.0	1.000000	NaN	strati

Create model using the standard fault displacement model

model = LS.GeologicalModel(origin, extent)
model.data = data
model.create_and_add_fault(
    "fault",
    10,
    nelements=1000,
    interpolator_type="PLI",
    buffer=0.5,
    major_axis=10,
    minor_axis=3,
    intermediate_axis=10,
)
model.create_and_add_foliation(
    "strati", nelements=1000, interpolator_type="PLI", faults=[model["fault"]]
)
model.update()
print(model['fault'].builder.fault_minor_axis, model['fault'].builder.up_to_date)

model['fault'].builder.fault_minor_axis = 6.0
print(model['fault'].builder.fault_minor_axis, model['fault'].builder.up_to_date)
model.update()
print(model['fault'].builder.fault_minor_axis, model['fault'].builder.up_to_date)

# import LoopStructural.visualisation as vis

# view = vis.Loop3DView(model)
# view.plot_surface(model.features[0], value=[0])
# view.plot_surface(model.features[1], value=5, paint_with=model.features[1])
# # view.add_vector_field(model["fault"][1], locations=model.regular_grid()[::100])

# view.display()

# ######################################################################
# # Define a fault displacement profile which
# # is a drag fault only on the footwall side.
# # In LoopStructural the displacement is defined by a function of the three
# # coordinates of the fault frame.
# # The fault profile in the fault surface field

# model['fault'].faultfunction.gx.plot()

# ######################################################################
# # The fault profile in the fault extent
# model['fault'].faultfunction.gy.plot()


# ######################################################################
# # The fault profile down dip is kept constant.
# # We will modify this profile so that the hanging wall is displaced by a constant value

# from LoopStructural.modelling.features.fault._fault_function import (
#     FaultDisplacement,
#     CubicFunction,
#     Ones,
# )

# fw = CubicFunction()
# fw.add_cstr(0, -1)
# fw.add_grad(0, 0)
# fw.add_cstr(-1, 0)
# fw.add_grad(-1, 0)
# fw.add_min(-1)
# hw = Ones()
# drag_fault = FaultDisplacement(hw=hw, fw=fw)

# drag_fault.gx.plot()
# drag_fault.gy.plot()
# drag_fault.gz.plot()

# model = LS.GeologicalModel(origin, extent)
# model.data = data
# model.create_and_add_fault(
#     "fault",
#     -1,
#     nelements=1000,
#     interpolator_type="PLI",
#     buffer=0.5,
#     major_axis=10,
#     minor_axis=6,
#     intermediate_axis=10,
#     faultfunction=drag_fault,
# )
# model.create_and_add_foliation(
#     "strati", nelements=1000, interpolator_type="PLI", faults=[model["fault"]]
# )


# view = vis.Loop3DView(model)
# model.bounding_box.nelements = 1e5
# view.plot_surface(model.features[0], value=[0])
# view.plot_surface(model['strati'], value=5)

# view.display()

3.0 <bound method FaultBuilder.up_to_date of <LoopStructural.modelling.features.builders._fault_builder.FaultBuilder object at 0x7f92289d7380>>
6.0 <bound method FaultBuilder.up_to_date of <LoopStructural.modelling.features.builders._fault_builder.FaultBuilder object at 0x7f92289d7380>>
6.0 <bound method FaultBuilder.up_to_date of <LoopStructural.modelling.features.builders._fault_builder.FaultBuilder object at 0x7f92289d7380>>