.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/howto_meshplotlib/plot_observation_points.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_howto_meshplotlib_plot_observation_points.py>`
        to download the full example code.

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

.. _sphx_glr_auto_examples_howto_meshplotlib_plot_observation_points.py:


How to plot data at observation points
======================================

.. sectionauthor:: Florian Zill (Helmholtz Centre for Environmental Research GmbH - UFZ)

In this example we plot the data values on observation points over all
timesteps. Since requested observation points don't necessarily coincide with
actual nodes of the mesh different interpolation options are available. See
:py:mod:`ogstools.meshlib.mesh_series.MeshSeries.probe` for more details.
Here we use a component transport example from the ogs benchmark gallery
(https://www.opengeosys.org/docs/benchmarks/hydro-component/elder/).

.. GENERATED FROM PYTHON SOURCE LINES 16-31

.. code-block:: Python



    import matplotlib.pyplot as plt
    import numpy as np

    from ogstools import examples, meshplotlib
    from ogstools.meshplotlib.utils import justified_labels
    from ogstools.propertylib import Scalar

    meshplotlib.setup.reset()
    mesh_series = examples.load_meshseries_CT_2D_XDMF()
    si = Scalar(
        data_name="Si", data_unit="", output_unit="%", output_name="Saturation"
    )








.. GENERATED FROM PYTHON SOURCE LINES 38-48

To read your own data as a mesh series you can do:

..  code-block:: python

  from ogstools.meshlib import MeshSeries
  mesh_series = MeshSeries("filepath/filename_pvd_or_xdmf")

You can also use a property from the available presets instead of needing to
create your own:
:ref:`sphx_glr_auto_examples_howto_propertylib_plot_propertylib.py`

.. GENERATED FROM PYTHON SOURCE LINES 50-51

Let's define 4 observation points and plot them on the mesh.

.. GENERATED FROM PYTHON SOURCE LINES 53-63

.. code-block:: Python

    points = np.asarray(
        [[x, 0, 60] for x in [0, 40, 80, 120]]
        + [[x, 0, 40] for x in [0, 40, 80, 120]]
    )
    fig = meshplotlib.plot(mesh_series.read(0), si)
    fig.axes[0].scatter(points[:, 0], points[:, 2], s=50, fc="none", ec="r", lw=3)
    for i, point in enumerate(points):
        fig.axes[0].annotate(str(i), (point[0], point[2] - 5), va="top")
    plt.rcdefaults()




.. image-sg:: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_001.png
   :alt: plot observation points
   :srcset: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 64-65

And now probe the points and the values over time:

.. GENERATED FROM PYTHON SOURCE LINES 67-72

.. code-block:: Python

    labels = [f"{i}: {label}" for i, label in enumerate(justified_labels(points))]
    fig = meshplotlib.plot_probe(
        mesh_series=mesh_series, points=points[:4], mesh_property=si,
        time_unit="a", labels=labels[:4]
    )



.. image-sg:: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_002.png
   :alt: plot observation points
   :srcset: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_002.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 73-78

You can also pass create your own matplotlib figure and pass the axes object.
Additionally, you can pass any keyword arguments which are known by
matplotlibs plot function to further customize the curves.
In this case `marker` and `linewidth` are not part of the API of `plot_probe`
but get processed correctly anyway.

.. GENERATED FROM PYTHON SOURCE LINES 80-88

.. code-block:: Python

    fig, axs = plt.subplots(nrows=2, figsize=[10, 5])
    meshplotlib.plot_probe(
        mesh_series, points[:4], si, time_unit="a", ax=axs[0], colors=["k"],
        labels=labels[:4], marker=".")
    meshplotlib.plot_probe(
        mesh_series, points[4:], si, time_unit="a", ax=axs[1], linestyles=["-"],
        labels=labels[4:], linewidth=1,
    )



.. image-sg:: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_003.png
   :alt: plot observation points
   :srcset: /auto_examples/howto_meshplotlib/images/sphx_glr_plot_observation_points_003.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_auto_examples_howto_meshplotlib_plot_observation_points.py:

.. only:: html

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

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

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

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

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