test_zmodel
===========

.. py:module:: test_zmodel

.. autoapi-nested-parse::

   Doublecheck codes read correctly the 3D velocity model of Zhu et al., 2015). + IASP91

   Issue reported from DK and GH.

   Potential issues:
   - Are the Vp and Vs values in the model changed during the scan?
     We can check this by plotting the vp and vs values along the same line and then we can see how it goes.
   - Some depth value might be missing in the model and this "squeezes" the migration,
     because we almost see that the depth scale is compressed.
   - The shape of the ray tracing itself is wrong for something that is probably related to velocity.
     This can result in the discontinuity red spots not being connected, the discontinuity being thick.



Attributes
----------

.. autoapisummary::

   test_zmodel.font
   test_zmodel.fig_size
   test_zmodel.min_lon
   test_zmodel.max_lon
   test_zmodel.min_lat
   test_zmodel.max_lat
   test_zmodel.work_dir
   test_zmodel.path_
   test_zmodel.path
   test_zmodel.path_zmodel_m60
   test_zmodel.parameters
   test_zmodel.longitudes
   test_zmodel.latitudes
   test_zmodel.depths
   test_zmodel.p_velocities
   test_zmodel.s_velocities
   test_zmodel.points_list
   test_zmodel.vp_values
   test_zmodel.vs_values
   test_zmodel.point
   test_zmodel.points
   test_zmodel.values_p
   test_zmodel.values_s
   test_zmodel.liner_interpolation_of_velocities_p
   test_zmodel.liner_interpolation_of_velocities_s
   test_zmodel.longitudes
   test_zmodel.latitudes
   test_zmodel.topos
   test_zmodel.thick_sediments
   test_zmodel.thick_upper
   test_zmodel.thick_lower
   test_zmodel.vp_sediments
   test_zmodel.vp_upper
   test_zmodel.vp_lower
   test_zmodel.vs_sediments
   test_zmodel.vs_upper
   test_zmodel.vs_lower
   test_zmodel.work_dir
   test_zmodel.path_epcrust
   test_zmodel.ln
   test_zmodel.lon
   test_zmodel.lat
   test_zmodel.ele
   test_zmodel.vp_sediments
   test_zmodel.vp_upper
   test_zmodel.vp_lower
   test_zmodel.vs_sediments
   test_zmodel.vs_upper
   test_zmodel.vs_lower
   test_zmodel.thick_sediments
   test_zmodel.thick_upper
   test_zmodel.thick_lower
   test_zmodel.points
   test_zmodel.p_velocities
   test_zmodel.s_velocities
   test_zmodel.z_0_
   test_zmodel.points
   test_zmodel.values_p
   test_zmodel.values_s
   test_zmodel.liner_interpolation_of_velocities_p_ep
   test_zmodel.liner_interpolation_of_velocities_s_ep
   test_zmodel.zmoho
   test_zmodel.minx
   test_zmodel.maxx
   test_zmodel.pasx
   test_zmodel.miny
   test_zmodel.maxy
   test_zmodel.pasy
   test_zmodel.x
   test_zmodel.y
   test_zmodel.z_
   test_zmodel.P_vel_3D_grid
   test_zmodel.S_vel_3D_grid
   test_zmodel.depths
   test_zmodel.pts
   test_zmodel.ax1


Functions
---------

.. autoapisummary::

   test_zmodel.get_iasp91


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

.. py:data:: font

.. py:data:: fig_size

.. py:data:: min_lon
   :value: 0


.. py:data:: max_lon
   :value: 32


.. py:data:: min_lat
   :value: 40


.. py:data:: max_lat
   :value: 55


.. py:data:: work_dir

.. py:data:: path_

.. py:data:: path
   :value: ''


.. py:data:: path_zmodel_m60
   :value: '/data/ZMODEL_M60/'


.. py:data:: parameters

.. py:data:: longitudes

.. py:data:: latitudes

.. py:data:: depths

.. py:data:: p_velocities

.. py:data:: s_velocities

.. py:data:: points_list
   :value: []


.. py:data:: vp_values
   :value: []


.. py:data:: vs_values
   :value: []


.. py:data:: point

.. py:data:: points

.. py:data:: values_p

.. py:data:: values_s

.. py:data:: liner_interpolation_of_velocities_p

.. py:data:: liner_interpolation_of_velocities_s

.. py:data:: longitudes
   :value: []


.. py:data:: latitudes
   :value: []


.. py:data:: topos
   :value: []


.. py:data:: thick_sediments
   :value: []


.. py:data:: thick_upper
   :value: []


.. py:data:: thick_lower
   :value: []


.. py:data:: vp_sediments
   :value: []


.. py:data:: vp_upper
   :value: []


.. py:data:: vp_lower
   :value: []


.. py:data:: vs_sediments
   :value: []


.. py:data:: vs_upper
   :value: []


.. py:data:: vs_lower
   :value: []


.. py:data:: work_dir

.. py:data:: path_epcrust

.. py:data:: ln

.. py:data:: lon

.. py:data:: lat

.. py:data:: ele

.. py:data:: vp_sediments

.. py:data:: vp_upper

.. py:data:: vp_lower

.. py:data:: vs_sediments

.. py:data:: vs_upper

.. py:data:: vs_lower

.. py:data:: thick_sediments

.. py:data:: thick_upper

.. py:data:: thick_lower

.. py:data:: points
   :value: []


.. py:data:: p_velocities
   :value: []


.. py:data:: s_velocities
   :value: []


.. py:data:: z_0_
   :value: -5.0


.. py:data:: points

.. py:data:: values_p

.. py:data:: values_s

.. py:data:: liner_interpolation_of_velocities_p_ep

.. py:data:: liner_interpolation_of_velocities_s_ep

.. py:function:: get_iasp91(x_, y, z, zmoho)

   Retrieves P-wave, S-wave velocities and depths
   from IASPEI91 global velocity model.

   :type x_: numpy.array
   :param x_: Numpy array of x values of the grid points.
   :type y_: numpy.array
   :param y_: Numpy array of y values of the grid points.
   :type z: numpy.array
   :param z: Numpy array of z values of the grid points.
   :type zmoho: int
   :param zmoho: Moho depth in km.

   :rtype: numpy.ndarrays
   :returns: Array of P-wave, S-wave velocities and their depths.


.. py:data:: zmoho
   :value: 35


.. py:data:: minx
   :value: -7.0


.. py:data:: maxx
   :value: 40.0


.. py:data:: pasx
   :value: 0.38


.. py:data:: miny
   :value: 38.0


.. py:data:: maxy
   :value: 60.0


.. py:data:: pasy
   :value: 0.27


.. py:data:: x

.. py:data:: y

.. py:data:: z_

.. py:data:: P_vel_3D_grid

.. py:data:: S_vel_3D_grid

.. py:data:: depths

.. py:data:: pts

.. py:data:: ax1