Reading and visualising MAS runs

First, load the required modules.

from psipy.model import MASOutput
import matplotlib.pyplot as plt

Load a set of MAS output files.

mas_path = '/Users/dstansby/github/psipy/data/helio'
model = MASOutput(mas_path)

Each MAS model contains a number of variables. The variable names can be accessed using the .variables attribute.

print(model.variables)

Out:

['p', 't', 'va', 'br', 'vp', 'jp', 'jr', 'bp', 'vr', 'bt', 'vt', 'rho', 'jt']

Plot a cut of the model at a constant radius

fig, ax = plt.subplots()
model['rho'].plot_radial_cut(0, ax=ax)

Plot a cut of the model at a constant longitude.

fig = plt.figure()
ax = plt.subplot(111, projection='polar')
model['rho'].plot_phi_cut(75, ax=ax)

Add the current sheet

fig = plt.figure()
ax = plt.subplot(111, projection='polar')
model['rho'].plot_phi_cut(75, ax=ax)
model['br'].contour_phi_cut(75, levels=[0], ax=ax, colors='white')
plt.show()