New techniques for calculating heat and particle source rates due to neutral beam injection in axisymmetric tokamaks