//--------------------------------------------------------------------- /// <summary> /// Grows an individual cohort for a year, incrementing its age by 1 /// and updating its biomass for annual growth and mortality. /// </summary> /// <param name="index"> /// The index of the cohort to grow; it must be between 0 and Count - 1. /// </param> /// <param name="site"> /// The site where the species' cohorts are located. /// </param> /// <returns> /// The index of the next younger cohort. Note this may be the same /// as the index passed in if that cohort dies due to senescence. /// </returns> public int GrowCohort(int index, ActiveSite site, bool annualTimestep) { Debug.Assert(0 <= index && index <= cohortData.Count); Debug.Assert(site != null); Cohort cohort = new Cohort(species, cohortData[index]); Debug.Assert(cohort.WoodBiomass <= cohort.WoodBiomass + cohort.LeafBiomass); Debug.Assert(cohort.LeafBiomass <= cohort.WoodBiomass + cohort.LeafBiomass); //if (isDebugEnabled) // log.DebugFormat(" grow cohort: {0}, {1} yrs, {2} Mg/ha", // cohort.Species.Name, cohort.Age, cohort.Biomass); // Check for senescence if (cohort.Age >= species.Longevity) { RemoveCohort(index, cohort, site, null); return(index); } if (annualTimestep) { cohort.IncrementAge(); } float[] biomassChange = Cohorts.BiomassCalculator.ComputeChange(cohort, site); Debug.Assert(-(cohort.WoodBiomass + cohort.LeafBiomass) <= biomassChange[0] + biomassChange[1]); // Cohort can't loss more biomass than it has //UI.WriteLine("B={0:0.00}, Age={1}, delta={2}", cohort.Biomass, cohort.Age, biomassChange); cohort.ChangeWoodBiomass(biomassChange[0]); cohort.ChangeLeafBiomass(biomassChange[1]); //if (isDebugEnabled) // log.DebugFormat(" biomass: change = {0}, cohort = {1}, site = {2}", // biomassChange, cohort.Biomass, siteBiomass); if (cohort.WoodBiomass + cohort.LeafBiomass > 0) { cohortData[index] = cohort.Data; return(index + 1); } else { RemoveCohort(index, cohort, site, null); return(index); } }
//--------------------------------------------------------------------- /// <summary> /// Grows an individual cohort for a year, incrementing its age by 1 /// and updating its biomass for annual growth and mortality. /// </summary> /// <param name="index"> /// The index of the cohort to grow; it must be between 0 and Count - 1. /// </param> /// <param name="site"> /// The site where the species' cohorts are located. /// </param> /// <returns> /// The index of the next younger cohort. Note this may be the same /// as the index passed in if that cohort dies due to senescence. /// </returns> public int GrowCohort(int index, ActiveSite site, bool annualTimestep) { Debug.Assert(0 <= index && index <= cohortData.Count); Debug.Assert(site != null); Cohort cohort = new Cohort(species, cohortData[index]); Debug.Assert(cohort.WoodBiomass <= cohort.WoodBiomass + cohort.LeafBiomass); Debug.Assert(cohort.LeafBiomass <= cohort.WoodBiomass + cohort.LeafBiomass); // Check for senescence if (cohort.Age >= species.Longevity) { RemoveCohort(index, cohort, site, null); return(index); } if (annualTimestep) { cohort.IncrementAge(); } float[] biomassChange = Cohorts.BiomassCalculator.ComputeChange(cohort, site); Debug.Assert(-(cohort.WoodBiomass + cohort.LeafBiomass) <= biomassChange[0] + biomassChange[1]); // Cohort can't loss more biomass than it has cohort.ChangeWoodBiomass(biomassChange[0]); cohort.ChangeLeafBiomass(biomassChange[1]); if (cohort.WoodBiomass + cohort.LeafBiomass > 0) { cohortData[index] = cohort.Data; return(index + 1); } else { RemoveCohort(index, cohort, site, null); return(index); } }