public void GetValues2B()
{
// This is a variation of GetValues2A. In this test the timeSeries is linked ID based to the
// top node of the upperRiver. The last upperRiver branch of the upperRiver is ID based connected
// to the top node of the lowerRiver. The last branch in the lowerRiver is linked to the trigger.
// The timeSeries provides data that changes over time. This is what makes this test different from
// GetValues2A, where everytning is the same for every time step.
try
{
ILinkableComponent timeSeries = new TimeSeriesComponent();
ILinkableComponent upperRiver = new RiverModelLC();
ILinkableComponent lowerRiver = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
timeSeries.Initialize(new Argument[0]);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] upperRiverArguments = new Argument[1];
upperRiverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
upperRiver.Initialize(upperRiverArguments);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] lowerRiverArguments = new Argument[1];
lowerRiverArguments[0] = new Argument("ModelID","lowerRiverModel",true,"argument");
lowerRiver.Initialize(lowerRiverArguments);
trigger.Initialize(new Argument[0]);
Assert.AreEqual("upperRiverModel",upperRiver.ModelID);
Assert.AreEqual("lowerRiverModel",lowerRiver.ModelID);
Link timeSeriesToUpperRiverLink = new Link();
timeSeriesToUpperRiverLink.ID = "timeSeriesToUpperRiverLink";
timeSeriesToUpperRiverLink.SourceComponent = timeSeries;
timeSeriesToUpperRiverLink.SourceElementSet = timeSeries.GetOutputExchangeItem(0).ElementSet; //last branch in the river
timeSeriesToUpperRiverLink.SourceQuantity = timeSeries.GetOutputExchangeItem(0).Quantity;
timeSeriesToUpperRiverLink.TargetComponent = upperRiver;
timeSeriesToUpperRiverLink.TargetElementSet = upperRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
timeSeriesToUpperRiverLink.TargetQuantity = upperRiver.GetInputExchangeItem(0).Quantity;
Link link = new Link();
link.ID = "RiverToRiverLink";
link.SourceComponent = upperRiver;
link.SourceElementSet = upperRiver.GetOutputExchangeItem(2).ElementSet; //last branch in the river
link.SourceQuantity = upperRiver.GetOutputExchangeItem(2).Quantity;
link.TargetComponent = lowerRiver;
link.TargetElementSet = lowerRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
link.TargetQuantity = lowerRiver.GetInputExchangeItem(0).Quantity;
Link triggerLink = new Link();
triggerLink.ID = "TriggerLink";
triggerLink.SourceComponent = lowerRiver;
triggerLink.SourceElementSet = lowerRiver.GetOutputExchangeItem(2).ElementSet;
triggerLink.SourceQuantity = lowerRiver.GetOutputExchangeItem(2).Quantity;
triggerLink.TargetComponent = trigger;
triggerLink.TargetElementSet = trigger.GetInputExchangeItem(0).ElementSet;
triggerLink.TargetQuantity = trigger.GetInputExchangeItem(0).Quantity;
timeSeries.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(link);
lowerRiver.AddLink(link);
lowerRiver.AddLink(triggerLink);
trigger.AddLink(triggerLink);
timeSeries.Prepare();
upperRiver.Prepare();
lowerRiver.Prepare();
trigger.Prepare();
double firstTriggerGetValuesTime = lowerRiver.TimeHorizon.Start.ModifiedJulianDay;
TimeStamp[] triggerTimes = new TimeStamp[2];
triggerTimes[0] = new TimeStamp(firstTriggerGetValuesTime + 12.5);
triggerTimes[1] = new TimeStamp(firstTriggerGetValuesTime + 16.2);
trigger.Run(triggerTimes);
double x1 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(0)).GetScalar(0);
Assert.AreEqual(315.0/32.0 + 13.0/64.0,x1);
double x2 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(1)).GetScalar(0);
Assert.AreEqual(315.0/32.0 + 17.0/64.0,x2);
upperRiver.Finish();
lowerRiver.Finish();
upperRiver.Dispose();
lowerRiver.Dispose();
}
catch (System.Exception e)
{
ExceptionHandler.WriteException(e);
throw (e);
}
}
// [Ignore ("This test fails when dt = 1 hour but works for dt = 0.5 day")]
public void GetValues2C()
{
// This test is a variation of GetValues 2B. What makes this test different is that the the two
// connected river models are not using the same timesteps.
// This test will test:
// - is the buffer working correcly with respect to interpolations and buffering
// - is the buffer working correctly with respect to clearing the buffer.
try
{
ILinkableComponent timeSeries = new TimeSeriesComponent();
ILinkableComponent upperRiver = new RiverModelLC();
ILinkableComponent lowerRiver = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
timeSeries.Initialize(new Argument[0]);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] upperRiverArguments = new Argument[2];
upperRiverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
upperRiverArguments[1] = new Argument("TimeStepLength","21600",true,"A time step length of 1 day");
upperRiver.Initialize(upperRiverArguments);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] lowerRiverArguments = new Argument[2];
lowerRiverArguments[0] = new Argument("ModelID","lowerRiverModel",true,"argument");
lowerRiverArguments[1] = new Argument("TimeStepLength","86400",true,"xx");
lowerRiver.Initialize(lowerRiverArguments);
trigger.Initialize(new Argument[0]);
Assert.AreEqual("upperRiverModel",upperRiver.ModelID);
Assert.AreEqual("lowerRiverModel",lowerRiver.ModelID);
Link timeSeriesToUpperRiverLink = new Link();
timeSeriesToUpperRiverLink.ID = "timeSeriesToUpperRiverLink";
timeSeriesToUpperRiverLink.SourceComponent = timeSeries;
timeSeriesToUpperRiverLink.SourceElementSet = timeSeries.GetOutputExchangeItem(0).ElementSet; //last branch in the river
timeSeriesToUpperRiverLink.SourceQuantity = timeSeries.GetOutputExchangeItem(0).Quantity;
timeSeriesToUpperRiverLink.TargetComponent = upperRiver;
timeSeriesToUpperRiverLink.TargetElementSet = upperRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
timeSeriesToUpperRiverLink.TargetQuantity = upperRiver.GetInputExchangeItem(0).Quantity;
Link link = new Link();
link.ID = "RiverToRiverLink";
link.SourceComponent = upperRiver;
link.SourceElementSet = upperRiver.GetOutputExchangeItem(2).ElementSet; //last branch in the river
link.SourceQuantity = upperRiver.GetOutputExchangeItem(2).Quantity;
link.TargetComponent = lowerRiver;
link.TargetElementSet = lowerRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
link.TargetQuantity = lowerRiver.GetInputExchangeItem(0).Quantity;
Link triggerLink = new Link();
triggerLink.ID = "TriggerLink";
triggerLink.SourceComponent = lowerRiver;
triggerLink.SourceElementSet = lowerRiver.GetOutputExchangeItem(2).ElementSet;
triggerLink.SourceQuantity = lowerRiver.GetOutputExchangeItem(2).Quantity;
triggerLink.TargetComponent = trigger;
triggerLink.TargetElementSet = trigger.GetInputExchangeItem(0).ElementSet;
triggerLink.TargetQuantity = trigger.GetInputExchangeItem(0).Quantity;
timeSeries.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(link);
lowerRiver.AddLink(link);
lowerRiver.AddLink(triggerLink);
trigger.AddLink(triggerLink);
timeSeries.Prepare();
upperRiver.Prepare();
lowerRiver.Prepare();
trigger.Prepare();
double firstTriggerGetValuesTime = lowerRiver.TimeHorizon.Start.ModifiedJulianDay;
TimeStamp[] triggerTimes = new TimeStamp[2];
triggerTimes[0] = new TimeStamp(firstTriggerGetValuesTime + 12.5);
triggerTimes[1] = new TimeStamp(firstTriggerGetValuesTime + 16.2);
trigger.Run(triggerTimes);
double x1 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(0)).GetScalar(0);
Assert.AreEqual(315.0/32.0 + 13.0/64.0,x1);
double x2 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(1)).GetScalar(0);
Assert.AreEqual(315.0/32.0 + 17.0/64.0,x2);
Assert.AreEqual(10,((RiverModelLC) upperRiver)._maxBufferSize); //Test that the buffer is cleared
upperRiver.Finish();
lowerRiver.Finish();
upperRiver.Dispose();
lowerRiver.Dispose();
}
catch (System.Exception e)
{
ExceptionHandler.WriteException(e);
throw (e);
}
}
public void GetValues2A()
{
// == Running with two instances of riverModelLC ==
//
// - The Two river are running with the same timestepping.
//
// - The link is ID Based link with flow from last branch of the source river to the top
// node of the target river.
//
// - The time argument in the GetValues from rive to river is of type ITimeSpan
//
//TODO: 1: The RiverModelEngine should change the inflow over time. As it is now the inflow is the same
// in all time steps. Another idea would be to have a output exchange item that hold the accumulated
// inflow, this could be useful when testing the manage state interface.
//
// 2: Make this test run with the two river using different timesteps and with the source river
// starting ealier that the target river.
//
// 3: In this test also events could be tested. Simply test if all the required events are
// thrown during the simulations.
try
{
ILinkableComponent upperRiver = new RiverModelLC();
ILinkableComponent lowerRiver = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] upperRiverArguments = new Argument[1];
upperRiverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
upperRiver.Initialize(upperRiverArguments);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] lowerRiverArguments = new Argument[1];
lowerRiverArguments[0] = new Argument("ModelID","lowerRiverModel",true,"argument");
lowerRiver.Initialize(lowerRiverArguments);
trigger.Initialize(new Argument[0]);
Assert.AreEqual("upperRiverModel",upperRiver.ModelID);
Assert.AreEqual("lowerRiverModel",lowerRiver.ModelID);
Link link = new Link();
link.ID = "RiverToRiverLink";
link.SourceComponent = upperRiver;
link.SourceElementSet = upperRiver.GetOutputExchangeItem(2).ElementSet; //last branch in the river
link.SourceQuantity = upperRiver.GetOutputExchangeItem(2).Quantity;
link.TargetComponent = lowerRiver;
link.TargetElementSet = lowerRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
link.TargetQuantity = lowerRiver.GetInputExchangeItem(0).Quantity;
Link triggerLink = new Link();
triggerLink.ID = "TriggerLink";
triggerLink.SourceComponent = lowerRiver;
triggerLink.SourceElementSet = lowerRiver.GetOutputExchangeItem(2).ElementSet;
triggerLink.SourceQuantity = lowerRiver.GetOutputExchangeItem(2).Quantity;
triggerLink.TargetComponent = trigger;
triggerLink.TargetElementSet = trigger.GetInputExchangeItem(0).ElementSet;
triggerLink.TargetQuantity = trigger.GetInputExchangeItem(0).Quantity;
upperRiver.AddLink(link);
lowerRiver.AddLink(link);
lowerRiver.AddLink(triggerLink);
trigger.AddLink(triggerLink);
upperRiver.Prepare();
lowerRiver.Prepare();
trigger.Prepare();
double firstTriggerGetValuesTime = lowerRiver.TimeHorizon.Start.ModifiedJulianDay;
TimeStamp[] triggerTimes = new TimeStamp[2];
triggerTimes[0] = new TimeStamp(firstTriggerGetValuesTime + 3);
triggerTimes[1] = new TimeStamp(firstTriggerGetValuesTime + 4.3);
trigger.Run(triggerTimes);
double x1 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(0)).GetScalar(0);
Assert.AreEqual(315.0/32.0,x1);
double x2 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(1)).GetScalar(0);
Assert.AreEqual(315.0/32.0,x2);
upperRiver.Finish();
lowerRiver.Finish();
upperRiver.Dispose();
lowerRiver.Dispose();
}
catch (System.Exception e)
{
ExceptionHandler.WriteException(e);
throw (e);
}
}
public void GetValues1B()
{
// Running with one instances of riverModelLC linked ID-based to trigger and to
// an instance of the TimeSeriesComponent.
try
{
ILinkableComponent timeSeries = new TimeSeriesComponent();
ILinkableComponent riverModelLC = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
timeSeries.Initialize(new Argument[0]);
Argument[] riverArguments = new Argument[2];
riverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
riverArguments[1] = new Argument("TimeStepLength","3600",true,"A time step length of 1 day");
riverModelLC.Initialize(riverArguments);
trigger.Initialize(new Argument[0]);
Link timeSeriesToRiverLink = new Link();
timeSeriesToRiverLink.ID = "timeSeriesToUpperRiverLink";
timeSeriesToRiverLink.SourceComponent = timeSeries;
timeSeriesToRiverLink.SourceElementSet = timeSeries.GetOutputExchangeItem(0).ElementSet; //last branch in the river
timeSeriesToRiverLink.SourceQuantity = timeSeries.GetOutputExchangeItem(0).Quantity;
timeSeriesToRiverLink.TargetComponent = riverModelLC;
timeSeriesToRiverLink.TargetElementSet = riverModelLC.GetInputExchangeItem(0).ElementSet; //first node in the river
timeSeriesToRiverLink.TargetQuantity = riverModelLC.GetInputExchangeItem(0).Quantity;
Link link = new Link(riverModelLC, riverModelLC.GetOutputExchangeItem(2).ElementSet,riverModelLC.GetOutputExchangeItem(2).Quantity,trigger,trigger.GetInputExchangeItem(0).ElementSet,trigger.GetInputExchangeItem(0).Quantity,"LinkID");
timeSeries.AddLink(timeSeriesToRiverLink);
riverModelLC.AddLink(timeSeriesToRiverLink);
riverModelLC.AddLink(link);
trigger.AddLink(link);
riverModelLC.Prepare();
double firstTriggerGetValuesTime = riverModelLC.TimeHorizon.Start.ModifiedJulianDay;
TimeStamp[] triggerTimes = new TimeStamp[2];
triggerTimes[0] = new TimeStamp(firstTriggerGetValuesTime + 12.1);
triggerTimes[1] = new TimeStamp(firstTriggerGetValuesTime + 16.7);
trigger.Run(triggerTimes);
double x1 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(0)).GetScalar(0);
Assert.AreEqual(35.0/4.0 + 13.0/8.0,x1);
double x2 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(1)).GetScalar(0);
Assert.AreEqual(35.0/4.0 + 17.0/8.0,x2);
riverModelLC.Finish();
riverModelLC.Dispose();
}
catch (System.Exception e)
{
ExceptionHandler.WriteException(e);
throw (e);
}
}
public void Validate()
{
ILinkableComponent upperRiver = new RiverModelLC();
ILinkableComponent lowerRiver = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] upperRiverArguments = new Argument[1];
upperRiverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
upperRiver.Initialize(upperRiverArguments);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] lowerRiverArguments = new Argument[1];
lowerRiverArguments[0] = new Argument("ModelID","lowerRiverModel",true,"argument");
lowerRiver.Initialize(lowerRiverArguments);
trigger.Initialize(new Argument[0]);
Assert.AreEqual("upperRiverModel",upperRiver.ModelID);
Assert.AreEqual("lowerRiverModel",lowerRiver.ModelID);
Dimension wrongDimension = new Dimension();
wrongDimension.SetPower(DimensionBase.Mass,1);
Quantity wrongQuantity = new Quantity(new Unit("dummy",0.0,0.0,"dummy"),"test","qid",global::OpenMI.Standard.ValueType.Vector,wrongDimension);
ElementSet wrongElementSet = new ElementSet("Wrong ElementSet","BadID",ElementType.XYPolyLine,new SpatialReference("no ref"));
Element element = new Element("dum Element");
element.AddVertex(new Vertex(4,5,0));
wrongElementSet.AddElement(element);
Link link = new Link();
link.ID = "RiverToRiverLink";
link.SourceComponent = upperRiver;
link.SourceElementSet = upperRiver.GetOutputExchangeItem(2).ElementSet; //last branch in the river
link.SourceQuantity = upperRiver.GetOutputExchangeItem(2).Quantity;
link.TargetComponent = lowerRiver;
link.TargetElementSet = lowerRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
link.TargetQuantity = wrongQuantity;
// link.AddDataOperation(upperRiver.GetOutputExchangeItem(6).GetDataOperation(0)); // bad data Operation
Link triggerLink = new Link();
triggerLink.ID = "TriggerLink";
triggerLink.SourceComponent = lowerRiver;
triggerLink.SourceElementSet = wrongElementSet;
triggerLink.SourceQuantity = wrongQuantity;
triggerLink.TargetComponent = trigger;
triggerLink.TargetElementSet = trigger.GetInputExchangeItem(0).ElementSet;
triggerLink.TargetQuantity = trigger.GetInputExchangeItem(0).Quantity;
upperRiver.AddLink(link);
lowerRiver.AddLink(link);
lowerRiver.AddLink(triggerLink);
trigger.AddLink(triggerLink);
bool isSilent = false;
if (!isSilent)
{
Console.WriteLine(lowerRiver.Validate());
foreach (string str in ((RiverModelLC) upperRiver).ValidationErrorMessages)
{
Console.WriteLine("Error upperRiver: " + str);
}
foreach (string str in ((RiverModelLC) lowerRiver).ValidationErrorMessages)
{
Console.WriteLine("Error lowerRiver: " + str);
}
foreach (string str in ((RiverModelLC) upperRiver).ValidationWarningMessages)
{
Console.WriteLine("Warning upperRiver: " + str);
}
foreach (string str in ((RiverModelLC) lowerRiver).ValidationWarningMessages)
{
Console.WriteLine("Warning lowerRiver: " + str);
}
}
Assert.AreEqual(0,((RiverModelLC) upperRiver).ValidationErrorMessages.Count);
Assert.AreEqual(4,((RiverModelLC) lowerRiver).ValidationErrorMessages.Count);
Assert.AreEqual(0,((RiverModelLC) upperRiver).ValidationWarningMessages.Count);
Assert.AreEqual(2,((RiverModelLC) lowerRiver).ValidationWarningMessages.Count);
}
public void InputExchangeItem()
{
ILinkableComponent riverModelLC = new RiverModelLC();
riverModelLC.Initialize(new Argument[0]);
IInputExchangeItem exchangeItem0 = riverModelLC.GetInputExchangeItem(0);
Assert.AreEqual("Node:0",exchangeItem0.ElementSet.ID);
Assert.AreEqual("Flow",exchangeItem0.Quantity.ID);
Assert.AreEqual("LiterPrSecond",exchangeItem0.Quantity.Unit.ID);
Assert.AreEqual(0.001,exchangeItem0.Quantity.Unit.ConversionFactorToSI);
Assert.AreEqual(0.0,exchangeItem0.Quantity.Unit.OffSetToSI);
Assert.AreEqual(3,exchangeItem0.Quantity.Dimension.GetPower(DimensionBase.Length));
Assert.AreEqual(-1,exchangeItem0.Quantity.Dimension.GetPower(DimensionBase.Time));
Assert.AreEqual(1,exchangeItem0.ElementSet.ElementCount);
Assert.AreEqual(ElementType.IDBased,exchangeItem0.ElementSet.ElementType);
IInputExchangeItem exchangeItem1 = riverModelLC.GetInputExchangeItem(1);
Assert.AreEqual("Node:1",exchangeItem1.ElementSet.ID);
Assert.AreEqual("Flow",exchangeItem1.Quantity.ID);
Assert.AreEqual(1,exchangeItem1.ElementSet.ElementCount);
Assert.AreEqual(ElementType.IDBased,exchangeItem1.ElementSet.ElementType);
IInputExchangeItem exchangeItem2 = riverModelLC.GetInputExchangeItem(2);
Assert.AreEqual("Node:2",exchangeItem2.ElementSet.ID);
Assert.AreEqual("Flow",exchangeItem2.Quantity.ID);
Assert.AreEqual(1,exchangeItem2.ElementSet.ElementCount);
Assert.AreEqual(ElementType.IDBased,exchangeItem2.ElementSet.ElementType);
IInputExchangeItem exchangeItem3 = riverModelLC.GetInputExchangeItem(3);
Assert.AreEqual("Node:3",exchangeItem3.ElementSet.ID);
Assert.AreEqual("Flow",exchangeItem3.Quantity.ID);
Assert.AreEqual(1,exchangeItem3.ElementSet.ElementCount);
Assert.AreEqual(ElementType.IDBased,exchangeItem3.ElementSet.ElementType);
IInputExchangeItem exchangeItem4 = riverModelLC.GetInputExchangeItem(4);
Assert.AreEqual("WholeRiver",exchangeItem4.ElementSet.ID);
Assert.AreEqual("Flow",exchangeItem4.Quantity.ID);
Assert.AreEqual(3,exchangeItem4.ElementSet.ElementCount);
Assert.AreEqual(ElementType.XYPolyLine,exchangeItem4.ElementSet.ElementType);
Assert.AreEqual(2,exchangeItem4.ElementSet.GetVertexCount(0));
Assert.AreEqual(2,exchangeItem4.ElementSet.GetVertexCount(1));
Assert.AreEqual(2,exchangeItem4.ElementSet.GetVertexCount(2));
Assert.AreEqual(3.0,exchangeItem4.ElementSet.GetXCoordinate(0,0));
Assert.AreEqual(5.0,exchangeItem4.ElementSet.GetXCoordinate(0,1));
Assert.AreEqual(5.0,exchangeItem4.ElementSet.GetXCoordinate(1,0));
Assert.AreEqual(8.0,exchangeItem4.ElementSet.GetXCoordinate(1,1));
Assert.AreEqual(8.0,exchangeItem4.ElementSet.GetXCoordinate(2,0));
Assert.AreEqual(8.0,exchangeItem4.ElementSet.GetXCoordinate(2,1));
Assert.AreEqual(9.0,exchangeItem4.ElementSet.GetYCoordinate(0,0));
Assert.AreEqual(7.0,exchangeItem4.ElementSet.GetYCoordinate(0,1));
Assert.AreEqual(7.0,exchangeItem4.ElementSet.GetYCoordinate(1,0));
Assert.AreEqual(7.0,exchangeItem4.ElementSet.GetYCoordinate(1,1));
Assert.AreEqual(7.0,exchangeItem4.ElementSet.GetYCoordinate(2,0));
Assert.AreEqual(3.0,exchangeItem4.ElementSet.GetYCoordinate(2,1));
}
public void RestoreState()
{
// This test is based on GetValues2C
try
{
ILinkableComponent timeSeries = new TimeSeriesComponent();
ILinkableComponent upperRiver = new RiverModelLC();
ILinkableComponent lowerRiver = new RiverModelLC();
HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger trigger = new HydroNumerics.OpenMI.Sdk.Wrapper.UnitTest.Trigger();
timeSeries.Initialize(new Argument[0]);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] upperRiverArguments = new Argument[2];
upperRiverArguments[0] = new Argument("ModelID","upperRiverModel",true,"argument");
upperRiverArguments[1] = new Argument("TimeStepLength","21600",true,"xx");
upperRiver.Initialize(upperRiverArguments);
// The ModelID is passes in ordet to make it easier to debug, otherwise you cannot se the difference between the two istances of RiverModelLC
Argument[] lowerRiverArguments = new Argument[2];
lowerRiverArguments[0] = new Argument("ModelID","lowerRiverModel",true,"argument");
lowerRiverArguments[1] = new Argument("TimeStepLength","86400",true,"xx");
lowerRiver.Initialize(lowerRiverArguments);
trigger.Initialize(new Argument[0]);
Assert.AreEqual("upperRiverModel",upperRiver.ModelID);
Assert.AreEqual("lowerRiverModel",lowerRiver.ModelID);
Link timeSeriesToUpperRiverLink = new Link();
timeSeriesToUpperRiverLink.ID = "timeSeriesToUpperRiverLink";
timeSeriesToUpperRiverLink.SourceComponent = timeSeries;
timeSeriesToUpperRiverLink.SourceElementSet = timeSeries.GetOutputExchangeItem(0).ElementSet; //last branch in the river
timeSeriesToUpperRiverLink.SourceQuantity = timeSeries.GetOutputExchangeItem(0).Quantity;
timeSeriesToUpperRiverLink.TargetComponent = upperRiver;
timeSeriesToUpperRiverLink.TargetElementSet = upperRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
timeSeriesToUpperRiverLink.TargetQuantity = upperRiver.GetInputExchangeItem(0).Quantity;
Link link = new Link();
link.ID = "RiverToRiverLink";
link.SourceComponent = upperRiver;
link.SourceElementSet = upperRiver.GetOutputExchangeItem(2).ElementSet; //last branch in the river
link.SourceQuantity = upperRiver.GetOutputExchangeItem(2).Quantity;
link.TargetComponent = lowerRiver;
link.TargetElementSet = lowerRiver.GetInputExchangeItem(0).ElementSet; //first node in the river
link.TargetQuantity = lowerRiver.GetInputExchangeItem(0).Quantity;
Link triggerLink = new Link();
triggerLink.ID = "TriggerLink";
triggerLink.SourceComponent = lowerRiver;
triggerLink.SourceElementSet = lowerRiver.GetOutputExchangeItem(2).ElementSet;
triggerLink.SourceQuantity = lowerRiver.GetOutputExchangeItem(2).Quantity;
triggerLink.TargetComponent = trigger;
triggerLink.TargetElementSet = trigger.GetInputExchangeItem(0).ElementSet;
triggerLink.TargetQuantity = trigger.GetInputExchangeItem(0).Quantity;
timeSeries.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(timeSeriesToUpperRiverLink);
upperRiver.AddLink(link);
lowerRiver.AddLink(link);
lowerRiver.AddLink(triggerLink);
trigger.AddLink(triggerLink);
timeSeries.Prepare();
upperRiver.Prepare();
lowerRiver.Prepare();
trigger.Prepare();
// double firstTriggerGetValuesTime = lowerRiver.TimeHorizon.Start.ModifiedJulianDay;
// TimeStamp[] triggerTimes = new TimeStamp[2];
// triggerTimes[0] = new TimeStamp(firstTriggerGetValuesTime + 12.5);
// triggerTimes[1] = new TimeStamp(firstTriggerGetValuesTime + 16.2);
double t = lowerRiver.TimeHorizon.Start.ModifiedJulianDay;
Assert.AreEqual(315.0/32.0 + 13.0/64.0, ((ScalarSet)lowerRiver.GetValues(new TimeStamp(t+12.5),"TriggerLink")).GetScalar(0));
string lowerRiverStateID = ((IManageState) lowerRiver).KeepCurrentState();
string upperRiverStateID = ((IManageState) upperRiver).KeepCurrentState();
Assert.AreEqual(315.0/32.0 + 17.0/64.0, ((ScalarSet)lowerRiver.GetValues(new TimeStamp(t+16.2),"TriggerLink")).GetScalar(0));
((IManageState) lowerRiver).RestoreState(lowerRiverStateID);
((IManageState) upperRiver).RestoreState(upperRiverStateID);
lowerRiver.GetValues(new TimeStamp(t + 14.0),"TriggerLink");
Assert.AreEqual(315.0/32.0 + 17.0/64.0, ((ScalarSet)lowerRiver.GetValues(new TimeStamp(t+16.2),"TriggerLink")).GetScalar(0));
// double x1 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(0)).GetScalar(0);
// Assert.AreEqual(315.0/32.0 + 13.0/64.0,x1);
//
// double x2 = ((IScalarSet)trigger.ResultsBuffer.GetValuesAt(1)).GetScalar(0);
// Assert.AreEqual(315.0/32.0 + 17.0/64.0,x2);
//
// Assert.AreEqual(2,((RiverModelLC) upperRiver)._maxBufferSize); //Test that the buffer is cleared
upperRiver.Finish();
lowerRiver.Finish();
upperRiver.Dispose();
lowerRiver.Dispose();
}
catch (System.Exception e)
{
ExceptionHandler.WriteException(e);
throw (e);
}
}