Output Number of Bonds at any Timestep

Hello,

I am using the BPM user package through LAMMPS and was interested in bond information, specifically how many bonds are present within my system at any given timestep. Is there a way to output this with a compute command?

My setup consists of a disk with 3759 particles. A velocity is applied and the particles are constantly shifting and breaking/forming bonds at each timestep.

Please see the examples in the examples/bpm/impact folder and then the LAMMPS manual for the documentation of the individual commands.