hermes-3
hermes-3 copied to clipboard
inconsistency between upper and lower target
On the lower target sin_alpha is missing:
https://github.com/bendudson/hermes-3/blob/cf5d2967622eb841e5c2245995f872b7fbece308/src/sheath_boundary.cxx#L250
There are also two different expressions for s_i used:
BoutReal s_i = (nesheath > 1e-5) ? nisheath / nesheath : 0.0; // Concentration
BoutReal s_i = clip(nisheath / floor(nesheath, 1e-10), 0, 1); // Concentration
I am currently writing a replacement for the sheath_boundary component, that will not only unify upper and lower target, but also FCI and non-FCI ...
This was fixed in 2faaa3c5