[gmx-users] thermostate

[hosein geraili]

Hi,
How can I get the temperature of a certain group from .edr or .trr, like the gmx energy, but the problem with gmx energy is, it just gives the temperature for the whole system.
Another question is, do you think it makes sense to use a thermostat with a coupling constant of 100 ps for a very long run (microsecond) on a protein to study its dynamic quantitatively?? Or is it better to use just the coupling on the solvent and not protein??
Also, Is it just GROMACS that have the ability to assign the thermostat to different groups with a different coupling constant?? I mean, it is a very hard decision, and it adds up to my doubt about not using the thermostat on protein and just on solvent when I see that AMBER MD engine doesn't even have the option to assign different coupling constant to different groups.
Best