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

.. only:: html

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

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

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

.. _sphx_glr_auto_examples_retrievals_plot_sonde_pbl.py:


Planetary Boundary Layer Height Profile via Radiosonde Retrievals
-----------------------------------------------------------------

This examples shoes how to estimate the planetary boundary layer height
via the Liu-Liang 2010 technique, a Bulk Richardson Number method,
and the Heffter technique.

Author: Joe O'Brien

.. GENERATED FROM PYTHON SOURCE LINES 11-65



.. image-sg:: /auto_examples/retrievals/images/sphx_glr_plot_sonde_pbl_001.png
   :alt: sgpsondewnpnC1.b1 on 20190101
   :srcset: /auto_examples/retrievals/images/sphx_glr_plot_sonde_pbl_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <matplotlib.legend.Legend object at 0x7f967fdc9e80>





|

.. code-block:: Python


    from arm_test_data import DATASETS

    import act

    # Read the SGP radiosonde data for an example
    filename_sonde = DATASETS.fetch('sgpsondewnpnC1.b1.20190101.053200.cdf')
    sonde_ds = act.io.arm.read_arm_netcdf(filename_sonde)

    # more explicit units are needed for Potential Temperature calculation
    sonde_ds['tdry'].attrs['units'] = 'degree_Celsius'

    # Estimate PBL Height via the Liu-Liang 2010 technique
    act.retrievals.calculate_pbl_liu_liang(sonde_ds)
    # Estimate PBL Height via the Bulk Richardson Number method
    act.retrievals.calculate_pbl_bulk_richardson(sonde_ds)
    # Estimate PBL Height via the Heffter technique
    act.retrievals.calculate_pbl_heffter(sonde_ds)

    # Calculate what pressure level corresponds to the PBL heights
    pbl_heffter = sonde_ds.where(
        abs(sonde_ds.alt_ss - sonde_ds.pblht_heffter) < 1, drop=True
    ).atm_pres_ss.data[0]
    pbl_liu = sonde_ds.where(
        abs(sonde_ds.alt_ss - sonde_ds.pblht_liu_liang) < 1, drop=True
    ).atm_pres_ss.data[0]
    pbl_rich = sonde_ds.where(
        abs(sonde_ds.alt_ss - sonde_ds.pblht_bulk_richardson_pt25) < 1, drop=True
    ).atm_pres_ss.data[0]

    # Display the Sonde PBL estimates over a Skew-T diagram
    display = act.plotting.SkewTDisplay(sonde_ds)

    display.plot_from_u_and_v(
        'u_wind',
        'v_wind',
        'pres',
        'tdry',
        'dp',
    )
    # narrowing in on lower atmosphere to highlight the PBL height estimates
    display.set_yrng([1000, 500])
    display.set_xrng([-30, 20])

    display.axes[0].axhline(
        pbl_heffter, linestyle="--", color="tab:orange", linewidth=1, label="Heffter PBL Height"
    )
    display.axes[0].axhline(
        pbl_liu, linestyle="--", color="tab:blue", linewidth=1, label="Liu & Liang PBL Height"
    )
    display.axes[0].axhline(
        pbl_rich, linestyle="--", color="tab:green", linewidth=1, label="Bulk Richardson PBL Height"
    )
    display.axes[0].legend(loc="upper right", fontsize=8)


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

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


.. _sphx_glr_download_auto_examples_retrievals_plot_sonde_pbl.py:

.. only:: html

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

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

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

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

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

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: plot_sonde_pbl.zip <plot_sonde_pbl.zip>`


.. only:: html

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

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