[lammps-users] fix_addforce

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

I don't. Once the fix is unfixed, anything it stores is gone. That
a system would gradually relax due to a constraint disappearing
does not seem unusual.

Steve