From df9ac5605011e3d03601dae38b5b4150c270c464 Mon Sep 17 00:00:00 2001 From: Aki Date: Fri, 31 Mar 2023 23:06:18 +0200 Subject: Take wind criterion into account when calculating SWI --- waterspout-radar.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/waterspout-radar.py b/waterspout-radar.py index c6d5ad5..ec58e64 100644 --- a/waterspout-radar.py +++ b/waterspout-radar.py @@ -35,7 +35,7 @@ def main(): el = calc.pressure_to_height_std(pressure) ccd = (el - lcl).to(units.ft) try: - swi = szilagyi.calculate_swi(dt, ccd) + swi = szilagyi.calculate_swi(dt, ccd, pred.at('wind_u', _L.H850)) # merge with wind_v except ValueError: swi = -10 print(pred.timestamp, dt, ccd, swi) -- cgit v1.1