WebSVN - public iLand - Diff - Rev 452 and 453 - /src/core/resourceunitspecies.cpp


Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/resourceunitspecies.cpp':
/** @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)
  */
#include "global.h"
#include "resourceunitspecies.h"

#include "species.h"
#include "resourceunit.h"


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);

    mRemovedGrowth = 0.;
    mLastYear = -1;
}


void ResourceUnitSpecies::calculate(const bool fromEstablishment)
{
    if (mLastYear == GlobalSettings::instance()->currentYear())
        return;
    mLastYear = GlobalSettings::instance()->currentYear();

    statistics().clear();
    if (mLAIfactor>0 || fromEstablishment==true) {
        mResponse.calculate();///< calculate environmental responses per species (vpd, temperature, ...)
        m3PG.calculate();///< production of NPP
    } 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.newSaplings() << mSapling.diedSaplings() << mSapling.recruitedSaplings();
        }
    //); // 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);
}
 

Subversion Repositories public iLand

(root)/src/core/resourceunitspecies.cpp @ 1222 - Rev 452 → 453

Rev 452		Rev 453
1	Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/resourceunitspecies.cpp':	1	Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/resourceunitspecies.cpp':
2	/** @class ResourceUnitSpecies	2	/** @class ResourceUnitSpecies
3	The class contains data available at ResourceUnit x Species scale.	3	The class contains data available at ResourceUnit x Species scale.
4	Data stored is either statistical (i.e. number of trees per species) or used	4	Data stored is either statistical (i.e. number of trees per species) or used
5	within the model (e.g. fraction of utilizable Radiation)	5	within the model (e.g. fraction of utilizable Radiation)
6	*/	6	*/
7	#include "global.h"	7	#include "global.h"
8	#include "resourceunitspecies.h"	8	#include "resourceunitspecies.h"
9		9
10	#include "species.h"	10	#include "species.h"
11	#include "resourceunit.h"	11	#include "resourceunit.h"
12		12
13		13
14	double ResourceUnitSpecies::leafArea() const	14	double ResourceUnitSpecies::leafArea() const
15	{	15	{
16	// Leaf area of the species:	16	// Leaf area of the species:
17	// total leaf area on the RU * fraction of leafarea	17	// total leaf area on the RU * fraction of leafarea
18	return mLAIfactor * ru()->leafAreaIndex();	18	return mLAIfactor * ru()->leafAreaIndex();
19	}	19	}
20		20
21	void ResourceUnitSpecies::setup(Species species, ResourceUnit ru)	21	void ResourceUnitSpecies::setup(Species species, ResourceUnit ru)
22	{	22	{
23	mSpecies = species;	23	mSpecies = species;
24	mRU = ru;	24	mRU = ru;
25	mResponse.setup(this);	25	mResponse.setup(this);
26	m3PG.setResponse(&mResponse);	26	m3PG.setResponse(&mResponse);
27	mEstablishment.setup(ru->climate(), this);	27	mEstablishment.setup(ru->climate(), this);
28	mSapling.setup(this);	28	mSapling.setup(this);
29	mStatistics.setResourceUnitSpecies(this);	29	mStatistics.setResourceUnitSpecies(this);
30	mStatisticsDead.setResourceUnitSpecies(this);	30	mStatisticsDead.setResourceUnitSpecies(this);
31	mStatisticsMgmt.setResourceUnitSpecies(this);	31	mStatisticsMgmt.setResourceUnitSpecies(this);
32		32
33	mRemovedGrowth = 0.;	33	mRemovedGrowth = 0.;
34	mLastYear = -1;	34	mLastYear = -1;
35	}	35	}
36		36
37		37
38	void ResourceUnitSpecies::calculate(const bool fromEstablishment)	38	void ResourceUnitSpecies::calculate(const bool fromEstablishment)
39	{	39	{
40	if (mLastYear == GlobalSettings::instance()->currentYear())	40	if (mLastYear == GlobalSettings::instance()->currentYear())
41	return;	41	return;
42	mLastYear = GlobalSettings::instance()->currentYear();	42	mLastYear = GlobalSettings::instance()->currentYear();
43		43
44	statistics().clear();	44	statistics().clear();
45	if (mLAIfactor>0 \|\| fromEstablishment==true) {	45	if (mLAIfactor>0 \|\| fromEstablishment==true) {
46	mResponse.calculate();///< calculate environmental responses per species (vpd, temperature, ...)	46	mResponse.calculate();///< calculate environmental responses per species (vpd, temperature, ...)
47	m3PG.calculate();///< production of NPP	47	m3PG.calculate();///< production of NPP
48	} else {	48	} else {
49	// if no LAI is present, then just clear the respones.	49	// if no LAI is present, then just clear the respones.
50	// note: subject to change when regeneration is added...	50	// note: subject to change when regeneration is added...
51	mResponse.clear();	51	mResponse.clear();
52	m3PG.clear();	52	m3PG.clear();
53	}	53	}
54	}	54	}
55		55
56		56
57	void ResourceUnitSpecies::updateGWL()	57	void ResourceUnitSpecies::updateGWL()
58	{	58	{
59	// removed growth is the running sum of all removed	59	// removed growth is the running sum of all removed
60	// tree volume. the current "GWL" therefore is current volume (standing) + mRemovedGrowth.	60	// tree volume. the current "GWL" therefore is current volume (standing) + mRemovedGrowth.
61	mRemovedGrowth+=statisticsDead().volume() + statisticsMgmt().volume();	61	mRemovedGrowth+=statisticsDead().volume() + statisticsMgmt().volume();
62	}	62	}
63		63
64	void ResourceUnitSpecies::calclulateEstablishment()	64	void ResourceUnitSpecies::calclulateEstablishment()
65	{	65	{
66	mEstablishment.calculate();	66	mEstablishment.calculate();
67	//DBGMODE(	67	//DBGMODE(
68	if (GlobalSettings::instance()->isDebugEnabled(GlobalSettings::dEstablishment)) {	68	if (GlobalSettings::instance()->isDebugEnabled(GlobalSettings::dEstablishment)) {
69	DebugList &out = GlobalSettings::instance()->debugList(ru()->index(), GlobalSettings::dEstablishment);	69	DebugList &out = GlobalSettings::instance()->debugList(ru()->index(), GlobalSettings::dEstablishment);
70	// establishment details	70	// establishment details
71	out << mSpecies->id() << ru()->index();	71	out << mSpecies->id() << ru()->index();
72	out << mEstablishment.avgSeedDensity();	72	out << mEstablishment.avgSeedDensity();
73	out << mEstablishment.TACAminTemp() << mEstablishment.TACAchill() << mEstablishment.TACAfrostFree() << mEstablishment.TACgdd();	73	out << mEstablishment.TACAminTemp() << mEstablishment.TACAchill() << mEstablishment.TACAfrostFree() << mEstablishment.TACgdd();
74	out << mEstablishment.TACAfrostDaysAfterBudBirst() << mEstablishment.abioticEnvironment();	74	out << mEstablishment.TACAfrostDaysAfterBudBirst() << mEstablishment.abioticEnvironment();
75	out << m3PG.fEnvYear() << mEstablishment.avgLIFValue() << mEstablishment.numberEstablished();	75	out << m3PG.fEnvYear() << mEstablishment.avgLIFValue() << mEstablishment.numberEstablished();
76	out << mSapling.livingSaplings() << mSapling.averageHeight() << mSapling.newSaplings() << mSapling.diedSaplings() << mSapling.recruitedSaplings();	76	out << mSapling.livingSaplings() << mSapling.averageHeight() << mSapling.newSaplings() << mSapling.diedSaplings() << mSapling.recruitedSaplings();
77	}	77	}
78	//); // DBGMODE()	78	//); // DBGMODE()
79		79
80		80
81	if ( logLevelDebug() )	81	if ( logLevelDebug() )
82	qDebug() << "establishment of RU" << mRU->index() << "species" << species()->id()	82	qDebug() << "establishment of RU" << mRU->index() << "species" << species()->id()
83	<< "seeds density:" << mEstablishment.avgSeedDensity()	83	<< "seeds density:" << mEstablishment.avgSeedDensity()
84	<< "abiotic environment:" << mEstablishment.abioticEnvironment()	84	<< "abiotic environment:" << mEstablishment.abioticEnvironment()
85	<< "f_env,yr:" << m3PG.fEnvYear()	85	<< "f_env,yr:" << m3PG.fEnvYear()
86	<< "N(established):" << mEstablishment.numberEstablished();	86	<< "N(established):" << mEstablishment.numberEstablished();
87		87
88	}	88	}
89		89
90	void ResourceUnitSpecies::calclulateSaplingGrowth()	90	void ResourceUnitSpecies::calclulateSaplingGrowth()
91	{	91	{
92	mSapling.calculateGrowth();	92	mSapling.calculateGrowth();
-		93	}
-		94
-		95	void ResourceUnitSpecies::visualGrid(Grid<float> &grid) const
-		96	{
-		97	mSapling.fillHeightGrid(grid);
93	}	98	}
94		99