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/>.

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

#include "productionout.h"

#include "debugtimer.h"

#include "model.h"

#include "resourceunit.h"

#include "species.h"

#include "speciesresponse.h"

ProductionOut::ProductionOut()

{

    setName("Production per month, species and resource unit", "production_month");

    setDescription("Details about the 3PG production submodule on monthly basis and for each species and resource unit.");

    columns() << OutputColumn::year() << OutputColumn::ru() << OutputColumn::id() << OutputColumn::species()

              << OutputColumn("month", "month of year", OutInteger)

              << OutputColumn("tempResponse", "monthly average of daily respose value temperature", OutDouble)

              << OutputColumn("waterResponse", "monthly average of daily respose value soil water", OutDouble)

              << OutputColumn("vpdResponse", "monthly vapour pressure deficit respose.", OutDouble)

              << OutputColumn("co2Response", "monthly response value for ambient co2.", OutDouble)

              << OutputColumn("nitrogenResponse", "yearly respose value nitrogen", OutDouble)

              << OutputColumn("radiation_m2", "global radiation PAR in MJ per m2 and month", OutDouble)

              << OutputColumn("utilizableRadiation_m2", "utilizable PAR in MJ per m2 and month (sum of daily rad*min(respVpd,respWater,respTemp))", OutDouble)

              << OutputColumn("GPP_kg_m2", "GPP (without Aging) in kg Biomass/m2", OutDouble);

 }

void ProductionOut::setup()

{

}

void ProductionOut::execute(const ResourceUnitSpecies *rus)

{

    const Production3PG &prod = rus->prod3PG();

    const SpeciesResponse *resp = prod.mResponse;

    for (int i=0;i<12;i++) {

        *this << currentYear() << rus->ru()->index() << rus->ru()->id() << rus->species()->id();

        *this << (i+1); // month

        // responses

        *this <<  resp->tempResponse()[i]

              << resp->soilWaterResponse()[i]

              << resp->vpdResponse()[i]

              << resp->co2Response()[i]

              << resp->nitrogenResponse()

              << resp->globalRadiation()[i]

              << prod.mUPAR[i]

              << prod.mGPP[i];

        writeRow();

    }

}

void ProductionOut::exec()

{

    DebugTimer t("ProductionOut");

    Model *m = GlobalSettings::instance()->model();

    foreach(ResourceUnit *ru, m->ruList()) {

        if (ru->id()==-1)

            continue; // do not include if out of project area

        foreach(const ResourceUnitSpecies *rus, ru->ruSpecies()) {

            execute(rus);

        }

    }

}