Consider stagnation case again

stella-bourdin · Sep 11, 2024 · 3c40bff · 3c40bff
1 parent 64e554d
commit 3c40bff
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/CPyS/theta.py b/CPyS/theta.py
@@ -20,8 +20,11 @@ def theta(x0=120, x1=130, y0=12, y1=10):
         0 for eastward, 90 for northward.
     """
     geodesic = pyproj.Geod(ellps="WGS84")
-    fwd_azimuth, back_azimuth, distance = geodesic.inv(x0, y0, x1, y1)
-    return -1 * (fwd_azimuth - 90) % 360
+    if (x1-x0, y1-y0) != (0,0):
+        fwd_azimuth, back_azimuth, distance = geodesic.inv(x0, y0, x1, y1)
+        return -1 * (fwd_azimuth - 90) % 360
+    else:
+        return np.nan
 
 
 def theta_track(lon, lat):