New paper (Turolla, Nobili 2013) will help with the calculations for the
code.