standstatistics.h

Go to the documentation of this file.

/********************************************************************************************
**    iLand - an individual based forest landscape and disturbance model
**    http://iland-model.org
**    Copyright (C) 2009-  Werner Rammer, Rupert Seidl
**
**    This program is free software: you can redistribute it and/or modify
**    it under the terms of the GNU General Public License as published by
**    the Free Software Foundation, either version 3 of the License, or
**    (at your option) any later version.
**
**    This program is distributed in the hope that it will be useful,
**    but WITHOUT ANY WARRANTY; without even the implied warranty of
**    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
**    GNU General Public License for more details.
**
**    You should have received a copy of the GNU General Public License
**    along with this program.  If not, see <http://www.gnu.org/licenses/>.
********************************************************************************************/
 
#ifndef STANDSTATISTICS_H
#define STANDSTATISTICS_H
class Tree;
struct TreeGrowthData;
class ResourceUnitSpecies;
class SaplingStat;
 
class StandStatistics
{
public:
    StandStatistics() { mRUS=0; clear();}
    void setResourceUnitSpecies(const ResourceUnitSpecies *rus) { mRUS=rus; }
 
    void add(const StandStatistics &stat); 
    void addAreaWeighted(const StandStatistics &stat, const double weight); 
    void add(const Tree *tree, const TreeGrowthData *tgd); 
    void addNPP(const TreeGrowthData *tgd); 
 
    void add(const SaplingStat *sapling); 
    void clear(); 
    void clearOnlyTrees(); 
    void calculate(); 
    void calculateAreaWeighted(); 
    // getters
    double count() const { return mCount; }
    double dbh_avg() const { return mAverageDbh; } 
    double height_avg() const { return mAverageHeight; } 
    double volume() const { return mSumVolume; } 
    double gwl() const { return mGWL;} 
    double basalArea() const { return mSumBasalArea; } 
    double leafAreaIndex() const { return mLeafAreaIndex; } 
    double leafAreaIndexSaplings() const { return mLAISaplings; } 
    double npp() const { return mNPP; } 
    double nppAbove() const { return mNPPabove; } 
    double nppSaplings() const { return mNPPsaplings; } 
    int cohortCount() const { return static_cast<int>(mCohortCount); } 
    int saplingCount() const { return static_cast<int>(mSaplingCount); } 
    double saplingAge() const { return mAverageSaplingAge; } 
    double saplingBasalArea() const { return mBasalAreaSaplings; } 
    // carbon/nitrogen cycle
    double cStem() const { return mCStem; }
    double nStem() const { return mNStem; }
    double cBranch() const { return mCBranch; }
    double nBranch() const { return mNBranch; }
    double cFoliage() const { return mCFoliage; }
    double nFoliage() const { return mNFoliage; }
    double cCoarseRoot() const { return mCCoarseRoot; }
    double nCoarseRoot() const { return mNCoarseRoot; }
    double cFineRoot() const { return mCFineRoot; }
    double nFineRoot() const { return mNFineRoot; }
    double cRegeneration() const { return mCRegeneration; }
    double nRegeneration() const { return mNRegeneration; }
    double totalCarbon() const { return mCStem + mCBranch + mCFoliage + mCFineRoot + mCCoarseRoot + mCRegeneration; }
 
private:
    inline void addBiomass(const double biomass, const double CNRatio, double *C, double *N);
    const ResourceUnitSpecies *mRUS; 
    double mCount;
    double mSumDbh;
    double mSumHeight;
    double mSumBasalArea;
    double mSumVolume;
    double mGWL;
    double mAverageDbh;
    double mAverageHeight;
    double mLeafAreaIndex;
    double mNPP;
    double mNPPabove;
    double mNPPsaplings; // carbon gain of saplings (kg Biomass)
    // regeneration layer
    double mCohortCount; 
    double mSaplingCount; 
    double mSumSaplingAge;
    double mAverageSaplingAge;
    double mLAISaplings; 
    double mBasalAreaSaplings; 
    // carbon and nitrogen pools
    double mCStem, mCFoliage, mCBranch, mCCoarseRoot, mCFineRoot;
    double mNStem, mNFoliage, mNBranch, mNCoarseRoot, mNFineRoot;
    double mCRegeneration, mNRegeneration;
};
 
 
class SystemStatistics
{
public:
    SystemStatistics() { reset(); }
    void reset() { treeCount=0; saplingCount=0; newSaplings=0;
                   tManagement = 0.; tApplyPattern=tReadPattern=tTreeGrowth=0.;
                   tSeedDistribution=tSapling=tEstablishment=tCarbonCycle=tWriteOutput=tTotalYear=0.; }
    void writeOutput();
    // the system counters
    int treeCount;
    int saplingCount;
    int newSaplings;
    // timings
    double tManagement;
    double tApplyPattern;
    double tReadPattern;
    double tTreeGrowth;
    double tSeedDistribution;
    double tSapling;
    double tEstablishment;
    double tCarbonCycle;
    double tWriteOutput;
    double tTotalYear;
 
};
 
#endif // STANDSTATISTICS_H