/// <summary>
/// Run ecological processes that operate across grid cells
/// </summary>
public void RunCrossGridCellEcology(ref uint dispersals, bool dispersalOnly, MadingleyModelInitialisation modelInitialisation)
{
// If we are running specific locations, then we do not run dispersal
if (SpecificLocations != true)
{
if (RunGridCellsInParallel)
{
// Loop through each grid cell, and run dispersal for each. Note that because cells essentially run independently, we do not need to thread-lock variables
// as they do not exchange information until they are all completed (they basically just build up delta arrays of cohorts to move). However, if cells
// should start to exchange information for dispersal, or all contribute to a single centralised variable, then thread-locked parallel variables would
// have to be used.
Parallel.For(0, _CellList.Count, ii =>
{
EcologyCrossGridCell TempMadingleyEcologyCrossGridCell = new EcologyCrossGridCell();
// Initialise cross grid cell ecology
TempMadingleyEcologyCrossGridCell.InitializeCrossGridCellEcology(_GlobalModelTimeStepUnit, DrawRandomly, modelInitialisation);
//Initialise the delta for dispersal lists for this grid cell
EcosystemModelGrid.DeltaFunctionalGroupDispersalArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCohortNumberDispersalArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCellToDisperseToArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint[]>();
EcosystemModelGrid.DeltaCellExitDirection[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCellEntryDirection[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
// We have looped through individal cells and calculated ecological processes for each. Now do this for cross grid cell processes
TempMadingleyEcologyCrossGridCell.RunCrossGridCellEcology(_CellList[ii], EcosystemModelGrid, dispersalOnly,
CohortFunctionalGroupDefinitions, StockFunctionalGroupDefinitions, CurrentMonth);
});
}
else
{
// Loop through each grid cell, and run dispersal for each.
// Note that currently dispersal is not parallelised, although it could be (though care would need to be taken to ensure that necessary variables are thread-locked
for (int ii = 0; ii < _CellList.Count; ii++)
{
//Initialise the delta for dispersal lists for this grid cell
EcosystemModelGrid.DeltaFunctionalGroupDispersalArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCohortNumberDispersalArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCellToDisperseToArray[_CellList[ii][0], _CellList[ii][1]] = new List<uint[]>();
EcosystemModelGrid.DeltaCellExitDirection[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
EcosystemModelGrid.DeltaCellEntryDirection[_CellList[ii][0], _CellList[ii][1]] = new List<uint>();
// We have looped through individal cells and calculated ecological processes for each. Now do this for cross grid cell processes
MadingleyEcologyCrossGridCell.RunCrossGridCellEcology(_CellList[ii], EcosystemModelGrid, dispersalOnly,
CohortFunctionalGroupDefinitions, StockFunctionalGroupDefinitions, CurrentMonth);
}
}
// Apply the changes in the delta arrays from dispersal
MadingleyEcologyCrossGridCell.UpdateCrossGridCellEcology(EcosystemModelGrid, ref dispersals, TrackCrossCellProcesses, CurrentTimeStep);
}
}
