deeplycloudy/iem_nexrad_wms_cartopy.py

## iem_nexrad_wms_cartopy.py
""" Overlay IEM NEXRAD composite on Cartopy for any date or time in the
long-running IEM archive. This works for any target projection, with the
reprojection handled by the Cartopy WMS functionality.

If you don't enter a time on an even 5 min boundary, you will get a blank map.
"""

iem_wms_nexrad = 'https://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r-t.cgi?'
layers = ["nexrad-n0r-wmst"]
wms_kw =  {"time":"2017-10-22T05:15:00Z", "transparent":'true'}

fig = plt.figure()
ax = fig.add_subplot(111, projection=ccrs.InterruptedGoodeHomolosine())
ax.coastlines()

ax.add_wms(iem_wms_nexrad, layers, wms_kwargs=wms_kw)

plt.show()
	""" Overlay IEM NEXRAD composite on Cartopy for any date or time in the
	long-running IEM archive. This works for any target projection, with the
	reprojection handled by the Cartopy WMS functionality.

	If you don't enter a time on an even 5 min boundary, you will get a blank map.
	"""

	iem_wms_nexrad = 'https://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r-t.cgi?'
	layers = ["nexrad-n0r-wmst"]
	wms_kw = {"time":"2017-10-22T05:15:00Z", "transparent":'true'}

	fig = plt.figure()
	ax = fig.add_subplot(111, projection=ccrs.InterruptedGoodeHomolosine())
	ax.coastlines()

	ax.add_wms(iem_wms_nexrad, layers, wms_kwargs=wms_kw)

	plt.show()