-
Notifications
You must be signed in to change notification settings - Fork 20
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Added species specific energy grids #113
base: master
Are you sure you want to change the base?
Conversation
Changed the EKEV, EBND, WE, and DE variables to be 2 dimensional with the first dimension being the species. Now each species can have their own energy grid range. This is in preperation for being able to simulate relativistic electrons
Started this pull request to keep track of the ongoing work. Initialization and restart files still need work, but otherwise this feature works. NOTE: Because we are using arbitrary energy limits we needed a new way to create the energy bins. As such, this breaks all tests. We can either update the tests, or we can code it such that the basic energy grids are the same as they were before and only new species have different energy grids. |
Update to the restart, initialization, and satellite routines to correctly handle the species dependent energy grids
Added an initialization file that includes the radiation belt electrons (1 MeV to 10 MeV) based on the work done by Xingzhi.
Added interpolation on intial read for the energy and pitch angles, but currently pitch angle interpolation does not work as we do not store the pitch angle bins in the initialization or restart files
Updated the initialization file to include a larger range of electron energies for both the default electron species and the rad belt electron species
Removed negatives from fluxes in the initialization file, also fixed an issue with the satellite outputs not having the correct energy grids
…energy_by_species
There have been several updates since I started this branch, these changes bring us to a working point with the current upstream changes
Changed the EKEV, EBND, WE, and DE variables to be 2 dimensional with the first dimension being the species. Now each species can have their own energy grid range. This is in preperation for being able to simulate relativistic electrons