Rev 540 | Rev 546 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 540 | Rev 541 | ||
---|---|---|---|
Line 6... | Line 6... | ||
6 | #include "expression.h"
|
6 | #include "expression.h"
|
7 | // for calculation of climate decomposition
|
7 | // for calculation of climate decomposition
|
8 | #include "resourceunit.h"
|
8 | #include "resourceunit.h"
|
9 | #include "watercycle.h"
|
9 | #include "watercycle.h"
|
10 | #include "climate.h"
|
10 | #include "climate.h"
|
- | 11 | #include "model.h"
|
|
11 | 12 | ||
12 | /** @class Snag
|
13 | /** @class Snag
|
13 | Snag deals with carbon / nitrogen fluxes from the forest until the reach soil pools.
|
14 | Snag deals with carbon / nitrogen fluxes from the forest until the reach soil pools.
|
14 | Snag lives on the level of the ResourceUnit; carbon fluxes from trees enter Snag, and parts of the biomass of snags
|
15 | Snag lives on the level of the ResourceUnit; carbon fluxes from trees enter Snag, and parts of the biomass of snags
|
15 | is subsequently forwarded to the soil sub model.
|
16 | is subsequently forwarded to the soil sub model.
|
Line 147... | Line 148... | ||
147 | /// calculation is done on the level of ResourceUnit because the water content per day is needed.
|
148 | /// calculation is done on the level of ResourceUnit because the water content per day is needed.
|
148 | double Snag::calculateClimateFactors() |
149 | double Snag::calculateClimateFactors() |
149 | {
|
150 | {
|
150 | double psi_kpa; |
151 | double psi_kpa; |
151 | double ft, fw; |
152 | double ft, fw; |
152 | const double min_kpa = -1500.; |
- | |
- | 153 | const double min_kpa = Model::settings().decayReferencePsi; |
|
153 | double f_sum = 0.; |
154 | double f_sum = 0.; |
154 | for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day) |
155 | for (const ClimateDay *day=mRU->climate()->begin(); day!=mRU->climate()->end(); ++day) |
155 | {
|
156 | {
|
156 | psi_kpa = mRU->waterCycle()->psi_kPa(day->day); |
157 | psi_kpa = mRU->waterCycle()->psi_kPa(day->day); |
157 | 158 |