public void TestEvaporation()
{
MockSoil soil = new MockSoil();
APSIM.Shared.Soils.Soil soilProperties = Setup();
APSIMReadySoil.Create(soilProperties);
MockClock clock = new MockClock();
clock.Today = new DateTime(2015, 6, 1);
MockWeather weather = new MockWeather();
weather.MaxT = 30;
weather.MinT = 10;
weather.Rain = 100;
weather.Radn = 25;
MockSurfaceOrganicMatter surfaceOrganicMatter = new MockSurfaceOrganicMatter();
surfaceOrganicMatter.Cover = 0.8;
List<ICanopy> canopies = new List<ICanopy>();
EvaporationModel evaporation = new EvaporationModel();
SetLink(soil, "properties", soilProperties);
SetLink(evaporation, "soil", soil);
SetLink(evaporation, "clock", clock);
SetLink(evaporation, "weather", weather);
SetLink(evaporation, "canopies", canopies);
SetLink(evaporation, "surfaceOrganicMatter", surfaceOrganicMatter);
// Empty profile.
soil.Water = MathUtilities.Multiply(soilProperties.Water.LL15, soilProperties.Water.Thickness);
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 3.00359));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 2.20072));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 1.57064));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.96006));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.75946));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.64851));
soil.Infiltration = 100;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 3.00359));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 2.20072));
}
public void TestEvaporation()
{
MockSoil soil = new MockSoil();
APSIM.Shared.Soils.Soil soilProperties = Setup();
APSIMReadySoil.Create(soilProperties);
MockClock clock = new MockClock();
clock.Today = new DateTime(2015, 6, 1);
MockWeather weather = new MockWeather();
weather.MaxT = 30;
weather.MinT = 10;
weather.Rain = 100;
weather.Radn = 25;
MockSurfaceOrganicMatter surfaceOrganicMatter = new MockSurfaceOrganicMatter();
surfaceOrganicMatter.Cover = 0.8;
List <ICanopy> canopies = new List <ICanopy>();
EvaporationModel evaporation = new EvaporationModel();
SetLink(soil, "properties", soilProperties);
SetLink(evaporation, "soil", soil);
SetLink(evaporation, "clock", clock);
SetLink(evaporation, "weather", weather);
SetLink(evaporation, "canopies", canopies);
SetLink(evaporation, "surfaceOrganicMatter", surfaceOrganicMatter);
// Empty profile.
soil.Water = MathUtilities.Multiply(soilProperties.Water.LL15, soilProperties.Water.Thickness);
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 3.00359));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 2.20072));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 1.57064));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.96006));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.75946));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 0.64851));
soil.Infiltration = 100;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 3.00359));
soil.Infiltration = 0;
evaporation.Calculate();
Assert.IsTrue(MathUtilities.FloatsAreEqual(evaporation.Es, 2.20072));
}