Rev 863 | Rev 936 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
1 | |||
187 | iland | 2 | #ifndef RESOURCEUNITSPECIES_H |
3 | #define RESOURCEUNITSPECIES_H |
||
115 | Werner | 4 | #include "production3pg.h" |
180 | werner | 5 | #include "standstatistics.h" |
193 | werner | 6 | #include "speciesresponse.h" |
440 | werner | 7 | #include "establishment.h" |
450 | werner | 8 | #include "sapling.h" |
453 | werner | 9 | #include "grid.h" |
475 | werner | 10 | #include "snag.h" |
111 | Werner | 11 | class Species; |
187 | iland | 12 | class ResourceUnit; |
111 | Werner | 13 | |
187 | iland | 14 | class ResourceUnitSpecies |
111 | Werner | 15 | { |
16 | public: |
||
521 | werner | 17 | ResourceUnitSpecies() : mLAIfactor(0.), mSpecies(0), mRU(0) {} |
468 | werner | 18 | ~ResourceUnitSpecies(); |
234 | werner | 19 | void setup(Species *species, ResourceUnit *ru); |
115 | Werner | 20 | |
449 | werner | 21 | // access |
209 | werner | 22 | const SpeciesResponse *speciesResponse() const { return &mResponse; } |
208 | werner | 23 | const Species *species() const { return mSpecies; } ///< return pointer to species |
24 | const ResourceUnit *ru() const { return mRU; } ///< return pointer to resource unit |
||
228 | werner | 25 | const Production3PG &prod3PG() const { return m3PG; } ///< the 3pg production model of this speies x resourceunit |
504 | werner | 26 | const Sapling &sapling() const { return mSapling; } ///< sapling growth submodel |
600 | werner | 27 | Sapling &changeSapling() { return mSapling; } ///< sapling growth submodel (non-const access) |
208 | werner | 28 | StandStatistics &statistics() { return mStatistics; } ///< statistics of this species on the resourceunit |
278 | werner | 29 | StandStatistics &statisticsDead() { return mStatisticsDead; } ///< statistics of died trees |
30 | StandStatistics &statisticsMgmt() { return mStatisticsMgmt; } ///< statistics of removed trees |
||
262 | werner | 31 | const StandStatistics &constStatistics() const { return mStatistics; } ///< const accessor |
32 | const StandStatistics &constStatisticsDead() const { return mStatisticsDead; } ///< const accessor |
||
278 | werner | 33 | const StandStatistics &constStatisticsMgmt() const { return mStatisticsMgmt; } ///< const accessor |
34 | |||
449 | werner | 35 | // actions |
277 | werner | 36 | void updateGWL(); |
37 | double removedVolume() const { return mRemovedGrowth; } ///< sum of volume with was remvoved because of death/management (m3) |
||
720 | werner | 38 | /// relative fraction of LAI of this species (0..1) (if total LAI on resource unit is >= 1, then the sum of all LAIfactors of all species = 1) |
39 | double LAIfactor() const { return mLAIfactor; } |
||
370 | werner | 40 | void setLAIfactor(const double newLAIfraction) { mLAIfactor=newLAIfraction; if (mLAIfactor<0 || mLAIfactor>1.00001) qDebug() << "invalid LAIfactor"<<mLAIfactor; } |
376 | werner | 41 | // properties |
42 | double leafArea() const; ///< total leaf area of the species on the RU (m2). |
||
115 | Werner | 43 | // action |
440 | werner | 44 | void calculate(const bool fromEstablishment=false); ///< calculate response for species, calculate actual 3PG production |
468 | werner | 45 | // establishment, sapling growth |
863 | werner | 46 | void calculateEstablishment(); ///< perform establishment calculations |
450 | werner | 47 | void calclulateSaplingGrowth(); ///< growth of saplings |
911 | werner | 48 | int addSapling(const QPoint &position, const float height=0.05f) { return mSapling.addSapling(position, height); } ///< add a saplings on a given position |
454 | werner | 49 | void clearSaplings(const QPoint &position) { mSapling.clearSaplings(position);} ///< clear saplings on a given position (after recruitment) |
462 | werner | 50 | bool hasSaplingAt(const QPoint &position) const { return mSapling.hasSapling(position); } ///< return true if a sapling of the current speices is present at 'position' |
453 | werner | 51 | // visualization/graphical output |
52 | void visualGrid(Grid<float> &grid) const; |
||
115 | Werner | 53 | |
111 | Werner | 54 | private: |
454 | werner | 55 | ResourceUnitSpecies(const ResourceUnitSpecies &); // no copy |
56 | ResourceUnitSpecies &operator=(const ResourceUnitSpecies &); // no copy |
||
376 | werner | 57 | double mLAIfactor; ///< relative amount of this species' LAI on this resource unit (0..1). Is calculated once a year. |
277 | werner | 58 | double mRemovedGrowth; ///< m3 volume of trees removed/managed (to calculate GWL) |
262 | werner | 59 | StandStatistics mStatistics; ///< statistics of a species on this resource unit |
60 | StandStatistics mStatisticsDead; ///< statistics of died trees (this year) of a species on this resource unit |
||
278 | werner | 61 | StandStatistics mStatisticsMgmt; ///< statistics of removed trees (this year) of a species on this resource unit |
234 | werner | 62 | Production3PG m3PG; ///< NPP prodution unit of this species |
63 | SpeciesResponse mResponse; ///< calculation and storage of species specific respones on this resource unit |
||
440 | werner | 64 | Establishment mEstablishment; ///< establishment for seedlings and sapling growth |
450 | werner | 65 | Sapling mSapling; ///< saplings storage/growth |
468 | werner | 66 | Species *mSpecies; ///< link to speices |
67 | ResourceUnit *mRU; ///< link to resource unit |
||
438 | werner | 68 | int mLastYear; |
111 | Werner | 69 | }; |
70 | |||
71 | #endif // RESSOURCEUNITSPECIES_H |