* Implement fluxvar specializations with energy on
* EnergyLocalResidual helps to add energy contributions efficiently
* The base implicit volVars now exist in an isothermal and a non-isothermal version
* Make 1p non-isothermal cc models run

* TODO
	* Fourier's law for box
	* Adjust volume variables / local residuals of other models