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

.. only:: html

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

        :ref:`Go to the end <sphx_glr_download_gallery_reproducing_ziolkowski2007.py>`
        to download the full example code

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

.. _sphx_glr_gallery_reproducing_ziolkowski2007.py:


Ziolkowski et al., 2007
=======================

Reproducing Figure 3 of Ziolkowski et al., 2007, Geophysics. This is a land
MTEM example.

**Reference**

- **Ziolkowski, A., B. Hobbs, and D. Wright,  2007**, Multitransient
  electromagnetic demonstration survey in France: Geophysics, 72, F197-F209;
  DOI: `10.1190/1.2735802 <https://doi.org/10.1190/1.2735802>`_.

.. GENERATED FROM PYTHON SOURCE LINES 15-20

.. code-block:: Python

    import empymod
    import numpy as np
    from copy import deepcopy as dc
    import matplotlib.pyplot as plt








.. GENERATED FROM PYTHON SOURCE LINES 21-23

Computation
-----------

.. GENERATED FROM PYTHON SOURCE LINES 23-46

.. code-block:: Python


    # Time
    t = np.linspace(0.001, 0.06, 101)

    # Target model
    inp2 = {'src': [0, 0, 0.001],
            'rec': [1000, 0, 0.001],
            'depth': [0, 500, 525],
            'res': [2e14, 20, 500, 20],
            'freqtime': t,
            'verb': 1}

    # HS model
    inp1 = dc(inp2)
    inp1['depth'] = inp2['depth'][0]
    inp1['res'] = inp2['res'][:2]

    # Compute responses
    sths = empymod.dipole(**inp1, signal=1)  # Step, HS
    sttg = empymod.dipole(**inp2, signal=1)  # " "   Target
    imhs = empymod.dipole(**inp1, signal=0, ft='fftlog')  # Impulse, HS
    imtg = empymod.dipole(**inp2, signal=0, ft='fftlog')  # " "      Target








.. GENERATED FROM PYTHON SOURCE LINES 47-49

Plot
----

.. GENERATED FROM PYTHON SOURCE LINES 49-81

.. code-block:: Python


    plt.figure(figsize=(9, 4))
    plt.subplots_adjust(wspace=.3)

    # Step response
    plt.subplot(121)
    plt.title('(a)')
    plt.plot(np.r_[0, 0, t], np.r_[0, sths[0], sths], 'k',
             label='Uniform half-space')
    plt.plot(np.r_[0, 0, t], np.r_[0, sttg[0], sttg], 'r',
             label='Hydrocarbon reservoir')
    plt.axis([-.02, 0.06, 0, 8e-9])
    plt.xlabel('Time (s)')
    plt.ylabel('Electric field amplitude (V/m/A-m)')
    plt.legend()

    # Impulse response
    plt.subplot(122)
    plt.title('(b)')

    # Normalize by max-response
    ntg = np.max(np.r_[imtg, imhs])

    plt.plot(np.r_[0, 0, t], np.r_[2, 0, imhs/ntg], 'k',
             label='Uniform half-space')
    plt.plot(np.r_[0, t], np.r_[0, imtg/ntg], 'r', label='Hydrocarbon reservoir')
    plt.axis([-.02, 0.06, 0, 1.02])
    plt.xlabel('Time (s)')
    plt.ylabel(r'Normalized derivative ($\Omega$/m$^2$/s)')
    plt.legend()
    plt.show()




.. image-sg:: /gallery/reproducing/images/sphx_glr_ziolkowski2007_001.png
   :alt: (a), (b)
   :srcset: /gallery/reproducing/images/sphx_glr_ziolkowski2007_001.png
   :class: sphx-glr-single-img





.. GENERATED FROM PYTHON SOURCE LINES 82-89

Original Figure
---------------

Figure 3 of Ziolkowski et al., 2007, Geophysics:

.. image:: ../../_static/figures/Ziolkowski2007.jpg


.. GENERATED FROM PYTHON SOURCE LINES 91-93

.. code-block:: Python


    empymod.Report()





.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <table style='border: 1.5px solid;'>
      <tr>
         <td style='text-align: center; font-weight: bold; font-size: 1.2em; border: 1px solid;' colspan='6'>Fri Mar 01 20:18:46 2024 UTC</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>OS</td>
        <td style='text-align: left; border: 1px solid;'>Linux</td>
        <td style='text-align: right; border: 1px solid;'>CPU(s)</td>
        <td style='text-align: left; border: 1px solid;'>2</td>
        <td style='text-align: right; border: 1px solid;'>Machine</td>
        <td style='text-align: left; border: 1px solid;'>x86_64</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>Architecture</td>
        <td style='text-align: left; border: 1px solid;'>64bit</td>
        <td style='text-align: right; border: 1px solid;'>RAM</td>
        <td style='text-align: left; border: 1px solid;'>7.5 GiB</td>
        <td style='text-align: right; border: 1px solid;'>Environment</td>
        <td style='text-align: left; border: 1px solid;'>Python</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>File system</td>
        <td style='text-align: left; border: 1px solid;'>ext4</td>
      </tr>
      <tr>
         <td style='text-align: center; border: 1px solid;' colspan='6'>Python 3.11.6 (main, Feb  1 2024, 16:47:41) [GCC 11.4.0]</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>numpy</td>
        <td style='text-align: left; border: 1px solid;'>1.26.4</td>
        <td style='text-align: right; border: 1px solid;'>scipy</td>
        <td style='text-align: left; border: 1px solid;'>1.12.0</td>
        <td style='text-align: right; border: 1px solid;'>numba</td>
        <td style='text-align: left; border: 1px solid;'>0.59.0</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>empymod</td>
        <td style='text-align: left; border: 1px solid;'>2.3.0</td>
        <td style='text-align: right; border: 1px solid;'>IPython</td>
        <td style='text-align: left; border: 1px solid;'>8.22.1</td>
        <td style='text-align: right; border: 1px solid;'>matplotlib</td>
        <td style='text-align: left; border: 1px solid;'>3.8.3</td>
      </tr>
    </table>
    </div>
    <br />
    <br />


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

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

**Estimated memory usage:**  10 MB


.. _sphx_glr_download_gallery_reproducing_ziolkowski2007.py:

.. only:: html

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

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

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

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

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


.. only:: html

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

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