public void Analyze(CelestialBody body, double altitude_mult, bool sunlight)
{
// shortcuts
CelestialBody sun = FlightGlobals.Bodies[0];
this.body = body;
altitude = body.Radius * altitude_mult;
landed = altitude <= double.Epsilon;
breathable = Sim.Breathable(body) && landed;
atmo_factor = Sim.AtmosphereFactor(body, 0.7071);
sun_dist = Sim.Apoapsis(Lib.PlanetarySystem(body)) - sun.Radius - body.Radius;
Vector3d sun_dir = (sun.position - body.position).normalized;
solar_flux = sunlight ? Sim.SolarFlux(sun_dist) * (landed ? atmo_factor : 1.0) : 0.0;
albedo_flux = sunlight ? Sim.AlbedoFlux(body, body.position + sun_dir * (body.Radius + altitude)) : 0.0;
body_flux = Sim.BodyFlux(body, altitude);
total_flux = solar_flux + albedo_flux + body_flux + Sim.BackgroundFlux();
temperature = !landed || !body.atmosphere ? Sim.BlackBodyTemperature(total_flux) : body.GetTemperature(0.0);
temp_diff = Sim.TempDiff(temperature, body, landed);
orbital_period = Sim.OrbitalPeriod(body, altitude);
shadow_period = Sim.ShadowPeriod(body, altitude);
shadow_time = shadow_period / orbital_period;
zerog = !landed && (!body.atmosphere || body.atmosphereDepth < altitude);
RadiationBody rb = Radiation.Info(body);
RadiationBody sun_rb = Radiation.Info(sun);
gamma_transparency = Sim.GammaTransparency(body, 0.0);
extern_rad = PreferencesStorm.Instance.ExternRadiation;
heliopause_rad = extern_rad + sun_rb.radiation_pause;
magnetopause_rad = heliopause_rad + rb.radiation_pause;
inner_rad = magnetopause_rad + rb.radiation_inner;
outer_rad = magnetopause_rad + rb.radiation_outer;
surface_rad = magnetopause_rad * gamma_transparency;
storm_rad = heliopause_rad + PreferencesStorm.Instance.StormRadiation * (solar_flux > double.Epsilon ? 1.0 : 0.0);
}
public void Analyze(CelestialBody body, double altitude_mult, Planner.SunlightState sunlight)
{
this.body = body;
CelestialBody mainSun;
Vector3d sun_dir;
solar_flux = Sim.SolarFluxAtBody(body, true, out mainSun, out sun_dir, out sun_dist);
altitude = body.Radius * altitude_mult;
landed = altitude <= double.Epsilon;
atmo_factor = Sim.AtmosphereFactor(body, 0.7071);
solar_flux = sunlight == Planner.SunlightState.Shadow ? 0.0 : solar_flux * (landed ? atmo_factor : 1.0);
breathable = Sim.Breathable(body) && landed;
albedo_flux = sunlight == Planner.SunlightState.Shadow ? 0.0 : Sim.AlbedoFlux(body, body.position + sun_dir * (body.Radius + altitude));
body_flux = Sim.BodyFlux(body, altitude);
total_flux = solar_flux + albedo_flux + body_flux + Sim.BackgroundFlux();
temperature = !landed || !body.atmosphere ? Sim.BlackBodyTemperature(total_flux) : body.GetTemperature(0.0);
temp_diff = Sim.TempDiff(temperature, body, landed);
orbital_period = Sim.OrbitalPeriod(body, altitude);
shadow_period = Sim.ShadowPeriod(body, altitude);
shadow_time = shadow_period / orbital_period;
zerog = !landed && (!body.atmosphere || body.atmosphereDepth < altitude);
RadiationBody rb = Radiation.Info(body);
RadiationBody sun_rb = Radiation.Info(mainSun); // TODO Kopernicus support: not sure if/how this work with multiple suns/stars
gamma_transparency = Sim.GammaTransparency(body, 0.0);
// add gamma radiation emitted by body and its sun
var gamma_radiation = Radiation.DistanceRadiation(rb.radiation_r0, altitude) / 3600.0;
#if DEBUG_RADIATION
Lib.Log("Planner/EA: " + body + " sun " + mainSun + " alt " + altitude + " sol flux " + solar_flux + " aalbedo flux " + albedo_flux + " body flux " + body_flux + " total flux " + total_flux);
Lib.Log("Planner/EA: body surface radiation " + Lib.HumanReadableRadiation(gamma_radiation, false));
#endif
var b = body;
while (b != null && b.orbit != null && b != mainSun)
{
if (b == b.referenceBody)
{
break;
}
var dist = b.orbit.semiMajorAxis;
b = b.referenceBody;
gamma_radiation += Radiation.DistanceRadiation(Radiation.Info(b).radiation_r0, dist) / 3600.0;
#if DEBUG_RADIATION
Lib.Log("Planner/EA: with gamma radiation from " + b + " " + Lib.HumanReadableRadiation(gamma_radiation, false));
Lib.Log("Planner/EA: semi major axis " + dist);
#endif
}
extern_rad = Settings.ExternRadiation / 3600.0;
heliopause_rad = gamma_radiation + extern_rad + sun_rb.radiation_pause;
magnetopause_rad = gamma_radiation + heliopause_rad + rb.radiation_pause;
inner_rad = gamma_radiation + magnetopause_rad + rb.radiation_inner;
outer_rad = gamma_radiation + magnetopause_rad + rb.radiation_outer;
surface_rad = magnetopause_rad * gamma_transparency + rb.radiation_surface / 3600.0;
storm_rad = heliopause_rad + PreferencesRadiation.Instance.StormRadiation * (solar_flux > double.Epsilon ? 1.0 : 0.0);
#if DEBUG_RADIATION
Lib.Log("Planner/EA: extern_rad " + Lib.HumanReadableRadiation(extern_rad, false));
Lib.Log("Planner/EA: heliopause_rad " + Lib.HumanReadableRadiation(heliopause_rad, false));
Lib.Log("Planner/EA: magnetopause_rad " + Lib.HumanReadableRadiation(magnetopause_rad, false));
Lib.Log("Planner/EA: inner_rad " + Lib.HumanReadableRadiation(inner_rad, false));
Lib.Log("Planner/EA: outer_rad " + Lib.HumanReadableRadiation(outer_rad, false));
Lib.Log("Planner/EA: surface_rad " + Lib.HumanReadableRadiation(surface_rad, false));
Lib.Log("Planner/EA: storm_rad " + Lib.HumanReadableRadiation(storm_rad, false));
#endif
}
