Subversion Repositories public iLand

#include "global.h"

#include "resourceunitspecies.h"

#include "species.h"

#include "resourceunit.h"

#include "model.h"

#include "snag.h"

#include "watercycle.h"

/** @class ResourceUnitSpecies

    The class contains data available at ResourceUnit x Species scale.

    Data stored is either statistical (i.e. number of trees per species) or used

    within the model (e.g. fraction of utilizable Radiation).

    Important submodules are:

    * 3PG production (Production3PG)

    * Establishment

    * Growth and Recruitment of Saplings

    * Snag dynamics

  */

ResourceUnitSpecies::~ResourceUnitSpecies()

{

    if (mSnag)

        delete mSnag;

    mSnag = 0;

}

double ResourceUnitSpecies::leafArea() const

{

    // Leaf area of the species:

    // total leaf area on the RU * fraction of leafarea

    return mLAIfactor * ru()->leafAreaIndex();

}

void ResourceUnitSpecies::setup(Species *species, ResourceUnit *ru)

{

    mSpecies = species;

    mRU = ru;

    mResponse.setup(this);

    m3PG.setResponse(&mResponse);

    mEstablishment.setup(ru->climate(), this);

    mSapling.setup(this);

    mStatistics.setResourceUnitSpecies(this);

    mStatisticsDead.setResourceUnitSpecies(this);

    mStatisticsMgmt.setResourceUnitSpecies(this);

    if (mSnag)

        delete mSnag;

    mSnag=0;

    if (Model::settings().carbonCycleEnabled) {

       mSnag = new Snag;

       mSnag->setup(mRU);

    }

    mRemovedGrowth = 0.;

    mLastYear = -1;

}

void ResourceUnitSpecies::calculate(const bool fromEstablishment)

{

    if (mLastYear == GlobalSettings::instance()->currentYear())

        return;

    statistics().clear();

    if (mLAIfactor>0 || fromEstablishment==true) {

        // execute the water calculation...

        if (fromEstablishment)

            const_cast<WaterCycle*>(mRU->waterCycle())->run(); // run the water sub model (only if this has not be done already)

        mResponse.calculate();// calculate environmental responses per species (vpd, temperature, ...)

        m3PG.calculate();// production of NPP

        mLastYear = GlobalSettings::instance()->currentYear(); // mark this year as processed

    } else {

        // if no LAI is present, then just clear the respones.

        // note: subject to change when regeneration is added...

        mResponse.clear();

        m3PG.clear();

    }

}

void ResourceUnitSpecies::updateGWL()

{

    // removed growth is the running sum of all removed

    // tree volume. the current "GWL" therefore is current volume (standing) + mRemovedGrowth.

    mRemovedGrowth+=statisticsDead().volume() + statisticsMgmt().volume();

}

void ResourceUnitSpecies::calclulateEstablishment()

{

    mEstablishment.calculate();

    //DBGMODE(

        if (GlobalSettings::instance()->isDebugEnabled(GlobalSettings::dEstablishment)) {

            DebugList &out = GlobalSettings::instance()->debugList(ru()->index(), GlobalSettings::dEstablishment);

            // establishment details

            out << mSpecies->id() << ru()->index();

            out << mEstablishment.avgSeedDensity();

            out << mEstablishment.TACAminTemp() << mEstablishment.TACAchill() << mEstablishment.TACAfrostFree() << mEstablishment.TACgdd();

            out << mEstablishment.TACAfrostDaysAfterBudBirst() << mEstablishment.abioticEnvironment();

            out << m3PG.fEnvYear() << mEstablishment.avgLIFValue() << mEstablishment.numberEstablished();

            out << mSapling.livingSaplings() << mSapling.averageHeight() << mSapling.averageAge() << mSapling.averageDeltaHPot() << mSapling.averageDeltaHRealized();

            out << mSapling.newSaplings() << mSapling.diedSaplings() << mSapling.recruitedSaplings() << mSpecies->saplingGrowthParameters().referenceRatio;

        }

    //); // DBGMODE()

    if ( logLevelDebug() )

        qDebug() << "establishment of RU" << mRU->index() << "species" << species()->id()

        << "seeds density:" << mEstablishment.avgSeedDensity()

        << "abiotic environment:" << mEstablishment.abioticEnvironment()

        << "f_env,yr:" << m3PG.fEnvYear()

        << "N(established):" << mEstablishment.numberEstablished();

}

void ResourceUnitSpecies::calclulateSaplingGrowth()

{

    mSapling.calculateGrowth();

}

void ResourceUnitSpecies::visualGrid(Grid<float> &grid) const

{

    mSapling.fillHeightGrid(grid);

}

void ResourceUnitSpecies::calculateSnagDynamics()

{

    if (!snag())

        return;

    snag()->processYear();

    if (GlobalSettings::instance()->isDebugEnabled(GlobalSettings::dSnagDynamics) && !snag()->isEmpty()) {

        DebugList &out = GlobalSettings::instance()->debugList(ru()->index(), GlobalSettings::dSnagDynamics);

        out << mSpecies->id() << ru()->index() << snag()->debugList(); // use private variables...

    }

}