Rev 331 | Rev 369 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 331 | Rev 367 | ||
|---|---|---|---|
| Line 167... | Line 167... | ||
| 167 | .arg(LAI) |
167 | .arg(LAI) |
| 168 | .arg(interception_fraction) |
168 | .arg(interception_fraction) |
| 169 | .arg(mLRI_modification) |
169 | .arg(mLRI_modification) |
| 170 | .arg(mStockedArea); |
170 | .arg(mStockedArea); |
| 171 | ); |
171 | ); |
| - | 172 | ||
| - | 173 | // calculate LAI fractions
|
|
| - | 174 | ResourceUnitSpecies *i; |
|
| - | 175 | QVector<ResourceUnitSpecies>::iterator iend = mRUSpecies.end(); |
|
| - | 176 | for (i=mRUSpecies.begin(); i!=iend; ++i) { |
|
| - | 177 | i->setLAIfactor(i->statistics().leafAreaIndex() / leafAreaIndex()); |
|
| - | 178 | }
|
|
| - | 179 | ||
| 172 | // soil water model - this determines soil water contents needed for response calculations
|
180 | // soil water model - this determines soil water contents needed for response calculations
|
| 173 | {
|
181 | {
|
| 174 | DebugTimer tw("water:run"); |
182 | DebugTimer tw("water:run"); |
| 175 | mWater->run(); |
183 | mWater->run(); |
| 176 | }
|
184 | }
|
| 177 | 185 | ||
| 178 | // invoke species specific calculation (3PG)
|
186 | // invoke species specific calculation (3PG)
|
| 179 | ResourceUnitSpecies *i; |
- | |
| 180 | QVector<ResourceUnitSpecies>::iterator iend = mRUSpecies.end(); |
- | |
| 181 | - | ||
| 182 | for (i=mRUSpecies.begin(); i!=iend; ++i) { |
187 | for (i=mRUSpecies.begin(); i!=iend; ++i) { |
| 183 | i->statistics().clear(); |
188 | i->statistics().clear(); |
| 184 | i->calculate(); |
189 | i->calculate(); |
| 185 | qDebug() << "species" << (*i).species()->id() << "raw_gpp_m2" << i->prod3PG().GPPperArea() << "area:" << productiveArea() << "gpp:" << productiveArea()*i->prod3PG().GPPperArea(); |
190 | qDebug() << "species" << (*i).species()->id() << "raw_gpp_m2" << i->prod3PG().GPPperArea() << "area:" << productiveArea() << "gpp:" << productiveArea()*i->prod3PG().GPPperArea(); |
| 186 | }
|
191 | }
|