athenapk
athenapk copied to clipboard
AthenaPK: a performance portable version of Athena++ built on Parthenon and Kokkos
This adds the precipitator pgen. (New PR to avoid issues with actions running on forks.) **This PR depends on:** * https://github.com/parthenon-hpc-lab/athenapk/pull/79 * https://github.com/parthenon-hpc-lab/athenapk/pull/80 * TODO: well-balanced reconstruction (for hydrostatic equilibrium)...
Python analysis/plotting example: ```python import matplotlib.pyplot as plt import numpy as np import sys sys.path.insert( 1, "/home/pgrete/src/athenapk/external/parthenon" + "/scripts/python/packages/parthenon_tools/parthenon_tools", ) try: import phdf except ModuleNotFoundError: print("Couldn't find module to read...
This adds a Markdown file with documentation on the convention for CIE cooling curves that is expected by AthenaPK, as well as documenting the approximations used in calculating the mean...
It would be useful to have tracer particles, following the [Parthenon example](https://github.com/parthenon-hpc-lab/parthenon/blob/88ece809dc56ae16d695db8d9051433bff01feb1/example/particle_tracers/particle_tracers.cpp>). We could even use the marker particle advection algorithm using the time-centered, face-centered velocities: https://ui.adsabs.harvard.edu/abs/1965PhFl....8.2182H/abstract
The cooling tables in the repository are inconsistent with the convention expected by AthenaPK (https://github.com/parthenon-hpc-lab/athenapk/pull/97). https://github.com/parthenon-hpc-lab/athenapk/tree/main/inputs/cooling_tables For temperatures about $10^{5.5}$ K where CIE is valid, this causes a 10-20% discrepancy...
New version of https://github.com/parthenon-hpc-lab/athenapk/pull/57.
New version of https://github.com/parthenon-hpc-lab/athenapk/pull/58. (Renamed the branch and didn't realize the PR would be closed automatically.)