diff --git a/lib/xmet/sounding.py b/lib/xmet/sounding.py index feef1af..6fd8fe1 100644 --- a/lib/xmet/sounding.py +++ b/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