If the user specifies starting_chemistry but runs in disk mode (return_arrays=False), the initial conditions aren't loaded correctly. Users should use abundLoadFile as a user parameter instead. See plot below:
Two fixes are possible:
- Use the provided abundances for starting_chemistry.
- Throw an error when starting_chemistry is specified but not return_array / return_makerates.