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

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

.. _sphx_glr_examples_interp_plot_interpolated_surface.py:


Interpolated with surface topo
==============================

Demonstrate interpolating two-dimensional model output when surface topography,
rather than bedrock topography, is present in the original data. The result is
the same as in the interpolated output example.

.. GENERATED FROM PYTHON SOURCE LINES 13-46



.. image-sg:: /examples/interp/images/sphx_glr_plot_interpolated_surface_001.png
   :alt: Interpolated output
   :srcset: /examples/interp/images/sphx_glr_plot_interpolated_surface_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    downloading https://raw.githubusercontent.com/juseg/hyoga-data/main/pism/pism.alps.vis.refined.nc...






|

.. code-block:: default


    import matplotlib.pyplot as plt
    import hyoga

    # initialize figure
    ax = plt.subplot()

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

        # compute surface altitude and remove bedrock altitude
        ds['usurf'] = ds.hyoga.getvar('surface_altitude')
        ds = ds.drop('topg')

        # compute isostatic adjustment from a reference input topography
        ds = ds.hyoga.assign_isostasy(hyoga.open.example('pism.alps.in.boot.nc'))

        # perform the actual interpolation
        ds = ds.hyoga.interp(hyoga.open.example('pism.alps.vis.refined.nc'))

        # plot model output
        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')

        # add scale bar
        ds.hyoga.plot.scale_bar()

    # set title
    ax.set_title('Interpolated output')

    # show
    plt.show()


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

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


.. _sphx_glr_download_examples_interp_plot_interpolated_surface.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_surface.py <plot_interpolated_surface.py>`

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

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


.. only:: html

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

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