establishment.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 ESTABLISHMENT_H
#define ESTABLISHMENT_H
#include <QtCore/QPoint>
class Climate;
class ResourceUnitSpecies;
 
class Establishment
{
public:
    Establishment();
    Establishment(const Climate *climate, const ResourceUnitSpecies *rus);
    void setup(const Climate *climate, const ResourceUnitSpecies *rus);
    void clear();
    void calculateAbioticEnvironment(); 
    void writeDebugOutputs();
    // some informations after execution
    double avgSeedDensity() const { return mPxDensity;} 
    double abioticEnvironment() const {return mPAbiotic; } 
    int numberEstablished() const { return mNumberEstablished; } 
    bool TACAminTemp() const { return mTACA_min_temp;} 
    bool TACAchill() const { return mTACA_chill;} 
    bool TACgdd() const { return mTACA_gdd;} 
    int growingDegreeDays() const { return mGDD; } 
    bool TACAfrostFree() const { return mTACA_frostfree;} 
    int TACAfrostDaysAfterBudBirst() const { return mTACA_frostAfterBuds; } 
    double avgLIFValue() const { return mLIFcount>0?mSumLIFvalue/double(mLIFcount):0.; } 
    double waterLimitation() const { return mWaterLimitation; } 
 
private:
    double mPAbiotic; 
    double calculateWaterLimitation();
    double calculateSOLDepthLimitation();
 
    const Climate *mClimate; 
    const ResourceUnitSpecies *mRUS; 
    // some statistics
    double mPxDensity;
    int mNumberEstablished; // number of established trees in the current year
    // TACA switches
    bool mTACA_min_temp; // minimum temperature threshold
    bool mTACA_chill;  // (total) chilling requirement
    bool mTACA_gdd;   // gdd-thresholds
    int mGDD; // growing degree days
    bool mTACA_frostfree; // frost free days in vegetation period
    int mTACA_frostAfterBuds; // frost days after bud birst
    double mSumLIFvalue;
    double mWaterLimitation; // scalar 0..1 signifying the drought limitation of establishment
    int mLIFcount;
 
};
 
#endif // ESTABLISHMENT_H