Reading and visualising PLUTO runs

First, load the required modules.

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

Load a set of PLUTO output files.

pluto_path = '/Users/dstansby/github/psipy/data/pluto'
model = PLUTOOutput(pluto_path)

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

print(model.variables)

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['Bx1'].contour_phi_cut(75, levels=[0], ax=ax, colors='white')
plt.show()