/// <summary>
/// Initializes the ecosystem model
/// </summary>
/// <param name="initialisation">An instance of the model initialisation class</param>
/// <param name="scenarioParameters">The parameters for the scenarios to run</param>
/// <param name="scenarioIndex">The index of the scenario being run</param>
/// <param name="outputFilesSuffix">The suffix to be applied to all outputs from this model run</param>
/// <param name="globalModelTimeStepUnit">The time step unit used in the model</param>
/// <param name="simulation">The index of the simulation being run</param>
public MadingleyModel(MadingleyModelInitialisation initialisation, ScenarioParameterInitialisation scenarioParameters, int scenarioIndex,
string outputFilesSuffix, string globalModelTimeStepUnit, int simulation)
{
// Assign the properties for this model run
AssignModelRunProperties(initialisation, scenarioParameters, scenarioIndex, outputFilesSuffix);
// Set up list of global diagnostics
SetUpGlobalDiagnosticsList();
// Set up the model grid
SetUpModelGrid(initialisation, scenarioParameters, scenarioIndex, simulation);
// Set up model outputs
SetUpOutputs(initialisation, simulation, scenarioIndex);
// Make the initial outputs
InitialOutputs(outputFilesSuffix, initialisation, CurrentMonth);
// Instance the array of process trackers
ProcessTrackers = new ProcessTracker[_CellList.Count];
// Temporary variables
Boolean varExists;
// Set up process trackers for each grid cell
for (int i = 0; i < _CellList.Count; i++)
{
ProcessTrackers[i] = new ProcessTracker(NumTimeSteps,
EcosystemModelGrid.Lats, EcosystemModelGrid.Lons,
_CellList,
initialisation.ProcessTrackingOutputs,
initialisation.TrackProcesses,
CohortFunctionalGroupDefinitions,
EcosystemModelGrid.GlobalMissingValue,
outputFilesSuffix,
initialisation.OutputPath, initialisation.ModelMassBins,
SpecificLocations, i, initialisation,
EcosystemModelGrid.GetEnviroLayer("Realm", 0, _CellList[i][0], _CellList[i][1], out varExists) == 2.0,
EcosystemModelGrid.LatCellSize,
EcosystemModelGrid.LonCellSize);
}
// Set up a cross cell process tracker
TrackCrossCellProcesses = new CrossCellProcessTracker(initialisation.TrackCrossCellProcesses, "DispersalData", initialisation.OutputPath, outputFilesSuffix);
//Set up a global process tracker
if (SpecificLocations) initialisation.TrackGlobalProcesses = false;
TrackGlobalProcesses = new GlobalProcessTracker(NumTimeSteps,
EcosystemModelGrid.Lats, EcosystemModelGrid.Lons,
_CellList,
initialisation.ProcessTrackingOutputs,
initialisation.TrackGlobalProcesses,
CohortFunctionalGroupDefinitions,
StockFunctionalGroupDefinitions,
EcosystemModelGrid.GlobalMissingValue,
outputFilesSuffix,
initialisation.OutputPath, initialisation.ModelMassBins,
SpecificLocations, initialisation,
EcosystemModelGrid.LatCellSize,
EcosystemModelGrid.LonCellSize);
//Set-up the instance of OutputModelState
WriteModelState = new OutputModelState(initialisation, outputFilesSuffix, simulation);
if (SpecificLocations) initialisation.RunRealm = "";
// Record the initial cohorts in the process trackers
RecordInitialCohorts();
// Initialise the class for cross-grid-cell ecology
MadingleyEcologyCrossGridCell = new EcologyCrossGridCell();
// Initialise the time step timer
TimeStepTimer = new StopWatch();
EcologyTimer = new StopWatch();
OutputTimer = new StopWatch();
// Set the global model time step unit
_GlobalModelTimeStepUnit = globalModelTimeStepUnit;
// Initialise the utility functions
Utilities = new UtilityFunctions();
// Initialise the climate change impacts class
ClimateChangeSimulator = new ClimateChange();
// Initialise the harvesting impacts class
HarvestingSimulator = new Harvesting(EcosystemModelGrid.Lats, EcosystemModelGrid.Lons, (float)EcosystemModelGrid.LatCellSize);
}
/// <summary>
/// Initializes the ecosystem model
/// </summary>
/// <param name="initialisation">An instance of the model initialisation class</param>
/// <param name="scenarioParameters">The parameters for the scenarios to run</param>
/// <param name="scenarioIndex">The index of the scenario being run</param>
/// <param name="outputFilesSuffix">The suffix to be applied to all outputs from this model run</param>
/// <param name="globalModelTimeStepUnit">The time step unit used in the model</param>
/// <param name="simulation">The index of the simulation being run</param>
public MadingleyModel(MadingleyModelInitialisation initialisation, ScenarioParameterInitialisation scenarioParameters, int scenarioIndex,
string outputFilesSuffix, string globalModelTimeStepUnit, int simulation)
{
// Assign the properties for this model run
AssignModelRunProperties(initialisation, scenarioParameters, scenarioIndex, outputFilesSuffix);
#endif
// Set up the model grid
#if true
this._CellList = initialisation.CellList.ToList();
SetUpModelGrid(initialisation);
EcosystemModelGrid.SetGridCells(gridCells, this._CellList);
#else
SetUpModelGrid(initialisation, scenarioParameters, scenarioIndex, simulation);
#endif
// Set up model outputs
SetUpOutputs(initialisation, simulation, scenarioIndex);
// Make the initial outputs
InitialOutputs(outputFilesSuffix, initialisation, CurrentMonth);
// Instance the array of process trackers
ProcessTrackers = new ProcessTracker[_CellList.Count];
// Temporary variables
Boolean varExists;
// Set up process trackers for each grid cell
for (int i = 0; i < _CellList.Count; i++)
{
ProcessTrackers[i] = new ProcessTracker(NumTimeSteps,
EcosystemModelGrid.Lats, EcosystemModelGrid.Lons,
_CellList,
initialisation.ProcessTrackingOutputs,
initialisation.TrackProcesses,
CohortFunctionalGroupDefinitions,
EcosystemModelGrid.GlobalMissingValue,
outputFilesSuffix,
initialisation.OutputPath, initialisation.ModelMassBins,
SpecificLocations, i, initialisation,
EcosystemModelGrid.GetEnviroLayer("Realm", 0, _CellList[i][0], _CellList[i][1], out varExists) == 2.0,
EcosystemModelGrid.LatCellSize,
EcosystemModelGrid.LonCellSize);
}
// Set up a cross cell process tracker
TrackCrossCellProcesses = new CrossCellProcessTracker(initialisation.TrackCrossCellProcesses, "DispersalData", initialisation.OutputPath, outputFilesSuffix);
//Set up a global process tracker
if (SpecificLocations) initialisation.TrackGlobalProcesses = false;
TrackGlobalProcesses = new GlobalProcessTracker(NumTimeSteps,
EcosystemModelGrid.Lats, EcosystemModelGrid.Lons,
_CellList,
initialisation.ProcessTrackingOutputs,
initialisation.TrackGlobalProcesses,
CohortFunctionalGroupDefinitions,
StockFunctionalGroupDefinitions,
EcosystemModelGrid.GlobalMissingValue,
outputFilesSuffix,
initialisation.OutputPath, initialisation.ModelMassBins,
SpecificLocations, initialisation,
EcosystemModelGrid.LatCellSize,
EcosystemModelGrid.LonCellSize);
//Set-up the instance of OutputModelState
WriteModelState = new OutputModelState(initialisation, outputFilesSuffix, simulation);
if (SpecificLocations) initialisation.RunRealm = "";
// Record the initial cohorts in the process trackers
RecordInitialCohorts();
// Initialise the time step timer
TimeStepTimer = new StopWatch();
EcologyTimer = new StopWatch();
OutputTimer = new StopWatch();
// Set the global model time step unit
_GlobalModelTimeStepUnit = globalModelTimeStepUnit;
}