.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "gallery/reproducing/hunziker2015.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_gallery_reproducing_hunziker2015.py: Hunziker et al., 2015, Geophysics ================================= Reproducing Figure 3 of the manual from `EMmod`. This example does, as such, not actually reproduce a figure of Hunziker et al., 2015, but of the manual that comes with the software accompanying the paper. With the software comes an example input file named `simplemod.scr`, and the corresponding result is shown in the manual of the code in Figure 3. If you are interested in reproducing the figures of the actual paper have a look at the notebooks in the repo `article-geo2017 `_. **Reference** - **Hunziker, J., J. Thorbecke, and E. Slob, 2015**, The electromagnetic response in a layered vertical transverse isotropic medium: A new look at an old problem: Geophysics, 80(1), F1–F18; DOI: `10.1190/geo2013-0411.1 `_; Software: `software.seg.org/2015/0001 `_. .. GENERATED FROM PYTHON SOURCE LINES 24-29 .. code-block:: Python import empymod import numpy as np import matplotlib.pyplot as plt .. GENERATED FROM PYTHON SOURCE LINES 30-34 Compute the data ---------------- Compute the electric field with the parameters defined in `simplemod.scr`. .. GENERATED FROM PYTHON SOURCE LINES 34-61 .. code-block:: Python # x- and y-offsets x = np.arange(4000)*7-1999.5*7 y = np.arange(1500)*10-749.5*10 # Create 2D arrays of them rx = np.repeat([x, ], np.size(y), axis=0) ry = np.repeat([y, ], np.size(x), axis=0) ry = ry.transpose() # Compute the electric field efield = empymod.dipole( src=[0, 0, 150], rec=[rx.ravel(), ry.ravel(), 200], depth=[0, 200, 1000, 1200], res=[2e14, 1/3, 1, 50, 1], aniso=[1, 1, np.sqrt(10), 1, 1], freqtime=0.5, epermH=[1, 80, 17, 2.1, 17], epermV=[1, 80, 17, 2.1, 17], mpermH=[1, 1, 1, 1, 1], mpermV=[1, 1, 1, 1, 1], ab=11, htarg={'pts_per_dec': -1}, ).reshape(np.shape(rx)) .. rst-class:: sphx-glr-script-out .. code-block:: none :: empymod END; runtime = 0:00:02.728431 :: 1 kernel call(s) .. GENERATED FROM PYTHON SOURCE LINES 62-64 Plot ---- .. GENERATED FROM PYTHON SOURCE LINES 64-92 .. code-block:: Python # Create a similar colormap as Hunziker et al., 2015. cmap = plt.cm.get_cmap("jet", 61) plt.figure(figsize=(9, 8)) # 1. Amplitude plt.subplot(211) plt.title('Amplitude (V/m)') plt.xlabel('Offset (km)') plt.ylabel('Offset (km)') plt.pcolormesh(x/1e3, y/1e3, np.log10(efield.amp()), cmap=cmap, vmin=-16, vmax=-7, shading='nearest') plt.colorbar() # 2. Phase plt.subplot(212) plt.title('Phase (°)') plt.xlabel('Offset (km)') plt.ylabel('Offset (km)') plt.pcolormesh(x/1e3, y/1e3, efield.pha(deg=False, unwrap=False, lag=True), cmap=cmap, vmin=-np.pi, vmax=np.pi, shading='nearest') plt.colorbar() plt.tight_layout() plt.show() .. image-sg:: /gallery/reproducing/images/sphx_glr_hunziker2015_001.png :alt: Amplitude (V/m), Phase (°) :srcset: /gallery/reproducing/images/sphx_glr_hunziker2015_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/empymod/checkouts/stable/examples/reproducing/hunziker2015.py:66: MatplotlibDeprecationWarning: The get_cmap function was deprecated in Matplotlib 3.7 and will be removed two minor releases later. Use ``matplotlib.colormaps[name]`` or ``matplotlib.colormaps.get_cmap(obj)`` instead. cmap = plt.cm.get_cmap("jet", 61) .. GENERATED FROM PYTHON SOURCE LINES 93-99 Original Figure --------------- Figure 3 of the manual of `EMmod`. .. image:: ../../_static/figures/Hunziker2015.png .. GENERATED FROM PYTHON SOURCE LINES 101-102 .. code-block:: Python empymod.Report() .. raw:: html
Fri Mar 01 20:18:37 2024 UTC
OS Linux CPU(s) 2 Machine x86_64
Architecture 64bit RAM 7.5 GiB Environment Python
File system ext4
Python 3.11.6 (main, Feb 1 2024, 16:47:41) [GCC 11.4.0]
numpy 1.26.4 scipy 1.12.0 numba 0.59.0
empymod 2.3.0 IPython 8.22.1 matplotlib 3.8.3


.. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 10.058 seconds) **Estimated memory usage:** 834 MB .. _sphx_glr_download_gallery_reproducing_hunziker2015.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: hunziker2015.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: hunziker2015.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_