Overlay IEM NEXRAD composite on Cartopy

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()