.. 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-19

.. code-block:: Python

    import empymod
    import numpy as np
    import matplotlib.pyplot as plt








.. GENERATED FROM PYTHON SOURCE LINES 20-22

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

.. GENERATED FROM PYTHON SOURCE LINES 22-47

.. code-block:: Python


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

    # Target model
    inp_tg = {
        '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
    inp_hs = inp_tg.copy()
    inp_hs['depth'] = inp_tg['depth'][0]
    inp_hs['res'] = inp_tg['res'][:2]

    # Compute responses
    sths = empymod.dipole(**inp_hs, signal=1)               # Step, Halfspace
    sttg = empymod.dipole(**inp_tg, signal=1)               # Step, Target
    imhs = empymod.dipole(**inp_hs, signal=0, ft='fftlog')  # Impulse, Halfspace
    imtg = empymod.dipole(**inp_tg, signal=0, ft='fftlog')  # Impulse, Target








.. GENERATED FROM PYTHON SOURCE LINES 48-50

Plot
----

.. GENERATED FROM PYTHON SOURCE LINES 50-76

.. code-block:: Python


    fig, (ax1, ax2) = plt.subplots(
            1, 2, figsize=(8, 4), sharex=True, constrained_layout=True)

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

    # Impulse response
    ax2.set_title('(b)')
    ax2.plot(np.r_[0, 0, t], np.r_[2, 0, imhs/imtg.max()],
             'k', label='Uniform half-space')
    ax2.plot(np.r_[0, t], np.r_[0, imtg/imtg.max()],
             'r', label='Hydrocarbon reservoir')
    ax2.axis([-.02, 0.06, 0, 1.02])
    ax2.set_xlabel('Time (s)')
    ax2.set_ylabel('Normalized derivative (Ω/m²/s)')
    ax2.legend()




.. 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 77-84

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

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

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


.. GENERATED FROM PYTHON SOURCE LINES 86-87

.. 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'>Sat Jan 17 19:16:46 2026 UTC</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>OS</td>
        <td style='text-align: left; border: 1px solid;'>Linux (Ubuntu 22.04)</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.6 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.12 (main, May  6 2025, 10:45:53) [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;'>2.3.5</td>
        <td style='text-align: right; border: 1px solid;'>scipy</td>
        <td style='text-align: left; border: 1px solid;'>1.17.0</td>
        <td style='text-align: right; border: 1px solid;'>numba</td>
        <td style='text-align: left; border: 1px solid;'>0.63.1</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>empymod</td>
        <td style='text-align: left; border: 1px solid;'>2.5.4.dev4+g4866dc83a</td>
        <td style='text-align: right; border: 1px solid;'>libdlf</td>
        <td style='text-align: left; border: 1px solid;'>0.3.0</td>
        <td style='text-align: right; border: 1px solid;'>IPython</td>
        <td style='text-align: left; border: 1px solid;'>9.9.0</td>
      </tr>
      <tr>
        <td style='text-align: right; border: 1px solid;'>matplotlib</td>
        <td style='text-align: left; border: 1px solid;'>3.10.8</td>
        <td style= border: 1px solid;'></td>
        <td style= border: 1px solid;'></td>
        <td style= border: 1px solid;'></td>
        <td style= border: 1px solid;'></td>
      </tr>
    </table>
    </div>
    <br />
    <br />


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

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

**Estimated memory usage:**  193 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>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: ziolkowski2007.zip <ziolkowski2007.zip>`


.. only:: html

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

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