Skip to content

Commit

Permalink
fix: error creating mesh for faults
Browse files Browse the repository at this point in the history
error caused by mesh origin/max
not being updated by fault builder.
  • Loading branch information
Lachlan Grose committed Dec 5, 2021
1 parent 009d2b8 commit 7efe910
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 6 deletions.
6 changes: 3 additions & 3 deletions LoopStructural/modelling/core/geological_model.py
Original file line number Diff line number Diff line change
Expand Up @@ -521,7 +521,7 @@ def data(self, data):
strike_dip_vector(
self._data.loc[mask, "strike"], self._data.loc[mask, "dip"]
)
* self._data.loc[mask, "polarity"][:, None]
* self._data.loc[mask, "polarity"].to_numpy()[:, None]
)
self._data.drop(["strike", "dip"], axis=1, inplace=True)

Expand Down Expand Up @@ -1459,10 +1459,10 @@ def create_and_add_fault(
# fault_frame_builder.origin = self.bounding_box[0,:]
# fault_frame_builder.maximum = self.bounding_box[1,:]
if "force_mesh_geometry" not in kwargs:

fault_frame_builder.set_mesh_geometry(
kwargs.get("fault_buffer", 0.2), 0
) # ,
# np.rad2deg(np.arccos(np.dot(fault_normal_vector[:2],np.array([0,1])))))
)
if "splay" in kwargs and "splayregion" in kwargs:
fault_frame_builder.add_splay(kwargs["splay"], kwargs["splayregion"])

Expand Down
7 changes: 4 additions & 3 deletions LoopStructural/modelling/fault/fault_builder.py
Original file line number Diff line number Diff line change
Expand Up @@ -31,9 +31,10 @@ def __init__(
"""

StructuralFrameBuilder.__init__(self, interpolator, interpolators, **kwargs)
self.origin = np.array([np.nan, np.nan, np.nan])
self.maximum = np.array([np.nan, np.nan, np.nan])
self.model = model
self.model = model
self.origin = self.model.bounding_box[0, :]
self.maximum = self.model.bounding_box[1, :]

# define a maximum area to mesh adding buffer to model
# buffer = .2
self.minimum_origin = self.model.bounding_box[
Expand Down

0 comments on commit 7efe910

Please sign in to comment.