Working version of the unsteady heat solver