//--------------------------------------------------------------------- private ExtensionMain[] LoadExtensions(ExtensionAndInitFile[] extensions) { ExtensionMain[] loadedExtensions = new ExtensionMain[extensions.Length]; foreach (int i in Indexes.Of(extensions)) { ExtensionAndInitFile extensionAndInitFile = extensions[i]; ui.WriteLine("Loading {0} extension ...", extensionAndInitFile.Info.Name); ExtensionMain loadedExtension = Loader.Load <ExtensionMain>(extensionAndInitFile.Info); loadedExtension.LoadParameters(extensionAndInitFile.InitFile, this); loadedExtensions[i] = loadedExtension; disturbAndOtherExtensions.Add(loadedExtension); } return(loadedExtensions); }
//--------------------------------------------------------------------- private PlugIn[] LoadAndInitPlugIns(PlugInAndInitFile[] plugIns) { PlugIn[] loadedPlugIns = new PlugIn[plugIns.Length]; foreach (int i in Indexes.Of(plugIns)) { PlugInAndInitFile plugInAndInitFile = plugIns[i]; UI.WriteLine("Loading {0} plug-in ...", plugInAndInitFile.Info.Name); PlugIn loadedPlugIn = Loader.Load <PlugIn>(plugInAndInitFile.Info); loadedPlugIn.Initialize(plugInAndInitFile.InitFile, this); loadedPlugIns[i] = loadedPlugIn; if (loadedPlugIn is PlugIns.I2PhaseInitialization) { plugInsWith2PhaseInit.Add((PlugIns.I2PhaseInitialization)loadedPlugIn); } if (loadedPlugIn is PlugIns.ICleanUp) { plugInsWithCleanUp.Add((PlugIns.ICleanUp)loadedPlugIn); } } return(loadedPlugIns); }
//--------------------------------------------------------------------- /// <summary> /// Runs a model scenario. /// </summary> public void Run(string scenarioPath, IUserInterface ui) { this.ui = ui; siteVarRegistry.Clear(); Scenario scenario = LoadScenario(scenarioPath); startTime = scenario.StartTime; endTime = scenario.EndTime; timeSinceStart = 0; currentTime = startTime; InitializeRandomNumGenerator(scenario.RandomNumberSeed); LoadSpecies(scenario.Species); LoadEcoregions(scenario.Ecoregions); ui.WriteLine("Initializing landscape from ecoregions map \"{0}\" ...", scenario.EcoregionsMap); Ecoregions.Map ecoregionsMap = new Ecoregions.Map(scenario.EcoregionsMap, ecoregions, rasterFactory); // -- ProcessMetadata(ecoregionsMap.Metadata, scenario); cellLength = scenario.CellLength.Value; cellArea = (float)((cellLength * cellLength) / 10000); ui.WriteLine("Cell length = {0} m, cell area = {1} ha", cellLength, cellArea); using (IInputGrid <bool> grid = ecoregionsMap.OpenAsInputGrid()) { ui.WriteLine("Map dimensions: {0} = {1:#,##0} cell{2}", grid.Dimensions, grid.Count, (grid.Count == 1 ? "" : "s")); // landscape = new Landscape(grid); landscape = landscapeFactory.CreateLandscape(grid); } ui.WriteLine("Sites: {0:#,##0} active ({1:p1}), {2:#,##0} inactive ({3:p1})", landscape.ActiveSiteCount, (landscape.Count > 0 ? ((double)landscape.ActiveSiteCount) / landscape.Count : 0), landscape.InactiveSiteCount, (landscape.Count > 0 ? ((double)landscape.InactiveSiteCount) / landscape.Count : 0)); ecoregionSiteVar = ecoregionsMap.CreateSiteVar(landscape); disturbAndOtherExtensions = new List <ExtensionMain>(); try { ui.WriteLine("Loading {0} extension ...", scenario.Succession.Info.Name); succession = Loader.Load <SuccessionMain>(scenario.Succession.Info); succession.LoadParameters(scenario.Succession.InitFile, this); succession.Initialize(); ExtensionMain[] disturbanceExtensions = LoadExtensions(scenario.Disturbances); InitExtensions(disturbanceExtensions); ExtensionMain[] otherExtensions = LoadExtensions(scenario.OtherExtensions); InitExtensions(otherExtensions); // Perform 2nd phase of initialization for non-succession extensions. foreach (ExtensionMain extension in disturbanceExtensions) { extension.InitializePhase2(); } foreach (ExtensionMain extension in otherExtensions) { extension.InitializePhase2(); } // Run output extensions for TimeSinceStart = 0 (time step 0) foreach (ExtensionMain extension in otherExtensions) { if (extension.Type.IsMemberOf("output")) { Run(extension); } } //******************// for Rob // Main time loop // //******************// for (currentTime++, timeSinceStart++; currentTime <= endTime; currentTime++, timeSinceStart++) { bool isSuccTimestep = IsExtensionTimestep(succession); List <ExtensionMain> distExtensionsToRun; distExtensionsToRun = GetExtensionsToRun(disturbanceExtensions); bool isDistTimestep = distExtensionsToRun.Count > 0; List <ExtensionMain> otherExtensionsToRun; otherExtensionsToRun = GetExtensionsToRun(otherExtensions); if (!isSuccTimestep && !isDistTimestep && otherExtensionsToRun.Count == 0) { continue; } ui.WriteLine("Current time: {0}", currentTime); if (isDistTimestep) { if (scenario.DisturbancesRandomOrder) { distExtensionsToRun = shuffle(distExtensionsToRun); } foreach (ExtensionMain distExtension in distExtensionsToRun) { Run(distExtension); } } if (isSuccTimestep || isDistTimestep) { Run(succession); } foreach (ExtensionMain otherExtension in otherExtensionsToRun) { Run(otherExtension); } } // main time loop } finally { foreach (ExtensionMain extension in disturbAndOtherExtensions) { extension.CleanUp(); } } ui.WriteLine("Model run is complete."); }
//--------------------------------------------------------------------- /// <summary> /// Runs a model scenario. /// </summary> public void Run(string scenarioPath) { siteVarRegistry.Clear(); Scenario scenario = LoadScenario(scenarioPath); startTime = scenario.StartTime; endTime = scenario.EndTime; timeSinceStart = 0; currentTime = startTime; InitializeRandomNumGenerator(scenario.RandomNumberSeed); LoadSpecies(scenario.Species); LoadEcoregions(scenario.Ecoregions); UI.WriteLine("Initializing landscape from ecoregions map \"{0}\" ...", scenario.EcoregionsMap); Ecoregions.Map ecoregionsMap = new Ecoregions.Map(scenario.EcoregionsMap, ecoregions, rasterDriverMgr); ProcessMetadata(ecoregionsMap.Metadata, scenario); using (IInputGrid <EcoregionCode> grid = ecoregionsMap.OpenAsInputGrid()) { UI.WriteLine("Map dimensions: {0} = {1:#,##0} cell{2}", grid.Dimensions, grid.Count, (grid.Count == 1 ? "" : "s")); landscape = new Landscape(grid, scenario.BlockSize.HasValue ? scenario.BlockSize.Value : 1); } UI.WriteLine("Sites: {0:#,##0} active ({1:p1}), {2:#,##0} inactive ({3:p1})", landscape.ActiveSiteCount, (landscape.Count > 0 ? ((double)landscape.ActiveSiteCount) / landscape.Count : 0), landscape.InactiveSiteCount, (landscape.Count > 0 ? ((double)landscape.InactiveSiteCount) / landscape.Count : 0)); ecoregionSiteVar = ecoregionsMap.CreateSiteVar(landscape); plugInsWith2PhaseInit = new List <PlugIns.I2PhaseInitialization>(); plugInsWithCleanUp = new List <PlugIns.ICleanUp>(); try { UI.WriteLine("Loading {0} plug-in ...", scenario.Succession.Info.Name); succession = Loader.Load <SuccessionPlugIn>(scenario.Succession.Info); succession.Initialize(scenario.Succession.InitFile, this); succession.InitializeSites(scenario.InitialCommunities, scenario.InitialCommunitiesMap); PlugIn[] disturbancePlugIns = LoadAndInitPlugIns(scenario.Disturbances); PlugIn[] otherPlugIns = LoadAndInitPlugIns(scenario.OtherPlugIns); // Perform 2nd phase of initialization for those plug-ins // that have it. foreach (PlugIns.I2PhaseInitialization plugIn in plugInsWith2PhaseInit) { plugIn.InitializePhase2(); } // Run output plug-ins for TimeSinceStart = 0 (time step 0) foreach (PlugIn plugIn in otherPlugIns) { if (plugIn.PlugInType.IsMemberOf("output")) { Run(plugIn); } } //******************// for Rob // Main time loop // //******************// for (currentTime++, timeSinceStart++; currentTime <= endTime; currentTime++, timeSinceStart++) { bool isSuccTimestep = IsPlugInTimestep(succession); List <PlugIn> distPlugInsToRun; distPlugInsToRun = GetPlugInsToRun(disturbancePlugIns); bool isDistTimestep = distPlugInsToRun.Count > 0; List <PlugIn> otherPlugInsToRun; otherPlugInsToRun = GetPlugInsToRun(otherPlugIns); if (!isSuccTimestep && !isDistTimestep && otherPlugInsToRun.Count == 0) { continue; } UI.WriteLine("Current time: {0}", currentTime); if (isDistTimestep) { if (scenario.DisturbancesRandomOrder) { Util.Random.Shuffle(distPlugInsToRun); } foreach (PlugIn distPlugIn in distPlugInsToRun) { Run(distPlugIn); } } if (isSuccTimestep || isDistTimestep) { Run(succession); } foreach (PlugIn otherPlugIn in otherPlugInsToRun) { Run(otherPlugIn); } } // main time loop } finally { foreach (PlugIns.ICleanUp plugIn in plugInsWithCleanUp) { plugIn.CleanUp(); } } }