Subversion Repositories public iLand

/** @class SpeciesResponse
    Environmental responses calculates production relevant responses for a tree species on resource unit level.
    SpeciesResponse combines data from different sources and converts information about the environment
    into responses of a species. The spatial level is the "ResourceUnit", where homogenetiy of environmental factors
    is assumed. The temporal aggregation depends on the factor, but usually, the daily environmental data is
    aggregated to monthly response values (which subsequently are used during 3PG production).
 Sources are:
    - vapour pressure deficit (dryness of atmosphere): directly from climate data (daily)
    - soil water status (dryness of soil): TODO (daily)
    - temperature: directly from climate data (daily)
    - phenology: @sa Phenology, combines several sources (quasi-monthly)
    - CO2: @sa SpeciesSet::co2Response() based on ambient CO2 level (climate data), nitrogen and soil water responses (yearly)
    - nitrogen: based on the amount of available nitrogen (yearly)
*/
#include "speciesresponse.h"

#include "resourceunit.h"
#include "species.h"
#include "resourceunitspecies.h"
#include "climate.h"
#include "model.h"

SpeciesResponse::SpeciesResponse()
{
    mSpecies=0;
    mRu=0;
}
void SpeciesResponse::clear()
{
    for (int i=0;i<12;i++)
        mVpdResponse[i]=mSoilWaterResponse[i]=mTempResponse[i]=mRadiation[i]=0.;

    mCO2Response=mNitrogenResponse=mSoilWaterResponseYear=0.;

}
void SpeciesResponse::setup(ResourceUnitSpecies *rus)
{
    mSpecies = rus->species();
    mRu = rus->ru();
    clear();
}

/// Main function that calculates monthly / annual species responses
void SpeciesResponse::calculate()
{
    const ClimateDay *begin, *end;
    double no_of_days;

    clear(); // reset values

    // calculate yearly responses
    mSoilWaterResponseYear = 0.5; // TODO
    double ambient_co2 = mRu->climate()->begin()->co2; // CO2 level of first day of year
    double nitrogen =  Model::settings().nitrogenAvailable;

    // Nitrogen response: a yearly value based on available nitrogen
    mNitrogenResponse = mSpecies->nitrogenResponse( nitrogen );

    // CO2 response: a yearly value based on atmospheric CO2 concentration and the response to nutrients and water in the soil.
    mCO2Response = mSpecies->speciesSet()->co2Response(ambient_co2,
                                                       mNitrogenResponse,
                                                       mSoilWaterResponseYear);

    const Phenology &pheno = mRu->climate()->phenology(mSpecies->phenologyClass());
    int veg_begin = pheno.vegetationPeriodStart();
    int veg_end = pheno.vegetationPeriodEnd();

    // calculate monthly responses
    int doy=0;
    for (int mon=0;mon<12;mon++) {
        mRu->climate()->monthRange(mon, &begin, &end);
        no_of_days = 0;
        for (const ClimateDay *day=begin; day!=end; ++day, no_of_days++) {
            // VPD response
            mVpdResponse[mon]+=mSpecies->vpdResponse( day->vpd );
            // Temperature Response
            mTempResponse[mon]+=mSpecies->temperatureResponse(day->temp_delayed);
            // radiation: only count days in vegetation period
            if (doy>=veg_begin && doy<=veg_end)
                mRadiation[mon] += day->radiation;
            // soil water: fake
            mSoilWaterResponse[mon] = mSoilWaterResponseYear;

            doy++;
        }
        mVpdResponse[mon] /= no_of_days; // vpd: average of month
        mTempResponse[mon] /= no_of_days; // temperature: average value of daily responses
    }
}