  This is an addendum to the paper by Cappellari (2008, MNRAS, 390, 71), which
presented a simple and efficient method to model the stellar kinematics of
axisymmetric stellar systems. The technique reproduces well the integral-field
kinematics of real galaxies. It allows for orbital anisotropy (three-integral
distribution function), multiple kinematic components, supermassive black holes
and dark matter. The paper described the derivation of the projected second
moments and we provided a reference software implementation. However only the
line-of-sight component was given in the paper. For completeness we provide
here all the six projected second moments, including radial velocities and
proper motions. We present a test against realistic N-body galaxy simulations.
