示例#1
0
            /// <summary>
            /// Update the 'sunsInfo' list and the 'mainSun', 'solarFluxTotal' variables.
            /// Uses discrete or analytic (for high timewarp speeds) evaluation methods based on the isAnalytic bool.
            /// Require the 'visibleBodies' variable to be set.
            /// </summary>
            // at the two highest timewarp speed, the number of sun visibility samples drop to the point that
            // the quantization error first became noticeable, and then exceed 100%, to solve this :
            // - we switch to an analytical estimation of the sunlight/shadow period
            // - atmo_factor become an average atmospheric absorption factor over the daylight period (not the whole day)
            public static void UpdateSunsInfo(VesselData vd, Vector3d vesselPosition)
            {
                Vessel v             = vd.Vessel;
                double lastSolarFlux = 0.0;

                vd.sunsInfo          = new List <SunInfo>(Sim.suns.Count);
                vd.solarFluxTotal    = 0.0;
                vd.rawSolarFluxTotal = 0.0;

                foreach (Sim.SunData sunData in Sim.suns)
                {
                    SunInfo sunInfo = new SunInfo(sunData);

                    if (vd.isAnalytic)
                    {
                        // get sun direction and distance
                        Lib.DirectionAndDistance(vesselPosition, sunInfo.sunData.body, out sunInfo.direction, out sunInfo.distance);
                        // analytical estimation of the portion of orbit that was in sunlight, current limitations :
                        // - the result is dependant on the vessel altitude at the time of evaluation,
                        //   consequently it gives inconsistent behavior with highly eccentric orbits
                        // - this totally ignore the orbit inclinaison, polar orbits will be treated as equatorial orbits
                        sunInfo.sunlightFactor = 1.0 - Sim.ShadowPeriod(v) / Sim.OrbitalPeriod(v);
                        // get atmospheric absorbtion
                        // for atmospheric bodies whose rotation period is less than 120 hours,
                        // determine analytic atmospheric absorption over a single body revolution instead
                        // of using a discrete value that would be unreliable at large timesteps :
                        if (vd.inAtmosphere)
                        {
                            sunInfo.atmoFactor = Sim.AtmosphereFactorAnalytic(v.mainBody, vesselPosition, sunInfo.direction);
                        }
                        else
                        {
                            sunInfo.atmoFactor = 1.0;
                        }
                    }
                    else
                    {
                        // determine if in sunlight, calculate sun direction and distance
                        sunInfo.sunlightFactor = Sim.IsBodyVisible(v, vesselPosition, sunData.body, vd.visibleBodies, out sunInfo.direction, out sunInfo.distance) ? 1.0 : 0.0;
                        // get atmospheric absorbtion
                        sunInfo.atmoFactor = Sim.AtmosphereFactor(v.mainBody, vesselPosition, sunInfo.direction);
                    }

                    // get resulting solar flux in W/m²
                    sunInfo.rawSolarFlux = sunInfo.sunData.SolarFlux(sunInfo.distance);
                    sunInfo.solarFlux    = sunInfo.rawSolarFlux * sunInfo.sunlightFactor * sunInfo.atmoFactor;
                    // increment total flux from all stars
                    vd.rawSolarFluxTotal += sunInfo.rawSolarFlux;
                    vd.solarFluxTotal    += sunInfo.solarFlux;
                    // add the star to the list
                    vd.sunsInfo.Add(sunInfo);
                    // the most powerful star will be our "default" sun. Uses raw flux before atmo / sunlight factor
                    if (sunInfo.rawSolarFlux > lastSolarFlux)
                    {
                        lastSolarFlux = sunInfo.rawSolarFlux;
                        vd.mainSun    = sunInfo;
                    }
                }

                vd.sunlightFactor = 0.0;
                foreach (SunInfo sunInfo in vd.sunsInfo)
                {
                    sunInfo.fluxProportion = sunInfo.rawSolarFlux / vd.rawSolarFluxTotal;
                    vd.sunlightFactor     += sunInfo.SunlightFactor * sunInfo.fluxProportion;
                }
                // avoid rounding errors
                if (vd.sunlightFactor > 0.99)
                {
                    vd.sunlightFactor = 1.0;
                }
            }
