fix: if downthrow not provided estimate abutting direction

Loop3D · Apr 28, 2022 · 5a7c8d5 · 5a7c8d5
1 parent bf272d1
commit 5a7c8d5
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 10 deletions.
diff --git a/LoopStructural/modelling/core/geological_model.py b/LoopStructural/modelling/core/geological_model.py
@@ -319,17 +319,25 @@ def from_processor(cls, processor):
                         model[edge[1]].splay[model[edge[0]].name] = region
                         splay = True
                 if splay == False:
+                    positive = None
+                    if (
+                        "downthrow_dir"
+                        in processor.stratigraphic_column["faults"][edge[0]]
+                    ):
+                        positive = (
+                            np.abs(
+                                processor.stratigraphic_column["faults"][edge[0]][
+                                    "downthrow_dir"
+                                ]
+                                - processor.stratigraphic_column["faults"][edge[1]][
+                                    "downthrow_dir"
+                                ]
+                            )
+                            < 90
+                        )
                     model[edge[1]].add_abutting_fault(
                         model[edge[0]],
-                        np.abs(
-                            processor.stratigraphic_column["faults"][edge[0]][
-                                "downthrow_dir"
-                            ]
-                            - processor.stratigraphic_column["faults"][edge[1]][
-                                "downthrow_dir"
-                            ]
-                        )
-                        < 90,
+                        positive=positive,
                     )
         for s in processor.stratigraphic_column.keys():
             if s != "faults":

diff --git a/LoopStructural/modelling/features/fault/_fault_segment.py b/LoopStructural/modelling/features/fault/_fault_segment.py
@@ -400,7 +400,6 @@ def add_abutting_fault(self, abutting_fault_feature, positive=None):
             positive = abut_value > 0
         # we want to crop the fault by the abutting fault so create a positive/neg region and include the fault centre and normal vector to help
         # outside of the fault interpolation support
-
         if positive:
             abutting_region = PositiveRegion(
                 abutting_fault_feature,