From 09ee041fb1773151f43b6246e75db148901b261e Mon Sep 17 00:00:00 2001 From: William Pringle Date: Mon, 5 Aug 2024 12:59:20 -0500 Subject: [PATCH] moving carq_ref to above holland B calculation and just using that to compute it instead of needing to take the nanmean --- stormevents/nhc/track.py | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) diff --git a/stormevents/nhc/track.py b/stormevents/nhc/track.py index 7bf6f31..87eb2e1 100644 --- a/stormevents/nhc/track.py +++ b/stormevents/nhc/track.py @@ -1270,16 +1270,17 @@ def clamp(n, minn, maxn): else: carq_forecast = carq_tracks[list(carq_tracks)[0]] + # Get CARQ from forecast hour 0 and isotach 34kt (i.e. the first item) + carq_ref = carq_forecast.loc[carq_forecast.forecast_hours == 0].iloc[0] + + # get the Holland B parameter for filling in central pressure relation = HollandBRelation() holland_b = relation.holland_b( - max_sustained_wind_speed=carq_forecast["max_sustained_wind_speed"], - background_pressure=carq_forecast["background_pressure"], - central_pressure=carq_forecast["central_pressure"], + max_sustained_wind_speed=carq_ref["max_sustained_wind_speed"], + background_pressure=carq_ref["background_pressure"], + central_pressure=carq_ref["central_pressure"], ) - holland_b[holland_b == numpy.inf] = numpy.nan - holland_b = numpy.nanmean(holland_b) - # Get CARQ from forecast hour 0 and isotach 34kt (i.e. the first item) carq_ref = carq_forecast.loc[carq_forecast.forecast_hours == 0].iloc[0]