public MeanWaterLevelInterpretation(Dam dam, double waterLevel)
{
this.minWaterLevel = dam.Reservoir.MinimumSupplyLevel;
this.maxFloodLevel = dam.Reservoir.MaximumFloodLevel;
this.capacity = dam.Reservoir.StorageCapacity;
this.waterLevel = waterLevel;
}
private Dictionary <int, Rainfall> ComputeWaterLevel(Dam dam)
{
Dictionary <int, Rainfall> result = new Dictionary <int, Rainfall>();
Dictionary <int, List <Rainfall> > stations = new Dictionary <int, List <Rainfall> >();
foreach (DamStation station in dam.Stations)
{
int id = station.Id;
List <Rainfall> rainFalls = new List <Rainfall>();
foreach (Forecast forecast in station.Forecasts)
{
Rainfall rainfall = new Rainfall(forecast.Timestamp, forecast.Rain);
rainFalls.Add(rainfall);
}
stations.Add(id, rainFalls);
}
foreach (int key in stations.Keys)
{
List <Rainfall> rainfalls = stations[key];
double rains = 0;
DateTime ts = new DateTime();
foreach (Rainfall rainfall in rainfalls)
{
rains += rainfall.Amount;
ts = rainfall.DateTime;
}
double meanPrecipitation = rains / rainfalls.Count;
Rainfall newRainfall = new Rainfall(ts, meanPrecipitation);
result.Add(key, newRainfall);
}
return(result);
}
public Dictionary <int, Rainfall> GetRainFall(int damId)
{
Dam dam = getDam(damId);
return(ComputeRainFall(dam));
}
public Dictionary <int, Rainfall> GetWaterLevel(int damId)
{
Dam dam = getDam(damId);
return(ComputeWaterLevel(dam)); // from mmm to m
}
public static void Main()
{
Dictionary <int, Dam> dams = new Dictionary <int, Dam>();
Location damLocation = new Location(16.9040145, 121.3094922, 12 /*11.75 */);
double dischargeCap = 30600; // mmm
double storageCapacity = 1080000000; // 1.8B mmm
double minSupply = 160; // masl
double maxSupply = 193; // masl
double observedDamWaterAmount = 108000000;
Reservoir reservoir = new Reservoir(storageCapacity, minSupply, maxSupply);
// STATIONS, FORECASTS
List <Forecast> forecasts = new List <Forecast>();
Forecast forecast = new Forecast(DateTime.Now, 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
forecast = new Forecast(DateTime.Now.AddHours(1), 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
Location stationLocation = new Location(11.3581, 120.7242, 26);
List <DamStation> stations = new List <DamStation>();
DamStation station = new DamStation(980001, "Pamalican (Amanpulo) 1", stationLocation, forecasts);
stations.Add(station);
station = new DamStation(980002, "Pamalican (Amanpulo) 2", stationLocation, forecasts);
stations.Add(station);
Dam dam = new Dam("Magat", damLocation, reservoir, dischargeCap, stations, observedDamWaterAmount);
dams.Add(1, dam);
damLocation = new Location(16.4587389, 120.7430438, 17);
dischargeCap = 30000; // mmm
storageCapacity = 121000000; // 121M mmm
minSupply = 740; // masl
maxSupply = 752; // masl
reservoir = new Reservoir(storageCapacity, minSupply, maxSupply);
observedDamWaterAmount = 12100000;
// STATIONS, FORECASTS
forecasts = new List <Forecast>();
forecast = new Forecast(DateTime.Now.AddHours(2), 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
forecast = new Forecast(DateTime.Now.AddHours(3), 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
stationLocation = new Location(11.3581, 120.7242, 26);
stations = new List <DamStation>();
station = new DamStation(980003, "Pamalican (Amanpulo) 1", stationLocation, forecasts);
stations.Add(station);
station = new DamStation(980004, "Pamalican (Amanpulo) 2", stationLocation, forecasts);
stations.Add(station);
dam = new Dam("Ambuklao", damLocation, reservoir, dischargeCap, stations, observedDamWaterAmount);
dams.Add(2, dam);
damLocation = new Location(16.3948684, 120.7279512, 17);
dischargeCap = 10521; // mmm
storageCapacity = 21000000; // 121M mmm
minSupply = 757; // masl
maxSupply = 566; // masl
reservoir = new Reservoir(storageCapacity, minSupply, maxSupply);
observedDamWaterAmount = 22000000;
// STATIONS, FORECASTS
forecasts = new List <Forecast>();
forecast = new Forecast(DateTime.Now.AddHours(4), 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
forecast = new Forecast(DateTime.Now.AddHours(5), 28.65, 36, 0, 1010.85, 0.2, 25.2, 44, 54, 34, 29, 34);
forecasts.Add(forecast);
stationLocation = new Location(11.3581, 120.7242, 26);
stations = new List <DamStation>();
station = new DamStation(980003, "Pamalican (Amanpulo) 1", stationLocation, forecasts);
stations.Add(station);
station = new DamStation(980004, "Pamalican (Amanpulo) 2", stationLocation, forecasts);
stations.Add(station);
dam = new Dam("Binga", damLocation, reservoir, dischargeCap, stations, observedDamWaterAmount);
dams.Add(3, dam);
MeanPrecipitationCalculator meanPrecipitation = new MeanPrecipitationCalculator(dams);
foreach (int damId in dams.Keys)
{
Console.WriteLine("\nDam Id: " + damId + " Name: " + dams[damId].Name);
Dictionary <int, Rainfall> meanRainfalls = meanPrecipitation.GetRainFall(damId);
foreach (int key in meanRainfalls.Keys)
{
Rainfall rainfall = meanRainfalls[key];
Console.WriteLine("Station Id: " + key
+ ", Timestamp: " + rainfall.DateTime + ", Mean Rainfall: " + rainfall.Amount
+ ", Interpretation: " + new MeanRainfallInterpretation(rainfall).Interpret()
);
}
Dictionary <int, Rainfall> meanWaterLevel = meanPrecipitation.GetWaterLevel(damId);
foreach (int key in meanWaterLevel.Keys)
{
Rainfall rainfall = meanWaterLevel[key];
double waterLevel = rainfall.Amount + dam.ObservedVolume;
Console.WriteLine("Station Id: " + key
+ ", Timestamp: " + rainfall.DateTime + ", Mean Water Level: " + waterLevel
+ ", Interpretation: " + new MeanWaterLevelInterpretation(
dams[damId], waterLevel).Interpret()
);
}
}
}