Temperature should be read only
Temperature is defined as:
myTemperature->values() = myEnergy->values()/(specificHeat*myDensity.values());
and "refreshed" as such here and there according to Disk::computeTemperature
.
In some boundary conditions, we set the outer most "layers" directly. That's when we do stuff like
if (first proc) {
temperature[rings(0)] = temperatur[ring(1)];
as long as we do the same thing as energy and density.
BUT from time to time, for real 3D specific BC, we set it to
real zTemperature = physic.adiabaticEoS->fullEnergyEoS->zBoundaryTemperature();
temperature[bot0] = zTemperature;
but even, density and energy are actually modified accordingly.
This is an indication that temperature does not need to be stored, read nor setGhost called on it.
@elena : am I missing something ?