Use of FullEnergyEoS and StarRadiation
I apologise since the use of FullEnergyEoS and StarRadiation was badly done in the legacy fargOCA version
We discussed a lot to improve (see fullEnergy issue #305 (closed) ) , however there is still something I do not understand, actually dicovered by @astr-mle
In physicExtractor:
if (cfg.get<bool>("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation", false)) {
DiskPhysic::AdiabaticEoS::FullEnergyEoS::Star
starCfg = {
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation.TStar", real(0)),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation.RStar", real(1)),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.SolidAccretion", real(0)),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.OpacityLaw", OL_BELL_LIN),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.OpacityLaw.Value", std::numeric_limits<real>::quiet_NaN()),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.DustToGas", 0.01),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation.Epsilon", 5e-7),
cfg.get("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation.ShadowAngle", real(0))
From this I understand that the code reads and uses the OpacityLaw for example
if in the config file I have:
{
Index 1.4
FullEnergyEoS true
{
OpacityLaw Constant
{
Value 0.01
}
}
Actually OpacityLaw is apparently not used, in fact:
Equation of State, adiabatic
Equation of State, full energy
* zBoundaryTemperature (user): 3
* zBoundaryTemperature (used): 0.000146253
* zBoundaryEnergy (used): 8.49867e-17
(from OUT.15222252 on /scratch/elena/TestMPI)
and by the way results are puzling
Instead using:
AdiabaticEoS true
{
Index 1.4
FullEnergyEoS true
{
StarRadiation true ------> We require also StarRadiation although we shouldn't to get into:
cfg.get<bool>("Disk.AdiabaticEoS.FullEnergyEoS.StarRadiation", false)
OpacityLaw Constant
{
Value 0.01
}
}
}
We have:
Equation of State, full energy
* zBoundaryTemperature (user): 3
* zBoundaryTemperature (used): 0.000146253
* zBoundaryEnergy (used): 8.49867e-17
Star Radiation:
* temperature: 0
* radius: 1
* solid accretion 0
* opacityLaw: CONSTANT
* constant kappa: 0.01
* dust to gas: 0.01
* epsilon: 5e-07
* fStar: 0
(from OUT.15222251 /scratch/elena/TestMPIfdr) And the code provides expected results
Would it be possible to have the same behaviour as the simulation in /scratch/elena/TestMPIfdr (submission script subm3.sh and config file: config.info) without requiring StarRadiation true ?
I hope it is clear