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

.. only:: html

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

        Click :ref:`here <sphx_glr_download_examples_interp_plot_interpolated_resolutions.py>`
        to download the full example code

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

.. _sphx_glr_examples_interp_plot_interpolated_resolutions.py:


Grid refinement
===============

Demonstrate two-dimensional model output interpolation onto grids with
increasing spatial resolution. The leftmost panel shows the original data with
a spatial resolution of 1 km, and other panels show interpolated results.

.. GENERATED FROM PYTHON SOURCE LINES 13-53



.. image-sg:: /examples/interp/images/sphx_glr_plot_interpolated_resolutions_001.png
   :alt: original, 500 m, 200 m, 100 m
   :srcset: /examples/interp/images/sphx_glr_plot_interpolated_resolutions_001.png
   :class: sphx-glr-single-img





.. code-block:: default


    import matplotlib.pyplot as plt
    import hyoga

    # initialize figure
    fig, axes = plt.subplots(ncols=4)
    resolutions = [500, 200, 100]

    # open 100m resolution topography for edits
    with hyoga.open.example('pism.alps.vis.refined.nc') as ds:
        topo = ds.hyoga.getvar('bedrock_altitude')
        topo = topo.sel(x=slice(402e3, 427e3))

    # open demo data
    with hyoga.open.example('pism.alps.out.2d.nc') as ds:
        ds = ds.hyoga.assign_isostasy(hyoga.open.example('pism.alps.in.boot.nc'))

        # plot original data
        ax = axes[0]
        ds.hyoga.plot.bedrock_altitude(ax=ax, cmap='Topographic', center=False)
        ds.hyoga.plot.surface_altitude_contours(ax=ax)
        ds.hyoga.plot.ice_margin(ax=ax, facecolor='w')
        ax.set_title('original')

        # plot interpolated results
        for ax, res in zip(axes[1:], resolutions):
            stride = int(res/100)
            interp = ds.hyoga.interp(topo[::stride, ::stride])
            interp.hyoga.plot.bedrock_altitude(
                ax=ax, cmap='Topographic', center=False)
            interp.hyoga.plot.surface_altitude_contours(ax=ax)
            interp.hyoga.plot.ice_margin(ax=ax, facecolor='w')
            ax.set_title(f'{res} m')

        # constrain original data axes limits
        axes[0].set_xlim(axes[1].get_xlim())
        axes[0].set_ylim(axes[1].get_ylim())

    # show
    plt.show()


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

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


.. _sphx_glr_download_examples_interp_plot_interpolated_resolutions.py:

.. only:: html

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


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

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

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

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


.. only:: html

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

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