Rev 376 | Rev 387 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 376 | Rev 381 | ||
|---|---|---|---|
| Line 84... | Line 84... | ||
| 84 | mOpacity = 1. - exp(- Model::settings().lightExtinctionCoefficientOpacity * mLeafArea / mStamp->crownArea()); |
84 | mOpacity = 1. - exp(- Model::settings().lightExtinctionCoefficientOpacity * mLeafArea / mStamp->crownArea()); |
| 85 | mNPPReserve = (1+species()->finerootFoliageRatio())*mFoliageMass; // initial value |
85 | mNPPReserve = (1+species()->finerootFoliageRatio())*mFoliageMass; // initial value |
| 86 | mDbhDelta = 0.1; // initial value: used in growth() to estimate diameter increment |
86 | mDbhDelta = 0.1; // initial value: used in growth() to estimate diameter increment |
| 87 | 87 | ||
| 88 | // initial value for tree aging...
|
88 | // initial value for tree aging...
|
| 89 | mRU->addTreeAging(mLeafArea,mSpecies->aging(mHeight, mAge)); |
- | |
| - | 89 | mRU->addTreeAging(mLeafArea,mSpecies->aging(mHeight, mAge,flag(Tree::TreeHasRealAge))); |
|
| 90 | }
|
90 | }
|
| 91 | 91 | ||
| 92 | //////////////////////////////////////////////////
|
92 | //////////////////////////////////////////////////
|
| 93 | //// Light functions (Pattern-stuff)
|
93 | //// Light functions (Pattern-stuff)
|
| 94 | //////////////////////////////////////////////////
|
94 | //////////////////////////////////////////////////
|
| Line 483... | Line 483... | ||
| 483 | double raw_gpp_per_area = mRU->resourceUnitSpecies(species()).prod3PG().GPPperArea(); |
483 | double raw_gpp_per_area = mRU->resourceUnitSpecies(species()).prod3PG().GPPperArea(); |
| 484 | // (2) GPP (without aging-effect) in kg Biomass / year
|
484 | // (2) GPP (without aging-effect) in kg Biomass / year
|
| 485 | double raw_gpp = raw_gpp_per_area * effective_area; |
485 | double raw_gpp = raw_gpp_per_area * effective_area; |
| 486 | 486 | ||
| 487 | // apply aging according to the state of the individuum
|
487 | // apply aging according to the state of the individuum
|
| 488 | const double aging_factor = mSpecies->aging(mHeight, mAge); |
- | |
| - | 488 | const double aging_factor = mSpecies->aging(mHeight, mAge, flag(Tree::TreeHasRealAge)); |
|
| 489 | mRU->addTreeAging(mLeafArea, aging_factor); |
489 | mRU->addTreeAging(mLeafArea, aging_factor); |
| 490 | double gpp = raw_gpp * aging_factor; // |
490 | double gpp = raw_gpp * aging_factor; // |
| 491 | d.NPP = gpp * 0.47; // respiration loss, cf. Waring et al 1998. |
491 | d.NPP = gpp * 0.47; // respiration loss, cf. Waring et al 1998. |
| 492 | 492 | ||
| 493 | //DBGMODE(
|
493 | //DBGMODE(
|