Commit ab64cad9 authored by Carl Schreck's avatar Carl Schreck

Added US_MSR_10M shaded relief

parent 25cebbca
......@@ -32,34 +32,36 @@ cmap = 'RdBu_r'
cjs.tstamp('Here we go!')
fig, axes = plt.subplots(nrows=1, ncols=1, figsize=(6.4, 6.4),
subplot_kw={'projection': cjs.conus_projection, })
for tif_name in ['SR_50M', 'MSR_50M', 'PRISMA_SR_50M', 'US_MSR']:
# for tif_name in ['US_MSR']:
fig, axes = plt.subplots(nrows=1, ncols=1, figsize=(6.4, 6.4),
subplot_kw={'projection': cjs.conus_projection, })
cjs.tstamp('Drawing scaled data')
ax = axes
helpers.draw_station_data(ax, station_path, levels=levels, cmap=cmap)
scaled_ds = xr.open_dataset(scaled_path)
scaled_data = scaled_ds[var_name].sel(time=time).drop('time')
plot = scaled_data.plot(ax=ax, transform=ccrs.PlateCarree(),
add_colorbar=True,
cmap=cmap, levels=levels, zorder=1,
vmin=-np.inf, vmax=np.inf, extend='both',
cbar_kwargs={'orientation': 'horizontal',
'label': '°C',
'pad': 0.05,
}
)
ax.set_title('c) Scaled data', loc='left')
ax.set_extent([-115.0, -100.0, 35.0, 42.5])
ax.outline_patch.set_visible(True)
ax.set_anchor('N')
cjs.tstamp('Drawing scaled data')
ax = axes
helpers.draw_station_data(ax, station_path, levels=levels, cmap=cmap)
scaled_ds = xr.open_dataset(scaled_path)
scaled_data = scaled_ds[var_name].sel(time=time).drop('time')
plot = scaled_data.plot(ax=ax, transform=ccrs.PlateCarree(),
add_colorbar=True,
cmap=cmap, levels=levels, zorder=1,
vmin=-np.inf, vmax=np.inf, extend='both',
cbar_kwargs={'orientation': 'horizontal',
'label': '°C',
'pad': 0.05,
}
)
ax.set_title('c) Scaled data', loc='left')
ax.set_extent([-120.0, -100.0, 35.0, 45.0])
ax.outline_patch.set_visible(True)
ax.set_anchor('N')
trans_gray = cjs.add_shaded_relief(ax)
trans_gray = cjs.add_shaded_relief(ax, tif_name=tif_name)
cjs.tstamp('Drawing figure')
plt.tight_layout()
plt.savefig('figures/temp_steps.png', dpi=300)
plt.show()
cjs.tstamp('Drawing figure')
plt.tight_layout()
plt.savefig(f'figures/{tif_name}.png', dpi=300)
plt.show()
cjs.tstamp('Thank you, come again.')
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