Implement a temperature() function in the volvars and use it
instead of volvars.fluidstate().xtemperature().
This reuqired the implementation of two template functions
to ensure that the correct temperature() function is called
depending on the number of energy equations considered.