Just pinged you some (alpha) code.
As for the calculation, its basically Q = m Cp dT
I enjoy doing this work and helping people out ....but in full disclosure I'll be using this same approach on my rig, so it will be nice to have worked out the bugs prior to implementation on my build. :mug: