Hi, Justin, Thanks for your response. By extracting temperature information out edr file, do you mean the "gmx enemat" command? Actually what I was up to is calculate the temperature development of a certain region instead of a group of atoms.