interpolate_moho_depths
=======================

.. py:module:: interpolate_moho_depths

.. autoapi-nested-parse::

   Script to create linear interpolation of the moho depths.
   And then get values with the surface.
   Meta aftes tis times mporoun na dior8w8oun mia mia...

   Location: Chavannes-pres-renens, CH
   Date: Feb 2023
   Author: Konstantinos Michailos



Attributes
----------

.. autoapisummary::

   interpolate_moho_depths.lon
   interpolate_moho_depths.lat
   interpolate_moho_depths.dep
   interpolate_moho_depths.ln
   interpolate_moho_depths.profile_A
   interpolate_moho_depths.f
   interpolate_moho_depths.znew


Functions
---------

.. autoapisummary::

   interpolate_moho_depths.interp_surf
   interpolate_moho_depths.create_profiles


Module Contents
---------------

.. py:data:: lon
   :value: []


.. py:data:: lat
   :value: []


.. py:data:: dep
   :value: []


.. py:data:: ln

.. py:function:: interp_surf(x_, y_, z_, plot=False, **kwargs)

   Gives the 2-D interpolated surface of exhumation rates.

   Use scipy interpolate to create a 2-D surface of a grid of
   observations.

   :type x_: list
   :param x_: list of coordinates (longitude)
   :type y_: list
   :param y_: list of coordinates (latitutes)
   :type z_: list
   :param z_: list of values to be interpolated
   :param kwargs: Any other arguments accepted by scipy.interpolate

   :returns: interpolated surface
   :rtype: function

   .. rubric:: Example        fig = plt.figure(figsize=(10, 6))
       ax = axes3d.Axes3D(fig)
       ax.plot_wireframe(x, y, z)
       ax.plot_surface(x, y, z, cmap=cm.viridis, alpha=0.2)
       ax.scatter3D(x, y, z, c='r')
       ax.set_zlim(.0, 12.0)
       ax.zaxis.set_major_locator(LinearLocator(10))
       ax.zaxis.set_major_formatter(FormatStrFormatter('%.02f'))
       plt.show()

   >>> f = interp_surf(x,y,z)
   >>> znew = f(xnew, ynew)


.. py:function:: create_profiles(profile_points)

.. py:data:: profile_A

.. py:data:: f

.. py:data:: znew