Mixed mode surface velocities, kinetic energies and spectral data from compressible 2D simulations of a red giant star

This repository holds the numerical data underlying the paper of De Vries et al. 2026 (in review) regarding the surface velocities and kinetic energies of mixed modes in 2D simulations of a red giant star.MESA_1D_RG.data and MESA_1D_RG.data.FGONG represent the 1D stellar model generated using MESA, in MESA and FGONG formats respectively. These were used as the initial conditions for the simulations, and used to perform the linear theory computation of the mixed modes. GYRE_2D_modes_RG_r0.02-0.90_nonint.tar.gz, GYRE_2D_modes_RG_r0.02-0.98_nonint.tar.gz, GYRE_2D_modes_RG_r0.0-1.0_nonint.tar.gz are tar files for the modes computed using linear theory in both truncated models and the full stellar model. All generated using GYRE. Both the summary file and the detail files for all modes with ell < 5 are included. Note that the wedge setup of the simulations means the values of ell are non-integer, reflected in these linearly computed modes. The simulations are 2D, so the azimuthal order m = 0 for all these modes.mode_outputs_vtot_redgiant.csv and mode_outputs_vtot_redgiant_0.98.csv are the csv files for the RGext90 and RGext98 simulations described in "Revealing mixed modes in compressible hydrodynamical simulations of red giant stars: surface amplitudes and mode kinetic energies, De Vries et al., 2026 ApJL." respectively. These csv files contain 5 columns per mode, in order: the value of the harmonic degree (ell), the frequency of the mode in the simulation (frequency), the frequency of the mode of the full stellar model used for extrapolation (extrpl_frequency), the mode kinetic energy (kin_energy) and the extrapolated surface velocity (surf_velocity).The frequencies are given in \muHz, the kinetic energies in ergs and the surface velocities in cm/s. The kinetic energies are computed solely over the truncated domain of the simulation. The surface velocities are obtained by extrapolating the modes in the simulation to the surface, assuming equal kinetic energy to the simulation in the simulation domain. The csv files are computed using the total velocity squared (v_tot^2 = v_r^2 + v_theta^2). spectrum_dtrec75d1_vr2_0.90_ell_max_50.0.h5, spectrum_dtrec75d1_vth2_0.90_ell_max_50.0.h5 are the spectra up to ell_max = 50 for the simulation RGext90. spectrum_dtrec75d1_vr2_0.98_ell_max_50.0.h5, spectrum_dtrec75d1_vth2_0.98_ell_max_50.0.h5 are the spectra up to ell_max = 50 for the simulation RGext98.The cadence of snapshots used to generate these spectra is dtrec = 750 seconds. Each of these spectra contain 4 datasets. Three of these datasets are 1-dimensional arrays: frequencies (key: frequency), radial grid points (key: radius), and harmonic degree (key: ells). The final dataset contains the spectrum itself (key: spectrum), either the spectrum of P[\hat{v}_r] = 2*|\hat{v}_r|^2 or P[\hat{v}_th] = 2*|\hat{v}_th|^2, indicated by vr2 or vth2 in the name of the spectra. This spectrum is 3-dimensional, with axes corresponding to in order: frequency, radius, harmonic degree.