Calculating the heat stored in a solid sub-domain

Hi,

I am solving a time-dependent heat transfer problem on a solid domain. The specific heat capacity (Cp) is a function of temperature. I would like to calculate the amount of heat stored in specific sub-domains as a function of time. Since Cp is not constant, I need to perform a temperature integration on Cp before I do a domain integration.

How can I do this in COMSOL? Any ideas/help would be greatly appreciated.

Thanks.

Ismail