Subversion Repositories public iLand

/********************************************************************************************

**    iLand - an individual based forest landscape and disturbance model

**    http://iland.boku.ac.at

**    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 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 [kg Biomass / m2]

    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