Note
Click here to download the full example code
Reading and visualising PLUTO runs
First, load the required modules.
from pathlib import Path
import matplotlib.pyplot as plt
from psipy.model import PLUTOOutput
Load a set of PLUTO output files.
pluto_path = 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()
Total running time of the script: ( 0 minutes 0.000 seconds)