psipy is a package for loading and visualising the output of PSI’s MAS model runs. This page provides some narrative documentation that should get you up and running with obtaining, loading, and visualising some model results.
The PSI MHDWeb pages give access to MAS model runs. The runs are indexed by Carrington rotation, and for each Carrington rotation there are generally a number of different runs, varying in the type model run and/or the boundary conditions.
To load data with psipy you need to manually download the files you are interested in to a directory on your computer.
psipy stores the output variables from a single MAS run in the
object. To create one of these, specify the directory which has all of the
.hdf files you want to load:
from psipy.model import MASOutput directory = '/path/to/files' mas_output = MASOutput(directory)
It is assmumed that the files have the filename structure
var is a variable name, and
is a three-digit zero-padded integer timestep.
To see which variables have been loaded, we can look at the
This will print a list of the variables that have been loaded. Each individual variable can then be accessed with square brackets, for example to get the radial magnetic field component:
br = mas_output['br']
The data stored in
Variable.data contains the values of the data as a normal
array, and in addition stores the coordinates of each data point.
MAS model outputs are defined on a 3D grid of points on a spherical grid. The
coordinate names are
'r', 'theta', 'phi'. The coordinate values along each
dimension can be accessed using the
r_coords, theta_coords, phi_coords
rvals = br.r_coords
Variable objects have a
Variable.sample_at_coords method, to take a sample of
the 3D data cube along a 1D trajectory. This is helpful for flying a ‘virtual
spacecraft’ through the model, in order to compare model results with in-situ
sample_at_coords requires arrays of longitude, latitude, and radial distance.
Given these coordinates, it uses linear interpolation to extract the values
of the variable at each of the coordinate points.
For an example of how all this works, see Comparing in-situ data to model output.