fix dip calculation

ucgmsim · Jul 2, 2024 · 8a1c48d · 8a1c48d
1 parent 4b3106b
commit 8a1c48d
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/nshmdb/fault.py b/nshmdb/fault.py
@@ -18,9 +18,9 @@
 
 import dataclasses
 from enum import Enum
+from typing import Optional
 
 import numpy as np
-from typing import Optional
 from qcore import coordinates, geo
 
 
@@ -110,6 +110,10 @@ def bottom_m(self) -> float:
         """
         return self.corners_nztm[-1, -1]
 
+    @property
+    def top_m(self) -> float:
+        return self.corners_nztm[0, -1]
+
     @property
     def width(self) -> float:
         """
@@ -193,7 +197,7 @@ def dip(self) -> float:
         float
             The dip angle of the fault.
         """
-        return np.degrees(np.arcsin(np.abs(self.bottom_m) / self.width_m))
+        return np.degrees(np.arcsin(np.abs(self.bottom_m - self.top_m) / self.width_m))
 
     def plane_coordinates_to_global_coordinates(
         self, plane_coordinates: np.ndarray