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

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download_api_gallery_climate_plot_09_manabe_wetherald_radiative_convective.py>`
        to download the full example code or to run this example in your browser via JupyterLite.

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

.. _sphx_glr_api_gallery_climate_plot_09_manabe_wetherald_radiative_convective.py:


Manabe and Wetherald's radiative-convective equilibrium (1967)
==============================================================

Syukuro Manabe and Richard Wetherald computed the temperature profile of a
single atmospheric column in which radiation and convection balance. Where
radiative equilibrium alone would be unstable, convection holds the
temperature at a fixed lapse rate. With relative humidity held fixed, doubling
CO2 warmed the surface by about 2.3 K. In this simplest version, the
troposphere passes through the effective temperature at the emission height,
and more CO2 lifts that height. The surface then warms by the lapse rate times
the lift. The model keeps the stratosphere at a fixed skin temperature, so it
misses the stratospheric cooling that Manabe and Wetherald also found.

.. GENERATED FROM PYTHON SOURCE LINES 15-40



.. image-sg:: /api/gallery/climate/images/sphx_glr_plot_09_manabe_wetherald_radiative_convective_001.png
   :alt: Radiative-convective profile, Surface warming 0.98 K
   :srcset: /api/gallery/climate/images/sphx_glr_plot_09_manabe_wetherald_radiative_convective_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    import matplotlib.pyplot as plt
    import numpy as np

    from earthsciencekit.climate import radiative_convective_temperature

    height = np.linspace(0, 25e3, 501)
    effective, lapse = 255.0, 6.5e-3
    profiles = []
    fig, (left, right) = plt.subplots(1, 2, figsize=(10, 4.5), sharey=True)
    for emission, label in [(5.0e3, "present"), (5.15e3, "more CO2: emission level 150 m higher")]:
        profiles.append(radiative_convective_temperature(height, effective, lapse, emission))
        assert np.isclose(radiative_convective_temperature(emission, effective, lapse, emission), effective)
        left.plot(profiles[-1], height / 1e3, label=label)
        left.plot(effective, emission / 1e3, "ko", ms=4)
    warming = profiles[1] - profiles[0]
    assert np.isclose(warming[0], lapse * 150)
    np.testing.assert_allclose(profiles[0][-1], effective / 2**0.25)
    right.plot(warming, height / 1e3)
    left.set(xlabel="temperature [K]", ylabel="height [km]", title="Radiative-convective profile")
    right.set(xlabel="change [K]", title=f"Surface warming {warming[0]:.2f} K")
    left.legend()
    fig.tight_layout()

    plt.show()


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

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


.. _sphx_glr_download_api_gallery_climate_plot_09_manabe_wetherald_radiative_convective.py:

.. only:: html

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

    .. container:: lite-badge

      .. image:: images/jupyterlite_badge_logo.svg
        :target: ../../../lite/lab/index.html?path=api/gallery/climate/plot_09_manabe_wetherald_radiative_convective.ipynb
        :alt: Launch JupyterLite
        :width: 150 px

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

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

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

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

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

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


.. only:: html

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

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