public static Mortality ( IInsect insect ) : void | ||
insect | IInsect | |
Résultat | void |
//--------------------------------------------------------------------- ///<summary> /// Run Biomass Insects extension at a particular timestep. ///</summary> public override void Run() { running = true; PlugIn.ModelCore.UI.WriteLine(" Processing landscape for Biomass Insect events ..."); SiteVars.SiteDefoliation.ActiveSiteValues = 0; foreach (IInsect insect in manyInsect) { //SiteVars.BiomassRemoved.ActiveSiteValues = 0; //SiteVars.InitialOutbreakProb.ActiveSiteValues = 0.0; if (insect.MortalityYear == PlugIn.ModelCore.CurrentTime) { Outbreak.Mortality(insect); } if (insect.MortalityYear != PlugIn.ModelCore.CurrentTime) { insect.LastBioRemoved = 0; } // Copy the data from current to last foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { double thisYearDefol = insect.ThisYearDefoliation[site]; insect.LastYearDefoliation[site] = thisYearDefol; SiteVars.SiteDefoliation[site] += (int)Math.Round(thisYearDefol * 100); if (thisYearDefol > 0) { SiteVars.TimeOfLastEvent[site] = PlugIn.ModelCore.CurrentTime - 1; SiteVars.InsectName[site] = insect.Name; } } insect.ThisYearDefoliation.ActiveSiteValues = 0.0; insect.ActiveOutbreak = false; insect.SingleOutbreakYear = false; //PlugIn.ModelCore.NormalDistribution.Mu = 0.0; //PlugIn.ModelCore.NormalDistribution.Sigma = 1.0; //double randomNum = PlugIn.ModelCore.NormalDistribution.NextDouble(); //randomNum = PlugIn.ModelCore.NormalDistribution.NextDouble(); double randomNum = Distribution.GenerateRandomNum(DistributionType.Normal, 0.0, 1.0); DistributionType distDuration = insect.DurationDistribution; double randomNumE = Distribution.GenerateRandomNum(distDuration, insect.DurationParameter1, insect.DurationParameter2); // First, has enough time passed since the last outbreak? This is calculated each year, // but only used in last year of an outbreak to generate next outbreak characteristics. double timeBetweenOutbreaks = insect.MeanTimeBetweenOutbreaks + (insect.StdDevTimeBetweenOutbreaks * randomNum); double duration = System.Math.Round(randomNumE); if (distDuration == DistributionType.Exponential || distDuration == DistributionType.Normal) { duration = System.Math.Max(duration + 1, 1); // Duration cannot be less than 1. Added to allow Normal distribution parameters. } else { duration = System.Math.Max(duration, 1); // Duration cannot be less than 1. } // Apply optional maximum - if not used maxDur = Inf if (duration > insect.MaxDuration) { duration = insect.MaxDuration; } double timeAfterDuration = timeBetweenOutbreaks - duration; // Users can parameterize model to have overlapping outbreaks, but then patches will not initialize correctly. // Do below to prevent overlapping outbreaks of same insect. This will affect the realized distribution of time between outbreaks somewhat. while (timeAfterDuration <= 1.0) { PlugIn.ModelCore.UI.WriteLine("CAUTION: Time Between Outbreaks TOO SMALL (< 1), you may need to adjust timing parameters: {0}.", timeAfterDuration); timeBetweenOutbreaks = timeBetweenOutbreaks + 1; timeAfterDuration = timeBetweenOutbreaks - duration; } //PlugIn.ModelCore.UI.WriteLine("Calculated time between = {0:0.00}. inputMeanTime={1:0.00}, inputStdTime={2:0.00}., timeAftDur={3:0.00}.", timeBetweenOutbreaks, insect.MeanTimeBetweenOutbreaks, insect.StdDevTimeBetweenOutbreaks, timeAfterDuration); // The logic below determines whether an outbreak is active. And sets a new outbreak duration and timeBetweenOutbreaks // for the next outbreak if is the last year of an outbreak. // The very first outbreak is set first. if (PlugIn.ModelCore.CurrentTime == 1) { PlugIn.ModelCore.UI.WriteLine(" Year 1: Setting initial start and stop times."); double randomNum1 = PlugIn.ModelCore.GenerateUniform(); //insect.OutbreakStartYear = (int) (timeBetweenOutbreaks / 2.0) + 1; insect.OutbreakStartYear = Math.Max(2, (int)(randomNum1 * timeBetweenOutbreaks + 1)); // New, try making 1st start year more random. 1st outbreak has to occur > year1 to for InitializeDefoliationPatches to work properly. insect.OutbreakStartYear = Math.Max(insect.OutbreakStartYear, insect.StartYear); // BRM - optional start year (first outbreak can't be before StartYear) insect.OutbreakStopYear = insect.OutbreakStartYear + (int)duration - 1; // PlugIn.ModelCore.UI.WriteLine(" {0} is not active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); } else if (insect.OutbreakStartYear <= PlugIn.ModelCore.CurrentTime && insect.OutbreakStopYear > PlugIn.ModelCore.CurrentTime) // first year of a multiyear outbreak { // PlugIn.ModelCore.UI.WriteLine(" An outbreak starts or continues. Start and stop time do not change."); insect.ActiveOutbreak = true; PlugIn.ModelCore.UI.WriteLine(" {0} is active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); insect.MortalityYear = PlugIn.ModelCore.CurrentTime + 1; insect.LastStartYear = insect.OutbreakStartYear; // Added here for Brian's log file. insect.LastStopYear = insect.OutbreakStopYear; // Added here for Brian's log file. } //Special case for single year outbreak. else if (insect.OutbreakStartYear == PlugIn.ModelCore.CurrentTime && insect.OutbreakStopYear == PlugIn.ModelCore.CurrentTime) { insect.ActiveOutbreak = true; PlugIn.ModelCore.UI.WriteLine(" {0} is active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); if (insect.OutbreakStartYear == insect.OutbreakStopYear) // shouldn't need this. JRF { insect.SingleOutbreakYear = true; } insect.MortalityYear = PlugIn.ModelCore.CurrentTime + 1; insect.LastStartYear = insect.OutbreakStartYear; // Added here for Brian's log file. insect.LastStopYear = insect.OutbreakStopYear; // Added here for Brian's log file. insect.OutbreakStartYear = PlugIn.ModelCore.CurrentTime + (int)timeBetweenOutbreaks; insect.OutbreakStopYear = insect.OutbreakStartYear + (int)duration - 1; } else if (insect.OutbreakStopYear <= PlugIn.ModelCore.CurrentTime && timeAfterDuration > PlugIn.ModelCore.CurrentTime - insect.OutbreakStopYear) { // PlugIn.ModelCore.UI.WriteLine(" In between outbreaks, reset start and stop times."); insect.ActiveOutbreak = true; PlugIn.ModelCore.UI.WriteLine(" {0} is active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); insect.MortalityYear = PlugIn.ModelCore.CurrentTime + 1; insect.LastStartYear = insect.OutbreakStartYear; // Added here for Brian's log file. insect.LastStopYear = insect.OutbreakStopYear; // Added here for Brian's log file. insect.OutbreakStartYear = PlugIn.ModelCore.CurrentTime + (int)timeBetweenOutbreaks; insect.OutbreakStopYear = insect.OutbreakStartYear + (int)duration - 1; } PlugIn.ModelCore.UI.WriteLine(" Insect Start Time = {0}, Stop Time = {1}.", insect.OutbreakStartYear, insect.OutbreakStopYear); // Now that logic determining when an outbreak will be active is done, tell model what to do when outbreak is occurring. if (insect.ActiveOutbreak) { // PlugIn.ModelCore.UI.WriteLine(" OutbreakStartYear={0}.", insect.OutbreakStartYear); if (insect.OutbreakStartYear == PlugIn.ModelCore.CurrentTime || insect.SingleOutbreakYear) { // Initialize neighborhoodGrowthReduction with patches Outbreak.InitializeDefoliationPatches(insect); } else { insect.NeighborhoodDefoliation.ActiveSiteValues = 0; } } // Now report on the previous year's defoliation, that which has been processed // through biomass succession. Calculations for logfile. double sumDefoliation = 0.0; double meanDefoliation = 0.0; int numSites0_33 = 0; int numSites33_66 = 0; int numSites66_100 = 0; int numInitialSites = 0; int numSitesActive = 0; // Just get a sum of all active sites to calculate mean defoliation accurately for log file. // ONly calculate for log file when outbreak or mortality is active <- Modified, JRF, add to log file each year. if (insect.ActiveOutbreak || insect.SingleOutbreakYear || (insect.LastStopYear + 1 >= PlugIn.ModelCore.CurrentTime) || (insect.LastBioRemoved > 0)) { foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { sumDefoliation += insect.LastYearDefoliation[site]; if (insect.LastYearDefoliation[site] > 0.0 && insect.LastYearDefoliation[site] <= 0.33) { numSites0_33++; } if (insect.LastYearDefoliation[site] > 0.33 && insect.LastYearDefoliation[site] <= 0.66) { numSites33_66++; } if (insect.LastYearDefoliation[site] > 0.66 && insect.LastYearDefoliation[site] <= 1.0) { numSites66_100++; } if (insect.Disturbed[site] && SiteVars.InitialOutbreakProb[site] > 0) { numInitialSites++; } numSitesActive++; } if (insect.OutbreakStartYear == PlugIn.ModelCore.CurrentTime) { insect.InitialSites = numInitialSites; } if (numSites0_33 + numSites33_66 + numSites66_100 > 0) { meanDefoliation = sumDefoliation / (double)numSitesActive; } //PlugIn.ModelCore.UI.WriteLine(" sumDefoliation={0}, numSites={1}.", sumDefoliation, numSites0_33 + numSites33_66 + numSites66_100); } int totalBioRemoved = 0; foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { totalBioRemoved += SiteVars.BiomassRemoved[site]; // kg across all defoliated sites } insect.LastBioRemoved = totalBioRemoved; //Assign variables for the logfile // PlugIn.ModelCore.UI.WriteLine(" totalBioRemoved={0}.", totalBioRemoved); // ONly add to log & output maps during outbreak <- Modified, JRF, add to log file each year. //if ((insect.ActiveOutbreak && insect.OutbreakStartYear < PlugIn.ModelCore.CurrentTime) || (meanDefoliation > 0) || (insect.LastBioRemoved > 0)) //{ /* Old code for log file * if (meanDefoliation > 0) * { * * log.Write("{0},{1},{2},{3},{4:0.0000},{5},{6},{7},{8},{9}", * PlugIn.ModelCore.CurrentTime-1, //0 * insect.Name, //1 * insect.LastStartYear, //2 * insect.LastStopYear, //3 * meanDefoliation, //4 * numSites0_33, //5 * numSites33_66, //6 * numSites66_100, //7 * insect.InitialSites, //8 * insect.LastBioRemoved //9 * ); * } * * else * { * log.Write("{0},{1},{2},{3},{4:0.0000},{5},{6},{7},{8},{9}", * PlugIn.ModelCore.CurrentTime - 1, //0 * insect.Name, //1 * 0, //2 * 0, //3 * meanDefoliation, //4 * numSites0_33, //5 * numSites33_66, //6 * numSites66_100, //7 * 0, //8 * insect.LastBioRemoved //9 * ); * } * * * //foreach (IEcoregion ecoregion in Ecoregions.Dataset) * // log.Write(",{0}", 1); * * log.WriteLine(""); */ eventLog.Clear(); EventsLog el = new EventsLog(); el.Time = PlugIn.ModelCore.CurrentTime - 1; el.InsectName = insect.Name; el.MeanDefoliation = meanDefoliation; el.NumSitesDefoliated0_33 = numSites0_33; el.NumSitesDefoliated33_66 = numSites33_66; el.NumSitesDefoliated66_100 = numSites66_100; el.NumOutbreakInitialSites = insect.InitialSites; el.MortalityBiomass = insect.LastBioRemoved; if (insect.ActiveOutbreak) { el.StartYear = insect.OutbreakStartYear; el.StopYear = insect.OutbreakStopYear; } else { el.StartYear = insect.LastStartYear; el.StopYear = insect.LastStopYear; } eventLog.AddObject(el); eventLog.WriteToFile(); //----- Write Insect Defoliation/GrowthReduction maps -------- string path = MapNames.ReplaceTemplateVars(mapNameTemplate, insect.Name, PlugIn.ModelCore.CurrentTime - 1); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { pixel.MapCode.Value = (short)(insect.LastYearDefoliation[site] * 100.0); } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); } } //----- Write Initial Patch maps -------- string path2 = MapNames.ReplaceTemplateVars(mapNameTemplate, ("InitialPatchMap" + insect.Name), PlugIn.ModelCore.CurrentTime); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path2, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { if (insect.Disturbed[site]) { pixel.MapCode.Value = (short)(SiteVars.InitialOutbreakProb[site] * 100); } else { pixel.MapCode.Value = 0; } } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); //Zero out the InitialOutbreakProb after output maps are written. SiteVars.InitialOutbreakProb[site] = 0; } } //----- Write Biomass Reduction maps -------- string path3 = MapNames.ReplaceTemplateVars(mapNameTemplate, ("BiomassRemoved" + insect.Name), PlugIn.ModelCore.CurrentTime); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path3, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { if (SiteVars.BiomassRemoved[site] > 0) { pixel.MapCode.Value = Math.Max((short)1, (short)(SiteVars.BiomassRemoved[site] / 100)); // convert to Mg/ha } else { pixel.MapCode.Value = 0; } } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); //Zero out the BiomassRemoved after the last insect mortality event in a given year. SiteVars.BiomassRemoved[site] = 0; } } //} //insect.ThisYearDefoliation.ActiveSiteValues = 0.0; //reset this year to 0 for all sites, this was already done at the top of loop to initialize defoliation patchs, Outbreak.cs } }
//--------------------------------------------------------------------- ///<summary> /// Run the BDA extension at a particular timestep. ///</summary> public override void Run() { running = true; PlugIn.ModelCore.UI.WriteLine(" Processing landscape for Biomass Insect events ..."); foreach (IInsect insect in manyInsect) { //SiteVars.BiomassRemoved.ActiveSiteValues = 0; //SiteVars.InitialOutbreakProb.ActiveSiteValues = 0.0; insect.BiomassRemoved.ActiveSiteValues = 0; //insect.InitialOutbreakProb.ActiveSiteValues = 0.0; if (insect.MortalityYear == PlugIn.ModelCore.CurrentTime) { Outbreak.Mortality(insect); } // Copy the data from current to last foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { insect.LastYearDefoliation[site] = insect.ThisYearDefoliation[site]; } insect.ActiveOutbreak = false; PlugIn.ModelCore.NormalDistribution.Mu = 0.0; PlugIn.ModelCore.NormalDistribution.Sigma = 1.0; double randomNum = PlugIn.ModelCore.NormalDistribution.NextDouble(); PlugIn.ModelCore.ExponentialDistribution.Lambda = 1 / (insect.MeanDuration); // Mean = 1/lambda double randomNumE = PlugIn.ModelCore.ExponentialDistribution.NextDouble(); // First, has enough time passed since the last outbreak? double timeBetweenOutbreaks = insect.MeanTimeBetweenOutbreaks + (insect.StdDevTimeBetweenOutbreaks * randomNum); double duration = System.Math.Round(randomNumE + 1); double timeAfterDuration = timeBetweenOutbreaks - duration; //PlugIn.ModelCore.UI.WriteLine("Calculated time between = {0}. inputMeanTime={1}, inputStdTime={2}.", timeBetweenOutbreaks, insect.MeanTimeBetweenOutbreaks, insect.StdDevTimeBetweenOutbreaks); //PlugIn.ModelCore.UI.WriteLine("Calculated duration = {0}. inputMeanDura={1}, inputStdDura={2}.", duration, insect.MeanDuration, insect.StdDevDuration); //PlugIn.ModelCore.UI.WriteLine("Insect Start Time = {0}, Stop Time = {1}.", insect.OutbreakStartYear, insect.OutbreakStopYear); if (PlugIn.ModelCore.CurrentTime == 1) { PlugIn.ModelCore.UI.WriteLine(" Year 1: Setting initial start and stop times."); insect.OutbreakStartYear = (int)(timeBetweenOutbreaks / 2.0) + 1; insect.OutbreakStopYear = insect.OutbreakStartYear + (int)duration - 1; //PlugIn.ModelCore.UI.WriteLine(" {0} is not active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); } else if (insect.OutbreakStartYear <= PlugIn.ModelCore.CurrentTime && insect.OutbreakStopYear >= PlugIn.ModelCore.CurrentTime) { //PlugIn.ModelCore.UI.WriteLine(" An outbreak starts or continues. Start and stop time do not change."); insect.ActiveOutbreak = true; //PlugIn.ModelCore.UI.WriteLine(" {0} is active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); insect.MortalityYear = PlugIn.ModelCore.CurrentTime + 1; } if (insect.OutbreakStopYear <= PlugIn.ModelCore.CurrentTime && timeAfterDuration > PlugIn.ModelCore.CurrentTime - insect.OutbreakStopYear) { //PlugIn.ModelCore.UI.WriteLine(" In between outbreaks, reset start and stop times."); insect.ActiveOutbreak = true; //PlugIn.ModelCore.UI.WriteLine(" {0} is active. StartYear={1}, StopYear={2}, CurrentYear={3}.", insect.Name, insect.OutbreakStartYear, insect.OutbreakStopYear, PlugIn.ModelCore.CurrentTime); insect.MortalityYear = PlugIn.ModelCore.CurrentTime + 1; //insect.OutbreakStartYear = PlugIn.ModelCore.CurrentTime + (int) timeBetweenOutbreaks; //insect.OutbreakStopYear = insect.OutbreakStartYear + (int) duration; } //PlugIn.ModelCore.UI.WriteLine(" Insect Start Time = {0}, Stop Time = {1}.", insect.OutbreakStartYear, insect.OutbreakStopYear); if (insect.ActiveOutbreak) { // PlugIn.ModelCore.UI.WriteLine(" OutbreakStartYear={0}.", insect.OutbreakStartYear); if (insect.OutbreakStartYear == PlugIn.ModelCore.CurrentTime) { // Initialize neighborhoodGrowthReduction with patches Outbreak.InitializeDefoliationPatches(insect); } else { insect.NeighborhoodDefoliation.ActiveSiteValues = 0; } } // Now report on the previous year's defoliation, that which has been processed // through biomass succession. double sumDefoliation = 0.0; int numSites0_33 = 0; int numSites33_66 = 0; int numSites66_100 = 0; int numInitialSites = 0; foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { sumDefoliation += insect.LastYearDefoliation[site]; if (insect.LastYearDefoliation[site] > 0.0 && insect.LastYearDefoliation[site] <= 0.33) { numSites0_33++; } else if (insect.LastYearDefoliation[site] > 0.33 && insect.LastYearDefoliation[site] <= 0.66) { numSites33_66++; } else if (insect.LastYearDefoliation[site] > 0.66 && insect.LastYearDefoliation[site] <= 1.0) { numSites66_100++; } else if (insect.LastYearDefoliation[site] > 1.0) { PlugIn.ModelCore.UI.WriteLine("Defoliation exceeds 1.0"); } if (insect.Disturbed[site] && insect.InitialOutbreakProb[site] > 0) { numInitialSites++; } } if (insect.OutbreakStartYear == PlugIn.ModelCore.CurrentTime) { insect.InitialSites = numInitialSites; } double meanDefoliation = 0.0; if (numSites0_33 + numSites33_66 + numSites66_100 > 0) { meanDefoliation = sumDefoliation / (double)(numSites0_33 + numSites33_66 + numSites66_100); } //PlugIn.ModelCore.UI.WriteLine(" sumDefoliation={0}, numSites={1}.", sumDefoliation, numSites0_33 + numSites33_66 + numSites66_100); int totalBioRemoved = 0; foreach (ActiveSite site in PlugIn.ModelCore.Landscape) { //totalBioRemoved += SiteVars.BiomassRemoved[site]; totalBioRemoved += insect.BiomassRemoved[site]; } //PlugIn.ModelCore.UI.WriteLine(" totalBioRemoved={0}.", totalBioRemoved); // ONly add to log & output maps during outbreak if ((insect.ActiveOutbreak && insect.OutbreakStartYear < PlugIn.ModelCore.CurrentTime) || (meanDefoliation > 0) || (insect.LastBioRemoved > 0)) { if (insect.ActiveOutbreak) { log.Write("{0},{1},{2},{3},{4:0.00},{5},{6},{7},{8},{9}", PlugIn.ModelCore.CurrentTime - 1, //0 insect.Name, //1 insect.OutbreakStartYear, //2 insect.OutbreakStopYear, //3 meanDefoliation, //4 numSites0_33, //5 numSites33_66, //6 numSites66_100, //7 insect.InitialSites, //8 insect.LastBioRemoved //9 ); } else { log.Write("{0},{1},{2},{3},{4:0.00},{5},{6},{7},{8},{9}", PlugIn.ModelCore.CurrentTime - 1, //0 insect.Name, //1 insect.LastStartYear, //2 insect.LastStopYear, //3 meanDefoliation, //4 numSites0_33, //5 numSites33_66, //6 numSites66_100, //7 insect.InitialSites, //8 insect.LastBioRemoved //9 ); } //foreach (IEcoregion ecoregion in Ecoregions.Dataset) // log.Write(",{0}", 1); log.WriteLine(""); //----- Write Insect GrowthReduction maps -------- string path = MapNames.ReplaceTemplateVars(mapNameTemplate, insect.Name, PlugIn.ModelCore.CurrentTime - 1); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { pixel.MapCode.Value = (short)(insect.LastYearDefoliation[site] * 100.0); } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); } } //----- Write Initial Patch maps -------- string path2 = MapNames.ReplaceTemplateVars(mapNameTemplate, ("InitialPatchMap-" + insect.Name), PlugIn.ModelCore.CurrentTime - 1); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path2, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { if (insect.Disturbed[site]) { pixel.MapCode.Value = (short)(insect.InitialOutbreakProb[site] * 100); } else { pixel.MapCode.Value = 0; } } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); } } //----- Write Biomass Reduction maps -------- string path3 = MapNames.ReplaceTemplateVars(mapNameTemplate, ("BiomassRemoved-" + insect.Name), PlugIn.ModelCore.CurrentTime); using (IOutputRaster <ShortPixel> outputRaster = modelCore.CreateRaster <ShortPixel>(path3, modelCore.Landscape.Dimensions)) { ShortPixel pixel = outputRaster.BufferPixel; foreach (Site site in PlugIn.ModelCore.Landscape.AllSites) { if (site.IsActive) { pixel.MapCode.Value = (short)(insect.BiomassRemoved[site] / 100); // convert to Mg/ha } else { // Inactive site pixel.MapCode.Value = 0; } outputRaster.WriteBufferPixel(); } } } insect.ThisYearDefoliation.ActiveSiteValues = 0.0; //reset this year to 0 for all sites insect.LastBioRemoved = totalBioRemoved; if (insect.OutbreakStopYear <= PlugIn.ModelCore.CurrentTime && timeAfterDuration > PlugIn.ModelCore.CurrentTime - insect.OutbreakStopYear) { insect.LastStartYear = insect.OutbreakStartYear; insect.LastStopYear = insect.OutbreakStopYear; insect.OutbreakStartYear = PlugIn.ModelCore.CurrentTime + (int)timeBetweenOutbreaks; insect.OutbreakStopYear = insect.OutbreakStartYear + (int)duration; //PlugIn.ModelCore.UI.WriteLine(" Insect Start Time = {0}, Stop Time = {1}.", insect.OutbreakStartYear, insect.OutbreakStopYear); } } }