Revert default none stn status & fixstyling

oceanmodeling · Aug 23, 2024 · 39b2962 · 39b2962
1 parent 475fd6e
commit 39b2962
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 9 deletions.
diff --git a/stormevents/nhc/track.py b/stormevents/nhc/track.py
@@ -1134,7 +1134,7 @@ def __compute_velocity(data: DataFrame) -> DataFrame:
                 .astype("timedelta64[s]")
                 .astype(float)
             )
-            with numpy.errstate(invalid='ignore'):
+            with numpy.errstate(invalid="ignore"):
                 speeds = pandas.Series(distances / abs(intervals), index=indices)
             bearings = pandas.Series(inverse_azimuths % 360, index=indices)
             # use forward azimuths for negative intervals

diff --git a/stormevents/stormevent.py b/stormevents/stormevent.py
@@ -432,7 +432,7 @@ def coops_product_within_region(
         region: Polygon,
         start_date: datetime = None,
         end_date: datetime = None,
-        status: StationStatus = StationStatus.ACTIVE,
+        status: StationStatus = None,
         datum: COOPS_TidalDatum = None,
         units: COOPS_Units = None,
         time_zone: COOPS_TimeZone = None,
@@ -491,10 +491,9 @@ def coops_product_within_region(
 
         stations = gpd.GeoDataFrame()
         if not region.is_empty:
-            stations = get_coops_stations(
-                region=region, metadata_source='main'
-            )
-            stations = stations[stations.status == status.value]
+            stations = get_coops_stations(region=region, metadata_source="main")
+            if status is not None:
+                stations = stations[stations.status == status.value]
 
         if len(stations) > 0:
             stations_data = []
@@ -508,12 +507,10 @@ def coops_product_within_region(
                     datum=datum,
                 ).data
 
-
                 station_data["nos_id"] = station.Index
                 station_data.set_index(["nos_id", station_data.index], inplace=True)
                 station_data = station_data.to_xarray()
 
-
                 if len(station_data["t"]) > 0:
                     station_data = station_data.assign_coords(
                         nws_id=("nos_id", [station.nws_id]),
@@ -535,7 +532,6 @@ def coops_product_within_region(
                 coords={"t": None, "nos_id": None, "nws_id": None, "x": None, "y": None}
             )
 
-
         return stations_ds
 
     def __repr__(self) -> str: