[FLASH-USERS] How does FLASH3 calculate temperatures for checkpoint and plot dumps?

[Dave]

Hi,

When examining the checkpoint file from a simple test run of a galaxy, 
I'm finding that the pres and dens values are what I expect, but the 
temp values seem to be way off. In an (unrealistic) supposedly 
isothermal test, the outer parts of the galaxy are ~10^7 K, as they 
should be, but there is a sharp discontinuity in the centre, where the 
temperature drops to ~10^7 K, according to the temp scalar - even though 
the pres and dens values are what they should be.

I'm trying to track down where in Flash the temperature is calculated 
from internal energy etc for output, to determine if this is just being 
dumped out wrong, or if the internal energy itself is wrong. Any idea 
where that goes on?

Cheers,
-David