WebSVN - public iLand - Diff - Rev 211 and 213


Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/speciesset.cpp':
#include <QtCore>
#include <QtSql>
#include "global.h"
#include "xmlhelper.h"
#include "speciesset.h"
#include "species.h"

SpeciesSet::SpeciesSet()
{
    mSetupQuery = 0;
}

SpeciesSet::~SpeciesSet()
{
   clear();
}

void SpeciesSet::clear()
{
    qDeleteAll(mSpecies.values());
    mSpecies.clear();
    mActiveSpecies.clear();
}

const Species *SpeciesSet::species(const int &index)
{
    foreach(Species *s, mSpecies)
        if (s->index() == index)
            return s;
    return NULL;
}

/** loads active species from a database table and creates/setups the species.
    The function uses the global database-connection.
  */
int SpeciesSet::setup()
{
    // setup phenology
    setupPhenology();

    const XmlHelper &xml = GlobalSettings::instance()->settings();
    QString tableName = xml.value("model.species.source", "species");
    QString readerFile = xml.value("model.species.reader", "reader.bin");
    readerFile = GlobalSettings::instance()->path(readerFile, "lip");
    mReaderStamp.load(readerFile);

    QSqlQuery query(GlobalSettings::instance()->dbin());
    mSetupQuery = &query;
    QString sql = QString("select * from %1").arg(tableName);
    query.exec(sql);
    clear();
    qDebug() << "attempting to load a species set from" << tableName;
    while (query.next()) {
        if (var("active").toInt()==0)
            continue;

        Species *s = new Species(this); // create
        // call setup routine (which calls SpeciesSet::var() to retrieve values
        s->setup();

        mSpecies.insert(s->id(), s); // store
        if (s->active())
            mActiveSpecies.append(s);
    } // while query.next()
    qDebug() << "loaded" << mSpecies.count() << "active species:";
    qDebug() << mSpecies.keys();

    mSetupQuery = 0;

    // setup nitrogen response
    XmlHelper resp(xml.node("model.species.nitrogenResponseClasses"));
    if (!resp.isValid())
        throw IException("model.species.nitrogenResponseClasses not present!");
    mNitrogen_1a = resp.valueDouble("class_1_a");
    mNitrogen_1b = resp.valueDouble("class_1_b");
    mNitrogen_2a = resp.valueDouble("class_2_a");
    mNitrogen_2b = resp.valueDouble("class_2_b");
    mNitrogen_3a = resp.valueDouble("class_3_a");
    mNitrogen_3b = resp.valueDouble("class_3_b");
    if (mNitrogen_1a*mNitrogen_1b*mNitrogen_2a*mNitrogen_2b*mNitrogen_3a*mNitrogen_3b == 0)
        throw IException("at least one parameter of model.species.nitrogenResponseClasses is not valid (value=0)!");

    // setup CO2 response
    XmlHelper co2(xml.node("model.species.CO2Response"));
    mCO2base = co2.valueDouble("baseConcentration");
    mCO2comp = co2.valueDouble("compensationPoint");
    mCO2beta0 = co2.valueDouble("beta0");
    mCO2p0 = co2.valueDouble("p0");
    if (mCO2base*mCO2comp*(mCO2base-mCO2comp)*mCO2beta0*mCO2p0==0)
        throw IException("at least one parameter of model.species.CO2Response is not valid!");

    return mSpecies.count();

}

void SpeciesSet::setupPhenology()
{
    mPhenology.clear();
    mPhenology.push_back(Phenology()); // id=0
    XmlHelper xml(GlobalSettings::instance()->settings().node("model.species.phenology"));
    int i=0;
    do {
        QDomElement n = xml.node(QString("type[%1]").arg(i));
        if (n.isNull())
            break;
        i++;
        int id;
        id = n.attribute("id", "-1").toInt();
        if (id<0) throw IException(QString("Error setting up phenology: id invalid\ndump: %1").arg(xml.dump("")));
 
        xml.setCurrentNode(n);
        Phenology item( id,
                        xml.valueDouble(".vpdMin",0.5), // use relative access to node (".x")
                        xml.valueDouble(".vpdMax", 5),
                        xml.valueDouble(".minDayLength",10),
                        xml.valueDouble(".maxDayLength",11),
                        xml.valueDouble(".tempMin", 2),
                        xml.valueDouble(".tempMax", 9) );

        mPhenology.push_back(item);
    } while(true);

}

