Commit fc85bee1 authored by Carl Schreck's avatar Carl Schreck
Browse files

Final products for August (patched non-zero infinities)

parent e5ee400d
......@@ -20,8 +20,8 @@ from datetime import datetime
# These are some parameters that could be useful to have up top
year = 2019
month = 8
first_day = 1
last_day = 7
first_day = 22
last_day = 30
time_range = slice(f'{year}-{month:02}-{first_day:02}',
f'{year}-{month:02}-{last_day:02}')
......@@ -61,6 +61,8 @@ for shape in shapereader.Reader(shp_path).records():
anom_data = 100. * total_data / clim_data
# Take care of infinities when the observed prcp is nonzero
anom_data = anom_data.where((total_data == 0) | (clim_data != 0), 600)
levels = [5, 10, 25, 50, 75, 100, 125, 150, 200, 300, 500]
plot = anom_data.plot(levels=levels, ax=ax, transform=ccrs.PlateCarree(),
cmap='BrBG', zorder=2,
......
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