PyPLUTO: a data analysis Python package for the PLUTO code
arXiv:2501.09748 · doi:10.21105/joss.08448
Abstract
In recent years, numerical simulations have become indispensable for addressing complex astrophysical problems. The MagnetoHydroDynamics (MHD) framework represents a key tool for investigating the dynamical evolution of astrophysical plasmas, which are described as a set of partial differential equations that enforce the conservation of mass, momentum, and energy, along with Maxwell's equations for the evolution of the electromagnetic fields. Due to the high nonlinearity of the MHD equations (regardless of their specifications, e.g., classical/relativistic or ideal/resistive), a general analytical solution is precluded, making the numerical approach crucial. Numerical simulations usually end up producing large sets of data files, and their scientific analysis leans on dedicated software designed for data visualization. However, in order to encompass all of the code output features, specialized tools focusing on the numerical code may represent a more versatile and built-in tool. Here, we present PyPLUTO, a Python package tailored for efficient loading, manipulation, and visualization of outputs produced with the PLUTO code (Mignone et al., 2007; Mignone et al., 2012). PyPLUTO uses memory mapping to optimize data loading and provides general routines for data manipulation and visualization. PyPLUTO also supports the particle modules of the PLUTO code, enabling users to load and visualize particles, such as cosmic rays (Mignone et al., 2018), Lagrangian (Vaidya et al., 2018), or dust (Mignone et al., 2019) particles, from hybrid simulations. A dedicated Graphical User Interface simplifies the generation of single-subplot figures, making PyPLUTO a powerful yet user-friendly toolkit for astrophysical data analysis.
7 pages, 3 figures. Accepted for publication in JOSS
References in corpus (14)
- Array Programming with NumPy
- PLUTO: a Numerical Code for Computational Astrophysics
- Resistive relativistic MHD simulations of astrophysical jets
- A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
- A -order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
- Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating
- Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions
- The polarization of the synchrotron radiation from a recollimated jet: application to high-energy BL Lacs
- Three-temperature radiation hydrodynamics with PLUTO: Thermal and kinematic signatures of accreting protoplanets
- Relativistic reconnection with effective resistivity: I. Dynamics and reconnection rate
- Magnetic dissipation in short gamma-ray burst jets. I. Resistive relativistic MHD evolution in a model environment
- A comparison of approximate non-linear Riemann solvers for Relativistic MHD
- Stellar-wind feedback and magnetic fields around young compact star clusters: 3D MHD simulations
- Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks