The magnetic skeleton finder is functionality that is designed to find and analyze magnetic skeleton structures using a trilinear method as described in Haynes et al. (2010) In this notebook we will briefly cover how the magnetic skeleton finder utilizes the method described in the paper in order to describe the structures of magnetic field.
import numpy as np
from plasmapy.analysis.magnetic_skeleton import magnetic_skeleton_find
--------------------------------------------------------------------------- ModuleNotFoundError Traceback (most recent call last) Input In [6], in <cell line: 2>() 1 import numpy as np ----> 2 from plasmapy.analysis.magnetic_skeleton import ( 3 magnetic_skeleton_find 4 ) ModuleNotFoundError: No module named 'plasmapy'
def example():
x_arr = np.linspace(-5, 5, 60)
y_arr = np.linspace(-5, 5, 60)
z_arr = np.linspace(-5, 5, 60)
x, y, z = np.meshgrid(
x_arr,
y_arr,
z_arr,
indexing="ij",
)
u, v, w = vspace_func_3(x, y, z)
nullpoints, spines, fan, separators = magnetic_skeleton_find(
x_arr, y_arr, z_arr, u, v, w
)
example()