static internal double minLOSCalc(bodydef body)
{
var sat = GraphWindow.numSats;
if (sat >= 3)
{
return((body.eqr * occmod(body)) / (Math.Cos(0.5 * (2f * Math.PI / sat))) - body.eqr);
}
return(0);
}
static internal double occmod(bodydef body)
{
if (!GraphWindow.occlusionModifiers)
{
return(1);
}
if (body.atm > 0)
{
return(GraphWindow.atmOcclusion);
}
else
{
return(GraphWindow.vacOcclusion);
}
}
static double burnCalc(orbitdef s, orbitdef c, bodydef b)
{
double sta = 0;
double cta = 0;
if (c.Ap == s.Ap)
{
cta = 180;
}
double sr = s.SMA * (1 - Math.Pow(s.e, 2)) / (1 + (s.e * Math.Cos(sta)));
double sdv = Math.Sqrt(b.GM * ((2 / sr) - (1 / s.SMA)));
double cr = c.SMA * (1 - Math.Pow(c.e, 2)) / (1 + (c.e * Math.Cos(cta)));
double cdv = Math.Sqrt(b.GM * ((2 / sr) - (1 / c.SMA)));
return(Math.Round(100 * Math.Abs(sdv - cdv)) / 100);
}
public static void Update()
{
body = new bodydef(PlanetSelection.selectedBody);
var sataltitude = GraphWindow.orbitAltitude;
if (sataltitude > body.SOIAlt())
{
disabled = true;
}
if (body.geoAlt > 0 && body.geoSMA < body.SOI)
{
synchrorbit = body.geoAlt.ToString("N") + " m";
}
else
{
synchrorbit = "n/a";
}
minLOS = minLOSCalc(body);
if (minLOS < body.SOIAlt() && minLOS > 0)
{
losorbit = minLOS.ToString("N") + " m";
losOrbitWarning = (minLOS < body.atm);
constellationWarning = sataltitude * 1.01 < minLOS;
}
else
{
losorbit = "n/a";
}
double satcount = GraphWindow.numSats;
double satratio = ((satcount + 1) / satcount);
satelliteorbit = new orbitdef(sataltitude, sataltitude, body);
carrierorbit = new orbitdef(sataltitude, sataltitude, body);
if (satcount > 0 && sataltitude > 0)
{
carrierorbit.modifyAp(newMAfromT(satelliteorbit.T * satratio, body));
if (carrierorbit.Ap > body.SOIAlt() || GraphWindow.flipOrbit == true)
{
satratio = ((satcount - 1) / satcount);
carrierorbit = new orbitdef(sataltitude, sataltitude, body);
carrierorbit.modifyPe(newMAfromT(satelliteorbit.T * satratio, body));
}
period = satelliteorbit.oph;
if (!periodEntry)
{
periodHour = satelliteorbit.op_p(orbitdef.timePos.hours).ToString();
periodMin = satelliteorbit.op_p(orbitdef.timePos.min).ToString();
periodSec = satelliteorbit.op_p(orbitdef.timePos.sec).ToString();
}
/* if (carrierorbit.Ap > 0) */
carrierAp = carrierorbit.Ap.ToString("N1") + " m";
carrierApWarning = (carrierorbit.Ap < body.atm || carrierorbit.Ap < 0);
/*if (carrierorbit.Pe > 0)*/
carrierPe = carrierorbit.Pe.ToString("N1") + " m";
carrierPeWarning = (carrierorbit.Pe < body.atm);
carrierPEUrgent = (carrierorbit.Pe <= 0);
carrierT = carrierorbit.oph;
dBurnDV = burnCalc(satelliteorbit, carrierorbit, body);
burnDV = dBurnDV.ToString("N2") + " m/s";
}
if (period == "")
{
header[0] = "";
header[1] = "";
}
else
{
string s = body.body.displayName.Substring(0, body.body.displayName.Length - 2);
header[0] = s;
header[0] += " " + satcount + "-satellite constellation";
header[1] = "Ap " + carrierAp;
header[1] += " Pe " + carrierPe;
header[1] += " Δv " + burnDV;
}
}
static public double newMAfromT(double T, bodydef body)
{
return(2 * Cbrt((body.GM * Math.Pow(T, 2f)) / 39.4784176)); // 39.4784176 = 4π^2
}
public orbitdef(double Ap, double Pe, bodydef body)
{
this.Ap = Ap;
this.Pe = Pe;
this.body = body;
}