public void CompareStreamAndLakes() { var StreamNetwork = NetworkBuilder.CreateBranch(10); var LakeNetwork = NetworkBuilder.CreateConnectedLakes(10); Model Streams = new Model(); Streams._waterBodies.AddRange(StreamNetwork.Cast <IWaterBody>()); Model Lakes = new Model(); Lakes._waterBodies.AddRange(LakeNetwork.Cast <IWaterBody>()); SinkSourceBoundary b1 = new SinkSourceBoundary(100); StreamNetwork.First().Sources.Add(b1); LakeNetwork.First().Sources.Add(b1); Stopwatch SW = new Stopwatch(); Stopwatch SW2 = new Stopwatch(); DateTime Start = new DateTime(2000, 1, 1); DateTime End = new DateTime(2000, 1, 10); Streams.SetState("Initial", Start, new WaterPacket(1)); Lakes.SetState("Initial", Start, new WaterPacket(1)); SW.Start(); Streams.MoveInTime(End, TimeSpan.FromHours(5)); SW.Stop(); SW2.Start(); Lakes.MoveInTime(End, TimeSpan.FromHours(5)); SW2.Stop(); TimespanSeries TS1 = StreamNetwork.Last().Output.Items.First() as TimespanSeries; TimespanSeries TS2 = LakeNetwork.Last().Output.Items.First() as TimespanSeries; for (int i = 0; i < TS1.Items.Count; i++) { Assert.AreEqual(TS1.Items[i].Value, TS2.Items[i].Value, 0.000001); Assert.AreEqual(TS1.Items[i].StartTime, TS2.Items[i].StartTime); } }
public void TracerTest() { int count = 3; double length = 10870; DateTime Start = new DateTime(2000, 1, 1); List <Lake> lakes = NetworkBuilder.CreateConnectedLakes(count); foreach (Lake L in lakes) { L.SurfaceArea = XYPolygon.GetSquare(length / count); L.Depth = 1; L.SetState("Initial", Start, new WaterPacket(L.Volume)); } Chemical c = ChemicalFactory.Instance.GetChemical(ChemicalNames.Cl); SinkSourceBoundary fb = new SinkSourceBoundary(10870.0 / (8.49 * 3600)); fb.WaterSample = new WaterPacket(1); lakes.First().Sources.Add(fb); WaterPacket plug = new WaterPacket(1); plug.AddChemical(c, 10000); lakes.First().AddWaterPacket(Start, Start.AddHours(1), plug.DeepClone()); lakes.First().Output.LogChemicalConcentration(c); lakes.Last().Output.LogChemicalConcentration(c); Stream us = new Stream("us", 1, 1, 1); us.Sources.Add(fb); Stream s = new Stream("s", 1000, 1, 1); Lake L2 = new Lake("L2", 870); Stream s1 = new Stream("s1", 9000, 1, 1); s.AddWaterPacket(Start, Start.AddSeconds(1), plug.DeepClone()); s1.Output.LogChemicalConcentration(c); us.AddDownStreamWaterBody(s); s.AddDownStreamWaterBody(L2); L2.AddDownStreamWaterBody(s1); Model m = new Model(); m._waterBodies.AddRange(lakes.Cast <IWaterBody>()); m._waterBodies.Add((IWaterBody)us); m._waterBodies.Add((IWaterBody)s); m._waterBodies.Add((IWaterBody)L2); m._waterBodies.Add((IWaterBody)s1); m.SetState("Initial", Start, new WaterPacket(1)); m.MoveInTime(Start.AddHours(15), TimeSpan.FromMinutes(1)); lakes.Last().Output.Save(@"C:\temp\LastLake.xts"); s1.Output.Save(@"C:\temp\Stream.xts"); int n = 15; List <IWaterBody> wbs = NetworkBuilder.CreateCombo(n, 10870 / n / 2.0); foreach (IWaterBody wb in wbs) { wb.SetState("Initial", Start, new WaterPacket(wb.Volume)); } wbs.First().AddWaterPacket(Start, Start.AddHours(1), plug.DeepClone()); us.AddDownStreamWaterBody(wbs.First()); us.RestoreState("Initial"); m._waterBodies.Clear(); m._waterBodies.Add(us); m._waterBodies.AddRange(wbs); m.SetState("Initial", Start, new WaterPacket(1)); ((Stream)wbs.Last()).Output.LogChemicalConcentration(c); m.MoveInTime(Start.AddHours(15), TimeSpan.FromMinutes(1)); ((Stream)wbs.Last()).Output.Save(@"C:\temp\Stream.xts"); }