climate.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 CLIMATE_H
#define CLIMATE_H
 
#include <QtSql>
#include "phenology.h"
struct ClimateDay
{
    int year; // year
    int month; // month (1..12)
    int dayOfMonth; // day of the month (1..31)
    double temperature; // average day  degree C (of the light hours)
    double min_temperature; // minimum temperature of the day
    double max_temperature; // maximum temperature of the day
    double mean_temp() const { return (min_temperature + max_temperature) / 2.; } // mean temperature
    double temp_delayed; // temperature delayed (after Maekela, 2008) for response calculations
    double preciptitation; // sum of day [mm]
    double radiation; // sum of day (MJ/m2)
    double vpd; // average of day [kPa] = [0.1 mbar] (1 bar = 100kPa)
    static double co2; // ambient CO2 content in ppm
    QString toString() const { return QString("%1.%2.%3").arg(dayOfMonth).arg(month).arg(year); }
    bool isValid() const  { return (year>=0); }
    int id() const { return year*10000 + month*100 + dayOfMonth; }
};
 
class Sun
{
public:
    void setup(const double latitude_rad);
    QString dump();
    const double &daylength(const int day) const { return mDaylength_h[day]; }
    int longestDay() const { return mDayWithMaxLength; }
    bool northernHemishere() const { return mDayWithMaxLength<300; }
    int dayShorter10_5hrs() const { return mDayWith10_5hrs; }
    int dayShorter14_5hrs() const { return mDayWith14_5hrs; }
private:
    double mLatitude; 
    int mDayWithMaxLength; 
    double mDaylength_h[366]; 
    int mDayWith10_5hrs; // last day of year with a day length > 10.5 hours (see Establishment)
    int mDayWith14_5hrs; // last doy with at least 14.5 hours of day length
};
 
class Climate
{
public:
    Climate();
    void setup(); 
    bool isSetup() const { return mIsSetup; }
    const QString &name() const { return mName; } 
    // activity
    void nextYear();
    // access to climate data
    const ClimateDay *dayOfYear(const int dayofyear) const { return mBegin + dayofyear;} 
    const ClimateDay *day(const int month, const int day) const; 
    int whichDayOfYear(const ClimateDay *climate) const {return climate-mBegin; } 
    void monthRange(const int month, const ClimateDay **rBegin, const ClimateDay **rEnd) const;
    double days(const int month) const; 
    int daysOfYear() const; 
    const ClimateDay *begin() const { return mBegin; } 
    const ClimateDay *end() const { return mEnd; } 
    void toDate(const int yearday, int *rDay=0, int *rMonth=0, int *rYear=0) const; 
    //
    double totalRadiation() const { return mAnnualRadiation; } 
    const double* precipitationMonth() const { return mPrecipitationMonth; }
    double meanAnnualTemperature() const { return mMeanAnnualTemperature; }
    double annualPrecipitation() const { double r=0.; for (int i=0;i<12;++i) r+=mPrecipitationMonth[i]; return r;}
    const double *temperatureMonth() const { return mTemperatureMonth; }
    int climateDataYear() const { return mBegin->year; }
    // access to other subsystems
    int phenologyGroupCount() const { return mPhenology.count(); }
    const Phenology &phenology(const int phenologyGroup) const; 
    const Sun &sun() const { return mSun; } 
    double daylength_h(const int doy) const { return sun().daylength(doy); } 
 
private:
    bool mIsSetup;
    bool mDoRandomSampling; 
    bool mTMaxAvailable; 
    QString mName;
    Sun mSun; 
    void load(); 
    void setupPhenology(); 
    void climateCalculations(const ClimateDay &lastDay); 
    ClimateDay mInvalidDay;
    int mLoadYears; // count of years to load ahead
    int mCurrentYear; // current year (relative)
    int mMinYear; // lowest year in store (relative)
    int mMaxYear;  // highest year in store (relative)
    double mTemperatureShift; // add this to daily temp
    double mPrecipitationShift; // multiply prec with that
    ClimateDay *mBegin; // pointer to the first day of the current year
    ClimateDay *mEnd; // pointer to the last day of the current year (+1)
    QVector<ClimateDay> mStore; 
    QVector<int> mDayIndices; 
    QSqlQuery mClimateQuery; 
    QList<Phenology> mPhenology; 
    QVector<int> mRandomYearList; 
    int mRandomListIndex; 
    double mAnnualRadiation;  
    double mPrecipitationMonth[12]; 
    double mTemperatureMonth[12]; 
    double mMeanAnnualTemperature; 
    static QVector<int> sampled_years; 
};
 
#endif // CLIMATE_H