public void BufferFlip(PlanetData PD)
{
Logger("BufferFlip! cycle = " + CellUpdater.run);
List <KWSCellMap <WeatherCell> > temp = PD.LiveMap;
PD.LiveMap = PD.BufferMap;
PD.BufferMap = temp;
KWSCellMap <SoilCell> temp1 = PD.LiveSoilMap;
PD.LiveSoilMap = PD.BufferSoilMap;
PD.BufferSoilMap = temp1;
List <KWSCellMap <WeatherCell> > temp2 = PD.LiveStratoMap;
PD.LiveStratoMap = PD.BufferStratoMap;
PD.BufferStratoMap = temp2;
PD.updateTime = currentTime; //update the update time
WeatherDatabase.PlanetaryData[PD.index] = PD;
CellUpdater.run++;
currentTime = 0;
MapOverlay.refreshCellColours();
//if(CellUpdater.run % 50 == 0) WeatherDatabase.SavePlanetaryData(cellindex, HighLogic.SaveFolder);
}
private static KWSCellMap <SoilCell> InitSoilCalcs(PlanetData PD)
{
KWSCellMap <SoilCell> tempMap = new KWSCellMap <SoilCell>(PD.gridLevel);
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
SoilCell wCell = tempMap[cell];
wCell.temperature = GetInitTemperature(PD, 0, cell);
tempMap[cell] = wCell;
}
return(tempMap);
}
public static void PrintStat(PlanetData PD) //TODO: make for multi-threading, this method will be inherently safe (doesn't write, and only uses objects not being changed)
{
int layerCount = PD.LiveMap.Count;
int stratoCount = PD.LiveStratoMap.Count;
KWSCellMap <WeatherCell> cellMap = new KWSCellMap <WeatherCell>(PD.gridLevel);
float avgValue = 0;
float maxValue = 0;
float minValue = 0;
uint maxIndex = 0;
uint minIndex = 0;
// Cell thisCell = new Cell();
string[] values = { "temperature", "pressure", "relativeHumidity", "CCN" }; //, "cloud.dropletSize", "cloud.rainyDuration", "windvector.y"}; //the values we want to get from the cell
// check if Stat folder exists, create it
if (!Directory.Exists(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat"))
{
Directory.CreateDirectory(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat");
}
//write to Statistics Log file
Debug.Log("Writing Statistics file");
int num = Directory.GetFiles(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat/").Length;
using (System.IO.StreamWriter file =
new System.IO.StreamWriter(@KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat/Stat" + num + "cy" + CellUpdater.run + ".txt"))
{
// int layer = 0;
String STS = " | ";
file.WriteLine("Body: " + PD.body.bodyName + STS + "Update Cycle: " + CellUpdater.run + STS + "AvgProcessTime (μs): " + WeatherSimulator.AvgCycleTime);
file.WriteLine();
// print: Total cells, other general data about the collection
file.WriteLine("Number of cells: " + Cell.CountAtLevel(PD.gridLevel));
file.WriteLine();
foreach (string s in values) // each valid s prints a section
{
file.WriteLine(s);
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = stratoCount - 1; i >= 0; i--)
{
cellMap = PD.LiveStratoMap[i];
if (cellMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
{
avgValue = cellMap.Average(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
minValue = cellMap.Min(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
maxValue = cellMap.Max(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue)
{
minIndex = cell.Index;
}
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue)
{
maxIndex = cell.Index;
}
}
}
else
{
avgValue = 0; minValue = 0; maxValue = 0; minIndex = 0; maxIndex = 0;
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "St_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
if (cellMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
{
avgValue = cellMap.Average(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
minValue = cellMap.Min(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
maxValue = cellMap.Max(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue)
{
minIndex = cell.Index;
}
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue)
{
maxIndex = cell.Index;
}
}
}
else
{
avgValue = 0; minValue = 0; maxValue = 0; minIndex = 0; maxIndex = 0;
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
// Soil layer
if (s == "temperature")
{
avgValue = PD.LiveSoilMap.Average(cell => cell.Value.temperature);
minValue = PD.LiveSoilMap.Min(cell => cell.Value.temperature);
maxValue = PD.LiveSoilMap.Max(cell => cell.Value.temperature);
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if (PD.LiveSoilMap[cell].temperature == minValue)
{
minIndex = cell.Index;
}
if (PD.LiveSoilMap[cell].temperature == maxValue)
{
maxIndex = cell.Index;
}
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Soil", STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
/*
* if (PD.LiveSoilMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
* {
* avgValue = PD.LiveSoilMap.Average(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
* minValue = PD.LiveSoilMap.Min(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
* maxValue = PD.LiveSoilMap.Max(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
* foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
* {
* if ((float)PD.LiveSoilMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue) { minIndex = cell.Index; }
* if ((float)PD.LiveSoilMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue) { maxIndex = cell.Index; }
* }
* file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Soil", STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
* }
*/
file.WriteLine(); // closing section
}
file.WriteLine("cloud.dropletSize");
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
avgValue = cellMap.Average(cell => cell.Value.cloud.dropletSize);
minValue = cellMap.Min(cell => cell.Value.cloud.dropletSize);
maxValue = cellMap.Max(cell => cell.Value.cloud.dropletSize);
minIndex = cellMap.Where(cell => cell.Value.cloud.dropletSize == minValue).First().Key.Index;
maxIndex = cellMap.Where(cell => cell.Value.cloud.dropletSize == maxValue).First().Key.Index;
/*
* foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
* {
* if (cellMap[cell].cloud.dropletSize == minValue) { minIndex = cell.Index; }
* if (cellMap[cell].cloud.dropletSize == maxValue) { maxIndex = cell.Index; }
* }
*/
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
file.WriteLine(); // closing section
file.WriteLine("cloud.rainyDuration");
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
avgValue = cellMap.Average(cell => (float)cell.Value.cloud.rainyDuration);
minValue = cellMap.Min(cell => (float)cell.Value.cloud.rainyDuration);
maxValue = cellMap.Max(cell => (float)cell.Value.cloud.rainyDuration);
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if (cellMap[cell].cloud.rainyDuration == minValue)
{
minIndex = cell.Index;
}
if (cellMap[cell].cloud.rainyDuration == maxValue)
{
maxIndex = cell.Index;
}
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
file.WriteLine(); // closing section
file.WriteLine("windVector.y");
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
avgValue = cellMap.Average(cell => (float)cell.Value.windVector.y);
minValue = cellMap.Min(cell => (float)cell.Value.windVector.y);
maxValue = cellMap.Max(cell => (float)cell.Value.windVector.y);
minIndex = cellMap.Where(cell => cell.Value.windVector.y == minValue).First().Key.Index;
maxIndex = cellMap.Where(cell => cell.Value.windVector.y == maxValue).First().Key.Index;
/*
* foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
* {
* if (cellMap[cell].windVector.y == minValue) { minIndex = cell.Index; }
* if (cellMap[cell].windVector.y == maxValue) { maxIndex = cell.Index; }
* }
*/
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
file.WriteLine(); // closing section
file.WriteLine("Horizontal Wind");
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
avgValue = cellMap.Average(cell => (float)Math.Sqrt(cell.Value.windVector.x * cell.Value.windVector.x + cell.Value.windVector.z * cell.Value.windVector.z));
minValue = cellMap.Min(cell => (float)Math.Sqrt(cell.Value.windVector.x * cell.Value.windVector.x + cell.Value.windVector.z * cell.Value.windVector.z));
maxValue = cellMap.Max(cell => (float)Math.Sqrt(cell.Value.windVector.x * cell.Value.windVector.x + cell.Value.windVector.z * cell.Value.windVector.z));
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
float WsH = (float)Math.Sqrt(cellMap[cell].windVector.x * cellMap[cell].windVector.x + cellMap[cell].windVector.z * cellMap[cell].windVector.z);
if (WsH == minValue)
{
minIndex = cell.Index;
}
if (WsH == maxValue)
{
maxIndex = cell.Index;
}
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
file.WriteLine(); // closing section
/* Section to print cell positional data
* file.WriteLine("cell "+ " .Lat.. .Long.." + " L" + " neigh .Lat.. .Long.. Dist. Dir(°) " + " neigh .Lat.. .Long.. Dist. Dir(°) "
+ " neigh .Lat.. .Long.. Dist. Dir(°) " + " neigh .Lat.. .Long.. Dist. Dir(°) "
+ " neigh .Lat.. .Long.. Dist. Dir(°) " + " neigh .Lat.. .Long.. Dist. Dir(°) ");
+ Vector3d North = new Vector3d(0.0, 1.0, 0.0);
+ Vector3d Up = new Vector3d(0.0, 1.0, 0.0);
+
+ foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
+ {
+ file.Write("{0,5}", cell.Index);
+ file.Write(" ");
+ //file.Write(cell.Position);
+ file.Write(String.Format("{0:+00.00;-00.00}", WeatherFunctions.GetCellLatitude(cell)));
+ file.Write(" ");
+ file.Write(String.Format("{0:+000.00;-000.00}", WeatherFunctions.GetCellLongitude(cell)));
+ file.Write(" ");
+ file.Write(cell.Level);
+ file.Write(" ");
+ foreach (Cell neighbor in cell.GetNeighbors(PD.gridLevel))
+ {
+ file.Write("{0,5}", neighbor.Index);
+ file.Write(" ");
+ file.Write(String.Format("{0:+00.00;-00.00}", WeatherFunctions.GetCellLatitude(neighbor)));
+ file.Write(" ");
+ file.Write(String.Format("{0:+000.00;-000.00}", WeatherFunctions.GetCellLongitude(neighbor)));
+ file.Write(" ");
+ file.Write(String.Format("{0:00000}", WeatherFunctions.GetDistanceBetweenCells(PD.index, PD.LiveSoilMap[cell].centroid, PD.LiveSoilMap[neighbor].centroid, 0)));
+ // file.Write("{0,5:N1}", WeatherFunctions.GetCentroidDirection(PD.LiveSoilMap[cell], PD.LiveSoilMap[neighbor]) * Mathf.Rad2Deg);
+ float Dir = WeatherFunctions.GetDirectionBetweenCells(PD.LiveSoilMap[cell].centroid, PD.LiveSoilMap[neighbor].centroid);
+ file.Write("{0,7:N2}", Dir * Mathf.Rad2Deg + (Dir < 0 ? 360 : 0));
+ // float DeltaLon = WeatherFunctions.GetCellLongitude(neighbor) - WeatherFunctions.GetCellLongitude(cell);
+ // file.Write("{0,5:N0}", Mathf.Rad2Deg*Math.Atan2((DeltaLon > 180 ? DeltaLon - 360 : DeltaLon < -180 ? DeltaLon + 360 : DeltaLon) * Math.Cos(WeatherFunctions.GetCellLatitude(cell)*Mathf.Deg2Rad),(WeatherFunctions.GetCellLatitude(neighbor) - WeatherFunctions.GetCellLatitude(cell))));
+ file.Write(" ");
+ }
+
+ file.WriteLine();
+ }
*/
}
}
public static void PrintStat(PlanetData PD) //TODO: make for multi-threading, this method will be inherently safe (doesn't write, and only uses objects not being changed)
{
int layerCount = PD.LiveMap.Count;
int stratoCount = PD.LiveStratoMap.Count;
KWSCellMap <WeatherCell> cellMap = new KWSCellMap <WeatherCell>(PD.gridLevel);
float avgValue = 0;
float maxValue = 0;
float minValue = 0;
uint maxIndex = 0;
uint minIndex = 0;
// Cell thisCell = new Cell();
string[] values = { "temperature", "pressure", "relativeHumidity", "CCN", "cloud.dropletSize", "cloud.rainyDuration", "windvector.y" }; //the values we want to get from the cell
// check if Stat folder exists, create it
if (!Directory.Exists(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat"))
{
Directory.CreateDirectory(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat");
}
//write to Statistics Log file
Debug.Log("Writing Statistics file");
int num = Directory.GetFiles(KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat/").Length;
using (System.IO.StreamWriter file =
new System.IO.StreamWriter(@KSPUtil.ApplicationRootPath + "/GameData/KerbalWeatherSystems/Stat/Stat" + num + "cy" + CellUpdater.run + ".txt"))
{
// int layer = 0;
String STS = " | ";
file.WriteLine("Body: " + PD.body.bodyName + STS + "Update Cycle: " + CellUpdater.run + STS + "AvgProcessTime (μs): " + WeatherSimulator.AvgCycleTime);
file.WriteLine();
// print: Total cells, other general data about the collection
file.WriteLine("Number of cells: " + Cell.CountAtLevel(PD.gridLevel));
file.WriteLine();
foreach (string s in values) // each valid s prints a section
{
file.WriteLine(s);
file.WriteLine("layer" + STS + " average " + STS + " minimum " + " cell " + STS + " maximum " + " cell " + STS);
for (int i = stratoCount - 1; i >= 0; i--)
{
cellMap = PD.LiveStratoMap[i];
if (cellMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
{
avgValue = cellMap.Average(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
minValue = cellMap.Min(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
maxValue = cellMap.Max(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue)
{
minIndex = cell.Index;
}
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue)
{
maxIndex = cell.Index;
}
}
}
else
{
avgValue = 0; minValue = 0; maxValue = 0; minIndex = 0; maxIndex = 0;
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "St_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
for (int i = layerCount - 1; i >= 0; i--)
{
cellMap = PD.LiveMap[i];
if (cellMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
{
avgValue = cellMap.Average(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
minValue = cellMap.Min(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
maxValue = cellMap.Max(cell => cell.Value.GetType().GetProperty(s).Name.Equals(s) ? (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null) : 0f);
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue)
{
minIndex = cell.Index;
}
if ((float)cellMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue)
{
maxIndex = cell.Index;
}
}
}
else
{
avgValue = 0; minValue = 0; maxValue = 0; minIndex = 0; maxIndex = 0;
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Tp_" + i, STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
// Soil layer
if (PD.LiveSoilMap[new Cell(0)].GetType().GetProperty(s) != null && cellMap[new Cell(0)].GetType().GetProperty(s).Name.Equals(s))
{
avgValue = PD.LiveSoilMap.Average(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
minValue = PD.LiveSoilMap.Min(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
maxValue = PD.LiveSoilMap.Max(cell => (float)cell.Value.GetType().GetProperty(s).GetValue(cell.Value, null));
foreach (Cell cell in Cell.AtLevel(PD.gridLevel))
{
if ((float)PD.LiveSoilMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == minValue)
{
minIndex = cell.Index;
}
if ((float)PD.LiveSoilMap[cell].GetType().GetProperty(s).GetValue(cellMap[cell], null) == maxValue)
{
maxIndex = cell.Index;
}
}
file.WriteLine("{0,-4} {1,5} {2,12:N4} {3,5} {4,12:N4} {5,6} {6,5} {7,12:N4} {8,6} {9,5}", "Soil", STS, avgValue, STS, minValue, minIndex, STS, maxValue, maxIndex, STS);
}
file.WriteLine(); // closing section
// avgValue = cellMap.Average(cell => cell.Value.cloud.GetType().GetProperty(s).GetValue(cell.Value, null));
}
}
}