/// <summary>
/// Utility routines for manipulating the DM_Pool type. AddDMPool adds the
/// contents of two pools together
/// </summary>
/// <param name="Pool1"></param>
/// <param name="Pool2"></param>
/// <returns></returns>
protected GrazType.DM_Pool AddDMPool(GrazType.DM_Pool Pool1, GrazType.DM_Pool Pool2)
{
int N = (int)GrazType.TOMElement.N;
int P = (int)GrazType.TOMElement.P;
int S = (int)GrazType.TOMElement.S;
GrazType.DM_Pool Result = new GrazType.DM_Pool();
Result.DM = Pool1.DM + Pool2.DM;
Result.Nu[N] = Pool1.Nu[N] + Pool2.Nu[N];
Result.Nu[P] = Pool1.Nu[P] + Pool2.Nu[P];
Result.Nu[S] = Pool1.Nu[S] + Pool2.Nu[S];
Result.AshAlk = Pool1.AshAlk + Pool2.AshAlk;
return Result;
}
/// <summary>
/// MultiplyDMPool scales the contents of a pool
/// </summary>
/// <param name="Src"></param>
/// <param name="X"></param>
/// <returns></returns>
protected GrazType.DM_Pool MultiplyDMPool(GrazType.DM_Pool Src, double X)
{
int N = (int)GrazType.TOMElement.N;
int P = (int)GrazType.TOMElement.P;
int S = (int)GrazType.TOMElement.S;
GrazType.DM_Pool Result = new GrazType.DM_Pool();
Result.DM = Src.DM * X;
Result.Nu[N] = Src.Nu[N] * X;
Result.Nu[P] = Src.Nu[P] * X;
Result.Nu[S] = Src.Nu[S] * X;
Result.AshAlk = Src.AshAlk * X;
return Result;
}
/// <summary>
/// Populate the dry matter pool
/// </summary>
/// <param name="Model"></param>
/// <param name="iCode"></param>
/// <param name="bUseYoung"></param>
/// <param name="bUseAll"></param>
/// <param name="bUseTag"></param>
/// <param name="aValue"></param>
/// <returns></returns>
public static bool PopulateDMPoolValue(TStockList Model, int iCode, bool bUseYoung, bool bUseAll, bool bUseTag, ref TDMPoolHead[] aValue)
{
int iNoPasses;
TAnimalGroup aGroup;
GrazType.DM_Pool Pool = new GrazType.DM_Pool();
GrazType.DM_Pool TotalPool = new GrazType.DM_Pool();
int iDenom;
int iPass, Idx;
bool Result = true;
for (int i = 0; i < aValue.Length; i++)
aValue[i] = new TDMPoolHead();
if (bUseTag)
iNoPasses = aValue.Length;
else
iNoPasses = 1;
for (iPass = 1; iPass <= iNoPasses; iPass++)
{
GrazType.ZeroDMPool(ref TotalPool);
iDenom = 0;
for (Idx = 1; Idx <= Model.Count(); Idx++)
{
if (!bUseTag || (Model.getTag(Idx) == iPass))
{
if (!bUseYoung)
aGroup = Model.At(Idx);
else
aGroup = Model.At(Idx).Young;
GrazType.ZeroDMPool(ref Pool);
if (aGroup != null)
{
int N = (int)GrazType.TOMElement.N;
int P = (int)GrazType.TOMElement.P;
int S = (int)GrazType.TOMElement.S;
switch (iCode)
{
case StockProps.prpINTAKE:
Pool.DM = aGroup.AnimalState.DM_Intake.Total;
Pool.Nu[N] = aGroup.AnimalState.CP_Intake.Total / GrazType.N2Protein;
Pool.Nu[P] = aGroup.AnimalState.Phos_Intake.Total;
Pool.Nu[S] = aGroup.AnimalState.Sulf_Intake.Total;
Pool.AshAlk = aGroup.AnimalState.PaddockIntake.AshAlkalinity * aGroup.AnimalState.PaddockIntake.Biomass
+ aGroup.AnimalState.SuppIntake.AshAlkalinity * aGroup.AnimalState.SuppIntake.Biomass;
break;
case StockProps.prpINTAKE_PAST:
Pool.DM = aGroup.AnimalState.DM_Intake.Herbage;
Pool.Nu[N] = aGroup.AnimalState.CP_Intake.Herbage / GrazType.N2Protein;
Pool.Nu[P] = aGroup.AnimalState.Phos_Intake.Herbage;
Pool.Nu[S] = aGroup.AnimalState.Sulf_Intake.Herbage;
Pool.AshAlk = aGroup.AnimalState.PaddockIntake.AshAlkalinity * aGroup.AnimalState.PaddockIntake.Biomass;
break;
case StockProps.prpINTAKE_SUPP:
Pool.DM = aGroup.AnimalState.DM_Intake.Supp;
Pool.Nu[N] = aGroup.AnimalState.CP_Intake.Supp / GrazType.N2Protein;
Pool.Nu[P] = aGroup.AnimalState.Phos_Intake.Supp;
Pool.Nu[S] = aGroup.AnimalState.Sulf_Intake.Supp;
Pool.AshAlk = aGroup.AnimalState.SuppIntake.AshAlkalinity * aGroup.AnimalState.SuppIntake.Biomass;
break;
case StockProps.prpFAECES:
GrazType.AddDMPool(aGroup.AnimalState.OrgFaeces, Pool);
GrazType.AddDMPool(aGroup.AnimalState.InOrgFaeces, Pool);
break;
case StockProps.prpINORG_FAECES:
GrazType.AddDMPool(aGroup.AnimalState.InOrgFaeces, Pool);
break;
default:
Result = false;
break;
}
}
if (!bUseTag && !bUseAll)
{
DMPool2Value(Pool, aValue[Idx - 1], (iCode == StockProps.prpINORG_FAECES));
}
else if (aGroup != null)
{
GrazType.AddDMPool(GrazType.MultiplyDMPool(Pool, aGroup.NoAnimals), TotalPool);
iDenom = iDenom + aGroup.NoAnimals;
}
}
} //_ loop over animal groups _
if (bUseTag || bUseAll)
{
if (iDenom > 0)
TotalPool = GrazType.PoolFraction(TotalPool, 1.0 / iDenom);
if (bUseAll)
{
DMPool2Value(TotalPool, aValue[0], (iCode == StockProps.prpINORG_FAECES));
}
else
{
DMPool2Value(TotalPool, aValue[iPass - 1], (iCode == StockProps.prpINORG_FAECES));
}
}
}
return Result;
}
