diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index 7c2f13a74a2dc7c3b463ee049da56fa940b84120..0f23e8bb213c3046d73a555b76c790cf18ab1ce6 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -1127,11 +1127,13 @@ def best_fit_altitude(amv_spd, amv_dir, amv_alt, amv_lat, amv_lon, fcst_spd, fcs # check contraints good_constraint = 1 + # Check for overall good agreement: want vdiff less than 4 m/s vdiff = np.sqrt((amv_uwind - sat_wind_best_fit_u) ** 2 + (amv_vwind - sat_wind_best_fit_v) ** 2) if vdiff > 4.0: good_constraint = 0 flag = 2 + # Check for a substantial secondary local minimum or if local minimum is too broad mm = np.where(fcst_alt > (sat_wind_best_fit_alt + constraint_half_width))[0] nn = np.where(fcst_alt < (sat_wind_best_fit_alt - constraint_half_width))[0]