Commit 123429f6 authored by Carl Schreck's avatar Carl Schreck

Switching to cjs.conus_projection and adding shaded relief.

parent 1e0b217a
......@@ -64,16 +64,9 @@ anom_data = anom_data.where(total_ds['prcp'][0].notnull())
cjs.tstamp('Draw map')
ax = plt.axes(projection=ccrs.AlbersEqualArea(central_longitude=360-95,
standard_parallels=(40,60)))
plt.gca().outline_patch.set_visible(False)
ax.set_extent([-118, -75, 23.5, 50.5])
shp_name = 'ne_10m_admin_1_states_provinces_lakes'
cjs.add_map(ax, shp_name, edgecolor='none', facecolor='silver', zorder=1,
usa_only=True)
cjs.add_map(ax, shp_name, edgecolor=(0.2, 0.2, 0.2), facecolor='none', zorder=3,
usa_only=True)
ax = plt.axes(projection=cjs.conus_projection)
cjs.add_conus_map(ax=ax)
cjs.add_shaded_relief(ax=ax)
cjs.tstamp('Draw data')
levels = [5, 10, 25, 50, 75, 100, 125, 150, 200, 300, 500]
......
......@@ -66,18 +66,8 @@ for first_day, last_day in zip(first_days, last_days):
anom_data = anom_data.where(total_ds['prcp'][0].notnull())
cjs.tstamp('Draw map')
ax = plt.axes(projection=ccrs.AlbersEqualArea(central_longitude=360 - 95,
standard_parallels=(40, 60)))
plt.gca().outline_patch.set_visible(False)
ax.set_extent([-118, -75, 23.5, 50.5])
shp_name = 'ne_10m_admin_1_states_provinces_lakes'
cjs.add_map(ax, shp_name, edgecolor='none', facecolor='silver', zorder=1,
usa_only=True)
cjs.add_map(ax, shp_name, edgecolor=(0.2, 0.2, 0.2), facecolor='none',
zorder=3,
usa_only=True)
ax = plt.axes(projection=cjs.conus_projection)
cjs.add_conus_map(ax=ax)
cjs.tstamp('Draw data')
levels = [5, 10, 25, 50, 75, 100, 125, 150, 200, 300, 500]
......
......@@ -58,16 +58,9 @@ cjs.tstamp('Anomalize')
anom_data = (total_data - clim_data) * 9 / 5
cjs.tstamp('Draw map')
ax = plt.axes(projection=ccrs.AlbersEqualArea(central_longitude=360-95,
standard_parallels=(40,60)))
plt.gca().outline_patch.set_visible(False)
ax.set_extent([-118, -75, 23.5, 50.5])
shp_name = 'ne_10m_admin_1_states_provinces_lakes'
cjs.add_map(ax, shp_name, edgecolor='none', facecolor='silver', zorder=1,
usa_only=True)
cjs.add_map(ax, shp_name, edgecolor=(0.2, 0.2, 0.2), facecolor='none', zorder=3,
usa_only=True)
ax = plt.axes(projection=cjs.conus_projection)
cjs.add_conus_map(ax=ax)
cjs.add_shaded_relief(ax=ax)
cjs.tstamp('Draw data')
levels = list(np.arange(-5,6)*3)
......
......@@ -62,17 +62,8 @@ for delta_days in [1, 5]:
anom_data = (total_data - clim_data) * 9 / 5
cjs.tstamp('Draw map')
ax = plt.axes(projection=ccrs.AlbersEqualArea(central_longitude=360-95,
standard_parallels=(40,60)))
plt.gca().outline_patch.set_visible(False)
ax.set_extent([-118, -75, 23.5, 50.5])
shp_name = 'ne_10m_admin_1_states_provinces_lakes'
cjs.add_map(ax, shp_name, edgecolor='none', facecolor='silver', zorder=1,
usa_only=True)
cjs.add_map(ax, shp_name, edgecolor=(0.2, 0.2, 0.2), facecolor='none',
zorder=3,
usa_only=True)
ax = plt.axes(projection=cjs.conus_projection)
cjs.add_conus_map(ax=ax)
cjs.tstamp('Draw data')
levels = list(np.arange(-5,6)*3)
......
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