import os import pyart import artview #from artview.view import view import matplotlib.pyplot as plt %matplotlib inline fdir = '/Users/guy/data/nexrad/KMPX/nex2/20130513' fname = 'KMPX20130513_185823_V06' fpath = os.path.join(fdir, fname) print(fpath) r = pyart.io.read(fpath) r.fields.keys() disp = pyart.graph.RadarDisplay(r) fig = plt.figure() ax = fig.add_subplot(111) disp.plot('reflectivity', 0, vmin=-32, vmax=64., cmap='pyart_NWSRef') disp.plot_range_rings([50, 100]) disp.plot_cross_hair(5.) ax.set_xlim(-150., 150) ax.set_ylim(-150., 150) artview.view.view(r)