The GALPROP program (http://galprop.stanford.edu) computes cosmic-ray propagation numerically. To obtain manifestly accurate solutions an explicit method is now available to replace the Crank-Nicolson operator-splitting method. This can be implemented in the public version of GALPROP with the update provided here.
Instructions: replace the four files propel.cc, propel_diagnostics, Galdef.cc and Galdef.h and rebuild the executable.
Add solution_method and solution_convergence to the galdef...