示例#2
0
            /// <summary>
            /// Update the 'sunsInfo' list and the 'mainSun', 'solarFluxTotal' variables.
            /// Uses discrete or analytic (for high timewarp speeds) evaluation methods based on the isAnalytic bool.
            /// Require the 'visibleBodies' variable to be set.
            /// </summary>
            // at the two highest timewarp speed, the number of sun visibility samples drop to the point that
            // the quantization error first became noticeable, and then exceed 100%, to solve this:
            // - we switch to an analytical estimation of the sunlight/shadow period
            // - atmo_factor become an average atmospheric absorption factor over the daylight period (not the whole day)
            public static void UpdateSunsInfo(VesselData vd, Vector3d vesselPosition, double elapsedSeconds)
            {
                Vessel v             = vd.Vessel;
                double lastSolarFlux = 0.0;

                vd.sunsInfo          = new List <SunInfo>(Sim.suns.Count);
                vd.solarFluxTotal    = 0.0;
                vd.rawSolarFluxTotal = 0.0;

                foreach (Sim.SunData sunData in Sim.suns)
                {
                    SunInfo sunInfo = new SunInfo(sunData);

                    if (vd.isAnalytic)
                    {
                        // get sun direction and distance
                        Lib.DirectionAndDistance(vesselPosition, sunInfo.sunData.body, out sunInfo.direction, out sunInfo.distance);
                        // analytical estimation of the portion of orbit that was in sunlight.
                        // it has some limitations, see the comments on Sim.ShadowPeriod

                        if (Settings.UseSamplingSunFactor)
                        {
                            // sampling estimation of the portion of orbit that is in sunlight
                            // until we will calculate again
                            sunInfo.sunlightFactor = Sim.SampleSunFactor(v, elapsedSeconds);
                        }

                        else
                        {
                            // analytical estimation of the portion of orbit that was in sunlight.
                            // it has some limitations, see the comments on Sim.ShadowPeriod
                            sunInfo.sunlightFactor = 1.0 - Sim.ShadowPeriod(v) / Sim.OrbitalPeriod(v);
                        }


                        // get atmospheric absorbtion
                        // for atmospheric bodies whose rotation period is less than 120 hours,
                        // determine analytic atmospheric absorption over a single body revolution instead
                        // of using a discrete value that would be unreliable at large timesteps :
                        if (vd.inAtmosphere)
                        {
                            sunInfo.atmoFactor = Sim.AtmosphereFactorAnalytic(v.mainBody, vesselPosition, sunInfo.direction);
                        }
                        else
                        {
                            sunInfo.atmoFactor = 1.0;
                        }
                    }
                    else
                    {
                        // determine if in sunlight, calculate sun direction and distance
                        sunInfo.sunlightFactor = Sim.IsBodyVisible(v, vesselPosition, sunData.body, vd.visibleBodies, out sunInfo.direction, out sunInfo.distance) ? 1.0 : 0.0;
                        // get atmospheric absorbtion
                        sunInfo.atmoFactor = Sim.AtmosphereFactor(v.mainBody, vesselPosition, sunInfo.direction);
                    }

                    // get resulting solar flux in W/m²
                    sunInfo.rawSolarFlux = sunInfo.sunData.SolarFlux(sunInfo.distance);
                    sunInfo.solarFlux    = sunInfo.rawSolarFlux * sunInfo.sunlightFactor * sunInfo.atmoFactor;
                    // increment total flux from all stars
                    vd.rawSolarFluxTotal += sunInfo.rawSolarFlux;
                    vd.solarFluxTotal    += sunInfo.solarFlux;
                    // add the star to the list
                    vd.sunsInfo.Add(sunInfo);
                    // the most powerful star will be our "default" sun. Uses raw flux before atmo / sunlight factor
                    if (sunInfo.rawSolarFlux > lastSolarFlux)
                    {
                        lastSolarFlux = sunInfo.rawSolarFlux;
                        vd.mainSun    = sunInfo;
                    }
                }

                vd.sunlightFactor = 0.0;
                foreach (SunInfo sunInfo in vd.sunsInfo)
                {
                    sunInfo.fluxProportion = sunInfo.rawSolarFlux / vd.rawSolarFluxTotal;
                    vd.sunlightFactor     += sunInfo.SunlightFactor * sunInfo.fluxProportion;
                }
                // avoid rounding errors
                if (vd.sunlightFactor > 0.99)
                {
                    vd.sunlightFactor = 1.0;
                }
            }