soil.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 SOIL_H
#define SOIL_H
 
#include "snag.h"
struct SoilParams; // forward
class ResourceUnit; // forward
class SoilInputOut; // forward
 
class Soil
{
public:
    // lifecycle
    Soil(ResourceUnit *ru=0);
    void setInitialState(const CNPool &young_labile_kg_ha, const CNPool &young_refractory_kg_ha, const CNPair &SOM_kg_ha, double young_labile_aboveground_frac, double young_refractory_aboveground_frac);
 
    // actions
    void setSoilInput(const CNPool &labile_input_kg_ha, const CNPool &refractory_input_kg_ha, double labile_aboveground_C, double refractory_aboveground_C); 
    void setClimateFactor(const double climate_factor_re) { mRE = climate_factor_re; } 
    void newYear(); 
    void calculateYear(); 
 
    void disturbance(double DWDfrac, double litterFrac, double soilFrac);
    void disturbanceBiomass(double DWD_kg_ha, double litter_kg_ha, double soil_kg_ha);
 
    // access
    const CNPool &youngLabile() const { return mYL;} 
    const CNPool &youngRefractory() const { return mYR;} 
    const CNPair &oldOrganicMatter() const { return mSOM;} 
    double youngLabileAbovegroundFraction() const { return mYLaboveground_frac; } 
    double youngRefractoryAbovegroundFraction() const { return mYRaboveground_frac; } 
    double availableNitrogen() const { return mAvailableNitrogen; } 
    double totalCarbon() const; 
 
    const CNPair &fluxToAtmosphere() const { return mTotalToAtmosphere; } 
    const CNPair &fluxToDisturbance() const { return mTotalToDisturbance; } 
 
    QList<QVariant> debugList(); 
private:
    ResourceUnit *mRU; 
    void fetchParameters(); 
    static SoilParams *mParams; // static container for parameters
    // variables
    double mRE; 
    double mAvailableNitrogen; 
    double mAvailableNitrogenFromLabile; 
    double mAvailableNitrogenFromRefractory; 
    double mKyl; 
    double mKyr; 
    double mKo; 
    double mH; 
 
    CNPool mInputLab; 
    CNPool mInputRef; 
    // state variables
    CNPool mYL; 
    CNPool mYR; 
    CNPair mSOM; 
    // aboveground fraction (=fraction of the pool content for which the source is AG biomass (leafs, branches, stems)
    double mYLaboveground_frac; 
    double mYRaboveground_frac; 
 
    CNPair mTotalToDisturbance; 
    CNPair mTotalToAtmosphere; 
 
    static double mNitrogenDeposition; 
    friend class Snapshot;
    friend class SoilInputOut;
};
 
#endif // SOIL_H