/** retrieves variables from the datasource available during the setup of species.
  */
QVariant SpeciesSet::var(const QString& varName)
{
    Q_ASSERT(mSetupQuery!=0);

    int idx = mSetupQuery->record().indexOf(varName);
    if (idx>=0)
        return mSetupQuery->value(idx);
    throw IException(QString("SpeciesSet: variable not set: %1").arg(varName));
    //throw IException(QString("load species parameter: field %1 not found!").arg(varName));
    // lookup in defaults
    //qDebug() << "variable" << varName << "not found - using default.";
    //return GlobalSettings::instance()->settingDefaultValue(varName);
}

inline double SpeciesSet::nitrogenResponse(const double &availableNitrogen, const double &NA, const double &NB) const
{
    if (availableNitrogen<=NB)
        return 0;
    double x = 1. - exp(NA * (availableNitrogen-NB));
    return x;
}

/// calculate nitrogen response for a given amount of available nitrogen and a respone class
/// for fractional values, the response value is interpolated between the fixedly defined classes (1,2,3)
double SpeciesSet::nitrogenResponse(const double availableNitrogen, const double &responseClass) const
{
    double value1, value2, value3;
    if (responseClass>2.) {
        if (responseClass==3.)
            return nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);
        else {
            // interpolate between 2 and 3
            value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
            value3 = nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);
            return value2 + (responseClass-2)*(value3-value2);
        }
    }
    if (responseClass==2)
        return nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
    if (responseClass==1)
        return nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);
    // last ressort: interpolate between 1 and 2
    value1 = nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);
    value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
    return value1 + (responseClass-1)*(value2-value1);
}

/** calculation for the CO2 response for the ambientCO2 for the water- and nitrogen responses given.
    The calculation follows Friedlingsstein 1995 (see also links to equations in code)
*/
double SpeciesSet::co2Response(const double ambientCO2, const double nitrogenResponse, const double soilWaterResponse) const
{
    if (nitrogenResponse==0)
        return 0.;

    double co2_water = 2. - soilWaterResponse;
    double beta = mCO2beta0 * co2_water * nitrogenResponse;

    double r =1. +  M_LN2 * beta; // NPP increase for a doubling of atmospheric CO2 (Eq. 17)

    // fertilization function (cf. Farquhar, 1980) based on Michaelis-Menten expressions
    double deltaC = mCO2base - mCO2comp;
    double K2 = ((2*mCO2base - mCO2comp) - r*deltaC ) / ((r-1.)*deltaC*(2*mCO2base - mCO2comp)); // Eq. 16
    double K1 = (1. + K2*deltaC) / deltaC;

    double response = mCO2p0 * K1*(ambientCO2 - mCO2comp) / (1 + K2*(ambientCO2-mCO2comp)); // Eq. 16
    return response;

}

/** return the phenology of the group... */
const Phenology &SpeciesSet::phenology(const int phenologyGroup)
{
    const Phenology &p = mPhenology.at(phenologyGroup);
    if (p.id() == phenologyGroup)
        return p;
    // search...
    for (int i=0;i<mPhenology.count();i++)
        if (mPhenology[i].id()==phenologyGroup)
            return mPhenology[i];
    throw IException(QString("Error at SpeciesSEt::phenology(): invalid group: %1").arg(phenologyGroup));
}
 

