Initial implementation of bulk shear

lib/xmet/sounding.py

@@ -229,6 +229,23 @@ class Sounding(DatabaseTable):
             90 - deg(math.atan(vd / ud))
         )
 
+    def bulk_shear(self):
+        levels = len(self.samples) - 1
+
+        shear_speed = 0.0
+        shear_dir   = 0.0
+
+        for level in range(0, levels):
+            shear = self.shear(level)
+
+            shear_speed += shear[0]
+            shear_dir   += shear[1]
+
+        return (
+            shear_speed / levels,
+            shear_dir   / levels
+        )
+
 def between(n, a, b):
     return n > a and n < b