.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/howto_plot/plot_animation.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_howto_plot_plot_animation.py: How to create Animations ======================== .. sectionauthor:: Florian Zill (Helmholtz Centre for Environmental Research GmbH - UFZ) To demonstrate the creation of an animated plot 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 13-21 .. code-block:: Python import matplotlib.pyplot as plt import numpy as np import ogstools as ot from ogstools import examples mesh_series = examples.load_meshseries_CT_2D_XDMF() .. GENERATED FROM PYTHON SOURCE LINES 22-28 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") .. GENERATED FROM PYTHON SOURCE LINES 30-34 You can choose which timesteps to render by passing either an int array corresponding to the indices, or a float array corresponding to the timevalues to render. If a requested timevalue is not part of the timeseries it will be interpolated. In this case every second frame will be interpolated. .. GENERATED FROM PYTHON SOURCE LINES 36-40 .. code-block:: Python timevalues = np.linspace( mesh_series.timevalues[0], mesh_series.timevalues[-1], num=25 ) .. GENERATED FROM PYTHON SOURCE LINES 41-47 Now, let's animate the saturation solution. A timescale at the top indicates existing timesteps and the position of the current timevalue. Note that rendering many frames in conjunction with large meshes might take a really long time. We can pass a ``plot_func`` which can apply custom formatting and / or plotting. To modify the domain, we can use the transform method of MeshSeries. .. GENERATED FROM PYTHON SOURCE LINES 47-59 .. code-block:: Python def mesh_func(mesh: ot.Mesh) -> ot.Mesh: "Clip the left half of the mesh." return mesh.clip("-x", [0, 0, 0]) def plot_func(ax: plt.Axes, timevalue: float) -> None: "Add the time to the title." ax.set_title(f"{timevalue/(365.25*86400):.1f} yrs") .. GENERATED FROM PYTHON SOURCE LINES 60-64 .. code-block:: Python anim = mesh_series.transform(mesh_func).animate( ot.variables.saturation, timevalues, plot_func, vmin=0, vmax=100, dpi=50 ) .. container:: sphx-glr-animation .. raw:: html
.. rst-class:: sphx-glr-script-out .. code-block:: none 0%| | 0/25 [00:00` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_animation.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_animation.zip `