sam tygier
-
2011-10-24
Menu
▾
▴
#23 SORT not recorded by CHAMBR
The manual claims that when a particle is stopped by by CHAMBR its IEX is set to -4 and SORT to its path length.
In chmbre.f, IEX is set by a call to KSTOP, but SORT is not set. I attach a patch that sets the SORT of the particle to SAR which seems to hold the path length coordinate.
This is not perfect as I think SAR is the path length at the end of the integration step rather than the point that it crosses the boundary. Also I imagine in actual use the location of the loss is of more interest than the path length.
