oceanmodeling / StormEvents

Python interfaces for observational data surrounding named storm events, born from @jreniel's ADCIRCpy

Home Page:https://stormevents.readthedocs.io

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Forecast data for active storm is filtered by `end_date`

SorooshMani-NOAA opened this issue · comments

While one can still use unfiltered_data instead of data this issue results in not being able to get forecast windswath and isotacks for an active storm

The data returns

return self.unfiltered_data.loc[
(self.unfiltered_data["datetime"] >= self.start_date)
& (self.unfiltered_data["datetime"] <= self.end_date)

The issue is when creating a StormEvent object, the end_date (which is not available for active storms) is taken from VortexTrack object created ad-hoc by VortexTrack.from_storm_name(...); this track object is based b-deck and doesn't have the latest forecast advisory date. So in the VortexTrack.end_date the unfiltered_data.iloc[-1] returns the latest advisory up to now (best track).

Also it might be better to use

data_start = min(self.unfiltered_data['datetime'].values)
data_end = max(self.unfiltered_data['datetime'].values)

instead of calling iloc with -1 and 0

As a workaround (at least in my specific use case) one can start from the VortexTrack object. This is what I was trying to do that failed:

storm = StormEvent('fiona', 2022)
track = storm.track(file_deck='a')
track.wind_swaths(wind_speed=34)

which resulted in up-to-now wind swath only. Instead I can replace it with:

track = VortexTrack.from_storm_name('fiona', 2022, file_deck='a', advisories=['OFCL'])
track.wind_swaths(wind_speed=34)