Rev 211	Rev 213
1	Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/speciesset.cpp':	1	Redirecting to URL 'https://iland.boku.ac.at/svn/iland/tags/release_1.0/src/core/speciesset.cpp':
2	#include <QtCore>	2	#include <QtCore>
3	#include <QtSql>	3	#include <QtSql>
4	#include "global.h"	4	#include "global.h"
5	#include "xmlhelper.h"	5	#include "xmlhelper.h"
6	#include "speciesset.h"	6	#include "speciesset.h"
7	#include "species.h"	7	#include "species.h"
8		8
9	SpeciesSet::SpeciesSet()	9	SpeciesSet::SpeciesSet()
10	{	10	{
11	mSetupQuery = 0;	11	mSetupQuery = 0;
12	}	12	}
13		13
14	SpeciesSet::~SpeciesSet()	14	SpeciesSet::~SpeciesSet()
15	{	15	{
16	clear();	16	clear();
17	}	17	}
18		18
19	void SpeciesSet::clear()	19	void SpeciesSet::clear()
20	{	20	{
21	qDeleteAll(mSpecies.values());	21	qDeleteAll(mSpecies.values());
22	mSpecies.clear();	22	mSpecies.clear();
23	mActiveSpecies.clear();	23	mActiveSpecies.clear();
24	}	24	}
25		25
26	const Species *SpeciesSet::species(const int &index)	26	const Species *SpeciesSet::species(const int &index)
27	{	27	{
28	foreach(Species *s, mSpecies)	28	foreach(Species *s, mSpecies)
29	if (s->index() == index)	29	if (s->index() == index)
30	return s;	30	return s;
31	return NULL;	31	return NULL;
32	}	32	}
33		33
34	/** loads active species from a database table and creates/setups the species.	34	/** loads active species from a database table and creates/setups the species.
35	The function uses the global database-connection.	35	The function uses the global database-connection.
36	*/	36	*/
37	int SpeciesSet::setup()	37	int SpeciesSet::setup()
38	{	38	{
39	// setup phenology	39	// setup phenology
40	setupPhenology();	40	setupPhenology();
41		41
42	const XmlHelper &xml = GlobalSettings::instance()->settings();	42	const XmlHelper &xml = GlobalSettings::instance()->settings();
43	QString tableName = xml.value("model.species.source", "species");	43	QString tableName = xml.value("model.species.source", "species");
44	QString readerFile = xml.value("model.species.reader", "reader.bin");	44	QString readerFile = xml.value("model.species.reader", "reader.bin");
45	readerFile = GlobalSettings::instance()->path(readerFile, "lip");	45	readerFile = GlobalSettings::instance()->path(readerFile, "lip");
46	mReaderStamp.load(readerFile);	46	mReaderStamp.load(readerFile);
47		47
48	QSqlQuery query(GlobalSettings::instance()->dbin());	48	QSqlQuery query(GlobalSettings::instance()->dbin());
49	mSetupQuery = &query;	49	mSetupQuery = &query;
50	QString sql = QString("select * from %1").arg(tableName);	50	QString sql = QString("select * from %1").arg(tableName);
51	query.exec(sql);	51	query.exec(sql);
52	clear();	52	clear();
53	qDebug() << "attempting to load a species set from" << tableName;	53	qDebug() << "attempting to load a species set from" << tableName;
54	while (query.next()) {	54	while (query.next()) {
55	if (var("active").toInt()==0)	55	if (var("active").toInt()==0)
56	continue;	56	continue;
57		57
58	Species *s = new Species(this); // create	58	Species *s = new Species(this); // create
59	// call setup routine (which calls SpeciesSet::var() to retrieve values	59	// call setup routine (which calls SpeciesSet::var() to retrieve values
60	s->setup();	60	s->setup();
61		61
62	mSpecies.insert(s->id(), s); // store	62	mSpecies.insert(s->id(), s); // store
63	if (s->active())	63	if (s->active())
64	mActiveSpecies.append(s);	64	mActiveSpecies.append(s);
65	} // while query.next()	65	} // while query.next()
66	qDebug() << "loaded" << mSpecies.count() << "active species:";	66	qDebug() << "loaded" << mSpecies.count() << "active species:";
67	qDebug() << mSpecies.keys();	67	qDebug() << mSpecies.keys();
68		68
69	mSetupQuery = 0;	69	mSetupQuery = 0;
70		70
71	// setup nitrogen response	71	// setup nitrogen response
72	XmlHelper resp(xml.node("model.species.nitrogenResponseClasses"));	72	XmlHelper resp(xml.node("model.species.nitrogenResponseClasses"));
73	if (!resp.isValid())	73	if (!resp.isValid())
74	throw IException("model.species.nitrogenResponseClasses not present!");	74	throw IException("model.species.nitrogenResponseClasses not present!");
75	mNitrogen_1a = resp.valueDouble("class_1_a");	75	mNitrogen_1a = resp.valueDouble("class_1_a");
76	mNitrogen_1b = resp.valueDouble("class_1_b");	76	mNitrogen_1b = resp.valueDouble("class_1_b");
77	mNitrogen_2a = resp.valueDouble("class_2_a");	77	mNitrogen_2a = resp.valueDouble("class_2_a");
78	mNitrogen_2b = resp.valueDouble("class_2_b");	78	mNitrogen_2b = resp.valueDouble("class_2_b");
79	mNitrogen_3a = resp.valueDouble("class_3_a");	79	mNitrogen_3a = resp.valueDouble("class_3_a");
80	mNitrogen_3b = resp.valueDouble("class_3_b");	80	mNitrogen_3b = resp.valueDouble("class_3_b");
81	if (mNitrogen_1amNitrogen_1bmNitrogen_2amNitrogen_2bmNitrogen_3a*mNitrogen_3b == 0)	81	if (mNitrogen_1amNitrogen_1bmNitrogen_2amNitrogen_2bmNitrogen_3a*mNitrogen_3b == 0)
82	throw IException("at least one parameter of model.species.nitrogenResponseClasses is not valid (value=0)!");	82	throw IException("at least one parameter of model.species.nitrogenResponseClasses is not valid (value=0)!");
83		83
84	// setup CO2 response	84	// setup CO2 response
85	XmlHelper co2(xml.node("model.species.CO2Response"));	85	XmlHelper co2(xml.node("model.species.CO2Response"));
86	mCO2base = co2.valueDouble("baseConcentration");	86	mCO2base = co2.valueDouble("baseConcentration");
87	mCO2comp = co2.valueDouble("compensationPoint");	87	mCO2comp = co2.valueDouble("compensationPoint");
88	mCO2beta0 = co2.valueDouble("beta0");	88	mCO2beta0 = co2.valueDouble("beta0");
89	mCO2p0 = co2.valueDouble("p0");	89	mCO2p0 = co2.valueDouble("p0");
90	if (mCO2basemCO2comp(mCO2base-mCO2comp)mCO2beta0mCO2p0==0)	90	if (mCO2basemCO2comp(mCO2base-mCO2comp)mCO2beta0mCO2p0==0)
91	throw IException("at least one parameter of model.species.CO2Response is not valid!");	91	throw IException("at least one parameter of model.species.CO2Response is not valid!");
92		92
93	return mSpecies.count();	93	return mSpecies.count();
94		94
95	}	95	}
96		96
97	void SpeciesSet::setupPhenology()	97	void SpeciesSet::setupPhenology()
98	{	98	{
99	mPhenology.clear();	99	mPhenology.clear();
100	mPhenology.push_back(Phenology()); // id=0	100	mPhenology.push_back(Phenology()); // id=0
101	XmlHelper xml(GlobalSettings::instance()->settings().node("model.species.phenology"));	101	XmlHelper xml(GlobalSettings::instance()->settings().node("model.species.phenology"));
102	int i=0;	102	int i=0;
103	do {	103	do {
104	QDomElement n = xml.node(QString("type[%1]").arg(i));	104	QDomElement n = xml.node(QString("type[%1]").arg(i));
105	if (n.isNull())	105	if (n.isNull())
106	break;	106	break;
107	i++;	107	i++;
108	int id;	108	int id;
109	id = n.attribute("id", "-1").toInt();	109	id = n.attribute("id", "-1").toInt();
110	if (id<0) throw IException(QString("Error setting up phenology: id invalid\ndump: %1").arg(xml.dump("")));	110	if (id<0) throw IException(QString("Error setting up phenology: id invalid\ndump: %1").arg(xml.dump("")));
111	xml.setCurrentNode(QString("type[%1]").arg(i));	-
-		111	xml.setCurrentNode(n);
112	Phenology item( id,	112	Phenology item( id,
113	xml.valueDouble(".vpdMin",0.5), // use relative access to node (".x")	113	xml.valueDouble(".vpdMin",0.5), // use relative access to node (".x")
114	xml.valueDouble(".vpdMax", 5),	114	xml.valueDouble(".vpdMax", 5),
115	xml.valueDouble(".minDayLength",10),	115	xml.valueDouble(".minDayLength",10),
116	xml.valueDouble(".maxDayLength",11),	116	xml.valueDouble(".maxDayLength",11),
117	xml.valueDouble(".tempMin", 2),	117	xml.valueDouble(".tempMin", 2),
118	xml.valueDouble(".tempMax", 9) );	118	xml.valueDouble(".tempMax", 9) );
119		119
120	mPhenology.push_back(item);	120	mPhenology.push_back(item);
121	} while(true);	121	} while(true);
122		122
123	}	123	}
124		124
125	/** retrieves variables from the datasource available during the setup of species.	125	/** retrieves variables from the datasource available during the setup of species.
126	*/	126	*/
127	QVariant SpeciesSet::var(const QString& varName)	127	QVariant SpeciesSet::var(const QString& varName)
128	{	128	{
129	Q_ASSERT(mSetupQuery!=0);	129	Q_ASSERT(mSetupQuery!=0);
130		130
131	int idx = mSetupQuery->record().indexOf(varName);	131	int idx = mSetupQuery->record().indexOf(varName);
132	if (idx>=0)	132	if (idx>=0)
133	return mSetupQuery->value(idx);	133	return mSetupQuery->value(idx);
134	throw IException(QString("SpeciesSet: variable not set: %1").arg(varName));	134	throw IException(QString("SpeciesSet: variable not set: %1").arg(varName));
135	//throw IException(QString("load species parameter: field %1 not found!").arg(varName));	135	//throw IException(QString("load species parameter: field %1 not found!").arg(varName));
136	// lookup in defaults	136	// lookup in defaults
137	//qDebug() << "variable" << varName << "not found - using default.";	137	//qDebug() << "variable" << varName << "not found - using default.";
138	//return GlobalSettings::instance()->settingDefaultValue(varName);	138	//return GlobalSettings::instance()->settingDefaultValue(varName);
139	}	139	}
140		140
141	inline double SpeciesSet::nitrogenResponse(const double &availableNitrogen, const double &NA, const double &NB) const	141	inline double SpeciesSet::nitrogenResponse(const double &availableNitrogen, const double &NA, const double &NB) const
142	{	142	{
143	if (availableNitrogen<=NB)	143	if (availableNitrogen<=NB)
144	return 0;	144	return 0;
145	double x = 1. - exp(NA * (availableNitrogen-NB));	145	double x = 1. - exp(NA * (availableNitrogen-NB));
146	return x;	146	return x;
147	}	147	}
148		148
149	/// calculate nitrogen response for a given amount of available nitrogen and a respone class	149	/// calculate nitrogen response for a given amount of available nitrogen and a respone class
150	/// for fractional values, the response value is interpolated between the fixedly defined classes (1,2,3)	150	/// for fractional values, the response value is interpolated between the fixedly defined classes (1,2,3)
151	double SpeciesSet::nitrogenResponse(const double availableNitrogen, const double &responseClass) const	151	double SpeciesSet::nitrogenResponse(const double availableNitrogen, const double &responseClass) const
152	{	152	{
153	double value1, value2, value3;	153	double value1, value2, value3;
154	if (responseClass>2.) {	154	if (responseClass>2.) {
155	if (responseClass==3.)	155	if (responseClass==3.)
156	return nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);	156	return nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);
157	else {	157	else {
158	// interpolate between 2 and 3	158	// interpolate between 2 and 3
159	value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);	159	value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
160	value3 = nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);	160	value3 = nitrogenResponse(availableNitrogen, mNitrogen_3a, mNitrogen_3b);
161	return value2 + (responseClass-2)*(value3-value2);	161	return value2 + (responseClass-2)*(value3-value2);
162	}	162	}
163	}	163	}
164	if (responseClass==2)	164	if (responseClass==2)
165	return nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);	165	return nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
166	if (responseClass==1)	166	if (responseClass==1)
167	return nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);	167	return nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);
168	// last ressort: interpolate between 1 and 2	168	// last ressort: interpolate between 1 and 2
169	value1 = nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);	169	value1 = nitrogenResponse(availableNitrogen, mNitrogen_1a, mNitrogen_1b);
170	value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);	170	value2 = nitrogenResponse(availableNitrogen, mNitrogen_2a, mNitrogen_2b);
171	return value1 + (responseClass-1)*(value2-value1);	171	return value1 + (responseClass-1)*(value2-value1);
172	}	172	}
173		173
174	/** calculation for the CO2 response for the ambientCO2 for the water- and nitrogen responses given.	174	/** calculation for the CO2 response for the ambientCO2 for the water- and nitrogen responses given.
175	The calculation follows Friedlingsstein 1995 (see also links to equations in code)	175	The calculation follows Friedlingsstein 1995 (see also links to equations in code)
176	*/	176	*/
177	double SpeciesSet::co2Response(const double ambientCO2, const double nitrogenResponse, const double soilWaterResponse) const	177	double SpeciesSet::co2Response(const double ambientCO2, const double nitrogenResponse, const double soilWaterResponse) const
178	{	178	{
179	if (nitrogenResponse==0)	179	if (nitrogenResponse==0)
180	return 0.;	180	return 0.;
181		181
182	double co2_water = 2. - soilWaterResponse;	182	double co2_water = 2. - soilWaterResponse;
183	double beta = mCO2beta0 * co2_water * nitrogenResponse;	183	double beta = mCO2beta0 * co2_water * nitrogenResponse;
184		184
185	double r =1. + M_LN2 * beta; // NPP increase for a doubling of atmospheric CO2 (Eq. 17)	185	double r =1. + M_LN2 * beta; // NPP increase for a doubling of atmospheric CO2 (Eq. 17)
186		186
187	// fertilization function (cf. Farquhar, 1980) based on Michaelis-Menten expressions	187	// fertilization function (cf. Farquhar, 1980) based on Michaelis-Menten expressions
188	double deltaC = mCO2base - mCO2comp;	188	double deltaC = mCO2base - mCO2comp;
189	double K2 = ((2mCO2base - mCO2comp) - rdeltaC ) / ((r-1.)deltaC(2*mCO2base - mCO2comp)); // Eq. 16	189	double K2 = ((2mCO2base - mCO2comp) - rdeltaC ) / ((r-1.)deltaC(2*mCO2base - mCO2comp)); // Eq. 16
190	double K1 = (1. + K2*deltaC) / deltaC;	190	double K1 = (1. + K2*deltaC) / deltaC;
191		191
192	double response = mCO2p0 * K1(ambientCO2 - mCO2comp) / (1 + K2(ambientCO2-mCO2comp)); // Eq. 16	192	double response = mCO2p0 * K1(ambientCO2 - mCO2comp) / (1 + K2(ambientCO2-mCO2comp)); // Eq. 16
193	return response;	193	return response;
194		194
195	}	195	}
196		196
197	/** return the phenology of the group... */	197	/** return the phenology of the group... */
198	const Phenology &SpeciesSet::phenology(const int phenologyGroup)	198	const Phenology &SpeciesSet::phenology(const int phenologyGroup)
199	{	199	{
200	const Phenology &p = mPhenology.at(phenologyGroup);	200	const Phenology &p = mPhenology.at(phenologyGroup);
201	if (p.id() == phenologyGroup)	201	if (p.id() == phenologyGroup)
202	return p;	202	return p;
203	// search...	203	// search...
204	for (int i=0;i<mPhenology.count();i++)	204	for (int i=0;i<mPhenology.count();i++)
205	if (mPhenology[i].id()==phenologyGroup)	205	if (mPhenology[i].id()==phenologyGroup)
206	return mPhenology[i];	206	return mPhenology[i];
207	throw IException(QString("Error at SpeciesSEt::phenology(): invalid group: %1").arg(phenologyGroup));	207	throw IException(QString("Error at SpeciesSEt::phenology(): invalid group: %1").arg(phenologyGroup));
208	}	208	}
209		209

Subversion Repositories public iLand

(root)/src/core/speciesset.cpp @ 1222 - Rev 211 → 213