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 "waterout.h"

#include "output.h"

#include "model.h"

#include "watercycle.h"

#include "resourceunit.h"

#include "climate.h"

WaterOut::WaterOut()

{

    setName("Water output", "water");

    setDescription("Annual water cycle output on resource unit/landscape unit.\n" \

                   "The output includes annual averages of precipitation, evapotranspiration, water excess, " \

                   "snow cover, and radiation input. The spatial resolution is landscape averages and/or resource unit level (i.e. 100m pixels). " \

                   "Landscape level averages are indicated by -1 for the 'ru' and 'index' columns.\n\n" \

                   "You can specify a 'condition' to limit output execution to specific years (variable 'year'). " \

                   "The 'conditionRU' can be used to suppress resource-unit-level details; eg. specifying 'in(year,100,200,300)' limits output on reosurce unit level to the years 100,200,300 " \

                   "(leaving 'conditionRU' blank enables details per default).");

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

              << OutputColumn("stocked_area", "area (ha/ha) which is stocked (covered by crowns, absorbing radiation)", OutDouble)

              << OutputColumn("stockable_area", "area (ha/ha) which is stockable (and within the project area)", OutDouble)

              << OutputColumn("precipitation_mm", "Annual precipitation sum (mm)", OutDouble)

              << OutputColumn("et_mm", "Evapotranspiration (mm)", OutDouble)

              << OutputColumn("excess_mm", "annual sum of water loss due to lateral outflow/groundwater flow (mm)", OutDouble)

              << OutputColumn("snowcover_days", "days with snowcover >0mm", OutInteger)

              << OutputColumn("total_radiation", "total incoming radiation over the year (MJ/m2), sum of data in climate input)", OutDouble)

              << OutputColumn("radiation_snowcover", "sum of radiation input (MJ/m2) for days with snow cover", OutInteger);

}

void WaterOut::exec()

{

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

    // global condition

    if (!mCondition.isEmpty() && mCondition.calculate(GlobalSettings::instance()->currentYear())==0.)

        return;

    bool ru_level = true;

    // switch off details if this is indicated in the conditionRU option

    if (!mConditionDetails.isEmpty() && mConditionDetails.calculate(GlobalSettings::instance()->currentYear())==0.)

        ru_level = false;

    double ru_count = 0.;

    int snow_days = 0;

    double et=0., excess=0., rad=0., snow_rad=0., p=0.;

    double stockable=0., stocked=0.;

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

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

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

        const WaterCycle *wc = ru->waterCycle();

        if (ru_level) {

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

            *this << ru->stockedArea()/cRUArea << ru->stockableArea()/cRUArea;

            *this << ru->climate()->annualPrecipitation();

            *this << wc->mTotalET << wc->mTotalExcess;

            *this << wc->mSnowDays;

            *this << ru->climate()->totalRadiation() << wc->mSnowRad;

            writeRow();

        }

        ++ru_count;

        stockable+=ru->stockableArea(); stocked+=ru->stockedArea();

        p+=ru->climate()->annualPrecipitation();

        et+=wc->mTotalET; excess+=wc->mTotalExcess; snow_days+=wc->mSnowDays;

        rad+=ru->climate()->totalRadiation();

        snow_rad+=wc->mSnowRad;

    }

    // write landscape sums

    if (ru_count==0.)

        return;

    *this << currentYear() << -1 << -1; // codes -1/-1 for landscape level

    *this << stocked/ru_count/cRUArea << stockable/ru_count/cRUArea;

    *this << p / ru_count; // mean precip

    *this << et / ru_count;

    *this << excess / ru_count;

    *this << snow_days / ru_count;

    *this << rad / ru_count << snow_rad / ru_count;

    writeRow();

}

void WaterOut::setup()

{

    // use a condition for to control execuation for the current year

    QString condition = settings().value(".condition", "");

    mCondition.setExpression(condition);

    condition = settings().value(".conditionRU", "");

    mConditionDetails.setExpression(condition);

}