FORCE_MIXING and QM TEMPERATURE

[Pietro Vidossich]

Dear All,

 I have a problem using the FORCE_MIXING in a QMMM run. The problem is that
when the QM subsystem changes its temperature makes jumps. I traced back this
to the fact that
md_ener%nfree_qm
(used to compute the QM TEMPERATURE in SUBROUTINE compute_conserved_quantity)
is not updated when rebuilding the QM/MM calculation.
I would like to fix the problem, but not much clear to me where in the code is
more convenient/appropriate to update md_ener%nfree_qm. Any suggestion?

Many thanks, pietro