Plot a cross section from a PPI volume#

An example which extracts a cross section at two azimuth angles from a volume of PPI scans and plots both cross sections.

xsapr-sg 225.0 Deg. 2011-05-20T06:42:20Z  Equivalent reflectivity factor, xsapr-sg 270.0 Deg. 2011-05-20T06:42:22Z  Equivalent reflectivity factor
print(__doc__)

# Author: Jonathan J. Helmus (jhelmus@anl.gov)
# License: BSD 3 clause

import matplotlib.pyplot as plt
import pyart
from pyart.testing import get_test_data

# Read the data, a cfradial file
filename = get_test_data('swx_20120520_0641.nc')
radar = pyart.io.read(filename)

# Create a cross section at 225 and 270 degrees azimuth
xsect = pyart.util.cross_section_ppi(radar, [225, 270])

# Set the colorbar label
colorbar_label = 'Equivalent \n reflectivity factor \n (dBZ)'

display = pyart.graph.RadarDisplay(xsect)
fig = plt.figure()
ax1 = fig.add_subplot(211)
display.plot('reflectivity_horizontal', 0, vmin=-32, vmax=64., colorbar_label=colorbar_label)
plt.ylim(0, 15)
ax2 = fig.add_subplot(212)
display.plot('reflectivity_horizontal', 1, vmin=-32, vmax=64., colorbar_label=colorbar_label)
plt.ylim(0, 15)

plt.tight_layout()
plt.show()

Total running time of the script: ( 0 minutes 1.725 seconds)

Gallery generated by Sphinx-Gallery