/// <summary>Add in reset management events to the APSIM spec for the specified paddock.</summary>
/// <param name="paddock">The paddock.</param>
/// <param name="simulation">The simulation to add the management operations to.</param>
/// <exception cref="System.Exception">Cannot find soil water reset date</exception>
/// <exception cref="Exception">Cannot find soil water reset date</exception>
private static void AddResetDatesToManagement(Paddock paddock, APSIMSpecification simulation)
{
// Reset
if (paddock.SoilWaterSampleDate != DateTime.MinValue)
{
if (paddock.SoilNitrogenSampleDate == DateTime.MinValue)
{
paddock.SoilNitrogenSampleDate = paddock.SoilWaterSampleDate;
}
Sow sowing = YieldProphetUtility.GetCropBeingSown(paddock.Management);
if (sowing != null && sowing.Date != DateTime.MinValue)
{
List <Management> resetActions = new List <Management>();
// reset at sowing if the sample dates are after sowing.
if (paddock.SoilWaterSampleDate > sowing.Date)
{
resetActions.Add(new ResetWater()
{
Date = sowing.Date
});
resetActions.Add(new ResetSurfaceOrganicMatter()
{
Date = sowing.Date
});
}
if (paddock.SoilNitrogenSampleDate > sowing.Date)
{
resetActions.Add(new ResetNitrogen()
{
Date = sowing.Date
});
}
// reset on the sample dates.
resetActions.Add(new ResetWater()
{
Date = paddock.SoilWaterSampleDate
});
resetActions.Add(new ResetSurfaceOrganicMatter()
{
Date = paddock.SoilWaterSampleDate
});
resetActions.Add(new ResetNitrogen()
{
Date = paddock.SoilNitrogenSampleDate
});
simulation.Management.InsertRange(0, resetActions);
}
}
}
/// <summary>Converts the old Yield Prophet XML to new XML format capable of deserialisation</summary>
/// <param name="yieldProphetXML">The old Yield Prophet XML</param>
/// <returns>The new Yield Prophet XML</returns>
public static YieldProphet YieldProphetFromXML(XmlNode node, string baseFolder)
{
List <Paddock> simulations = new List <Paddock>();
List <XmlNode> paddocks = XmlUtilities.ChildNodes(node, "Paddock");
for (int p = 0; p < paddocks.Count; p++)
{
try
{
Paddock paddock = CreateSimulationSpec(paddocks[p], baseFolder);
simulations.Add(paddock);
}
catch (Exception err)
{
string name = XmlUtilities.Value(paddocks[p], "Name");
throw new Exception(err.Message + "\r\nPaddock name: " + name);
}
}
YieldProphet simulationsSpec = new YieldProphet
{
Paddock = simulations
};
// Some top level simulation metadata.
string reportDescription = XmlUtilities.Value(node, "ReportDescription");
if (reportDescription != "")
{
simulationsSpec.ReportName = reportDescription;
}
string reportType = XmlUtilities.Value(node, "ReportType");
if (reportType != string.Empty)
{
throw new Exception("Not sure what to do with a ReportType of " + reportType);
}
return(simulationsSpec);
}
/// <summary>Converts the paddock XML.</summary>
/// <param name="paddock">The paddock.</param>
/// <exception cref="System.Exception">Bad paddock name: + name</exception>
private static Paddock CreateSimulationSpec(XmlNode paddock, string baseFolder)
{
Paddock simulation = new Paddock();
string name = XmlUtilities.NameAttr(paddock);
string delimiter = "^";
int posCaret = name.IndexOf(delimiter);
if (posCaret == -1)
{
throw new Exception("Bad paddock name: " + name);
}
string remainder = StringUtilities.SplitOffAfterDelimiter(ref name, delimiter);
string growerName;
string paddockName = StringUtilities.SplitOffAfterDelimiter(ref remainder, delimiter);
if (paddockName == string.Empty)
{
growerName = name;
paddockName = remainder;
}
else
{
growerName = remainder;
}
// Give the paddock a name.
string fullName = string.Format("{0}^{1}^{2}", GetDate(paddock, "SowDateFull").Year, growerName, paddockName);
simulation.Name = fullName;
// Set the report date.
simulation.NowDate = GetDate(paddock.ParentNode, "TodayDateFull");
if (simulation.NowDate == DateTime.MinValue)
{
simulation.NowDate = DateTime.Now;
}
// Store any rainfall data in the simulation.
string rainfallSource = GetString(paddock, "RainfallSource");
if (rainfallSource != "Weather station")
{
string rainFileName = GetString(paddock, "RainfallFilename");
if (rainFileName != string.Empty)
{
string fullFileName = Path.Combine(baseFolder, rainFileName);
if (!File.Exists(fullFileName))
{
throw new Exception("Cannot find file: " + fullFileName);
}
simulation.ObservedData = ApsimTextFile.ToTable(fullFileName);
simulation.ObservedData.TableName = rainfallSource;
if (simulation.ObservedData.Rows.Count == 0)
{
simulation.ObservedData = null;
}
}
}
// Set the reset dates
simulation.SoilWaterSampleDate = GetDate(paddock, "ResetDateFull");
simulation.SoilNitrogenSampleDate = GetDate(paddock, "SoilNitrogenSampleDateFull");
simulation.StationNumber = GetInteger(paddock, "StationNumber");
simulation.StationName = GetString(paddock, "StationName");
// Create a sowing management
Sow sowing = new Sow();
simulation.Management.Add(sowing);
sowing.Date = GetDate(paddock, "SowDateFull");
sowing.EmergenceDate = GetDate(paddock, "EmergenceDateFull");
sowing.Crop = GetString(paddock, "Crop");
sowing.Cultivar = GetString(paddock, "Cultivar");
sowing.SkipRow = GetString(paddock, "SkipRow");
sowing.SowingDensity = GetInteger(paddock, "SowingDensity");
sowing.MaxRootDepth = GetInteger(paddock, "MaxRootDepth") * 10; // cm to mm
sowing.BedWidth = GetInteger(paddock, "BedWidth");
sowing.BedRowSpacing = GetDouble(paddock, "BedRowSpacing");
// Make sure we have a stubbletype
simulation.StubbleType = GetString(paddock, "StubbleType");
if (simulation.StubbleType == string.Empty || simulation.StubbleType == "None")
{
simulation.StubbleType = "Wheat";
}
simulation.StubbleMass = GetDouble(paddock, "StubbleMass");
simulation.Slope = GetDouble(paddock, "Slope");
simulation.SlopeLength = GetDouble(paddock, "SlopeLength");
simulation.UseEC = GetBoolean(paddock, "UseEC");
// Fertilise nodes.
List <XmlNode> fertiliserNodes = XmlUtilities.ChildNodes(paddock, "Fertilise");
for (int f = 0; f < fertiliserNodes.Count; f++)
{
Fertilise fertilise = new Fertilise();
simulation.Management.Add(fertilise);
fertilise.Date = GetDate(fertiliserNodes[f], "FertDateFull");
fertilise.Amount = GetDouble(fertiliserNodes[f], "FertAmount");
}
// Irrigate nodes.
List <XmlNode> irrigateNodes = XmlUtilities.ChildNodes(paddock, "Irrigate");
for (int i = 0; i < irrigateNodes.Count; i++)
{
Irrigate irrigate = new Irrigate();
simulation.Management.Add(irrigate);
irrigate.Date = GetDate(irrigateNodes[i], "IrrigateDateFull");
irrigate.Amount = GetDouble(irrigateNodes[i], "IrrigateAmount");
irrigate.Efficiency = GetDouble(irrigateNodes[i], "IrrigateEfficiency");
}
// Tillage nodes.
foreach (XmlNode tillageNode in XmlUtilities.ChildNodes(paddock, "Tillage"))
{
Tillage tillage = new Tillage();
simulation.Management.Add(tillage);
tillage.Date = GetDate(tillageNode, "TillageDateFull");
string disturbance = GetString(tillageNode, "Disturbance");
if (disturbance == "Low")
{
tillage.Disturbance = Tillage.DisturbanceEnum.Low;
}
else if (disturbance == "Medium")
{
tillage.Disturbance = Tillage.DisturbanceEnum.Medium;
}
else
{
tillage.Disturbance = Tillage.DisturbanceEnum.High;
}
}
// Stubble removed nodes.
foreach (XmlNode stubbleRemovedNode in XmlUtilities.ChildNodes(paddock, "StubbleRemoved"))
{
StubbleRemoved stubbleRemoved = new StubbleRemoved();
simulation.Management.Add(stubbleRemoved);
stubbleRemoved.Date = GetDate(stubbleRemovedNode, "StubbleRemovedDateFull");
stubbleRemoved.Percent = GetDouble(stubbleRemovedNode, "StubbleRemovedAmount");
}
// Look for a soil node.
XmlNode soilNode = XmlUtilities.FindByType(paddock, "Soil");
if (soilNode != null)
{
string testValue = XmlUtilities.Value(soilNode, "Water/Layer/Thickness");
if (testValue != string.Empty)
{
// old format.
soilNode = ConvertSoilNode.Upgrade(soilNode);
}
}
// Fix up soil sample variables.
Sample sample1 = new Sample
{
Name = "Sample1",
Thickness = GetArray(paddock, "Sample1Thickness"),
SW = GetArray(paddock, "SW")
};
if (sample1.SW == null)
{
// Really old way of doing samples - they are stored under soil.
List <XmlNode> sampleNodes = XmlUtilities.ChildNodes(soilNode, "Sample");
if (sampleNodes.Count > 0)
{
sample1 = XmlUtilities.Deserialise(sampleNodes[0], typeof(Sample)) as Sample;
}
}
else
{
sample1.NO3 = GetArray(paddock, "NO3");
sample1.NH4 = GetArray(paddock, "NH4");
}
Sample sample2 = null;
double[] sample2Thickness = GetArray(paddock, "Sample2Thickness");
if (sample2Thickness == null)
{
// Really old way of doing samples - they are stored under soil.
List <XmlNode> sampleNodes = XmlUtilities.ChildNodes(soilNode, "Sample");
if (sampleNodes.Count > 1)
{
sample2 = XmlUtilities.Deserialise(sampleNodes[1], typeof(Sample)) as Sample;
}
}
else
{
sample2 = sample1;
if (!MathUtilities.AreEqual(sample2Thickness, sample1.Thickness))
{
sample2 = new Sample
{
Name = "Sample2"
};
}
sample2.OC = GetArray(paddock, "OC");
sample2.EC = GetArray(paddock, "EC");
sample2.PH = GetArray(paddock, "PH");
sample2.CL = GetArray(paddock, "CL");
sample2.OCUnits = SoilOrganicMatter.OCUnitsEnum.WalkleyBlack;
sample2.PHUnits = Analysis.PHUnitsEnum.CaCl2;
}
// Make sure we have NH4 values.
if (sample1.NH4 == null && sample1.NO3 != null)
{
string[] defaultValues = StringUtilities.CreateStringArray("0.1", sample1.NO3.Length);
sample1.NH4 = MathUtilities.StringsToDoubles(defaultValues);
}
RemoveNullFieldsFromSample(sample1);
if (sample2 != null)
{
RemoveNullFieldsFromSample(sample2);
}
// Fix up <WaterFormat>
string waterFormatString = GetString(paddock, "WaterFormat");
if (waterFormatString.Contains("Gravimetric"))
{
sample1.SWUnits = Sample.SWUnitsEnum.Gravimetric;
}
else
{
sample1.SWUnits = Sample.SWUnitsEnum.Volumetric;
}
if (MathUtilities.ValuesInArray(sample1.Thickness))
{
simulation.Samples.Add(sample1);
}
if (sample2 != null && MathUtilities.ValuesInArray(sample2.Thickness) && sample2 != sample1)
{
simulation.Samples.Add(sample2);
}
// Check for InitTotalWater & InitTotalNitrogen
simulation.InitTotalWater = GetDouble(paddock, "InitTotalWater");
simulation.InitTotalNitrogen = GetDouble(paddock, "InitTotalNitrogen");
// Check to see if we need to convert the soil structure.
simulation.SoilPath = GetString(paddock, "SoilName");
if (soilNode != null)
{
// See if there is a 'SWUnits' value. If found then copy it into
// <WaterFormat>
string waterFormat = XmlUtilities.Value(paddock, "WaterFormat");
if (waterFormat == string.Empty)
{
int sampleNumber = 0;
foreach (XmlNode soilSample in XmlUtilities.ChildNodes(soilNode, "Sample"))
{
string swUnits = XmlUtilities.Value(soilSample, "SWUnits");
if (swUnits != string.Empty)
{
XmlUtilities.SetValue(paddock, "WaterFormat", swUnits);
}
// Also make sure we don't have 2 samples with the same name.
string sampleName = "Sample" + (sampleNumber + 1).ToString();
XmlUtilities.SetAttribute(soilSample, "name", sampleName);
sampleNumber++;
}
}
simulation.Soil = SoilUtilities.FromXML(soilNode.OuterXml);
if (simulation.Samples != null)
{
simulation.Soil.Samples = simulation.Samples;
}
}
return(simulation);
}
/// <summary>Fills the auto-calculated fields.</summary>
/// <param name="paddock">The paddock.</param>
/// <param name="observedData">The observed data.</param>
/// <param name="weatherData">The weather data.</param>
public static void FillInCalculatedFields(Paddock paddock, DataTable observedData, string workingFolder)
{
IEnumerable<Tillage> tillages = paddock.Management.OfType<Tillage>();
if (tillages.Count() > 0)
paddock.StubbleIncorporatedPercent = YieldProphetUtility.CalculateAverageTillagePercent(tillages);
DateTime lastRainfallDate = GetLastRainfallDate(observedData);
if (lastRainfallDate != DateTime.MinValue)
paddock.DateOfLastRainfallEntry = lastRainfallDate.ToString("dd/MM/yyyy");
string[] metFiles = Directory.GetFiles(workingFolder, "*.met");
if (metFiles.Length > 0)
{
string firstMetFile = Path.Combine(workingFolder, metFiles[0]);
ApsimTextFile textFile = new ApsimTextFile();
textFile.Open(firstMetFile);
DataTable data = textFile.ToTable();
textFile.Close();
paddock.RainfallSinceSoilWaterSampleDate = SumTableAfterDate(data, "Rain", paddock.SoilWaterSampleDate);
if (data.Rows.Count > 0)
{
DataRow lastweatherRow = data.Rows[data.Rows.Count - 1];
paddock.LastClimateDate = DataTableUtilities.GetDateFromRow(lastweatherRow);
}
}
}
/// <summary>Converts the paddock XML.</summary>
/// <param name="paddock">The paddock.</param>
/// <exception cref="System.Exception">Bad paddock name: + name</exception>
private static Paddock CreateSimulationSpec(XmlNode paddock, string baseFolder)
{
Paddock simulation = new Paddock();
string name = XmlUtilities.NameAttr(paddock);
string delimiter = "^";
int posCaret = name.IndexOf(delimiter);
if (posCaret == -1)
throw new Exception("Bad paddock name: " + name);
string remainder = StringUtilities.SplitOffAfterDelimiter(ref name, delimiter);
string growerName;
string paddockName = StringUtilities.SplitOffAfterDelimiter(ref remainder, delimiter);
if (paddockName == string.Empty)
{
growerName = name;
paddockName = remainder;
}
else
growerName = remainder;
simulation.StartSeasonDate = GetDate(paddock, "StartSeasonDateFull");
// Give the paddock a name.
string fullName = string.Format("{0}^{1}^{2}", simulation.StartSeasonDate.Year, growerName, paddockName);
simulation.Name = fullName;
// Set the report date.
simulation.NowDate = GetDate(paddock.ParentNode, "TodayDateFull");
if (simulation.NowDate == DateTime.MinValue)
simulation.NowDate = DateTime.Now;
// Store any rainfall data in the simulation.
simulation.RainfallSource = GetString(paddock, "RainfallSource");
if (simulation.RainfallSource != "Weather station")
{
string rainFileName = GetString(paddock, "RainfallFilename");
if (rainFileName != string.Empty)
{
string fullFileName = Path.Combine(baseFolder, rainFileName);
if (!File.Exists(fullFileName))
throw new Exception("Cannot find file: " + fullFileName);
simulation.ObservedData = ApsimTextFile.ToTable(fullFileName);
if (simulation.ObservedData.Rows.Count == 0)
simulation.ObservedData = null;
}
}
// Set the reset dates
simulation.SoilWaterSampleDate = GetDate(paddock, "ResetDateFull");
simulation.SoilNitrogenSampleDate = GetDate(paddock, "SoilNitrogenSampleDateFull");
simulation.StationNumber = GetInteger(paddock, "StationNumber");
simulation.StationName = GetString(paddock, "StationName");
// Create a sowing management
Sow sowing = new Sow();
simulation.Management.Add(sowing);
sowing.Date = GetDate(paddock, "SowDateFull");
sowing.EmergenceDate = GetDate(paddock, "EmergenceDateFull");
sowing.Crop = GetString(paddock, "Crop");
sowing.Cultivar = GetString(paddock, "Cultivar");
sowing.SkipRow = GetString(paddock, "SkipRow");
sowing.SowingDensity = GetInteger(paddock, "SowingDensity");
sowing.MaxRootDepth = GetInteger(paddock, "MaxRootDepth") * 10; // cm to mm
sowing.BedWidth = GetInteger(paddock, "BedWidth");
sowing.BedRowSpacing = GetDouble(paddock, "BedRowSpacing");
// Make sure we have a stubbletype
simulation.StubbleType = GetString(paddock, "StubbleType");
if (simulation.StubbleType == string.Empty || simulation.StubbleType == "None")
simulation.StubbleType = "Wheat";
simulation.StubbleMass = GetDouble(paddock, "StubbleMass");
simulation.Slope = GetDouble(paddock, "Slope");
simulation.SlopeLength = GetDouble(paddock, "SlopeLength");
simulation.UseEC = GetBoolean(paddock, "UseEC");
// Fertilise nodes.
List<XmlNode> fertiliserNodes = XmlUtilities.ChildNodes(paddock, "Fertilise");
for (int f = 0; f < fertiliserNodes.Count; f++)
{
Fertilise fertilise = new Fertilise();
simulation.Management.Add(fertilise);
fertilise.Date = GetDate(fertiliserNodes[f], "FertDateFull");
fertilise.Amount = GetDouble(fertiliserNodes[f], "FertAmount");
fertilise.Scenario = GetBoolean(fertiliserNodes[f], "Scenario");
}
// Irrigate nodes.
List<XmlNode> irrigateNodes = XmlUtilities.ChildNodes(paddock, "Irrigate");
for (int i = 0; i < irrigateNodes.Count; i++)
{
Irrigate irrigate = new Irrigate();
simulation.Management.Add(irrigate);
irrigate.Date = GetDate(irrigateNodes[i], "IrrigateDateFull");
irrigate.Amount = GetDouble(irrigateNodes[i], "IrrigateAmount");
irrigate.Efficiency = GetDouble(irrigateNodes[i], "IrrigateEfficiency");
irrigate.Scenario = GetBoolean(irrigateNodes[i], "Scenario");
}
// Tillage nodes.
foreach (XmlNode tillageNode in XmlUtilities.ChildNodes(paddock, "Tillage"))
{
Tillage tillage = new Tillage();
simulation.Management.Add(tillage);
tillage.Date = GetDate(tillageNode, "TillageDateFull");
string disturbance = GetString(tillageNode, "Disturbance");
if (disturbance == "Low")
tillage.Disturbance = Tillage.DisturbanceEnum.Low;
else if (disturbance == "Medium")
tillage.Disturbance = Tillage.DisturbanceEnum.Medium;
else
tillage.Disturbance = Tillage.DisturbanceEnum.High;
tillage.Scenario = GetBoolean(tillageNode, "Scenario");
}
// Stubble removed nodes.
foreach (XmlNode stubbleRemovedNode in XmlUtilities.ChildNodes(paddock, "StubbleRemoved"))
{
StubbleRemoved stubbleRemoved = new StubbleRemoved();
simulation.Management.Add(stubbleRemoved);
stubbleRemoved.Date = GetDate(stubbleRemovedNode, "StubbleRemovedDateFull");
stubbleRemoved.Percent = GetDouble(stubbleRemovedNode, "StubbleRemovedAmount");
}
// Look for a soil node.
XmlNode soilNode = XmlUtilities.FindByType(paddock, "Soil");
if (soilNode != null)
{
string testValue = XmlUtilities.Value(soilNode, "Water/Layer/Thickness");
if (testValue != string.Empty)
{
// old format.
soilNode = ConvertSoilNode.Upgrade(soilNode);
}
}
// Fix up soil sample variables.
Sample sample1 = new Sample();
sample1.Name = "Sample1";
sample1.Thickness = GetArray(paddock, "Sample1Thickness");
sample1.SW = GetArray(paddock, "SW");
if (sample1.SW == null)
{
// Really old way of doing samples - they are stored under soil.
List<XmlNode> sampleNodes = XmlUtilities.ChildNodes(soilNode, "Sample");
if (sampleNodes.Count > 0)
sample1 = XmlUtilities.Deserialise(sampleNodes[0], typeof(Sample)) as Sample;
}
else
{
sample1.NO3 = GetArray(paddock, "NO3");
sample1.NH4 = GetArray(paddock, "NH4");
}
Sample sample2 = null;
double[] sample2Thickness = GetArray(paddock, "Sample2Thickness");
if (sample2Thickness == null)
{
// Really old way of doing samples - they are stored under soil.
List<XmlNode> sampleNodes = XmlUtilities.ChildNodes(soilNode, "Sample");
if (sampleNodes.Count > 1)
sample2 = XmlUtilities.Deserialise(sampleNodes[1], typeof(Sample)) as Sample;
}
else
{
sample2 = sample1;
if (!MathUtilities.AreEqual(sample2Thickness, sample1.Thickness))
{
sample2 = new Sample();
sample2.Name = "Sample2";
}
sample2.OC = GetArray(paddock, "OC");
sample2.EC = GetArray(paddock, "EC");
sample2.PH = GetArray(paddock, "PH");
sample2.CL = GetArray(paddock, "CL");
sample2.OCUnits = SoilOrganicMatter.OCUnitsEnum.WalkleyBlack;
sample2.PHUnits = Analysis.PHUnitsEnum.CaCl2;
}
// Make sure we have NH4 values.
if (sample1.NH4 == null && sample1.NO3 != null)
{
string[] defaultValues = StringUtilities.CreateStringArray("0.1", sample1.NO3.Length);
sample1.NH4 = MathUtilities.StringsToDoubles(defaultValues);
}
RemoveNullFieldsFromSample(sample1);
if (sample2 != null)
RemoveNullFieldsFromSample(sample2);
// Fix up <WaterFormat>
string waterFormatString = GetString(paddock, "WaterFormat");
if (waterFormatString.Contains("Gravimetric"))
sample1.SWUnits = Sample.SWUnitsEnum.Gravimetric;
else
sample1.SWUnits = Sample.SWUnitsEnum.Volumetric;
if (MathUtilities.ValuesInArray(sample1.Thickness))
simulation.Samples.Add(sample1);
if (sample2 != null && MathUtilities.ValuesInArray(sample2.Thickness) && sample2 != sample1)
simulation.Samples.Add(sample2);
// Check for InitTotalWater & InitTotalNitrogen
simulation.InitTotalWater = GetDouble(paddock, "InitTotalWater");
simulation.InitTotalNitrogen = GetDouble(paddock, "InitTotalNitrogen");
// Check to see if we need to convert the soil structure.
simulation.SoilPath = GetString(paddock, "SoilName");
if (soilNode != null)
{
// See if there is a 'SWUnits' value. If found then copy it into
// <WaterFormat>
string waterFormat = XmlUtilities.Value(paddock, "WaterFormat");
if (waterFormat == string.Empty)
{
int sampleNumber = 0;
foreach (XmlNode soilSample in XmlUtilities.ChildNodes(soilNode, "Sample"))
{
string swUnits = XmlUtilities.Value(soilSample, "SWUnits");
if (swUnits != string.Empty)
XmlUtilities.SetValue(paddock, "WaterFormat", swUnits);
// Also make sure we don't have 2 samples with the same name.
string sampleName = "Sample" + (sampleNumber + 1).ToString();
XmlUtilities.SetAttribute(soilSample, "name", sampleName);
sampleNumber++;
}
}
simulation.Soil = SoilUtilities.FromXML(soilNode.OuterXml);
if (simulation.Samples != null)
simulation.Soil.Samples = simulation.Samples;
}
return simulation;
}
/// <summary>Creates a base APSIM simulation spec for the Yield Prophet spec.</summary>
/// <param name="yieldProphet">The yield prophet specification.</param>
/// <returns>The created APSIM simulation spec.</returns>
private static APSIMSpecification CreateBaseSimulation(Paddock paddock)
{
Paddock copyOfPaddock = paddock; // XmlUtilities.Clone(paddock) as JobsService.Paddock;
copyOfPaddock.ObservedData = paddock.ObservedData;
APSIMSpecification shortSimulation = new APSIMSpecification();
shortSimulation.Name = "Base";
shortSimulation.WeatherFileName = shortSimulation.Name + ".met";
// Start date of simulation should be the earliest of ResetDate, SowDate and StartSeasonDate
Sow sow = YieldProphetUtility.GetCropBeingSown(paddock.Management);
if (sow == null)
{
throw new Exception("No sowing specified for paddock: " + paddock.Name);
}
if (sow.Date == DateTime.MinValue)
{
throw new Exception("No sowing DATE specified for paddock: " + paddock.Name);
}
if (sow.Crop == null || sow.Crop == "" || sow.Crop == "None")
{
throw new Exception("No sowing CROP specified for paddock: " + paddock.Name);
}
shortSimulation.StartDate = DateTime.MaxValue;
if (paddock.SoilWaterSampleDate != DateTime.MinValue &&
paddock.SoilWaterSampleDate < shortSimulation.StartDate)
{
shortSimulation.StartDate = paddock.SoilWaterSampleDate;
}
if (paddock.SoilNitrogenSampleDate != DateTime.MinValue &&
paddock.SoilNitrogenSampleDate < shortSimulation.StartDate)
{
shortSimulation.StartDate = paddock.SoilNitrogenSampleDate;
}
if (sow != null && sow.Date < shortSimulation.StartDate && sow.Date != DateTime.MinValue)
{
shortSimulation.StartDate = sow.Date;
}
if (paddock.StartSeasonDate < shortSimulation.StartDate)
{
shortSimulation.StartDate = paddock.StartSeasonDate;
}
if (paddock.LongtermStartYear == 0)
{
shortSimulation.LongtermStartYear = 1957;
}
else
{
shortSimulation.LongtermStartYear = paddock.LongtermStartYear;
}
if (paddock.RunType == Paddock.RunTypeEnum.SingleSeason)
{
shortSimulation.EndDate = copyOfPaddock.NowDate.AddDays(-1);
if ((shortSimulation.EndDate - shortSimulation.StartDate).Days > 360)
{
shortSimulation.EndDate = shortSimulation.StartDate.AddDays(360);
}
}
else if (paddock.RunType == Paddock.RunTypeEnum.LongTermPatched)
{
shortSimulation.EndDate = shortSimulation.StartDate.AddDays(360);
}
else if (paddock.RunType == Paddock.RunTypeEnum.LongTerm)
{
// Sow opp reports.
shortSimulation.StartDate = new DateTime(shortSimulation.LongtermStartYear, 1, 1);
shortSimulation.EndDate = copyOfPaddock.NowDate.AddDays(-1);
}
shortSimulation.NowDate = copyOfPaddock.NowDate.AddDays(-1);
if (shortSimulation.NowDate == DateTime.MinValue)
{
shortSimulation.NowDate = DateTime.Now;
}
shortSimulation.DailyOutput = paddock.DailyOutput;
shortSimulation.MonthlyOutput = paddock.MonthlyOutput;
shortSimulation.YearlyOutput = paddock.YearlyOutput;
shortSimulation.ObservedData = copyOfPaddock.ObservedData;
shortSimulation.Soil = copyOfPaddock.Soil;
shortSimulation.SoilPath = copyOfPaddock.SoilPath;
shortSimulation.Samples = new List <APSIM.Shared.Soils.Sample>();
shortSimulation.Samples.AddRange(copyOfPaddock.Samples);
shortSimulation.InitTotalWater = copyOfPaddock.InitTotalWater;
shortSimulation.InitTotalNitrogen = copyOfPaddock.InitTotalNitrogen;
shortSimulation.StationNumber = copyOfPaddock.StationNumber;
shortSimulation.StubbleMass = copyOfPaddock.StubbleMass;
shortSimulation.StubbleType = copyOfPaddock.StubbleType;
shortSimulation.Management = new List <Management>();
shortSimulation.Management.AddRange(copyOfPaddock.Management);
shortSimulation.UseEC = paddock.UseEC;
shortSimulation.WriteDepthFile = false;
if (paddock.RunType == Paddock.RunTypeEnum.LongTermPatched)
{
shortSimulation.RunType = APSIMSpecification.RunTypeEnum.LongTermPatched;
}
else
{
shortSimulation.RunType = APSIMSpecification.RunTypeEnum.Normal;
}
shortSimulation.DecileDate = paddock.StartSeasonDate;
shortSimulation.NUnlimited = paddock.NUnlimited;
shortSimulation.NUnlimitedFromToday = paddock.NUnlimitedFromToday;
shortSimulation.WriteDepthFile = paddock.WriteDepthFile;
shortSimulation.Next10DaysDry = paddock.Next10DaysDry;
AddResetDatesToManagement(copyOfPaddock, shortSimulation);
// Do a stable sort on management actions.
shortSimulation.Management = shortSimulation.Management.OrderBy(m => m.Date).ToList();
if (paddock.RunType == Paddock.RunTypeEnum.LongTerm)
{
foreach (Management man in shortSimulation.Management)
{
man.IsEveryYear = true;
}
}
return(shortSimulation);
}