Skip to content

Commit 3d96d79

Browse files
remove from fast data type
1 parent ee152ec commit 3d96d79

File tree

1 file changed

+13
-17
lines changed

1 file changed

+13
-17
lines changed

climada/hazard/tc_tracks.py

Lines changed: 13 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -1644,7 +1644,7 @@ def from_FAST(cls, folder_name: str):
16441644
"""
16451645

16461646
LOGGER.info("Reading %s files.", len(get_file_names(folder_name)))
1647-
data: list = []
1647+
data = []
16481648
for file in get_file_names(folder_name):
16491649
if Path(file).suffix != ".nc":
16501650
continue
@@ -1653,32 +1653,28 @@ def from_FAST(cls, folder_name: str):
16531653
for i in dataset.n_trk:
16541654

16551655
# Select track
1656-
track: xr.Dataset = dataset.sel(n_trk=i, year=year)
1656+
track = dataset.sel(n_trk=i, year=year)
16571657
# chunk dataset at first NaN value
1658-
lon: np.ndarray = track.lon_trks.data
1659-
last_valid_index: int = np.where(np.isfinite(lon))[0][-1]
1660-
track: xr.Dataset = track.isel(
1661-
time=slice(0, last_valid_index + 1)
1662-
)
1658+
lon = track.lon_trks.data
1659+
last_valid_index = np.where(np.isfinite(lon))[0][-1]
1660+
track = track.isel(time=slice(0, last_valid_index + 1))
16631661
# Select lat, lon
1664-
lat: np.ndarray = track.lat_trks.data
1665-
lon: np.ndarray = track.lon_trks.data
1662+
lat = track.lat_trks.data
1663+
lon = track.lon_trks.data
16661664
# Convert lon from 0-360 to -180 - 180
1667-
lon: np.ndarray = ((lon + 180) % 360) - 180
1665+
lon = ((lon + 180) % 360) - 180
16681666
# Convert time to pandas Datetime "yyyy.mm.dd"
16691667
reference_time = (
16701668
f"{track.tc_years.item()}-{int(track.tc_month.item())}-01"
16711669
)
1672-
time: np.datetime64 = pd.to_datetime(
1670+
time = pd.to_datetime(
16731671
track.time.data, unit="s", origin=reference_time
16741672
).astype("datetime64[s]")
16751673
# Define variables
1676-
ms_to_kn: float = 1.943844
1677-
max_wind_kn: np.ndarray = track.vmax_trks.data * ms_to_kn
1678-
env_pressure: float = BASIN_ENV_PRESSURE[
1679-
track.tc_basins.data.item()
1680-
]
1681-
cen_pres: np.ndarray = _estimate_pressure(
1674+
ms_to_kn = 1.943844
1675+
max_wind_kn = track.vmax_trks.data * ms_to_kn
1676+
env_pressure = BASIN_ENV_PRESSURE[track.tc_basins.data.item()]
1677+
cen_pres = _estimate_pressure(
16821678
np.full(lat.shape, np.nan),
16831679
lat,
16841680
lon,

0 commit comments

Comments
 (0)