How to extract force when the center of mass of molecules is put in the box in GROMACS? To extract the force as in the title I used the gmx trjconv command like below.

gmx trjconv -f alad.trr -s alad.tpr -n alad.ndx -pbc mol -novel -force -o alad.gro


However, the force was not extracted in alad.gro file.
Of course, I got the coordinates and the force at equal intervals when getting the alad.trr file.

Another way, I used gmx traj command to get force.

gmx traj -f alad.trr -s alad.tpr -n alad.ndx -pbc -mol -ox coordi.xvg -of force.xvg


However, the force obtained from force.xvg didn't look like the value I wanted.

So, how do I get the force of molecules put in a single periodic box?
I solved the problem. As gmx trjconv -pbc mol option only make molecules translate, the cooridinates obtained from alad.gro can correspond to forces obtained from gmx traj -f alad.trr -s alad.tpr -n alad.ndx -of force-dump.xvg.