diff --git a/GMW/gmw.py b/GMW/gmw.py index d4e24d9..408eff5 100644 --- a/GMW/gmw.py +++ b/GMW/gmw.py @@ -66,7 +66,7 @@ def get_mangrove_locations_from_tiles(gmwdir: str, tilenames: list[str]) -> np.n return np.vstack([latitude, longitude]).T -# TODO: actually return BoundingBox rather than Polygon +# TODO: actually return BoundingBox rather than Polygon, once debugged def get_tiles(tilenames: list[str]) -> list[BoundingBox]: """Return a list of Bounding Boxes representing the 1x1 degree tiles containing mangroves.""" return [Polygon(get_tile_corners(tn)) for tn in tilenames] diff --git a/windspeedcomparison.py b/windspeedcomparison.py index 2038cee..1187efd 100644 --- a/windspeedcomparison.py +++ b/windspeedcomparison.py @@ -1,6 +1,6 @@ """ -This script performs a comparison between maximum wind speeds from different storms as reported in IBTrACS and ERA5 -datasets. +This script performs a comparison between maximum wind speeds from different storms as +reported in IBTrACS and ERA5 datasets. """ import cdsapi @@ -15,8 +15,11 @@ from Spherical.functions import addtolon -mps2kts = 1.94384 +mps2kts = 1.94384 # unit conversion ratio + +# where to save the data from each storm savepath = "/Users/fcseidl/EarthLab-local/BioExtremes/WindComparisonBasins.csv" + np.random.seed(0) # IBTrACS data where max wind and basin are recorded