float   TSolarSystem::CalcBodyTemperature(TSpaceBody &CBody, double time)
{
      vec3d  bodyPos = GetBodyBarycentricPos(&CBody, time);
      float  Temperature = 0.0f;
      float  TempExo = 0.0f;

      for (int i=0; i<NSunsBarys; ++i)
      {
          TSpaceBody &CSun = Body[i];
          if (CSun.Type != PT_SUN) continue;
          double sunDist = (GetBodyBarycentricPos(&CSun, time) - bodyPos).length() / km;
          double solAngRad = CSun.Surface.Radius / sunDist;
          if (solAngRad >= 1.0) solAngRad = 1.0;
          Temperature += CSun.StarData.GetTemperature() * sqrt(0.5 * asin(solAngRad));
          TempExo += 1500.0f * pow(sqrt(CSun.StarData.GetLuminosityUV() / 0.072f) / float(sunDist*km/AU), 0.21f);
      }

      // For slowly rotating planets
      //if (CBody.Surface.RotPeriod > 10.0f) Temperature *= 1.189207115f;

      Temperature = Temperature * pow(1.0f - CBody.Albedo, 0.25f) + PlanetMinTemperature;
      if (CBody.Atmosphere.Greenhouse != NO_DATA)
          Temperature += CBody.Atmosphere.Greenhouse;

      CBody.Atmosphere.TempExo = max2f(TempExo + PlanetMinTemperature, Temperature);
      CBody.Temperature = Temperature;

      return CBody.Temperature;
}

         /*Calculates a Stars total output
          */
         public static double getStarFlux(Star s)
         {
             double fluxM = Math.Pow(s.temperature, 4) * (Helper.sigma); //sigma = 5.6704E-008;
             double area = 4 * Math.PI * Math.Pow(s.radius * 1000, 2);
             return fluxM * area;

         }
         /*Gets a stars output at a given distance
          */
         public static double getStellarConstant(double starFlux, double distance)
         {
             double area = 4 * Math.PI * (distance * distance);
             return starFlux / area;

         }

         public static double getPlanetTemp(double starConstant, double abledo, double greenhouse)
         {
             double fluxAbsorbed = (starConstant / 4) * (1 - abledo);

             return Math.Pow((fluxAbsorbed / Helper.sigma), .25) + greenhouse + Helper.CMB; //CMB = 2.725
         }