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

.. only:: html

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

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

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

.. _sphx_glr_examples_12_compare_theis_quasi_steady.py:


Quasi steady convergence
========================

The quasi steady is reached, when the radial shape of the drawdown in not
changing anymore.

.. GENERATED FROM PYTHON SOURCE LINES 8-42



.. image-sg:: /examples/images/sphx_glr_12_compare_theis_quasi_steady_001.png
   :alt: 12 compare theis quasi steady
   :srcset: /examples/images/sphx_glr_12_compare_theis_quasi_steady_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    import numpy as np
    from matplotlib import pyplot as plt

    from anaflow import theis, thiem

    time = [10, 100, 1000]
    rad = np.geomspace(0.1, 10)
    r_ref = 10.0

    head_ref = theis(
        time,
        np.full_like(rad, r_ref),
        storage=1e-3,
        transmissivity=1e-4,
        rate=-1e-4,
    )
    head1 = theis(time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4) - head_ref
    head2 = theis(time, rad, storage=1e-3, transmissivity=1e-4, rate=-1e-4, r_bound=r_ref)
    head3 = thiem(rad, r_ref, transmissivity=1e-4, rate=-1e-4)

    for i, step in enumerate(time):
        label_1 = "Theis quasi steady" if i == 0 else None
        label_2 = "Theis bounded" if i == 0 else None
        plt.plot(rad, head1[i], label=label_1, color="C" + str(i), linestyle="--")
        plt.plot(rad, head2[i], label=label_2, color="C" + str(i))

    plt.plot(rad, head3, label="Thiem", color="k", linestyle=":")

    plt.xlabel("r in [m]")
    plt.ylabel("h in [m]")
    plt.legend()
    plt.tight_layout()
    plt.show()


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

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


.. _sphx_glr_download_examples_12_compare_theis_quasi_steady.py:

.. only:: html

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

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

      :download:`Download Jupyter notebook: 12_compare_theis_quasi_steady.ipynb <12_compare_theis_quasi_steady.ipynb>`

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

      :download:`Download Python source code: 12_compare_theis_quasi_steady.py <12_compare_theis_quasi_steady.py>`

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

      :download:`Download zipped: 12_compare_theis_quasi_steady.zip <12_compare_theis_quasi_steady.zip>`


.. only:: html

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

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