Hello!

I’m trying to simulate a packing in where a box filled with spherical particles with gran/hooke potential is isobarically compressed such that density is monotinically increasing, and average contact number/atom increases over time.

I can monitor the change in density and average number of contacts/atom through thermo_style, but I wish to implement an exit condition which monitors the change of these two variables between every N timesteps.

(To clarify: Exit condition should be that if the change of the two scalar values is smaller than a set tolerance, the script should exit.)

Is this possible?