private void ReadFile(string inputFile) { TimeBufferer output = new TimeInterpolator(); TimeSet timeSet = output.TTimeSet; IElementSet elementSet; StreamReader reader = new StreamReader(inputFile); // Read quantity and imply unit string line = GetNextLine(reader).Trim(' ', '"'); Quantity quantity = new Quantity(line); if (line.Equals("Flow", StringComparison.InvariantCultureIgnoreCase)) { quantity.Unit = new Unit(PredefinedUnits.CubicMeterPerSecond); } else if (line.Equals("WaterLevel", StringComparison.InvariantCultureIgnoreCase)) { quantity.Unit = new Unit(PredefinedUnits.Meter); } else { quantity.Unit = new Unit("Unspecified unit"); } output.ValueDefinition = quantity; // Read elementset line = GetNextLine(reader); string[] elements = line.Split(';'); string elementTypeString = elements[0].Trim('"'); if (elementTypeString.Equals("IdBased", StringComparison.InvariantCultureIgnoreCase)) { var idelementSet = new ElementSet(inputFile + "-" + quantity.Caption, inputFile, ElementType.IdBased, ""); for (int i = 1; i < elements.Length; i++) { idelementSet.AddElement(new Element(elements[i].Trim('"'))); } elementSet = idelementSet; } else if (elementTypeString.Equals("Points", StringComparison.InvariantCultureIgnoreCase)) { var pelementSet = new ElementSet(inputFile + "-" + quantity.Caption, inputFile, ElementType.Point, ""); for (int i = 1; i < elements.Length; i++) { string[] coordinates = elements[i].Trim('"', '(', ')').Split(','); if (coordinates.Length < 2) { throw new InvalidDataException("Invalid file format: only one coordinate for point: " + inputFile); } Element element = new Element(); double x = Double.Parse(coordinates[0], NumberFormatInfo.InvariantInfo); double y = Double.Parse(coordinates[1], NumberFormatInfo.InvariantInfo); Coordinate elmtCoor = new Coordinate(x, y); element.Vertices = new Coordinate[] { elmtCoor }; pelementSet.AddElement(element); } elementSet = pelementSet; } else if (elementTypeString.Equals("RegularGrid", StringComparison.InvariantCultureIgnoreCase)) { Spatial2DRegularGrid grid = new Spatial2DRegularGrid(); grid.IsNodeBased = false; string[] parts = elements[1].Trim('"', '(', ')').Split(','); grid.X0 = Double.Parse(parts[0], NumberFormatInfo.InvariantInfo); grid.Y0 = Double.Parse(parts[1], NumberFormatInfo.InvariantInfo); grid.Dx = Double.Parse(parts[2], NumberFormatInfo.InvariantInfo); grid.Dy = Double.Parse(parts[3], NumberFormatInfo.InvariantInfo); grid.XCount = Int32.Parse(parts[4], NumberFormatInfo.InvariantInfo); grid.YCount = Int32.Parse(parts[5], NumberFormatInfo.InvariantInfo); grid.Orientation = (parts.Length < 5) ? 0 : Int32.Parse(parts[6], NumberFormatInfo.InvariantInfo); elementSet = new Spatial2DGridWrapper(grid); } else { throw new InvalidDataException("Invalid file format: Element type not reckognized: " + inputFile); } output.SpatialDefinition = elementSet; // Read times and values. First element is time, following are element values. ITimeSpaceValueSet valueSet = output.Values; while ((line = GetNextLine(reader)) != null) { string[] fileValues = line.Split(';'); if (fileValues.Length - 1 != elementSet.ElementCount) { throw new InvalidDataException("Number of data does not match number of elements on line: \n" + line); } DateTime timestamp = DateTime.ParseExact(fileValues[0].Trim('"'), "yyyy-MM-dd HH:mm", CultureInfo.InvariantCulture); Time time = new Time(timestamp); double[] locationValues = new double[fileValues.Length - 1]; for (int i = 1; i < fileValues.Length; i++) { locationValues[i - 1] = Double.Parse(fileValues[i].Trim('"'), NumberFormatInfo.InvariantInfo); } timeSet.Times.Add(time); valueSet.Values2D.Add(locationValues); } reader.Close(); timeSet.SetTimeHorizonFromTimes(); _outputs.Add(output); }
public void CouplingGwRiver() { /// bit 1: Decides whether the timeInterpolator or grid-to-line adaptor comes first /// bit 2: when true, using a 16x16 gw grid (instead of 2x2) /// bit 3: when true, bi-directinal: adds a link from gwModel to river, with the gw-level for (int runNumber = 0; runNumber < 8; runNumber++) { //if (runNumber != 3) // continue; Console.Out.WriteLine("runNumber: " + runNumber); // Create trigger inputs Input queryDischargeItem = CreateDischargeInput(); Input queryVolume = CreateVolumeInput(); // Create models LinkableEngine riverModel = CreateRiverModel(); ITimeSpaceComponent gwModel = CreateGwModel(); // Add arguments and initialize IDictionary <string, IArgument> gwArgs = gwModel.Arguments.Dictionary(); // Increasing model grid size (otherwise GW model runs full too fast) gwArgs.UpdateValue("dx", 400.0); gwArgs.UpdateValue("dy", 400.0); gwArgs.UpdateValue("x0", 200.0); gwArgs.UpdateValue("y0", 200.0); int gwGridSize = 2 * 2; if ((runNumber & 2) == 2) // set 16 x 16 grid { gwArgs.UpdateValue("dx", 50.0); gwArgs.UpdateValue("dy", 50.0); gwArgs.UpdateValue("XCount", 16); gwArgs.UpdateValue("ycount", 16); gwGridSize = 16 * 16; } gwModel.Initialize(); IDictionary <string, IArgument> riverArgs = riverModel.Arguments.Dictionary(); // Increasing model grid size (otherwise GW model runs full too fast) riverArgs.UpdateValue("xyscale", 100.0); riverModel.Initialize(); // Connect triggering inputs ITimeSpaceOutput flowOnBranch = UTHelper.FindOutputItem(riverModel, "Branch:2:Flow"); TimeInterpolator flowOnBranch2 = new TimeInterpolator(flowOnBranch); flowOnBranch.AddAdaptedOutput(flowOnBranch2); flowOnBranch2.AddConsumer(queryDischargeItem); ITimeSpaceOutput storageInGw = UTHelper.FindOutputItem(gwModel, "Grid.Storage"); TimeInterpolator storageInGw2 = new TimeInterpolator(storageInGw); storageInGw.AddAdaptedOutput(storageInGw2); storageInGw2.AddConsumer(queryVolume); //========== Couple leakage items ========== // put leakage from river into ground water model { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel, "WholeRiver:Leakage"); ITimeSpaceInput gwInflowInput = UTHelper.FindInputItem(gwModel, "Grid.Inflow"); // Two adaptors are added: Time buffer and line-to-grid adaptor // they can be added in any order (though time buffer first will use less memory) if ((runNumber & 1) == 1) { // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); riverLeakageOutputGrid.AddConsumer(gwInflowInput); } else { // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput, gwInflowInput.ElementSet()); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutputGrid); // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutputGrid); riverLeakageOutputGrid.AddAdaptedOutput(riverLeakageOutput2); riverLeakageOutput2.AddConsumer(gwInflowInput); } } //========== Couple ground water level items ========== if ((runNumber & 4) == 4) { // put ground water level from ground water model into river ITimeSpaceInput riverGwleveInput = UTHelper.FindInputItem(riverModel, "WholeRiver:GroundWaterLevel"); ITimeSpaceOutput gwLevelOutput = UTHelper.FindOutputItem(gwModel, "Grid.gwLevel"); // Two adaptors are added: Time buffer and grid-to-line adaptor // they can be added in any order (though time buffer last will use less memory) if ((runNumber & 1) == 1) { // Time interpolator var gwLevelOutput2 = new TimeExtrapolator(gwLevelOutput); gwLevelOutput.AddAdaptedOutput(gwLevelOutput2); // Element mapper from polyline to polygon, weighted sum version var gwLevelOutputLine = new ElementMapperAdaptedOutput(new Identifier("ElementMapper700"), gwLevelOutput2, riverGwleveInput.ElementSet()); gwLevelOutput2.AddAdaptedOutput(gwLevelOutputLine); gwLevelOutputLine.AddConsumer(riverGwleveInput); } else { // Element mapper from polyline to polygon, weighted sum version var gwLevelOutputLine = new ElementMapperAdaptedOutput(new Identifier("ElementMapper700"), gwLevelOutput, riverGwleveInput.ElementSet()); gwLevelOutput.AddAdaptedOutput(gwLevelOutputLine); // Time interpolator var gwLevelOutput2 = new TimeExtrapolator(gwLevelOutputLine); gwLevelOutputLine.AddAdaptedOutput(gwLevelOutput2); gwLevelOutput2.AddConsumer(riverGwleveInput); } } //========== Run ========== // Validate riverModel.Validate(); Assert.IsTrue(riverModel.Status == LinkableComponentStatus.Valid); gwModel.Validate(); Assert.IsTrue(gwModel.Status == LinkableComponentStatus.Valid); // Prepare riverModel.Prepare(); Assert.IsTrue(riverModel.Status == LinkableComponentStatus.Updated); gwModel.Prepare(); Assert.IsTrue(gwModel.Status == LinkableComponentStatus.Updated); // specify query times double triggerTime0 = riverModel.CurrentTime.StampAsModifiedJulianDay; double triggerTime1 = triggerTime0 + 1; double triggerTime2 = triggerTime0 + 2; double triggerTime3 = triggerTime0 + 12.1; double triggerTime4 = triggerTime0 + 16.7; /// Properties of the river, without gw-level input /// Inflow into each node from rainfall runoff is 10 L/s /// Inflow to node 1: 10 L/s - leaking 5 L/s on branch 1 /// Inflow to node 2: 10 + 5 L/s - leaking 15/2 L/s on branch 2 /// Inflow to node 3: 10 + 15/2 L/s - leaking 35/4 L/s on branch 3 /// Total leakage 5+15/2+35/4 = (20+30+35)/4 = 85/4 L/s /// /// Number of seconds in a day: 60*60*24 = 86400 // check initial values Assert.AreEqual(1, ValueSet.GetElementCount(flowOnBranch.Values), "#values for " + flowOnBranch.Id); Assert.AreEqual(7.0, (double)flowOnBranch.Values.GetValue(0, 0), "Value[0] as property"); Assert.AreEqual(gwGridSize, ValueSet.GetElementCount(storageInGw.Values), "#values for " + storageInGw.Id); Assert.AreEqual(0, SumTimeStep(storageInGw.Values, 0)); // get values for specified query times, 1 days // Totally leaking: 86400 * 85/4 = 1.836e6 // For the bi-directional coupling: // the entire first day the river uses extrapolated values from the // gwModel, which gives a gwLevel of -10, hence same value as for the uni-directional queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime1); ITimeSpaceValueSet valuesV = storageInGw2.GetValues(queryDischargeItem); ITimeSpaceValueSet valuesQ = flowOnBranch2.GetValues(queryDischargeItem); Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(1.836e6, SumTimeStep(valuesV, 0), 1e-4); // Print out, to load in a plotting program for verification //StringBuilder b = new StringBuilder(); //foreach (double val in valuesV.GetElementValuesForTime(0)) // b.AppendLine(val.ToString(NumberFormatInfo.InvariantInfo)); //Console.Out.WriteLine(b.ToString()); // get values for specified query times, 2 days // 2 * 86400 * 85/4 = 3.672e6 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime2); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); if ((runNumber & 4) != 4) // unidirectional { Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(3.672e6, SumTimeStep(valuesV, 0), 1e-4); } else if ((runNumber & 2) != 2) // bi-directional 2x2 grid { Assert.AreEqual(8.843648, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(3.66390879366e6, SumTimeStep(valuesV, 0), 1e-4); } else // bi-directional 16x16 grid { Assert.AreEqual(9.65307, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(3.59397465219e6, SumTimeStep(valuesV, 0), 1e-4); } // get values for specified query times, 12.1 days // 12.1 * 86400 * 85/4 = 2.22156e7 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime3); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); if ((runNumber & 4) != 4) // unidirectional { Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(2.22156e7, SumTimeStep(valuesV, 0), 1e-4); } else if ((runNumber & 2) != 2) // bi-directional 2x2 grid { Assert.AreEqual(9.87828, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(2.16704019338e7, SumTimeStep(valuesV, 0), 1e-4); } else // bi-directional 16x16 grid { Assert.AreEqual(18.546999, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(1.722400002557e7, SumTimeStep(valuesV, 0), 1e-4); } // get values for specified query times, 16.7 days // 16.7 * 86400 * 85/4 = 3.06612e7 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime4); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); if ((runNumber & 4) != 4) // unidirectional { Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(3.06612e7, SumTimeStep(valuesV, 0), 1e-4); } else if ((runNumber & 2) != 2) // bi-directional 2x2 grid { Assert.AreEqual(10.255535, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(2.9595872035072e7, SumTimeStep(valuesV, 0), 1e-4); } else // bi-directional 16x16 grid { Assert.AreEqual(20.98699, (double)valuesQ.GetValue(0, 0), 1e-4); Assert.AreEqual(2.12991179998e7, SumTimeStep(valuesV, 0), 1e-4); } } }
public void CouplingGwRiver2() { /// runNumber 0: Using MultiInput /// runNumber 1: Using MultiInputAdaptor /// runNumber 2: Using MultiInputAdaptorFactory for (int runNumber = 0; runNumber < 3; runNumber++) { Console.Out.WriteLine("runNumber: " + runNumber); // Create trigger inputs Input queryDischargeItem = CreateDischargeInput(); Input queryVolume = CreateVolumeInput(); // Create models LinkableEngine riverModel = CreateRiverModel(); LinkableEngine riverModel2 = CreateRiverModel(); ITimeSpaceComponent gwModel = CreateGwModel(); // Add arguments and initialize IDictionary <string, IArgument> gwArgs = gwModel.Arguments.Dictionary(); // Increasing model grid size (otherwise GW model runs full too fast) gwArgs.UpdateValue("dx", 50.0); gwArgs.UpdateValue("dy", 50.0); gwArgs.UpdateValue("x0", 0.0); gwArgs.UpdateValue("y0", 200.0); gwArgs.UpdateValue("XCount", 24); gwArgs.UpdateValue("ycount", 16); if (runNumber == 0) { gwArgs.UpdateValue("UseMultiInput", true); } gwModel.Initialize(); int gwGridSize = 24 * 16; IDictionary <string, IArgument> riverArgs = riverModel.Arguments.Dictionary(); // Increasing model grid size (otherwise GW model runs full too fast) riverArgs.UpdateValue("xyscale", 100.0); riverModel.Initialize(); IDictionary <string, IArgument> river2Args = riverModel2.Arguments.Dictionary(); // Increasing model grid size (otherwise GW model runs full too fast) river2Args.UpdateValue("xyscale", 100.0); // Move river2 sligthly away from river1 river2Args.UpdateValue("xoffset", -220.0); river2Args.UpdateValue("yoffset", 180.0); riverModel2.Initialize(); // Connect triggering inputs ITimeSpaceOutput flowOnBranch = UTHelper.FindOutputItem(riverModel, "Branch:2:Flow"); TimeInterpolator flowOnBranch2 = new TimeInterpolator(flowOnBranch); flowOnBranch.AddAdaptedOutput(flowOnBranch2); flowOnBranch2.AddConsumer(queryDischargeItem); ITimeSpaceOutput storageInGw = UTHelper.FindOutputItem(gwModel, "Grid.Storage"); TimeInterpolator storageInGw2 = new TimeInterpolator(storageInGw); storageInGw.AddAdaptedOutput(storageInGw2); storageInGw2.AddConsumer(queryVolume); //========== Couple leakage items ========== ITimeSpaceInput gwInflowInput = UTHelper.FindInputItem(gwModel, "Grid.Inflow"); //========== IBaseMultiInput linking ========== if (runNumber == 0) { /// Example of adding up two outputs into one input, by the use of /// an IBaseMultiInput implementation Assert.IsTrue(gwInflowInput is IBaseMultiInput); Assert.IsTrue(gwInflowInput is ITimeSpaceMultiInput); // put leakage from river1 into ground water model { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel, "WholeRiver:Leakage"); // Two adaptors are added: Time buffer and line-to-grid adaptor // they can be added in any order (though time buffer first will use less memory) // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: No special action riverLeakageOutputGrid.AddConsumer(gwInflowInput); } // put leakage from river2 into ground water model { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel2, "WholeRiver:Leakage"); // Two adaptors are added: Time buffer and line-to-grid adaptor // they can be added in any order (though time buffer first will use less memory) // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: No special action riverLeakageOutputGrid.AddConsumer(gwInflowInput); } } //========== MultiInputAdaptor linking ========== if (runNumber == 1) { /// Example of adding up two outputs into one input, by the use of /// a MultiInputAdaptor class // Note !!!: Creating a MultiInputAdaptor MultiInputAdaptor sourceAdder = new MultiInputAdaptor("SomeId") { SpatialDefinition = gwInflowInput.SpatialDefinition }; // put leakage from river1 into ground water model // Two adaptors are added: Time buffer and line-to-grid adaptor { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel, "WholeRiver:Leakage"); // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: Adding to the list of adaptees sourceAdder.Adaptees.Add(riverLeakageOutputGrid); riverLeakageOutputGrid.AddAdaptedOutput(sourceAdder); } // put leakage from river2 into ground water model // Two adaptors are added: Time buffer and line-to-grid adaptor { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel2, "WholeRiver:Leakage"); // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: Adding to the list of adaptees sourceAdder.Adaptees.Add(riverLeakageOutputGrid); riverLeakageOutputGrid.AddAdaptedOutput(sourceAdder); } // Note !!!: Connect the gwInflowInput and the multiInputAdaptor sourceAdder.AddConsumer(gwInflowInput); } //========== MultiInputAdaptorFactory linking ========== if (runNumber == 2) { /// Example of adding up two outputs into one input, by the use of /// an MultiInputAdaptorFactory implementation var factory = new MultiInputAdaptorFactory(gwModel); // put leakage from river1 into ground water model // Two adaptors are added: Time buffer and line-to-grid adaptor { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel, "WholeRiver:Leakage"); // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: Creating a new AdaptedOutput and adding it IIdentifiable[] identifiables = factory.GetAvailableAdaptedOutputIds(riverLeakageOutputGrid, gwInflowInput); IBaseAdaptedOutput myOutput = factory.CreateAdaptedOutput(identifiables[0], riverLeakageOutputGrid, gwInflowInput); myOutput.AddConsumer(gwInflowInput); } // put leakage from river2 into ground water model // Two adaptors are added: Time buffer and line-to-grid adaptor { ITimeSpaceOutput riverLeakageOutput = UTHelper.FindOutputItem(riverModel2, "WholeRiver:Leakage"); // Time interpolator TimeInterpolator riverLeakageOutput2 = new TimeInterpolator(riverLeakageOutput); riverLeakageOutput.AddAdaptedOutput(riverLeakageOutput2); // Element mapper from polyline to polygon, weighted sum version ElementMapperAdaptedOutput riverLeakageOutputGrid = new ElementMapperAdaptedOutput(new Identifier("ElementMapper501"), riverLeakageOutput2, gwInflowInput.ElementSet()); riverLeakageOutput2.AddAdaptedOutput(riverLeakageOutputGrid); // Note !!!: Creating a new AdaptedOutput and adding it IIdentifiable[] identifiables = factory.GetAvailableAdaptedOutputIds(riverLeakageOutputGrid, gwInflowInput); IBaseAdaptedOutput myOutput = factory.CreateAdaptedOutput(identifiables[0], riverLeakageOutputGrid, gwInflowInput); myOutput.AddConsumer(gwInflowInput); } } //========== Run ========== // Validate riverModel.Validate(); Assert.IsTrue(riverModel.Status == LinkableComponentStatus.Valid); riverModel2.Validate(); Assert.IsTrue(riverModel2.Status == LinkableComponentStatus.Valid); gwModel.Validate(); Assert.IsTrue(gwModel.Status == LinkableComponentStatus.Valid); // Prepare riverModel.Prepare(); Assert.IsTrue(riverModel.Status == LinkableComponentStatus.Updated); riverModel2.Prepare(); Assert.IsTrue(riverModel2.Status == LinkableComponentStatus.Updated); gwModel.Prepare(); Assert.IsTrue(gwModel.Status == LinkableComponentStatus.Updated); // specify query times double triggerTime0 = riverModel.CurrentTime.StampAsModifiedJulianDay; double triggerTime1 = triggerTime0 + 1; double triggerTime2 = triggerTime0 + 2; double triggerTime3 = triggerTime0 + 12.1; double triggerTime4 = triggerTime0 + 16.7; /// Properties of the river, without gw-level input /// Inflow into each node from rainfall runoff is 10 L/s /// Inflow to node 1: 10 L/s - leaking 5 L/s on branch 1 /// Inflow to node 2: 10 + 5 L/s - leaking 15/2 L/s on branch 2 /// Inflow to node 3: 10 + 15/2 L/s - leaking 35/4 L/s on branch 3 /// Total leakage 5+15/2+35/4 = (20+30+35)/4 = 85/4 L/s /// /// Number of seconds in a day: 60*60*24 = 86400 // check initial values Assert.AreEqual(1, ValueSet.GetElementCount(flowOnBranch.Values), "#values for " + flowOnBranch.Id); Assert.AreEqual(7.0, (double)flowOnBranch.Values.GetValue(0, 0), "Value[0] as property"); Assert.AreEqual(gwGridSize, ValueSet.GetElementCount(storageInGw.Values), "#values for " + storageInGw.Id); Assert.AreEqual(0, SumTimeStep(storageInGw.Values, 0)); // get values for specified query times, 1 days // Totally leaking: 86400 * 85/4 = 1.836e6 // For the bi-directional coupling: // the entire first day the river uses extrapolated values from the // gwModel, which gives a gwLevel of -10, hence same value as for the uni-directional queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime1); ITimeSpaceValueSet valuesV = storageInGw2.GetValues(queryDischargeItem); ITimeSpaceValueSet valuesQ = flowOnBranch2.GetValues(queryDischargeItem); Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(2 * 1.836e6, SumTimeStep(valuesV, 0), 1e-4); // Print out, to load in a plotting program for verification StringBuilder b = new StringBuilder(); IList valV = valuesV.GetElementValuesForTime(0); int ivalvV = 0; for (int i = 0; i < 16; i++) { for (int j = 0; j < 24; j++) { b.Append(((double)valV[ivalvV++]).ToString(NumberFormatInfo.InvariantInfo)); b.Append(" "); } b.AppendLine(); } //Console.Out.WriteLine(b.ToString()); // get values for specified query times, 2 days // 2 * 86400 * 85/4 = 3.672e6 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime2); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(2 * 3.672e6, SumTimeStep(valuesV, 0), 1e-4); // get values for specified query times, 12.1 days // 12.1 * 86400 * 85/4 = 2.22156e7 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime3); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(2 * 2.22156e7, SumTimeStep(valuesV, 0), 1e-4); // get values for specified query times, 16.7 days // 16.7 * 86400 * 85/4 = 3.06612e7 queryDischargeItem.TimeSet.SetSingleTimeStamp(triggerTime4); valuesV = storageInGw2.GetValues(queryDischargeItem); valuesQ = flowOnBranch2.GetValues(queryDischargeItem); Assert.AreEqual(35.0 / 4.0, (double)valuesQ.GetValue(0, 0)); Assert.AreEqual(2 * 3.06612e7, SumTimeStep(valuesV, 0), 1e-4); } }
public void GetValues() { IQuantity riverModelOutflowQuantity = simpleRiver.Outputs[0].ValueDefinition as IQuantity; IElementSet riverModelNodeElement = ((ITimeSpaceOutput)simpleRiver.Outputs[0]).ElementSet(); Input queryItem1; var dischargeQuantity = new Quantity(new Unit(PredefinedUnits.CubicMeterPerSecond), null, "Discharge"); var waterlevelQuantity = new Quantity(new Unit(PredefinedUnits.Meter), null, "Water Level"); var idBasedElementSetA = new ElementSet(null, "ElmSet-A", ElementType.IdBased); idBasedElementSetA.AddElement(new Element("elm-1")); queryItem1 = new Input("discharge, to be retrieved from some output item", dischargeQuantity, idBasedElementSetA); queryItem1.TimeSet = new TimeSet(); // Connect the queryItem with the first output item using a time interpolator TimeBufferer buffer = new TimeInterpolator((ITimeSpaceOutput)simpleRiver.Outputs[0]); simpleRiver.Outputs[0].AddAdaptedOutput(buffer); buffer.AddConsumer(queryItem1); simpleRiver.Validate(); simpleRiver.Prepare(); //IListener eventListener = new EventListener(); //for (int i = 0; i < eventListener.GetAcceptedEventTypeCount(); i++) //{ // for (int n = 0; n < simpleRiver.GetPublishedEventTypeCount(); n++) // { // if (eventListener.GetAcceptedEventType(i) == simpleRiver.GetPublishedEventType(n)) // { // simpleRiver.Subscribe(eventListener, eventListener.GetAcceptedEventType(i)); // } // } //} double startTime = simpleRiver.TimeExtent.TimeHorizon.StampAsModifiedJulianDay; ITime[] times = new ITime[5]; times[0] = new Time(startTime + 1); times[1] = new Time(startTime + 2); times[2] = new Time(startTime + 3); times[3] = new Time(startTime + 4); times[4] = new Time(startTime + 5); ITimeSpaceValueSet values; queryItem1.TimeSet.SetSingleTimeStamp(times[0].StampAsModifiedJulianDay); values = ((ITimeSpaceOutput)queryItem1.Provider).GetValues(queryItem1); Assert.AreEqual(0.55, values.GetValue(0, 0)); queryItem1.TimeSet.SetSingleTimeStamp(times[1].StampAsModifiedJulianDay); values = ((ITimeSpaceOutput)queryItem1.Provider).GetValues(queryItem1); Assert.AreEqual(1.1, values.GetValue(0, 0)); queryItem1.TimeSet.SetSingleTimeStamp(times[2].StampAsModifiedJulianDay); values = ((ITimeSpaceOutput)queryItem1.Provider).GetValues(queryItem1); Assert.AreEqual(1.6, values.GetValue(0, 0)); queryItem1.TimeSet.SetSingleTimeStamp(times[3].StampAsModifiedJulianDay); values = ((ITimeSpaceOutput)queryItem1.Provider).GetValues(queryItem1); Assert.AreEqual(2.1, values.GetValue(0, 0)); queryItem1.TimeSet.SetSingleTimeStamp(times[4].StampAsModifiedJulianDay); values = ((ITimeSpaceOutput)queryItem1.Provider).GetValues(queryItem1); Assert.AreEqual(2.6, values.GetValue(0, 0)); }