Commit 5582fac2 authored by Carl Schreck's avatar Carl Schreck

Added SS category to individual storm output

parent 27a747de
......@@ -18,7 +18,7 @@ cjs.tstamp('Here we go!')
# These are some parameters that could be useful to have up top
basins = ['NA', 'EP']
parser = argparse.ArgumentParser()
parser.add_argument('--year', default=datetime.today().year)
parser.add_argument('--year', default=2019)
parser.add_argument('--month', default=(datetime.today().month - 1))
args = parser.parse_args()
end_year = int(args.year)
......@@ -98,10 +98,13 @@ for basin in basins:
& ibt.check_basin(basin))
any_eligible = np.any(is_eligible, axis=1)
for storm in ibt.ds.storm[any_eligible]:
print(f'name: {ibt.ds.name[storm].item().decode().title()}')
print(f'id: {ibt.ds.sid[storm].item().decode().title()}')
max_wind = ibt.ds.usa_wind[storm,is_eligible[storm,:]].max().values
min_pres = ibt.ds.usa_pres[storm,is_eligible[storm,:]].min().values
nature = ibt.ds.nature[
storm,max_wind==ibt.ds.usa_wind[storm,:]][0].values
print(f'id: {ibt.ds.sid[storm].item().decode().title()}')
print(f'name: {ibt.ds.name[storm].item().decode().title()}')
print(f'category: {ibt.ss_cat(max_wind,nature)}')
print(f'max wind: {max_wind} kt / {max_wind * 1.15} mph')
print(f'min pressure: {min_pres}')
min_time_ind = np.where(is_eligible[storm,:])[0].min()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment