.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/plotting/plot_rhi_two_panel.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_plotting_plot_rhi_two_panel.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_plotting_plot_rhi_two_panel.py:


===========================
Create a two panel RHI plot
===========================

An example which creates a two panel RHI plot of a cfradial file.  The fields
included in the two panels are reflectivity and doppler velocity.

.. GENERATED FROM PYTHON SOURCE LINES 10-54



.. image-sg:: /examples/plotting/images/sphx_glr_plot_rhi_two_panel_001.png
   :alt: Reflectivity and Velocity   Azimuth: 37.299°  2011-05-24 01:56:04.0 UTC
   :srcset: /examples/plotting/images/sphx_glr_plot_rhi_two_panel_001.png
   :class: sphx-glr-single-img





.. code-block:: Python

    print(__doc__)

    # Author: Max Grover (mgrover@anl.gov)
    # License: BSD 3 clause

    import matplotlib.pyplot as plt
    import numpy as np

    import pyart
    from pyart.testing import get_test_data

    # Read the data and create the display object
    filename = get_test_data("sgpxsaprrhicmacI5.c0.20110524.015604_NC4.nc")
    radar = pyart.io.read_cfradial(filename)
    display = pyart.graph.RadarDisplay(radar)

    # Fields to plot and ranges
    fields_to_plot = ["reflectivity_horizontal", "mean_doppler_velocity"]
    ranges = [(-20, 20), (-17.0, 17.0)]
    cmaps = ["pyart_HomeyerRainbow", "pyart_balance"]

    # Plot the data
    nplots = len(fields_to_plot)
    plt.figure(figsize=[5 * nplots, 4])

    # Plot each field
    for plot_num in range(nplots):
        field = fields_to_plot[plot_num]
        vmin, vmax = ranges[plot_num]
        cmap = cmaps[plot_num]

        plt.subplot(1, nplots, plot_num + 1)
        display.plot(field, 0, vmin=vmin, vmax=vmax, title_flag=False, cmap=cmap)
        display.set_limits(ylim=[0, 17])

    # Grab the fixed angle and time from the first sweep
    fixed_angle = radar.fixed_angle["data"][0]
    time = radar.time["units"][13:]

    # Add the metadata to the title
    plt.suptitle(
        f"Reflectivity and Velocity \n Azimuth: {np.around(fixed_angle, 3)}\u00B0 {time} UTC"
    )
    plt.show()


.. rst-class:: sphx-glr-timing

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


.. _sphx_glr_download_examples_plotting_plot_rhi_two_panel.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: plot_rhi_two_panel.ipynb <plot_rhi_two_panel.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: plot_rhi_two_panel.py <plot_rhi_two_panel.py>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_