Steve,
I was able to get the required output now. But, once you remove the force using unfix, the potential energy contribution due to this force will decrease gradually but not suddenly. I want to know if there is any way we can track the gradual decrease of PE contribution of the force. I hope you understand what i mean…
Thanks
Ajit
If I add these commands to bench/in.lj
fix 5 all addforce 0.0 0.0 1.0
fix_modify 5 energy yes
thermo_style custom temp pe ke etotal press f_5
then I get output. Are you using the most current version
of LAMMPS?
Steve