[gmx-developers] zero temperature II

[Michael Brunsteiner]

dear all,
In a previous mail I mentioned a problem
with zero temperature output. I seem to have 
found what the problem is and this looks like
a bug, or at least an undocumented limitation.
When I do md with a system containing only frozen atoms 
and (not-frozen) water both ener.edr and md.log say 
that the kinetic energy and the temperature are
zero. The thermostate seems to work though, because
my system does not "explode".
If I "un-freeze" the frozen atoms the
output is back to normal. The fact that frozen atoms
are formally provided with an infinite mass
might somehow cause this problem.
Any ideas for a work-around ?
Where in the code is E_kin calculated ?
cheers,
Michael


__________________________________
Do you Yahoo!?
Yahoo! Small Business $15K Web Design Giveaway 
http://promotions.yahoo.com/design_giveaway/