FixAzimuthCalculation

cgre-aachen · Oct 17, 2024 · 023f64b · 023f64b
1 parent 01bce2e
commit 023f64b
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/gemgis/postprocessing.py b/gemgis/postprocessing.py
@@ -566,9 +566,12 @@ def calculate_dip_and_azimuth_from_mesh(
 
     # Calculating the azimuths
     azimuths = [
-        np.rad2deg(np.arctan(mesh["Normals"][i][0] / mesh["Normals"][i][1])) + 180
+        np.rad2deg(np.arctan2(mesh["Normals"][i][0], mesh["Normals"][i][1]))
         for i in range(len(mesh["Normals"]))
     ]
+
+    # Shifting values
+    azimuths[azimuths < 0] += 360
 
     # Assigning dips and azimuths to scalars
     mesh["Dips [°]"] = dips