From 0d767268d704a37d4f46a8e0c8e43eac5d0ea709 Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Fri, 16 Oct 2020 12:02:13 -0500 Subject: [PATCH] snapshot... --- modules/amv/intercompare.py | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index 412c682f..afe68d6a 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -1125,13 +1125,15 @@ def best_fit_altitude(amv_spd, amv_dir, amv_alt, amv_lat, amv_lon, fcst_spd, fcs if a2 < sat_wind_best_fit_alt: lev_below = imin lev_above = imin + 1 - Prop = (sat_wind_best_fit_alt - a2) / (a3 - a2) + prop = (sat_wind_best_fit_alt - a2) / (a3 - a2) else: + lev_below = imin - 1 + lev_above = imin + prop = (sat_wind_best_fit_alt - a1) / (a2 - a1) pass - sat_wind_best_fit_u = fcst_uwind[lev_below] * (1.0 - Prop) + fcst_uwind[lev_above] * Prop - sat_wind_best_fit_v = fcst_vwind[lev_below] * (1.0 - Prop) + fcst_vwind[lev_above] * Prop - + sat_wind_best_fit_u = fcst_uwind[lev_below] * (1.0 - prop) + fcst_uwind[lev_above] * prop + sat_wind_best_fit_v = fcst_vwind[lev_below] * (1.0 - prop) + fcst_vwind[lev_above] * prop return bf_tup -- GitLab