Rev 226 |
Rev 434 |
Go to most recent revision |
Blame |
Compare with Previous |
Last modification |
View Log
| RSS feed
#ifndef PHENOLOGY_H
#define PHENOLOGY_H
class Climate;
class Phenology
{
public:
Phenology() {mClimate=0; mId=0; mMinVpd=mMaxVpd=mMinDayLength=mMaxDayLength=mMinTemp=mMaxTemp=0.; mDayStart=0;mDayEnd=365;}
Phenology(const int id, const Climate* climate, const double minVpd, const double maxVpd,
const double minDayLength, const double maxDayLength,
const double minTemp, const double maxTemp): mId(id), mClimate(climate), mMinVpd(minVpd), mMaxVpd(maxVpd),
mMinDayLength(minDayLength), mMaxDayLength(maxDayLength), mMinTemp(minTemp), mMaxTemp(maxTemp) {}
// access
const int id() const { return mId; }
/// calculate the phenology for the current year
void calculate();
/// get result of phenology calcualtion for this year (a pointer to a array of 12 values between 0..1: 0: no days with foliage)
const double *month() const { return mPhenoFraction; }
int vegetationPeriodLength() const { return mDayEnd - mDayStart; } ///< length of vegetation period in days, returs 365 for evergreens
int vegetationPeriodStart() const { return mDayStart; } ///< day of year when vegeation period starts
int vegetationPeriodEnd() const { return mDayEnd; } ///< day of year when vegeation period stops
private:
int mId; ///< identifier of this Phenology group
const Climate *mClimate; ///< link to relevant climate source
double mMinVpd; ///< minimum vpd [kPa]
double mMaxVpd; ///< maximum vpd [kPa]
double mMinDayLength; ///< minimum daylength [hours]
double mMaxDayLength; ///< maximum daylength [hours]
double mMinTemp; ///< minimum temperature [°]
double mMaxTemp; ///< maximum temperature [°]
double mPhenoFraction[12]; ///< fraction [0..1] of month i [0..11] to are inside the vegetation period, i.e. have leafs
int mDayStart; ///< start of vegetation period (in day of year)
int mDayEnd; ///< end of vegetation period (in days of year, 1.1. = 0)
};
#endif // PHENOLOGY_H