Thanks Axel,
It does work now, for the temperature control.
I have other question for output of the X, Y and Vx ,Vy of COM of every rigid body. The system include such bodies about one thousand.
Is there any simple method to handle it, except writing f_[][] from 1 to 1000 in the line of thermo_style custom.
Thanks again.
------------------ Original ------------------