Good afternoon,
I found a bug in the QuasiLocalMeasures thorn, in the computation of the ADM momentum. In the source code file qlm_analyse.F90, in line 361 there is a plus sign instead of a minus sign, contradicting the ADM momentum formula: P^i_{ADM} = \frac{1}{8\pi}\lim_{r\to\infty} \oint_S (K^i_j - \delta^i_j K)dS^j [1]
As expected for slices where K \neq 0 one was getting a wrong result. The correction is easy, just replace the plus sign by a minus sign. In any case you can find a patch file in attachment. Should I in any case file a ticket or a Pull Request somewhere regarding this?
Best regards,
Marco Brito.
[1] M. Alcubierre, Introduction to 3+1 Numerical Relativity, International Series of Monographs on Physics (Oxford University Press, Oxford, 2008).