private void OnDoSoilWaterMovement(object sender, EventArgs e)
{
// Calculate lateral flow.
lateralFlowModel.Calculate();
if (LateralFlow.Length > 0)
{
Water = MathUtilities.Subtract(Water, LateralFlow);
}
// Calculate runoff.
Runoff = runoffModel.Value();
// Calculate infiltration.
Infiltration = PotentialInfiltration - Runoff;
Water[0] = Water[0] + Infiltration + Runon;
// Allow irrigation to infiltrate.
foreach (var irrigation in irrigations)
{
if (irrigation.Amount > 0)
{
int irrigationLayer = soil.LayerIndexOfDepth(Convert.ToInt32(irrigation.Depth, CultureInfo.InvariantCulture));
Water[irrigationLayer] += irrigation.Amount;
if (irrigationLayer == 0)
{
Infiltration += irrigation.Amount;
}
if (no3 != null)
{
no3.kgha[irrigationLayer] += irrigation.NO3;
}
if (nh4 != null)
{
nh4.kgha[irrigationLayer] += irrigation.NH4;
}
if (cl != null)
{
cl.kgha[irrigationLayer] += irrigation.CL;
}
}
}
// Saturated flow.
Flux = saturatedFlow.Values;
// Add backed up water to runoff.
Water[0] = Water[0] - saturatedFlow.backedUpSurface;
// Now reduce the infiltration amount by what backed up.
Infiltration = Infiltration - saturatedFlow.backedUpSurface;
// Turn the proportion of the infiltration that backed up into runoff.
Runoff = Runoff + saturatedFlow.backedUpSurface;
// Should go to pond if one exists.
// pond = Math.Min(Runoff, max_pond);
MoveDown(Water, Flux);
double[] no3Values = no3.kgha;
double[] ureaValues = urea.kgha;
// Calcualte solute movement down with water.
double[] no3Down = CalculateSoluteMovementDown(no3Values, Water, Flux, SoluteFluxEfficiency);
MoveDown(no3Values, no3Down);
double[] ureaDown = CalculateSoluteMovementDown(ureaValues, Water, Flux, SoluteFluxEfficiency);
MoveDown(ureaValues, ureaDown);
// Calculate evaporation and remove from top layer.
double es = evaporationModel.Calculate();
Water[0] = Water[0] - es;
// Calculate unsaturated flow of water and apply.
Flow = unsaturatedFlow.Values;
MoveUp(Water, Flow);
// Check for errors in water variables.
//CheckForErrors();
// Calculate water table depth.
waterTableModel.Calculate();
// Calculate and apply net solute movement.
double[] no3Up = CalculateNetSoluteMovement(no3Values, Water, Flow, SoluteFlowEfficiency);
MoveUp(no3Values, no3Up);
double[] ureaUp = CalculateNetSoluteMovement(ureaValues, Water, Flow, SoluteFlowEfficiency);
MoveUp(ureaValues, ureaUp);
// Update flow output variables.
FlowNO3 = MathUtilities.Subtract(no3Down, no3Up);
FlowUrea = MathUtilities.Subtract(ureaDown, ureaUp);
// Set solute state variables.
no3.SetKgHa(SoluteSetterType.Soil, no3Values);
urea.SetKgHa(SoluteSetterType.Soil, ureaValues);
// Now that we've finished moving water, calculate volumetric water
waterVolumetric = MathUtilities.Divide(Water, soil.Thickness);
}
private void OnDoSoilWaterMovement(object sender, EventArgs e)
{
// Calculate lateral flow.
LateralFlow = lateralFlowModel.Values;
MathUtilities.Subtract(Water, LateralFlow);
// Calculate runoff.
Runoff = runoffModel.Value();
// Calculate infiltration.
Infiltration = PotentialInfiltration - Runoff;
Water[0] = Water[0] + Infiltration;
// Allow irrigation to infiltrate.
if (!irrigation.WillRunoff)
{
int irrigationLayer = APSIM.Shared.APSoil.SoilUtilities.FindLayerIndex(properties, Convert.ToInt32(irrigation.Depth, CultureInfo.InvariantCulture));
Water[irrigationLayer] = irrigation.IrrigationApplied;
Infiltration += irrigation.IrrigationApplied;
// DeanH - haven't implemented solutes in irrigation water yet.
// NO3[irrigationLayer] = irrigation.NO3;
// NH4[irrigationLayer] = irrigation.NH4;
// CL[irrigationLayer] = irrigation.Cl;
}
// Saturated flow.
Flux = saturatedFlow.Values;
// Add backed up water to runoff.
Water[0] = Water[0] - saturatedFlow.backedUpSurface;
// Now reduce the infiltration amount by what backed up.
Infiltration = Infiltration - saturatedFlow.backedUpSurface;
// Turn the proportion of the infiltration that backed up into runoff.
Runoff = Runoff + saturatedFlow.backedUpSurface;
// Should go to pond if one exists.
// pond = Math.Min(Runoff, max_pond);
MoveDown(Water, Flux);
double[] NO3Values = soilNitrogen.CalculateNO3();
double[] NH4Values = soilNitrogen.CalculateNH4();
// Calcualte solute movement down with water.
double[] NO3Down = CalculateSoluteMovementDown(NO3Values, Water, Flux, SoluteFluxEfficiency);
double[] NH4Down = CalculateSoluteMovementDown(NH4Values, Water, Flux, SoluteFluxEfficiency);
MoveDown(NO3Values, NO3Down);
MoveDown(NH4Values, NH4Down);
double es = evaporationModel.Calculate();
Water[0] = Water[0] - es;
Flow = unsaturatedFlow.Values;
MoveUp(Water, Flow);
CheckForErrors();
double waterTableDepth = waterTableModel.Value();
double[] NO3Up = CalculateSoluteMovementUpDown(soilNitrogen.CalculateNO3(), Water, Flow, SoluteFlowEfficiency);
double[] NH4Up = CalculateSoluteMovementUpDown(soilNitrogen.CalculateNH4(), Water, Flow, SoluteFlowEfficiency);
MoveUp(NO3Values, NO3Up);
MoveUp(NH4Values, NH4Up);
// Set deltas
NO3.SetKgHa(SoluteSetterType.Soil, MathUtilities.Subtract(soilNitrogen.CalculateNO3(), NO3Values));
NH4.SetKgHa(SoluteSetterType.Soil, MathUtilities.Subtract(soilNitrogen.CalculateNH4(), NH4Values));
}