full energy configuration
@elena Ok, I see that it remains confusing (unfortunately it comes from the fact that I had 2 different pieces of code for fullenergy and star and I badly merged them at some point)
let's try to clarify (not sure):
The gas thermodynamics is considered in a non-equilibrium situation, the gas and the radiation having two different temperatures. The thermal energy density (Egas ) and the radiation energy density (Erad ) are governed by a set of two coupled equations solved by RadiationManager::solve
Now there are 2 optional terms in this set of equations:
- stellarrad = 0 if we do not want to take into account radiation coming from the star
- planetrad = 0 if we do not want to take into account radiation from delivery of solids onto a planetary core
When we ask for FullEnergyEoS true we need to have Tstar =0
so that:
real fStar = cfg().fStar();
-> fStar == 0 if Tstar == 0
if (fStar > 0) { // optimisation, fStar == 0 would not do anything
in this way stellarrad = 0
and by the way we need StarRadiation == true to have radiative energy
while StarRadiation should be true only to receive a temperature of the star different from zero and a stellar radius to compute fstar
tell me if it remains obscure ...