The gas_energy polargrid is defined in function main.c. Its initialization is done in function
InitGasEnergy(see Pframeforce.c), and is consistent with that of the disk surface density (
InitGasDensity, in Pframeforce.c) and temperature for a locally isothermal disk:
where the value of is defined by
ADIABATICINDEX, that of by
SIGMASLOPE, and by
FLARINGINDEX. Note that the rotation profile is the same either with an energy equation or with a locally isothermal equation of state (the initialization of the gas velocity is now implemented in
InitGasVelocities, in Pframeforce.c).