Subversion Repositories public iLand

#ifndef PRODUCTION3PG_H

#define PRODUCTION3PG_H

class SpeciesResponse;

class ProductionOut;

class Production3PG

{

public:

    Production3PG();

    void setResponse(const SpeciesResponse *response) { mResponse=response;}

    double calculate(); ///< return  year GPP/rad: kg Biomass/MJ PAR/m2

    void clear(); ///< clear production values

    double rootFraction() const { return mRootFraction; } /// fraction of biomass that should be distributed to roots

    double GPPperArea() const { return mGPPperArea; } ///<  GPP production (yearly) (kg Biomass) per m2 (effective area)

    double fEnvYear() const { return mEnvYear; } ///< f_env,yr: aggregate environmental factor [0..1}

private:

    inline double calculateUtilizablePAR(const int month) const;

    inline double calculateEpsilon(const int month) const;

    inline double abovegroundFraction() const; ///< calculate fraction of biomass

    const SpeciesResponse *mResponse; ///< species specific responses

    double mUPAR[12]; ///< utilizable radiation MJ/m2 and month

    double mGPP[12]; ///< monthly Gross Primary Production gC/MJ radiation

    double mRootFraction; ///< fraction of production that flows into roots

    double mGPPperArea; ///< kg GPP Biomass / m2 interception area

    double mEnvYear; ///< f_env,yr: factor that aggregates the environment for the species over the year (weighted with the radiation pattern)

    friend class ProductionOut;

};

#endif // PRODUCTION3PG_H