Dataset of ``Plasma Distribution Solver: A Model for Field-Aligned Plasma Profiles Based on Spatial Variation of Velocity Distribution Functions"

This dataset contains the plasma distribution data in the Jupiter–Io system, calculated from the Plasma Distribution Solver and used for figures in the paper “Plasma Distribution Solver: A model for field-aligned plasma profiles based on spatial variation of velocity distribution functions” by K. Saito et al. (2023). The contents of files ‘all_Case_1.csv’ and ‘all_Case_2.csv’ are as follows: Position along the magnetic field line (0 at the magnetic equator) [m] (column 1) Distance from the Jovian center [km] (column 2) Magnetic latitude [rad]([degree]) (column 3(4)) Magnetic flux density [T] (column 5) The initial condition of electrostatic potential [V] (column 6) The result of electrostatic potential [V] (column 7) Number density profiles [m-3] (columns 8-17) Charge density profiles obtained from the integration of velocity distribution functions [C m-3] (column 18) Charge density profiles obtained from Poisson’s equation [C m-3] (column 19) Convergence value (column 20) Particle flux density [m-2 s-1] (columns 21-30) Mean flow velocity parallel to the field line [m s-1] (columns 31-40) Plasma pressure perpendicular to the field line [Pa] (columns 41-50) Plasma pressure parallel to the field line [Pa] (columns 51-60) Plasma dynamic pressure [Pa] (columns 61-70) Perpendicular temperature [J] (columns 71-80) Parallel temperature [J] (columns 81-90) Alfvén speed considering the displacement current term in Ampère’s law [m s-1] (column 91) Alfvén speed per the speed of light (column 92) Ion inertial length using averaged mass [m] (column 93) Electron inertial length [m] (column 94) Ion Larmor radius using averaged mass [m] (column 95) Ion acoustic gyroradius using averaged mass [m] (column 96) Electron Larmor radius [m] (column 97) Current density [A m-2] (column 98) The Python codes ‘plot_all.py,’ ‘plot_plasma_beta_comparison.py,’ and ‘plot_Alfven_speed_comparison.py’ can plot Figures 5, 6, 7, and 9 of the paper using the above CSV files. The files ‘boundary_conditions_Case_1.csv’ and ‘boundary_conditions_Case_2.csv’ contain the boundary conditions for Cases 1 and 2. The zip files ‘probability_density_function_Case_1_H_Io.zip’ and ‘probability_density_function_Case_1_H_Jupiter_North.zip’ are zipped CSV files with the same name. The contents of these files are as follows: Magnetic latitude [degree] (column 1) Perpendicular velocity at the particle position [m s-1] (column 2) Parallel velocity at the particle position [m s-1] (column 3) Perpendicular velocity at the boundary [m s-1] (column 4) Parallel velocity at the boundary [m s-1] (column 5) Probability density function [s3 m-3] (column 6) Differential flux per number density [cm-2 s-1 sr-1 keV-1] (column 7) The Python code ‘plot_velocity_distribution_function.py’ can plot Figure 8 of the paper using this CSV file.